error while writing variable: zeta

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
shiyf
Posts: 51
Joined: Mon Sep 09, 2013 11:09 pm
Location: OUC

error while writing variable: zeta

#1 Unread post by shiyf »

hi,
i am a new mode guy,when i run my application first,there are some errors:
error while writing variable: zeta;

the log file:

Code: Select all

 Model Input Parameters:  ROMS/TOMS version 3.6  
                          Thursday - November 21, 2013 -  7:42:55 PM
 -----------------------------------------------------------------------------

 SHIYF

 Operating system : Linux
 CPU/hardware     : x86_64
 Compiler system  : ifort
 Compiler command : /lustre/inspur/mpi/openmpi-intel/bin/mpif90
 Compiler flags   : -heap-arrays -fp-model precise -ip -O3 -free -free -free

 Input Script  : ocean_shiyf.in

 SVN Root URL  : https://www.myroms.org/svn/src/trunk
 SVN Revision  : Unversioned directory

 Local Root    : /lustre/home/shiyongfangWKJ/roms/trunk
 Header Dir    : /lustre/home/shiyongfangWKJ/roms/projects/roms_shiyf
 Header file   : shiyf.h
 Analytical Dir: /lustre/home/shiyongfangWKJ/roms/projects/roms_shiyf

 Resolution, Grid 01: 0449x0160x016,  Parallel Nodes:   9,  Tiling: 003x003


 Physical Parameters, Grid: 01
 =============================

       1440  ntimes          Number of timesteps for 3-D equations.
    300.000  dt              Timestep size (s) for 3-D equations.
         30  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
          1  ERstr           Starting ensemble/perturbation run number.
          1  ERend           Ending ensemble/perturbation run number.
          0  nrrec           Number of restart records to read from disk.
          T  LcycleRST       Switch to recycle time-records in restart file.
        288  nRST            Number of timesteps between the writing of data
                               into restart fields.
          1  ninfo           Number of timesteps between print of information
                               to standard output.
          T  ldefout         Switch to create a new output NetCDF file(s).
         72  nHIS            Number of timesteps between the writing fields
                               into history file.
          1  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
         72  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
 0.0000E+00  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 0.0000E+00  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 5.0000E+00  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 1.0000E-06  Akt_bak(01)     Background vertical mixing coefficient (m2/s)
                               for tracer 01: temp
 1.0000E-06  Akt_bak(02)     Background vertical mixing coefficient (m2/s)
                               for tracer 02: salt
 1.0000E-05  Akv_bak         Background vertical mixing coefficient (m2/s)
                               for momentum.
 5.0000E-06  Akk_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent energy.
 5.0000E-06  Akp_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent generic statistical field.
      3.000  gls_p           GLS stability exponent.
      1.500  gls_m           GLS turbulent kinetic energy exponent.
     -1.000  gls_n           GLS turbulent length scale exponent.
 7.6000E-06  gls_Kmin        GLS minimum value of turbulent kinetic energy.
 1.0000E-12  gls_Pmin        GLS minimum value of dissipation.
 5.4770E-01  gls_cmu0        GLS stability coefficient.
 1.4400E+00  gls_c1          GLS shear production coefficient.
 1.9200E+00  gls_c2          GLS dissipation coefficient.
-4.0000E-01  gls_c3m         GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p         GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk        GLS constant Schmidt number for TKE.
 1.3000E+00  gls_sigp        GLS constant Schmidt number for PSI.
   1400.000  charnok_alpha   Charnok factor for Zos calculation.
      0.500  zos_hsig_alpha  Factor for Zos calculation using Hsig(Awave).
      0.250  sz_alpha        Factor for Wave dissipation surface tke flux .
    100.000  crgban_cw       Factor for Craig/Banner surface tke flux.
 3.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2           Quadratic bottom drag coefficient.
 2.0000E-02  Zob             Bottom roughness (m).
 2.0000E-02  Zos             Surface roughness (m).
          1  lmd_Jwt         Jerlov water type.
          1  Vtransform      S-coordinate transformation equation.
          1  Vstretching     S-coordinate stretching function.
 6.0000E+00  theta_s         S-coordinate surface control parameter.
 0.0000E+00  theta_b         S-coordinate bottom  control parameter.
     10.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
      0.000  dstart          Time-stamp assigned to model initialization (days).
       0.00  time_ref        Reference time for units attribute (yyyymmdd.dd)
 2.0000E+00  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 2.0000E+00  Tnudg(02)       Nudging/relaxation time scale (days)
                               for tracer 02: salt
 0.0000E+00  Znudg           Nudging/relaxation time scale (days)
                               for free-surface.
 0.0000E+00  M2nudg          Nudging/relaxation time scale (days)
                               for 2D momentum.
 0.0000E+00  M3nudg          Nudging/relaxation time scale (days)
                               for 3D momentum.
 0.0000E+00  obcfac          Factor between passive and active
                               open boundary conditions.
          F  VolCons(1)      NLM western  edge boundary volume conservation.
          T  VolCons(2)      NLM southern edge boundary volume conservation.
          T  VolCons(3)      NLM eastern  edge boundary volume conservation.
          F  VolCons(4)      NLM northern edge boundary volume conservation.
     14.000  T0              Background potential temperature (C) constant.
     35.000  S0              Background salinity (PSU) constant.
      1.000  gamma2          Slipperiness variable: free-slip (1.0) or 
                                                    no-slip (-1.0).
          T  Hout(idFsur)    Write out free-surface.
          T  Hout(idUbar)    Write out 2D U-momentum component.
          T  Hout(idVbar)    Write out 2D V-momentum component.
          T  Hout(idUvel)    Write out 3D U-momentum component.
          T  Hout(idVvel)    Write out 3D V-momentum component.
          T  Hout(idWvel)    Write out W-momentum component.
          T  Hout(idOvel)    Write out omega vertical velocity.
          T  Hout(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt

          T  Aout(idFsur)    Write out averaged free-surface.
          T  Aout(idUbar)    Write out averaged 2D U-momentum component.
          T  Aout(idVbar)    Write out averaged 2D V-momentum component.
          T  Aout(idUvel)    Write out averaged 3D U-momentum component.
          T  Aout(idVvel)    Write out averaged 3D V-momentum component.
          T  Aout(idWvel)    Write out averaged W-momentum component.
          T  Aout(idOvel)    Write out averaged omega vertical velocity.
          T  Aout(idTvar)    Write out averaged tracer 01: temp
          T  Aout(idTvar)    Write out averaged tracer 02: salt

 Output/Input Files:

             Output Restart File:  ocean_rst.nc
             Output History File:  ocean_his.nc
            Output Averages File:  ocean_avg.nc
                 Input Grid File:  shiyf_grd.nc
    Input Nonlinear Initial File:  shiyf_ini.nc
           Input Forcing File 01:  shiyf_frc.nc
             Input Boundary File:  shiyf_bry.nc

 Tile partition information for Grid 01:  0449x0160x0016  tiling: 003x003

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1      150        1       53   127200
        1      151      300        1       53   127200
        2      301      449        1       53   126352
        3        1      150       54      107   129600
        4      151      300       54      107   129600
        5      301      449       54      107   128736
        6        1      150      108      160   127200
        7      151      300      108      160   127200
        8      301      449      108      160   126352

 Tile minimum and maximum fractional grid coordinates:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50   150.50     0.50    53.50  RHO-points
        1   150.50   300.50     0.50    53.50  RHO-points
        2   300.50   449.50     0.50    53.50  RHO-points
        3     0.50   150.50    53.50   107.50  RHO-points
        4   150.50   300.50    53.50   107.50  RHO-points
        5   300.50   449.50    53.50   107.50  RHO-points
        6     0.50   150.50   107.50   160.50  RHO-points
        7   150.50   300.50   107.50   160.50  RHO-points
        8   300.50   449.50   107.50   160.50  RHO-points

        0     1.00   150.50     0.50    53.50    U-points
        1   150.50   300.50     0.50    53.50    U-points
        2   300.50   449.00     0.50    53.50    U-points
        3     1.00   150.50    53.50   107.50    U-points
        4   150.50   300.50    53.50   107.50    U-points
        5   300.50   449.00    53.50   107.50    U-points
        6     1.00   150.50   107.50   160.50    U-points
        7   150.50   300.50   107.50   160.50    U-points
        8   300.50   449.00   107.50   160.50    U-points

        0     0.50   150.50     1.00    53.50    V-points
        1   150.50   300.50     1.00    53.50    V-points
        2   300.50   449.50     1.00    53.50    V-points
        3     0.50   150.50    53.50   107.50    V-points
        4   150.50   300.50    53.50   107.50    V-points
        5   300.50   449.50    53.50   107.50    V-points
        6     0.50   150.50   107.50   160.00    V-points
        7   150.50   300.50   107.50   160.00    V-points
        8   300.50   449.50   107.50   160.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     324
               HaloSizeJ(1) =     132
                TileSide(1) =     156
                TileSize(1) =    9360


 Lateral Boundary Conditions: NLM
 ============================

 Variable               Grid    West Edge   South Edge  East Edge   North Edge
 ---------              ----    ----------  ----------  ----------  ----------

 zeta                     1     Closed      Clamped     Clamped     Closed

 ubar                     1     Closed      Clamped     Clamped     Closed

 vbar                     1     Closed      Clamped     Clamped     Closed

 u                        1     Closed      Clamped     Clamped     Closed

 v                        1     Closed      Clamped     Clamped     Closed

 temp                     1     Closed      Clamped     Clamped     Closed

 salt                     1     Closed      Clamped     Clamped     Closed

 tke                      1     Closed      Clamped     Clamped     Closed

 Activated C-preprocessing Options:

 SHIYF               SHIYF
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 AVERAGES            Writing out time-averaged nonlinear model fields.
 CURVGRID            Orthogonal curvilinear grid.
 DIFF_GRID           Horizontal diffusion coefficient scaled by grid size.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 GLS_MIXING          Generic Length-Scale turbulence closure.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 MIX_GEO_TS          Mixing of tracers along geopotential surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 NO_WRITE_GRID       Not Writing grid arrays into NetCDF ouput files.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 RADIATION_2D        Use tangential phase speed in radiation conditions.
 RST_SINGLE          Single precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SCORRECTION         Surface salinity flux correction.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 SPLINES             Conservative parabolic spline reconstruction.
 TS_U3HADVECTION     Third-order upstream horizontal advection of tracers.
 TS_C4VADVECTION     Fourth-order centered vertical advection of tracers.
 TS_DIF2             Harmonic mixing of tracers.
 UV_ADV              Advection of momentum.
 UV_COR              Coriolis term.
 UV_U3HADVECTION     Third-order upstream horizontal advection of 3D momentum.
 UV_C4VADVECTION     Fourth-order centered vertical advection of momentum.
 UV_LDRAG            Linear bottom stress.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 VISC_GRID           Horizontal viscosity coefficient scaled by grid size.

 Process Information:

 Node #  4 (pid=     963) is active.
 Node #  8 (pid=   29422) is active.
 Node #  3 (pid=     962) is active.
 Node #  0 (pid=   27658) is active.
 Node #  5 (pid=     964) is active.
 Node #  2 (pid=   27660) is active.
 Node #  1 (pid=   27659) is active.
 Node #  6 (pid=   29420) is active.
 Node #  7 (pid=   29421) is active.

 INITIAL: Configuring and initializing forward nonlinear model ...


 Vertical S-coordinate System: 

 level   S-coord     Cs-curve   Z   at hmin       at hc    half way     at hmax

    16   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    15  -0.0625000  -0.0019030       -0.891      -0.625      -5.506     -10.121
    14  -0.1250000  -0.0040767       -1.821      -1.250     -11.707     -21.593
    13  -0.1875000  -0.0068304       -2.831      -1.875     -19.395     -35.959
    12  -0.2500000  -0.0105560       -3.978      -2.500     -29.576     -55.174
    11  -0.3125000  -0.0157835       -5.335      -3.125     -43.610     -81.884
    10  -0.3750000  -0.0232566       -7.006      -3.750     -63.403    -119.801
     9  -0.4375000  -0.0340388       -9.140      -4.375     -91.684    -174.228
     8  -0.5000000  -0.0496640      -11.953      -5.000    -132.388    -252.823
     7  -0.5625000  -0.0723554      -15.755      -5.625    -191.217    -366.678
     6  -0.6250000  -0.1053416      -20.998      -6.250    -276.451    -531.904
     5  -0.6875000  -0.1533158      -28.339      -6.875    -400.130    -771.921
     4  -0.7500000  -0.2231040      -38.735      -7.500    -579.762   -1120.789
     3  -0.8125000  -0.3246355      -53.574      -8.125    -840.815   -1628.056
     2  -0.8750000  -0.4723564      -74.880      -8.750   -1220.344   -2365.809
     1  -0.9375000  -0.6872846     -105.595      -9.375   -1772.260   -3438.925
     0  -1.0000000  -1.0000000     -150.000     -10.000   -2575.000   -5000.000

 Time Splitting Weights: ndtfast =  30    nfast =  42

    Primary            Secondary            Accumulated to Current Step

  1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
  2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
  3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
  4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
  5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
  6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
  7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
  8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
  9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
 10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
 11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
 12 0.0043042707117221 0.0335513677379153-0.0022367614257356 0.4026337758325831
 13 0.0060106451121704 0.0334078920475245 0.0037738836864348 0.4360416678801076
 14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
 15 0.0099910761708920 0.0329439123123590 0.0216737068001213 0.5021931174029188
 16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429148
 17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276978
 18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780105
 19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976138
 21 0.0256737146312911 0.0297153720153352 0.1342125541712342 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524359 0.7189082721405746
 23 0.0316613792205220 0.0279045862015855 0.1945237931729578 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236653 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893981 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
 27 0.0432936737565710 0.0230872677537126 0.3506756411451924 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772280 0.8685254094681656
 29 0.0481720587108284 0.0201154903974257 0.4447073467880565 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256415 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325521 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348098 0.9524525304259084
 34 0.0531891349131379 0.0115726061288397 0.7060109510479478 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213239 0.9738247715198165
 36 0.0503158038019030 0.0080591141492892 0.8085423793232269 0.9818838856691057
 37 0.0473537721847153 0.0063819206892258 0.8558961515079423 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497611 0.9930692679747334
 39 0.0376397765791564 0.0033640675316746 0.9367177506289175 0.9964333355064080
 40 0.0305543017255206 0.0021094083123694 0.9672720523544381 0.9985427438187774
 41 0.0217382098544504 0.0010909315881854 0.9890102622088885 0.9996336754069628
 42 0.0109897377911118 0.0003663245930371 1.0000000000000004 0.9999999999999999

 ndtfast, nfast =   30  42   nfast/ndtfast =  1.40000

 Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

    1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

 Minimum X-grid spacing, DXmin =  2.76948572E+01 km
 Maximum X-grid spacing, DXmax =  5.55664201E+01 km
 Minimum Y-grid spacing, DYmin =  2.80117339E+01 km
 Maximum Y-grid spacing, DYmax =  5.55653933E+01 km
 Minimum Z-grid spacing, DZmin =  8.91413705E-01 m
 Maximum Z-grid spacing, DZmax =  1.56107504E+03 m

 Minimum barotropic Courant Number =  9.76308271E-03
 Maximum barotropic Courant Number =  9.47936129E-02
 Maximum Coriolis   Courant Number =  3.78274280E-02

 Horizontal mixing scaled by grid size, GRDMAX =  5.55658756E+01 km

 NLM: GET_STATE - Read state initial conditions,             t =     0 00:00:01
                   (File: shiyf_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  1.00000000E+37 Max = -1.00000000E+37)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min = -2.82544277E+00 Max =  3.01532629E+01)
                - salinity
                   (Min =  7.48295316E+00 Max =  4.41049689E+01)
    GET_2DFLD   - surface u-momentum stress,                 t =   345 00:00:00
                   (Rec=0012, Index=1, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -6.73725163E+02 Max =  6.74438722E+02)
    GET_2DFLD   - surface v-momentum stress,                 t =   345 00:00:00
                   (Rec=0012, Index=1, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.23820027E+02 Max =  3.91016607E+02)
    GET_2DFLD   - solar shortwave radiation flux,            t =   345 00:00:00
                   (Rec=0012, Index=1, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  7.66000516E-07 Max =  6.73535484E-05)
    GET_2DFLD   - surface net heat flux,                     t =   345 00:00:00
                   (Rec=0012, Index=1, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.22467813E-04 Max =  3.12575784E-05)
    GET_2DFLD   - surface net freswater flux, (E-P),         t =   345 00:00:00
                   (Rec=0012, Index=1, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.01846074E-07 Max =  1.15583962E-07)
    GET_2DFLD   - sea surface salinity climatology,          t =   345 00:00:00
                   (Rec=0012, Index=1, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.65167007E+01 Max =  4.11551018E+01)
    GET_NGFLD   - free-surface eastern boundary condition,   t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.83470606E-01 Max =  7.99989567E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.50685381E-01 Max =  7.99989567E+00)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.60444052E-02 Max =  1.37641049E-02)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.79761074E-03 Max =  1.07012369E-03)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.24812288E-03 Max =  9.84361104E-03)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.02250826E+00 Max =  9.09095033E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.66099335E-01 Max =  3.02597453E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.07493064E-01 Max =  6.03430119E-02)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -6.48846883E-02 Max =  6.44035122E-02)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.08044189E+00 Max =  1.25793587E+00)
    GET_NGFLD   - temperature eastern boundary condition,    t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.67032195E+00 Max =  2.78780000E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.33691664E+01 Max =  3.52194348E+01)
    GET_NGFLD   - temperature southern boundary condition,   t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.15385186E+00 Max =  2.97989341E+01)
    GET_NGFLD   - salinity southern boundary condition,      t =   345 00:00:00
                   (Rec=0012, Index=2, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.20250773E+01 Max =  3.59597283E+01)

 Maximum grid stiffness ratios:  rx0 =   1.483417E-01 (Beckmann and Haidvogel)
                                 rx1 =   8.076796E-01 (Haney)


 Initial basin volumes: TotVolume =  3.7465697418E+17 m3
                        MinVolume =  6.9680682010E+08 m3
                        MaxVolume =  4.8199230429E+12 m3
                          Max/Min =  6.9171582479E+03


 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00001440)

    GET_2DFLD   - surface u-momentum stress,                 t =    15 00:00:00
                   (Rec=0001, Index=2, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -5.63892333E+02 Max =  7.27163233E+02)
    GET_2DFLD   - surface v-momentum stress,                 t =    15 00:00:00
                   (Rec=0001, Index=2, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.21693244E+02 Max =  3.89433969E+02)
    GET_2DFLD   - solar shortwave radiation flux,            t =    15 00:00:00
                   (Rec=0001, Index=2, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.56437748E-06 Max =  7.13031920E-05)
    GET_2DFLD   - surface net heat flux,                     t =    15 00:00:00
                   (Rec=0001, Index=2, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.28904424E-04 Max =  3.03673525E-05)
    GET_2DFLD   - surface net freswater flux, (E-P),         t =    15 00:00:00
                   (Rec=0001, Index=2, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.56906815E-07 Max =  1.14237990E-07)
    GET_2DFLD   - sea surface salinity climatology,          t =    15 00:00:00
                   (Rec=0001, Index=2, File: shiyf_frc.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.77527008E+01 Max =  4.05689011E+01)
    GET_NGFLD   - free-surface eastern boundary condition,   t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.50845128E-01 Max =  8.00844646E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -3.30109051E-01 Max =  8.00844646E+00)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -6.11860132E-03 Max =  2.18489387E-02)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.53135821E-03 Max =  1.14542320E-03)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -1.23071690E-03 Max =  2.52131175E-02)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.37649144E-01 Max =  9.63868906E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -2.11111437E-01 Max =  2.42682754E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.15717638E-02 Max =  6.30581806E-02)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -6.39798238E-02 Max =  7.41423366E-02)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min = -9.81503760E-01 Max =  1.40600143E+00)
    GET_NGFLD   - temperature eastern boundary condition,    t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.67032195E+00 Max =  2.74140332E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.35107269E+01 Max =  3.50963225E+01)
    GET_NGFLD   - temperature southern boundary condition,   t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  1.15385186E+00 Max =  2.98525752E+01)
    GET_NGFLD   - salinity southern boundary condition,      t =    15 00:00:00
                   (Rec=0001, Index=1, File: shiyf_bry.nc)
                   (Tmin=         15.0000 Tmax=        345.0000)
                   (Min =  3.05486010E+01 Max =  3.58689762E+01)

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
          C => (i,j,k)       Cu            Cv            Cw         Max Speed

      0     0 00:00:01  1.567195E-06  2.153599E+04  2.153599E+04  3.828818E+17
          (166,001,05)  0.000000E+00  3.603793E-03  0.000000E+00  9.407347E-01
      DEF_HIS   - creating history file: ocean_his.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      DEF_AVG   - creating average file: ocean_avg.nc
      1     0 00:05:01 NaN           NaN           NaN           NaN          
          (149,003,16)  2.419516E+01  3.961077E+00  0.000000E+00  0.000000E+00

 Blowing-up: Saving latest model state into  RESTART file


 WRT_RST - error while writing variable: zeta
           into restart NetCDF file for time record:    1

 Elapsed CPU time (seconds):

 Node   #  2 CPU:       1.219
 Node   #  1 CPU:       1.242
 Node   #  0 CPU:       1.211

 ROMS/TOMS - Output NetCDF summary for Grid 01:
             number of time records written in HISTORY file = 00000001
             number of time records written in RESTART file = 00000001

 Analytical header files used:

     ROMS/Functionals/ana_btflux.h
     ROMS/Functionals/ana_hmixcoef.h
     ROMS/Functionals/ana_nudgcoef.h

 ROMS/TOMS - Output error ............ exit_flag:   3


 MAIN: Abnormal termination: BLOWUP.
 Node   #  7 CPU:       1.214
 Node   #  4 CPU:       1.231
 Node   #  8 CPU:       1.216
 Node   #  3 CPU:       1.212
 Node   #  6 CPU:       1.196
 Node   #  5 CPU:       1.202
i get some messages from roms discussion,the above error may be result from Vertical S-coordinate System which have a zero value,is it? could you help me?
How to do to avoid this kind of error?

User avatar
arango
Site Admin
Posts: 1349
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: error while writing variable: zeta

#2 Unread post by arango »

shiyf wrote:

Code: Select all

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
          C => (i,j,k)       Cu            Cv            Cw         Max Speed

      0     0 00:00:01  1.567195E-06  2.153599E+04  2.153599E+04  3.828818E+17
          (166,001,05)  0.000000E+00  3.603793E-03  0.000000E+00  9.407347E-01
      DEF_HIS   - creating history file: ocean_his.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      DEF_AVG   - creating average file: ocean_avg.nc
      1     0 00:05:01 NaN           NaN           NaN           NaN          
          (149,003,16)  2.419516E+01  3.961077E+00  0.000000E+00  0.000000E+00

 Blowing-up: Saving latest model state into  RESTART file


 WRT_RST - error while writing variable: zeta
           into restart NetCDF file for time record:    1
:idea: Please read all the errors message that ROMS is giving you. The information is all there :!: Many users fail at this. The model is blowing-up and it is very likely that it is a CFL violation (too large time-step) because ROMS blows-up right away and gives you NaNs. If you don't know what a CFL violation is, you can :arrow: Google it and you will be surprised how much information you can find by your own...

The error in the restart NetCDF file is secondary. The free-surface (zeta) has NaNs values because the model is blowing-up and the NetCDF library cannot handle such values.

shiyf
Posts: 51
Joined: Mon Sep 09, 2013 11:09 pm
Location: OUC

Re: error while writing variable: zeta

#3 Unread post by shiyf »

arango wrote:
shiyf wrote:

Code: Select all

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
          C => (i,j,k)       Cu            Cv            Cw         Max Speed

      0     0 00:00:01  1.567195E-06  2.153599E+04  2.153599E+04  3.828818E+17
          (166,001,05)  0.000000E+00  3.603793E-03  0.000000E+00  9.407347E-01
      DEF_HIS   - creating history file: ocean_his.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      DEF_AVG   - creating average file: ocean_avg.nc
      1     0 00:05:01 NaN           NaN           NaN           NaN          
          (149,003,16)  2.419516E+01  3.961077E+00  0.000000E+00  0.000000E+00

 Blowing-up: Saving latest model state into  RESTART file


 WRT_RST - error while writing variable: zeta
           into restart NetCDF file for time record:    1
:idea: Please read all the errors message that ROMS is giving you. The information is all there :!: Many users fail at this. The model is blowing-up and it is very likely that it is a CFL violation (too large time-step) because ROMS blows-up right away and gives you NaNs. If you don't know what a CFL violation is, you can :arrow: Google it and you will be surprised how much information you can find by your own...

The error in the restart NetCDF file is secondary. The free-surface (zeta) has NaNs values because the model is blowing-up and the NetCDF library cannot handle such values.
Thanks a lot!! the error have solved,but when the mode run a short period of time,then stop !! I can not find error message,could you help me!
Attachments
log.in
(143.86 KiB) Downloaded 241 times

Post Reply