Model blowing up

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
gloeg002

Model blowing up

#1 Unread post by gloeg002 »

I recently changed my grid to be more resolved and now my model is blowing up. The Courant number is fine, so I know it's not my time step. Here is the log file. Any ideas?

Code: Select all

 Resolution, Grid 01: 0302x0145x030,  Parallel Nodes:  16,  Tiling: 004x004


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

       4000  ntimes          Number of timesteps for 3-D equations.
    100.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.
          9  nSTA            Number of timesteps between the writing of data
                               the stations file.
          T  ldefout         Switch to create a new output NetCDF file(s).
          9  nHIS            Number of timesteps between the writing fields
                               into history file.
          1  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
          9  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
          1  ntsDIA          Starting timestep for the accumulation of output
                               time-averaged diagnostics data.
          9  nDIA            Number of timesteps between the writing of
                               time-averaged data into diagnostics 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.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  Vtransform      S-coordinate transformation equation.
          4  Vstretching     S-coordinate stretching function.
 3.0000E+00  theta_s         S-coordinate surface control parameter.
 0.0000E+00  theta_b         S-coordinate bottom  control parameter.
     25.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)
 0.0000E+00  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 0.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.
          F  VolCons(2)      NLM southern edge boundary volume conservation.
          F  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.
   1027.000  R0              Background density (kg/m3) used in linear Equation
                               of State.
 1.7000E-04  Tcoef           Thermal expansion coefficient (1/Celsius).
 0.0000E+00  Scoef           Saline contraction coefficient (1/PSU).
      1.000  gamma2          Slipperiness variable: free-slip (1.0) or 
                                                    no-slip (-1.0).
          T  Hout(idFsur)    Write out free-surface.
          T  Hout(idUvel)    Write out 3D U-momentum component.
          T  Hout(idVvel)    Write out 3D V-momentum component.





 Output/Input Files:

             Output Restart File:  /home/austinj/gloegel/roms-models/model-001/Output/model-001_rst.nc
             Output History File:  /home/austinj/gloegel/roms-models/model-001/Output/model-001_his.nc
            Output Averages File:  /home/austinj/gloegel/roms-models/model-001/Output/model-001_avg.nc
         Output Diagnostics File:  /home/austinj/gloegel/roms-models/model-001/Output/model-001_dia.nc
            Output Stations File:  /home/austinj/gloegel/roms-models/model-001/Output/model-001_sta.nc
                 Input Grid File:  /home/austinj/gloegel/roms-grids/LS_grid_2km.nc
          Station positions File:  /home/austinj/gloegel/roms-models/model-001/run-files/stations.in

 Tile partition information for Grid 01:  0302x0145x0030  tiling: 004x004

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       75        1       36    81000
        1       76      151        1       36    82080
        2      152      227        1       36    82080
        3      228      302        1       36    81000
        4        1       75       37       73    83250
        5       76      151       37       73    84360
        6      152      227       37       73    84360
        7      228      302       37       73    83250
        8        1       75       74      110    83250
        9       76      151       74      110    84360
       10      152      227       74      110    84360
       11      228      302       74      110    83250
       12        1       75      111      145    78750
       13       76      151      111      145    79800
       14      152      227      111      145    79800
       15      228      302      111      145    78750

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

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    75.50     0.50    36.50  RHO-points
        1    75.50   151.50     0.50    36.50  RHO-points
        2   151.50   227.50     0.50    36.50  RHO-points
        3   227.50   302.50     0.50    36.50  RHO-points
        4     0.50    75.50    36.50    73.50  RHO-points
        5    75.50   151.50    36.50    73.50  RHO-points
        6   151.50   227.50    36.50    73.50  RHO-points
        7   227.50   302.50    36.50    73.50  RHO-points
        8     0.50    75.50    73.50   110.50  RHO-points
        9    75.50   151.50    73.50   110.50  RHO-points
       10   151.50   227.50    73.50   110.50  RHO-points
       11   227.50   302.50    73.50   110.50  RHO-points
       12     0.50    75.50   110.50   145.50  RHO-points
       13    75.50   151.50   110.50   145.50  RHO-points
       14   151.50   227.50   110.50   145.50  RHO-points
       15   227.50   302.50   110.50   145.50  RHO-points

        0     1.00    75.50     0.50    36.50    U-points
        1    75.50   151.50     0.50    36.50    U-points
        2   151.50   227.50     0.50    36.50    U-points
        3   227.50   302.00     0.50    36.50    U-points
        4     1.00    75.50    36.50    73.50    U-points
        5    75.50   151.50    36.50    73.50    U-points
        6   151.50   227.50    36.50    73.50    U-points
        7   227.50   302.00    36.50    73.50    U-points
        8     1.00    75.50    73.50   110.50    U-points
        9    75.50   151.50    73.50   110.50    U-points
       10   151.50   227.50    73.50   110.50    U-points
       11   227.50   302.00    73.50   110.50    U-points
       12     1.00    75.50   110.50   145.50    U-points
       13    75.50   151.50   110.50   145.50    U-points
       14   151.50   227.50   110.50   145.50    U-points
       15   227.50   302.00   110.50   145.50    U-points

        0     0.50    75.50     1.00    36.50    V-points
        1    75.50   151.50     1.00    36.50    V-points
        2   151.50   227.50     1.00    36.50    V-points
        3   227.50   302.50     1.00    36.50    V-points
        4     0.50    75.50    36.50    73.50    V-points
        5    75.50   151.50    36.50    73.50    V-points
        6   151.50   227.50    36.50    73.50    V-points
        7   227.50   302.50    36.50    73.50    V-points
        8     0.50    75.50    73.50   110.50    V-points
        9    75.50   151.50    73.50   110.50    V-points
       10   151.50   227.50    73.50   110.50    V-points
       11   227.50   302.50    73.50   110.50    V-points
       12     0.50    75.50   110.50   145.00    V-points
       13    75.50   151.50   110.50   145.00    V-points
       14   151.50   227.50   110.50   145.00    V-points
       15   227.50   302.50   110.50   145.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     264
               HaloSizeJ(1) =     147
                TileSide(1) =      82
                TileSize(1) =    3526



 Stations Parameters, Grid: 01
 =============================

          3  Nstation        Number of stations to write out into stations file.
          T  Sout(idTvar)    Write out tracer 01: temp
          T  Sout(idUsms)    Write out surface U-momentum stress.
          T  Sout(idVsms)    Write out surface V-momentum stress.
 
             Flag and positions for station 0001:  0    80.0000   40.0000
             Flag and positions for station 0002:  0    80.0000   30.0000
             Flag and positions for station 0003:  0    80.0000   20.0000

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

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

 zeta                     1     Periodic    Closed      Periodic    Closed

 ubar                     1     Periodic    Closed      Periodic    Closed

 vbar                     1     Periodic    Closed      Periodic    Closed

 u                        1     Periodic    Closed      Periodic    Closed

 v                        1     Periodic    Closed      Periodic    Closed

 temp                     1     Periodic    Closed      Periodic    Closed

 salt                     1     Periodic    Closed      Periodic    Closed

 tke                      1     Periodic    Closed      Periodic    Closed

 Activated C-preprocessing Options:

 UPWELLING           Wind-Driven Upwelling/Downwelling over a Periodic Channel
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_INITIAL         Analytical initial conditions.
 ANA_SMFLUX          Analytical kinematic surface momentum flux.
 ANA_SSFLUX          Analytical kinematic surface salinity flux.
 ANA_STFLUX          Analytical kinematic surface temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 AVERAGES            Writing out time-averaged nonlinear model fields.
 DIAGNOSTICS_TS      Computing and writing tracer diagnostic terms.
 DIAGNOSTICS_UV      Computing and writing momentum diagnostic terms.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 MASKING             Land/Sea masking.
 MIX_S_TS            Mixing of tracers along constant S-surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 NONLINEAR           Nonlinear Model.
 !NONLIN_EOS         Linear Equation of State for seawater.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLVE3D             Solving 3D Primitive Equations.
 SPLINES             Conservative parabolic spline reconstruction.
 STATIONS            Writing out station data.
 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.

 Process Information:

 Node #  0 (pid=    1871) is active.
 Node #  3 (pid=    1874) is active.
 Node #  5 (pid=    1876) is active.
 Node #  1 (pid=    1872) is active.
 Node #  2 (pid=    1873) is active.
 Node #  4 (pid=    1875) is active.
 Node #  6 (pid=    1877) is active.
 Node #  7 (pid=    1878) is active.
 Node #  8 (pid=   16684) is active.
 Node #  9 (pid=   16685) is active.
 Node # 10 (pid=   16686) is active.
 Node # 11 (pid=   16687) is active.
 Node # 12 (pid=   16688) is active.
 Node # 13 (pid=   16689) is active.
 Node # 14 (pid=   16690) is active.
 Node # 15 (pid=   16691) is active.

 INITIAL: Configuring and initializing forward nonlinear model ...

 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f
 reading f

 Vertical S-coordinate System: 

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

    30   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    29  -0.0333333  -0.0005519       -0.032      -0.424      -0.840      -0.993
    28  -0.0666667  -0.0022130       -0.064      -0.861      -1.879      -2.407
    27  -0.1000000  -0.0050000       -0.096      -1.313      -3.120      -4.246
    26  -0.1333333  -0.0089408       -0.129      -1.778      -4.568      -6.522
    25  -0.1666667  -0.0140748       -0.161      -2.259      -6.230      -9.250
    24  -0.2000000  -0.0204535       -0.193      -2.756      -8.115     -12.448
    23  -0.2333333  -0.0281406       -0.225      -3.268     -10.234     -16.142
    22  -0.2666667  -0.0372130       -0.258      -3.798     -12.602     -20.360
    21  -0.3000000  -0.0477616       -0.290      -4.347     -15.234     -25.137
    20  -0.3333333  -0.0598920       -0.323      -4.915     -18.150     -30.512
    19  -0.3666667  -0.0737256       -0.355      -5.505     -21.371     -36.532
    18  -0.4000000  -0.0894007       -0.388      -6.118     -24.922     -43.249
    17  -0.4333333  -0.1070744       -0.421      -6.755     -28.832     -50.723
    16  -0.4666667  -0.1269234       -0.454      -7.420     -33.131     -59.020
    15  -0.5000000  -0.1491465       -0.487      -8.114     -37.856     -68.215
    14  -0.5333333  -0.1739660       -0.520      -8.841     -43.047     -78.392
    13  -0.5666667  -0.2016303       -0.553      -9.604     -48.747     -89.647
    12  -0.6000000  -0.2324164       -0.586     -10.405     -55.007    -102.083
    11  -0.6333333  -0.2666323       -0.619     -11.250     -61.882    -115.817
    10  -0.6666667  -0.3046205       -0.653     -12.141     -69.433    -130.979
     9  -0.7000000  -0.3467612       -0.686     -13.085     -77.729    -147.712
     8  -0.7333333  -0.3934761       -0.720     -14.085     -86.844    -166.176
     7  -0.7666667  -0.4452328       -0.754     -15.149     -96.863    -186.549
     6  -0.8000000  -0.5025493       -0.789     -16.282    -107.879    -209.026
     5  -0.8333333  -0.5659992       -0.823     -17.492    -119.995    -233.824
     4  -0.8666667  -0.6362176       -0.858     -18.786    -133.323    -261.183
     3  -0.9000000  -0.7139071       -0.893     -20.174    -147.991    -291.371
     2  -0.9333333  -0.7998455       -0.928     -21.665    -164.138    -324.680
     1  -0.9666667  -0.8948927       -0.964     -23.269    -181.917    -361.436
     0  -1.0000000  -1.0000000       -1.000     -25.000    -201.500    -402.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.00000000E+00 km
 Maximum X-grid spacing, DXmax =  2.00000000E+00 km
 Minimum Y-grid spacing, DYmin =  2.00000000E+00 km
 Maximum Y-grid spacing, DYmax =  2.00000000E+00 km
 Minimum Z-grid spacing, DZmin =  3.20725078E-02 m
 Maximum Z-grid spacing, DZmax =  4.05638507E+01 m

 Minimum barotropic Courant Number =  7.38241153E-03
 Maximum barotropic Courant Number =  1.48016891E-01
 Maximum Coriolis   Courant Number =  1.10000000E-02


 Maximum grid stiffness ratios:  rx0 =   9.682540E-01 (Beckmann and Haidvogel)
                                 rx1 =   2.974274E+01 (Haney)


 Initial basin volumes: TotVolume =  1.1818404000E+13 m3
                        MinVolume =  1.2829003124E+05 m3
                        MaxVolume =  1.6225540264E+08 m3
                          Max/Min =  1.2647545649E+03


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


   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:00  0.000000E+00  8.702628E+02  8.702628E+02  1.228302E+13
          (000,000,00)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
      DEF_HIS   - creating history file: /home/austinj/gloegel/roms-models/model-001/Output/model-001_his.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      DEF_AVG   - creating average file: /home/austinj/gloegel/roms-models/model-001/Output/model-001_avg.nc
      DEF_DIAGS - creating diagnostics file: /home/austinj/gloegel/roms-models/model-001/Output/model-001_dia.nc
      DEF_STATION - creating stations file: /home/austinj/gloegel/roms-models/model-001/Output/model-001_sta.nc
      1     0 00:01:40  2.294946E-04  8.702606E+02  8.702608E+02  1.228258E+13
          (113,086,01)  2.467573E-02  2.539287E-02  0.000000E+00  1.135981E+00
      2     0 00:03:20  7.700370E-04  8.702543E+02  8.702551E+02  1.228156E+13
          (111,096,01)  3.220059E-02  4.033535E-02  4.664799E+01  2.059685E+00
      3     0 00:05:00  1.530783E-03  8.702329E+02  8.702345E+02  1.228025E+13
          (111,096,01)  5.648131E-02  6.783622E-02  8.453411E+01  2.789785E+00
      4     0 00:06:40  2.470958E-03  8.701939E+02  8.701964E+02  1.227881E+13
          (111,096,01)  8.368497E-02  9.658389E-02  1.121651E+02  3.431938E+00
      5     0 00:08:20  3.557556E-03  8.701424E+02  8.701460E+02  1.227733E+13
          (111,096,01)  1.083242E-01  1.223930E-01  1.316511E+02  3.963897E+00
      6     0 00:10:00  4.755493E-03  8.700833E+02  8.700880E+02  1.227585E+13
          (111,096,01)  1.285067E-01  1.426415E-01  1.431047E+02  4.765708E+00
      7     0 00:11:40  6.035806E-03  8.700189E+02  8.700249E+02  1.227440E+13
          (111,096,01)  1.460593E-01  1.607027E-01  1.469665E+02  6.306291E+00
      8     0 00:13:20  7.936571E-03  8.699506E+02  8.699586E+02  1.227296E+13
          (111,096,01)  1.602873E-01  1.783562E-01  1.490090E+02  2.462400E+01

 Blowing-up: Saving latest model state into  RESTART file

 Node   #  6 CPU:      15.849
 Node   #  5 CPU:      16.202
 Node   #  1 CPU:      17.056
 Node   #  7 CPU:      16.050
 Node   #  2 CPU:      17.246
 Node   #  3 CPU:      17.195
 Node   #  4 CPU:      16.399
      WRT_RST   - wrote re-start fields (Index=1,1) into time record = 0000001

 Elapsed CPU time (seconds):

 Node   #  0 CPU:      15.668
 Node   # 14 CPU:      16.701
 Node   # 15 CPU:      16.095
 Node   #  9 CPU:      16.065
 Node   #  8 CPU:      15.173
 Node   # 11 CPU:      16.406
 Node   # 12 CPU:      17.201
 Node   # 13 CPU:      17.160
 Node   # 10 CPU:      15.873

 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
             number of time records written in STATION file = 00000001

 Analytical header files used:

     ROMS/Functionals/ana_btflux.h
     ROMS/Functionals/ana_initial.h
     ROMS/Functionals/ana_nudgcoef.h
     ROMS/Functionals/ana_smflux.h
     ROMS/Functionals/ana_stflux.h

 ROMS/TOMS: DONE... Saturday - June 29, 2013 -  8:24:44 PM

miguel.solano
Posts: 12
Joined: Mon Jan 30, 2012 8:02 pm
Location: University of Texas at Dallas

Re: Model blowing up

#2 Unread post by miguel.solano »

gloeg002:
I have roughly 1 year of experience with ROMS but I'm pretty sure your problem is with your vertical grid. As you can see the model is blowing up at your first level (111 , 96, 01) , but more importantly your maximum grid stiffness ratios are off the roof! This is specially important if you are doing an upwelling case. There is extensive literature on bathymetry smoothing indicating r-factor should be between 0.20 and 0.30 (rx0), although I find rx0<0.25 ideal. You should also watch out for your rx1 factor which the ROMS documentation warns users to be careful with high values. What kind of smoothing did you use?

Post Reply