Blow up problems for an Internal Tide case

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
jinguangzhen
Posts: 6
Joined: Tue Aug 05, 2014 5:31 pm
Location: Ocean University of China

Blow up problems for an Internal Tide case

#1 Unread post by jinguangzhen »

I am a new ROMS modeler and try to run a simple internal tide case in the South China Sea. I have prepared all the input files (.nc files, .h file and .in file) using the ROMS_tools. Followings are the details of the results.

Code: Select all

 SCS TIDAL MODEL
       720  ntimes   Total number of timesteps for 3D equations.
     36.00  dt       Timestep [sec] for 3D equations
         1  ndtfast  Number of 2D timesteps within each 3D step.
         1  ninfo    Number of timesteps between runtime diagnostics.

 6.000E+00  theta_s  S-coordinate surface control parameter.
 0.000E+00  theta_b  S-coordinate bottom control parameter.
 1.000E+01  Tcline   S-coordinate surface/bottom layer width used in
                                vertical coordinate stretching, meters.
          Grid File:  ROMS_FILES/roms_grd.nc
  Forcing Data File:  ROMS_FILES/roms_frc.nc

 WARNING: Unrecognized keyword: climatology  --> DISREGARDED.

 Initial State File:  ROMS_FILES/roms_ini.nc    Record:  1
       Restart File:  ROMS_FILES/roms_rst.nc    nrst =   720    rec/file:   -1
       History File:  ROMS_FILES/roms_his.nc  Create new: T  nwrt =  72  rec/file =  0
         1  ntsavg      Starting timestep for the accumulation of output
                                time-averaged data.
        72  navg        Number of timesteps between writing of time-averaged
                                data into averages file.
      Averages File:  ROMS_FILES/roms_avg.nc rec/file =  0

 Fields to be saved in history file: (T/F)
      T  write zeta  free-surface.
      T  write UBAR  2D U-momentum component.
      T  write VBAR  2D V-momentum component.
      T  write U     3D U-momentum component.
      T  write V     3D V-momentum component.
      T  write T( 1)  Tracer of index 1.
      T  write T( 2)  Tracer of index 2.

      F  write RHO   Density anomaly.
      F  write Omega Omega vertical velocity.
      T  write W     True vertical velocity.
      F  write Akv   Vertical viscosity.
      T  write Akt   Vertical diffusivity for temperature.
      F  write Aks   Vertical diffusivity for salinity.
      T  write Hbl   Depth of KPP-model boundary layer.
      T  write Hbbl   Depth of bottom planetary boundary layer.

      T  write Bostr Bottom Stress.
      T  write Wstress Wind Stress.
      T  write U-Wstress comp. U-Wind Stress.
      T  write V-Wstress comp. V-Wind Stress.
      T  write Shflx [W/m2] Surface net heat flux
      T  write Swflx [cm/day] Surface freshwater flux (E-P)
      T  write Shflx_rsw [W/m2] Short-wave surface radiation
      T  write Visc3d Horizontal diffusivity.

 WARNING: Unrecognized keyword: gls_history_fields  --> DISREGARDED.


 Fields to be saved in averages file: (T/F)
      T  write zeta  free-surface.
      T  write UBAR  2D U-momentum component.
      T  write VBAR  2D V-momentum component.
      T  write U     3D U-momentum component.
      T  write V     3D V-momentum component.
      T  write T( 1)  Tracer of index 1.
      T  write T( 2)  Tracer of index 2.

      F  write RHO   Density anomaly
      T  write Omega Omega vertical velocity.
      T  write W     True vertical velocity.
      F  write Akv   Vertical viscosity
      T  write Akt   Vertical diffusivity for temperature.
      F  write Aks   Vertical diffusivity for salinity.
      T  write Hbl   Depth of KPP-model boundary layer
      T  write Hbbl   Depth of the bottom planetary boundary layer

      T  write Bostr Bottom Stress.
      T  write Wstr Wind Stress.
      T  write U-Wstress comp. U-Wind Stress.
      T  write V-Wstress comp. V-Wind Stress.
      T  write Shflx [W/m2] Surface net heat flux.
      T  write Swflx [cm/day] Surface freshwater flux (E-P)
      T  write Shflx_rsw [W/m2] Short-wave surface radiation.
      T  write diff3d Horizontal diffusivity

 WARNING: Unrecognized keyword: gls_averages  --> DISREGARDED.

 1025.0000  rho0     Boussinesq approximation mean density, kg/m3.

 WARNING: Unrecognized keyword: lateral_visc  --> DISREGARDED.


 WARNING: Unrecognized keyword: tracer_diff2  --> DISREGARDED.

 0.000E+00  tnu4( 1)  Horizontal biharmonic mixing coefficient [m4/s]
                                for tracer  1.
 0.000E+00  tnu4( 2)  Horizontal biharmonic mixing coefficient [m4/s]
                                for tracer  2.

 WARNING: Unrecognized keyword: vertical_mixing  --> DISREGARDED.

 3.000E-04  rdrg     Linear bottom drag coefficient (m/si).
 0.000E+00  rdrg2    Quadratic bottom drag coefficient.
 0.000E+00  Zob      Bottom roughness for logarithmic law (m).
 1.000E-04  Cdb_min  Minimum bottom drag coefficient.
 1.000E-01  Cdb_max  Maximum bottom drag coefficient.

      1.00  gamma2   Slipperiness parameter: free-slip +1, or no-slip -1.

 WARNING: Unrecognized keyword: sponge  --> DISREGARDED.


 WARNING: Unrecognized keyword: nudg_cof  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagnostics  --> DISREGARDED.


 WARNING: Unrecognized keyword: diag_avg  --> DISREGARDED.


 WARNING: Unrecognized keyword: diag3D_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diag2D_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diag3D_average_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diag2D_average_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagnosticsM  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagM_avg  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagM_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagM_average_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagnostics_bio  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbio_avg  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbioFlux_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbioVSink_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbioGasExc_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbioFlux_average_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbioVSink_average_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: diagbioGasExc_average_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: biology  --> DISREGARDED.


 WARNING: Unrecognized keyword: sediments  --> DISREGARDED.


 WARNING: Unrecognized keyword: sediment_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: bbl_history_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: floats  --> DISREGARDED.


 WARNING: Unrecognized keyword: float_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: stations  --> DISREGARDED.


 WARNING: Unrecognized keyword: station_fields  --> DISREGARDED.


 WARNING: Unrecognized keyword: psource  --> DISREGARDED.


 WARNING: Unrecognized keyword: psource_ncfile  --> DISREGARDED.


 Activated C-preprocessing Options:

          REGIONAL
          SCS_LR
          OPENMP
          OBC_EAST
          OBC_WEST
          OBC_NORTH
          OBC_SOUTH
          CURVGRID
          SPHERICAL
          MASKING
          WET_DRY
          NEW_S_COORD
          SOLVE3D
          UV_COR
          UV_ADV
          SALINITY
          NONLIN_EOS
          SPLIT_EOS
          TS_HADV_RSUP3
          TS_HADV_C4
          TS_DIF4
          LMD_MIXING
          LMD_SKPP
          LMD_BKPP
          LMD_RIMIX
          LMD_CONVEC
          LMD_NONLOCAL
          ANA_DIURNAL_SW
          QCORRECTION
          SFLX_CORR
          ANA_BSFLUX
          ANA_BTFLUX
          OBC_M2CHARACT
          OBC_M3ORLANSKI
          OBC_TORLANSKI
          AVERAGES
          AVERAGES_K
          TS_MIX_GEO
          VAR_RHO_2D
          M2FILTER_POWER
          UV_VADV_SPLINES
          TS_MIX_IMP
          DIF_COEF_3D
          NTRA_T3DMIX
          TS_VADV_AKIMA
          LIMIT_BSTRESS
          ANA_BRY_WKB
          LMD_SKPP2005
          DBLEPREC
          Linux
          QUAD
          QuadZero
          GLOBAL_2D_ARRAY
          GLOBAL_1D_ARRAYXI
          GLOBAL_1D_ARRAYETA
          START_2D_ARRAY
          START_1D_ARRAYXI
          START_1D_ARRAYETA
          PRIVATE_1D_SCRATCH_ARRAY
          PRIVATE_2D_SCRATCH_ARRAY
          PRIVATE_1DXI_SCRATCH_ARRAY
          PRIVATE_1DETA_SCRATCH_ARRAY
          WESTERN_EDGE
          EASTERN_EDGE
          SOUTHERN_EDGE
          NORTHERN_EDGE
          MYID
          MPI_master_only
          ZEROTH_TILE
          SINGLE_TILE_MODE
          LF_AM_STEP
          FIRST_TIME_STEP
          FIRST_2D_STEP
          NOT_LAST_2D_STEP
          PUT_GRID_INTO_RESTART
          PUT_GRID_INTO_HISTORY
          PUT_GRID_INTO_AVERAGES
          NF_FTYPE
          nf_get_att_FTYPE
          nf_put_att_FTYPE
          nf_get_var1_FTYPE
          nf_put_var1_FTYPE
          nf_get_vara_FTYPE
          nf_put_vara_FTYPE
          nf_put_var_FTYPE
          NF_FOUT
          fast_indx_out


 NUMBER OF THREADS:  2 BLOCKING:  1 x  4.

 Spherical grid detected.


 hmin        hmax         grdmin         grdmax         Cu_min      Cu_max
   50.000000  5000.000000 .802237814E+04 .926106764E+04  0.12175196  1.34977021
   volume=1.230119167693202400000E+16   open_cross=1.502076713765440750122E+10

lonmin= 99.0000 lonmax=129.0000latmin=  0.0000 latmax= 30.0286

 Vertical S-coordinate System:

 level   S-coord     Cs-curve          at_hmin  over_slope     at_hmax

    20   0.0000000   0.0000000           0.000       0.000       0.000
    19  -0.0500000  -0.0002259          -0.426      -1.060      -1.626
    18  -0.1000000  -0.0009240          -0.872      -3.297      -5.609
    17  -0.1500000  -0.0021577          -1.340      -6.867     -12.264
    16  -0.2000000  -0.0040388          -1.835     -12.049     -22.150
    15  -0.2500000  -0.0067379          -2.364     -19.268     -36.117
    14  -0.3000000  -0.0104998          -2.937     -29.133     -55.388
    13  -0.3500000  -0.0156655          -3.569     -42.494     -81.664
    12  -0.4000000  -0.0227035          -4.279     -60.517    -117.283
    11  -0.4500000  -0.0322519          -5.094     -84.791    -165.429
    10  -0.5000000  -0.0451767          -6.049    -117.472    -230.422
     9  -0.5500000  -0.0626496          -7.194    -161.478    -318.112
     8  -0.6000000  -0.0862553          -8.594    -220.755    -436.403
     7  -0.6500000  -0.1181340         -10.339    -300.633    -595.978
     6  -0.7000000  -0.1611766         -12.549    -408.308    -811.260
     5  -0.7500000  -0.2192860         -15.387    -553.501   -1101.726
     4  -0.8000000  -0.2977313         -19.072    -749.331   -1493.669
     3  -0.8500000  -0.4036258         -23.901   -1013.510   -2022.584
     2  -0.9000000  -0.5465716         -30.274   -1369.949   -2736.385
     1  -0.9500000  -0.7395307         -38.730   -1850.923   -3699.754
     0  -1.0000000  -1.0000000         -50.000   -2500.000   -5000.000

 Time splitting: ndtfast =  1    nfast =  1


 Maximum grid stiffness ratios:   rx0 =0.2366677417  rx1 =  1.6033734673

      GET_INITIAL -- Processing data for time =   0.000     record =   1
      GET_SMFLUX -- Read surface momentum stresses for time =    345.0    
      GET_SMFLUX -- Read surface momentum stresses for time =    15.00    
      GET_STFLUX -- Read surface flux of tracer  1 for time =    345.0    
      GET_STFLUX -- Read surface flux of tracer  1 for time =    15.00    
      GET_SST    -- Read SST and dQdSST fields     for time =    345.0    
      GET_SST    -- Read SST and dQdSST fields     for time =    15.00    
      GET_STFLUX -- Read surface flux of tracer  2 for time =    345.0    
      GET_STFLUX -- Read surface flux of tracer  2 for time =    15.00    
      GET_SSS    -- Read SSS fields                for time =    345.0    
      GET_SSS    -- Read SSS fields                for time =    15.00    
      GET_SRFLUX -- Read solar shortwave radiation for time =    345.0    
      GET_SRFLUX -- Read solar shortwave radiation for time =    15.00    
      DEF_HIS/AVG - Created new netCDF file 'ROMS_FILES/roms_his.nc'.
      WRT_GRID -- wrote grid data into file 'ROMS_FILES/roms_his.nc'.
      WRT_HIS -- wrote history fields into time record =   1 /   1

 MAIN: started time-steping.

 STEP   time[DAYS] KINETIC_ENRG    POTEN_ENRG    TOTAL_ENRG    NET_VOLUME   trd
       0     0.00000 0.000000000E+00 3.5992966E+01 3.5992966E+01 1.2301192E+16  0
       1     0.00042 8.695225696E-08 4.2114368E+01 4.2114368E+01 1.2081542E+16  1
       2     0.00083 3.454807995E-07 4.2114369E+01 4.2114369E+01 1.2081542E+16  0
       3     0.00125 7.738019801E-07 4.2114365E+01 4.2114366E+01 1.2081542E+16  0
       4     0.00167 1.370372795E-06 4.2114357E+01 4.2114359E+01 1.2081542E+16  1
       5     0.00208 2.134435943E-06 4.2114345E+01 4.2114347E+01 1.2081542E+16  1
       6     0.00250 3.064890743E-06 4.2114329E+01 4.2114332E+01 1.2081542E+16  0
       7     0.00292 4.161216080E-06 4.2114310E+01 4.2114314E+01 1.2081542E+16  1
       8     0.00333 5.422368738E-06 4.2114288E+01 4.2114293E+01 1.2081542E+16  0
       9     0.00375 6.847865684E-06 4.2114263E+01 4.2114270E+01 1.2081542E+16  1
      10     0.00417 8.436592958E-06 4.2114237E+01 4.2114245E+01 1.2081542E+16  0
      11     0.00458 1.018805389E-05 4.2114209E+01 4.2114220E+01 1.2081542E+16  0
      12     0.00500 1.210105515E-05 4.2114182E+01 4.2114194E+01 1.2081543E+16  0
      13     0.00542 1.417518628E-05 4.2114155E+01 4.2114170E+01 1.2081543E+16  0
      14     0.00583 1.640937775E-05 4.2114131E+01 4.2114147E+01 1.2081543E+16  1
      15     0.00625 1.880319455E-05 4.2114108E+01 4.2114127E+01 1.2081543E+16  0
      16     0.00667 2.135558885E-05 4.2114090E+01 4.2114111E+01 1.2081543E+16  0
      17     0.00708 2.406608693E-05 4.2114077E+01 4.2114101E+01 1.2081543E+16  0
      18     0.00750 2.693369028E-05 4.2114070E+01 4.2114097E+01 1.2081544E+16  1
      19     0.00792 2.995794005E-05 4.2114070E+01 4.2114100E+01 1.2081544E+16  1
      20     0.00833 3.313792726E-05 4.2114079E+01 4.2114112E+01 1.2081544E+16  1
      21     0.00875 3.647321323E-05 4.2114099E+01 4.2114135E+01 1.2081544E+16  0
      22     0.00917 3.996301533E-05 4.2114130E+01 4.2114170E+01 1.2081545E+16  0
      23     0.00958 4.360699859E-05 4.2114174E+01 4.2114218E+01 1.2081545E+16  0
      24     0.01000 4.740449958E-05 4.2114233E+01 4.2114281E+01 1.2081545E+16  1
      25     0.01042 5.135530507E-05 4.2114309E+01 4.2114360E+01 1.2081545E+16  1
      26     0.01083 5.545894233E-05 4.2114403E+01 4.2114458E+01 1.2081546E+16  0
      27     0.01125 5.971538507E-05 4.2114516E+01 4.2114576E+01 1.2081546E+16  0
      28     0.01167 6.412435017E-05 4.2114652E+01 4.2114716E+01 1.2081546E+16  0
      29     0.01208 6.868605770E-05 4.2114811E+01 4.2114879E+01 1.2081547E+16  0
      30     0.01250 7.340040635E-05 4.2114995E+01 4.2115069E+01 1.2081547E+16  0
      31     0.01292 7.826796023E-05 4.2115207E+01 4.2115285E+01 1.2081547E+16  0
      32     0.01333 8.328886627E-05 4.2115449E+01 4.2115532E+01 1.2081548E+16  0
      33     0.01375 8.846400609E-05 4.2115722E+01 4.2115811E+01 1.2081548E+16  0
      34     0.01417 9.379388540E-05 4.2116030E+01 4.2116124E+01 1.2081549E+16  1
      35     0.01458 9.927971636E-05 4.2116374E+01 4.2116473E+01 1.2081549E+16  1
      36     0.01500 1.049223998E-04 4.2116757E+01 4.2116862E+01 1.2081549E+16  0
      37     0.01542 1.107235563E-04 4.2117182E+01 4.2117293E+01 1.2081550E+16  1
      38     0.01583 1.166845978E-04 4.2117651E+01 4.2117767E+01 1.2081550E+16  0
      39     0.01625 1.228075781E-04 4.2118167E+01 4.2118290E+01 1.2081551E+16  0
      40     0.01667 1.290944797E-04 4.2118733E+01 4.2118862E+01 1.2081551E+16  0
      41     0.01708 1.355478976E-04 4.2119353E+01 4.2119488E+01 1.2081552E+16  1
      42     0.01750 1.421705088E-04 4.2120029E+01 4.2120171E+01 1.2081552E+16  0
      43     0.01792 1.489655289E-04 4.2120766E+01 4.2120915E+01 1.2081553E+16  1
      44     0.01833 1.559364143E-04 4.2121567E+01 4.2121723E+01 1.2081553E+16  1
      45     0.01875 1.630871986E-04 4.2122437E+01 4.2122600E+01 1.2081554E+16  0
      46     0.01917 1.704222502E-04 4.2123379E+01 4.2123549E+01 1.2081554E+16  1
      47     0.01958 1.779465488E-04 4.2124399E+01 4.2124576E+01 1.2081555E+16  0
      48     0.02000 1.856655626E-04 4.2125501E+01 4.2125686E+01 1.2081555E+16  0
      49     0.02042 1.935855346E-04 4.2126691E+01 4.2126885E+01 1.2081556E+16  0
      50     0.02083 2.017133579E-04 4.2127976E+01 4.2128177E+01 1.2081556E+16  0
      51     0.02125 2.100565502E-04 4.2129361E+01 4.2129571E+01 1.2081557E+16  0
      52     0.02167 2.186236281E-04 4.2130853E+01 4.2131072E+01 1.2081558E+16  0
      53     0.02208 2.274240723E-04 4.2132461E+01 4.2132689E+01 1.2081558E+16  0
      54     0.02250 2.364684042E-04 4.2134193E+01 4.2134429E+01 1.2081559E+16  0
      55     0.02292 2.457683773E-04 4.2136057E+01 4.2136303E+01 1.2081559E+16  0
      56     0.02333 2.553370215E-04 4.2138065E+01 4.2138320E+01 1.2081560E+16  1
      57     0.02375 2.651889265E-04 4.2140227E+01 4.2140492E+01 1.2081561E+16  1
      58     0.02417 2.753402319E-04 4.2142555E+01 4.2142830E+01 1.2081561E+16  0
      59     0.02458 2.858089764E-04 4.2145063E+01 4.2145349E+01 1.2081562E+16  0
      60     0.02500 2.966153176E-04 4.2147766E+01 4.2148063E+01 1.2081563E+16  1
      61     0.02542 3.077817377E-04 4.2150681E+01 4.2150988E+01 1.2081563E+16  0
      62     0.02583 3.193333257E-04 4.2153825E+01 4.2154144E+01 1.2081564E+16  0
      63     0.02625 3.312980904E-04 4.2157219E+01 4.2157550E+01 1.2081565E+16  0
      64     0.02667 3.437073344E-04 4.2160886E+01 4.2161229E+01 1.2081566E+16  0
      65     0.02708 3.565960416E-04 4.2164850E+01 4.2165207E+01 1.2081566E+16  0
      66     0.02750 3.700033552E-04 4.2169141E+01 4.2169511E+01 1.2081567E+16  0
      67     0.02792 3.839730633E-04 4.2173788E+01 4.2174172E+01 1.2081568E+16  0
      68     0.02833 3.985541958E-04 4.2178827E+01 4.2179226E+01 1.2081569E+16  0
      69     0.02875 4.138018623E-04 4.2184298E+01 4.2184711E+01 1.2081569E+16  1
      70     0.02917 4.297769642E-04 4.2190242E+01 4.2190672E+01 1.2081570E+16  0
      71     0.02958 4.465484249E-04 4.2196711E+01 4.2197157E+01 1.2081571E+16  0
      72     0.03000 4.641933217E-04 4.2203757E+01 4.2204221E+01 1.2081572E+16  0
      WRT_HIS -- wrote history fields into time record =   2 /   2
      DEF_HIS/AVG - Created new netCDF file 'ROMS_FILES/roms_avg.nc'.
      WRT_GRID -- wrote grid data into file 'ROMS_FILES/roms_avg.nc'.
      WRT_AVG -- wrote averaged fields into time record =   1 /   1
      73     0.03042 4.827978930E-04 4.2211442E+01 4.2211925E+01 1.2081572E+16  0
      74     0.03083 5.024586419E-04 4.2219836E+01 4.2220338E+01 1.2081573E+16  0
      75     0.03125 5.232835938E-04 4.2229014E+01 4.2229538E+01 1.2081574E+16  0
      76     0.03167 5.453938302E-04 4.2239065E+01 4.2239610E+01 1.2081575E+16  1
      77     0.03208 5.689246833E-04 4.2250083E+01 4.2250652E+01 1.2081576E+16  0
      78     0.03250 5.940275724E-04 4.2262180E+01 4.2262774E+01 1.2081577E+16  1
      79     0.03292 6.208716873E-04 4.2275476E+01 4.2276097E+01 1.2081577E+16  0
      80     0.03333 6.496462995E-04 4.2290110E+01 4.2290760E+01 1.2081578E+16  0
      81     0.03375 6.805630045E-04 4.2306236E+01 4.2306917E+01 1.2081579E+16  1
      82     0.03417 7.138582369E-04 4.2324027E+01 4.2324741E+01 1.2081580E+16  0
      83     0.03458 7.497960085E-04 4.2343679E+01 4.2344428E+01 1.2081581E+16  1
      84     0.03500 7.886712222E-04 4.2365408E+01 4.2366197E+01 1.2081582E+16  0
      85     0.03542 8.308131533E-04 4.2389462E+01 4.2390293E+01 1.2081583E+16  1
      86     0.03583 8.765895649E-04 4.2416114E+01 4.2416991E+01 1.2081583E+16  0
      87     0.03625 9.264108019E-04 4.2445675E+01 4.2446601E+01 1.2081584E+16  0
      88     0.03667 9.807356368E-04 4.2478489E+01 4.2479470E+01 1.2081585E+16  0
      89     0.03708 1.040076696E-03 4.2514946E+01 4.2515986E+01 1.2081586E+16  0
      90     0.03750 1.105007172E-03 4.2555480E+01 4.2556585E+01 1.2081587E+16  0
      91     0.03792 1.176169028E-03 4.2600580E+01 4.2601757E+01 1.2081588E+16  0
      92     0.03833 1.254280645E-03 4.2650793E+01 4.2652047E+01 1.2081589E+16  0
      93     0.03875 1.340151205E-03 4.2706731E+01 4.2708071E+01 1.2081590E+16  0
      94     0.03917 1.434688293E-03 4.2769081E+01 4.2770516E+01 1.2081590E+16  0
      95     0.03958 1.538917691E-03 4.2838613E+01 4.2840152E+01 1.2081591E+16  0
      96     0.04000 1.653997351E-03 4.2916191E+01 4.2917845E+01 1.2081592E+16  0
      97     0.04042 1.781243653E-03 4.3002783E+01 4.3004564E+01 1.2081593E+16  0
      98     0.04083 1.922174439E-03 4.3099480E+01 4.3101402E+01 1.2081594E+16  1
      99     0.04125 2.078407958E-03 4.3207493E+01 4.3209572E+01 1.2081595E+16  0
     100     0.04167 2.252674177E-03 4.3328298E+01 4.3330550E+01 1.2081596E+16  1
     101     0.04208 2.441442250E-03 4.3462318E+01 4.3464760E+01 1.2081597E+16  0
     102     0.04250 2.703502594E-03 4.3600877E+01 4.3603581E+01 1.2081598E+16  1
     103     0.04292 3.049854473E-03-4.5130322E+03-4.5130291E+03 1.2081599E+16  0
                                         
                                         
 ======================================= 
 =                                     = 
 =   STEP2D:   ABNORMAL JOB END        = 
 =                 BLOW UP             = 
 =                                     = 
 ======================================= 
                                         
  VMAX (M/S) =********
  IMAX JMAX  =   335   302
  IINT IEXT  =   104     1
I think this is because of the wrong set of #define in the .h file, which is given here:

Code: Select all

! $Id: cppdefs.h 1560 2014-06-19 09:15:09Z gcambon $
!
!======================================================================
! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France
! The two other branches from UCLA (Shchepetkin et al) 
! and Rutgers University (Arango et al) are under MIT/X style license.
! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license.
! 
! ROMS_AGRIF website : http://www.romsagrif.org
!======================================================================
!
/*
   This is "cppdefs.h": MODEL CONFIGURATION FILE
   ==== == ============ ===== ============= ====
*/
#undef  BASIN           /* Basin Example */
#undef  CANYON_A        /* Canyon_A Example */
#undef  CANYON_B        /* Canyon_B Example */
#undef  EQUATOR         /* Equator Example  */
#undef  GRAV_ADJ        /* Graviational Adjustment Example */
#undef  INNERSHELF      /* Inner Shelf Example */
#undef  RIVER           /* River run-off Example */
#undef  OVERFLOW        /* Graviational/Overflow Example */
#undef  SEAMOUNT        /* Seamount Example */
#undef  SHELFRONT       /* Shelf Front Example */
#undef  SOLITON         /* Equatorial Rossby Wave Example */
#undef  UPWELLING       /* Upwelling Example */
#undef  VORTEX          /* Baroclinic Vortex Example */
#undef  INTERNAL        /* Internal Tide Example */
#undef  JET             /* Baroclinic Jet Example */
#undef  RIP             /* Rip Current Test Case */
#undef  SHOREFACE       /* Shoreface Test Case on a Planar Beach */
#undef  THACKER         /* Thacker wetting-drying Example */
#define REGIONAL        /* REGIONAL Applications */


#if defined REGIONAL
/*
!====================================================================
!               REGIONAL (realistic) Configurations
!==================================================================== 
!
!----------------------
! BASIC OPTIONS
!----------------------
!
*/
                      /* Configuration Name */
# define SCS_LR
                      /* Parallelization */
# define  OPENMP
# undef  MPI
                      /* Nesting */
# undef  AGRIF
# undef  AGRIF_2WAY
                      /* OA Coupling via OASIS (MPI) */
# undef  OA_COUPLING
                      /* Open Boundary Conditions */
# undef  TIDES
# define OBC_EAST
# define OBC_WEST
# define OBC_NORTH
# define OBC_SOUTH
                      /* Applications */
# undef  BIOLOGY
# undef  FLOATS
# undef  STATIONS
# undef  PASSIVE_TRACER
# undef  SEDIMENT
# undef  BBL
/*!
!-------------------------------------------------
! PRE-SELECTED OPTIONS
!
! ADVANCED OPTIONS ARE IN SET_GLOBAL_DEFINITIONS.H
!-------------------------------------------------
*/
                      /* Parallelization */
# ifdef MPI
#  undef  PARALLEL_FILES
# endif
# undef  AUTOTILING
# undef  ETALON_CHECK
                      /* Grid configuration */
# define CURVGRID
# define SPHERICAL
# define MASKING
# define  WET_DRY
# define  NEW_S_COORD
                      /* Model dynamics */
# define SOLVE3D
# define UV_COR
# define UV_ADV
                      /* Equation of State */
# define SALINITY
# define NONLIN_EOS
# define SPLIT_EOS
                      /* Lateral Explicit Momentum Mixing */
# undef  UV_VIS2
# ifdef UV_VIS2
#  define UV_MIX_S
#  define UV_VIS_SMAGO
# endif
                      /* Lateral Tracer Advection (default UP3) */
# define TS_HADV_RSUP3
# undef  TS_HADV_UP5
# undef  TS_HADV_C4
# undef  TS_HADV_WENO5
                      /* Lateral Explicit Tracer Mixing */
# undef TS_MIX_S
# ifdef TS_HADV_C4
#  define  TS_DIF2
#  undef   TS_DIF4
#  define  TS_DIF_SMAGO
#  define  TS_MIX_ISO
# endif
                      /* Sponge layers for UV and TS */
# undef SPONGE
                      /* Vertical Mixing */
# undef  BODYFORCE
# undef  BVF_MIXING
# define LMD_MIXING
# undef  GLS_MIXING
# ifdef LMD_MIXING
#  define LMD_SKPP
#  define LMD_BKPP
#  define LMD_RIMIX
#  define LMD_CONVEC
#  undef  LMD_DDMIX
#  define LMD_NONLOCAL
# endif
# ifdef GLS_MIXING
#  define GLS_KKL
#  undef  GLS_KOMEGA
#  undef  GLS_KEPSILON
#  undef  GLS_GEN
#  undef  KANTHA_CLAYSON
#  undef  CRAIG_BANNER
#  undef  CANUTO_A
#  undef  ZOS_HSIG
# endif
                      /* Surface Forcing */
# undef  BULK_FLUX
# ifdef BULK_FLUX
#  define BULK_FAIRALL
#  define BULK_LW
#  define BULK_EP
#  define BULK_SMFLUX
#  undef  SST_SKIN
#  undef  ANA_DIURNAL_SW
#  undef  ONLINE
# else
#  define QCORRECTION
#  define SFLX_CORR
#  define ANA_DIURNAL_SW
# endif
                      /* Lateral Forcing */
# undef CLIMATOLOGY
# ifdef CLIMATOLOGY
#  define ZCLIMATOLOGY
#  define M2CLIMATOLOGY
#  define M3CLIMATOLOGY
#  define TCLIMATOLOGY

#  define ZNUDGING
#  define M2NUDGING
#  define M3NUDGING
#  define TNUDGING
#  undef  ROBUST_DIAG
# endif

# undef  FRC_BRY
# ifdef FRC_BRY
#  define Z_FRC_BRY
#  define M2_FRC_BRY
#  define M3_FRC_BRY
#  define T_FRC_BRY
# endif
                      /* Bottom Forcing */
# define ANA_BSFLUX
# define ANA_BTFLUX
                      /* Point Sources - Rivers */
# undef PSOURCE
# undef PSOURCE_NCFILE
# ifdef PSOURCE_NCFILE                    
#   define PSOURCE_NCFILE_TS
# endif
                      /* Open Boundary Conditions */
# ifdef TIDES
#  define SSH_TIDES
#  define UV_TIDES
#  define TIDERAMP
#  define OBC_M2FLATHER
# else
#  undef  OBC_M2SPECIFIED
#  undef  OBC_M2FLATHER
#  define OBC_M2CHARACT
#  undef  OBC_M2ORLANSKI
#  ifdef  OBC_M2ORLANSKI
#   define OBC_VOLCONS
#  endif
# endif
# define OBC_M3ORLANSKI
# define OBC_TORLANSKI
# undef  OBC_M3SPECIFIED
# undef  OBC_TSPECIFIED
                      /* Input/Output & Diagnostics */
# define AVERAGES
# define AVERAGES_K
# undef  DIAGNOSTICS_TS
# undef  DIAGNOSTICS_UV
# ifdef DIAGNOSTICS_TS
#  undef DIAGNOSTICS_TS_ADV
#  undef DIAGNOSTICS_TS_MLD
# endif
/*
!           Applications:
!---------------------------------
! Biology, floats, Stations, 
! Passive tracer, Sediments, BBL
!---------------------------------
!
   Quasi-monotone lateral advection scheme (WENO5)
   for passive/biology/sediment tracers 
*/
# if defined PASSIVE_TRACER || defined BIOLOGY || defined SEDIMENT
#  undef BIO_HADV_WENO5
# endif
                      /*   Choice of Biology models   */
# ifdef BIOLOGY
#  undef  PISCES
#  undef BIO_NChlPZD
#  undef  BIO_N2ChlPZD2
#  define BIO_BioEBUS
                      /*   Biology options    */
#  ifdef PISCES
#   define DIURNAL_INPUT_SRFLX
#   define key_trc_pisces
#   define key_passivetrc
#  endif
#  ifdef BIO_NChlPZD
#   define  OXYGEN
#  endif
#  ifdef BIO_BioEBUS
#   define NITROUS_OXIDE
#  endif
                      /*   Biology diagnostics    */
#  undef DIAGNOSTICS_BIO
#  if defined DIAGNOSTICS_BIO && defined PISCES
#   define key_trc_diaadd
#   define key_trc_dia3d
#  endif
# endif
                      /*   Lagrangian floats model    */
# ifdef FLOATS
#  undef  FLOATS_GLOBAL_ATTRIBUTES
#  undef  IBM
#  undef  RANDOM_WALK
#  ifdef RANDOM_WALK
#   define DIEL_MIGRATION
#   define RANDOM_VERTICAL
#   define RANDOM_HORIZONTAL
#  endif
# endif
                      /*   Stations recording    */
# ifdef STATIONS
#  define ALL_SIGMA
# endif
                      /*   Sediment dynamics model     */
# ifdef SEDIMENT
#  define ANA_SEDIMENT
#  undef  BED_ARMOR
#  undef  ANA_SPFLUX
#  undef  ANA_BPFLUX
# endif
                      /*   Bottom Boundary Layer model     */
# ifdef BBL
#  define ANA_WWAVE
# endif
/*
!
!==========================================================
!              IDEALIZED CONFIGURATIONS
!==========================================================
!
*/
#elif defined BASIN
/*
!                       Basin Example
!                       ===== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define UV_MIX_S
# define SOLVE3D
# define TS_DIF2
# define TS_MIX_S
# define BODYFORCE
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_BTFLUX

#elif defined CANYON_A
/*
!                       First Canyon Example
!                       ===== ====== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define SOLVE3D
# define EW_PERIODIC
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_BTFLUX

#elif defined CANYON_B
/*
!                       Second Canyon Example
!                       ====== ====== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define SOLVE3D
# define EW_PERIODIC
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_BTFLUX
# define ANA_VMIX

#elif defined EQUATOR
/*
!                       Equator Example
!                       ======= =======
! Boccaletti, G., R.C. Pacanowski, G.H. Philander and A.V. Fedorov, 2004,
! The Thermal Structure of the Upper Ocean, J.Phys.Oceanogr., 34, 888-902.
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define UV_MIX_S
# define SOLVE3D
# define SALINITY
# define TS_DIF2
# define TS_MIX_S
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SRFLUX
# define ANA_SSFLUX
# define ANA_BTFLUX
# define ANA_BSFLUX
# define QCORRECTION
# define ANA_SST
# define LMD_SKPP /* problem with MPI in Xi direction */
# define LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC

#elif defined GRAV_ADJ
/*
!                       Gravitational Adjustment Example
!                       ============= ========== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define SOLVE3D
# define UV_ADV
# define UV_VIS2
# define UV_MIX_S
# undef  TS_HADV_WENO5
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_BTFLUX

#elif defined INNERSHELF
/*
!                       Inner Shelf Example
!                       ===== ===== =======
*/
# define ETALON_CHECK
# undef  OPENMP
# undef  MPI
# define INNERSHELF_EKMAN
# define INNERSHELF_APG
# define SOLVE3D
# define UV_COR
# define ANA_GRID
# define ANA_INITIAL
# define AVERAGES
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_STFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_SMFLUX
# define NS_PERIODIC
# define OBC_WEST
# define OBC_M2ORLANSKI
# define OBC_VOLCONS
# define OBC_TORLANSKI
# define OBC_M3ORLANSKI
# define SPONGE
# ifndef INNERSHELF_EKMAN
#  define UV_ADV
#  define SALINITY
#  define NONLIN_EOS
#  define SPLIT_EOS
#  define LMD_MIXING
#  ifdef LMD_MIXING
#   define LMD_SKPP
#   define LMD_BKPP
#   define LMD_RIMIX
#   define LMD_CONVEC
#  else
#   define GLS_MIXING
#   define GLS_KKL
#  endif
# endif

#elif defined INTERNAL
/*
!                       Internal Tide Example
!                       ======== ==== =======
!
! Di Lorenzo, E, W.R. Young and S.L. Smith, 2006, Numerical and anlytical estimates of M2
! tidal conversion at steep oceanic ridges, J. Phys. Oceanogr., 36, 1072-1084.  
*/
# define SCS_LR
# undef  ETALON_CHECK
# undef  OPENMP
# undef  MPI
# define SOLVE3D
# define UV_COR
# define UV_ADV
# define BODYTIDE
# define ANA_GRID
# undef  INTERNALSHELF
# define ANA_INITIAL
# define ANA_BTFLUX
# define ANA_SMFLUX
# define ANA_SRFLUX
# define ANA_STFLUX
# define ANA_VMIX

# define EW_PERIODIC
# define NS_PERIODIC

# undef  UV_VIS2
# undef  UV_MIX_GEO
# undef  TS_DIF2
# undef  TS_MIX_GEO
# undef  SPONGE
# undef  ANA_SSH
# undef  ANA_M2CLIMA
# undef  ANA_M3CLIMA
# undef  ANA_TCLIMA
# undef  ZCLIMATOLOGY
# undef  M2CLIMATOLOGY
# undef  M3CLIMATOLOGY
# undef  TCLIMATOLOGY
# undef  ZNUDGING
# undef  M2NUDGING
# undef  M3NUDGING
# undef  TNUDGING
# undef  OBC_EAST
# undef  OBC_WEST
# undef  OBC_M2CHARACT
# undef  OBC_M2FLATHER
# undef  OBC_TORLANSKI
# undef  OBC_M3ORLANSKI

#elif defined RIVER
/*
!                       River run-off test problem
!                       ==========================
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define SOLVE3D
# define UV_ADV
# define UV_COR
# define M2FILTER_FLAT
# define NONLIN_EOS
# define SPLIT_EOS
# define SALINITY
# define ANA_GRID
# define MASKING
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_BTFLUX
# define ANA_BSFLUX
# define LMD_MIXING
# define LMD_SKPP
# define LMD_BKPP
# define LMD_RIMIX
# define LMD_CONVEC
# define PSOURCE
# define ANA_PSOURCE
# define NS_PERIODIC
# define FLOATS
# ifdef FLOATS
#   define RANDOM_WALK
#   ifdef RANDOM_WALK
#      define DIEL_MIGRATION
#      define RANDOM_VERTICAL
#      define RANDOM_HORIZONTAL
#   endif
# endif

#elif defined SEAMOUNT
/*
!                       Seamount Example
!                       ======== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define SOLVE3D
# define SALINITY
# define NONLIN_EOS
# define SPLIT_EOS
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_BTFLUX
# define ANA_BSFLUX

# elif defined SHELFRONT
/*
!                       Shelf Front Example
!                       ===== ===== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define SOLVE3D
# define SALINITY
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_BTFLUX
# define ANA_BSFLUX
# define EW_PERIODIC

#elif defined SOLITON
/*
!                       Equatorial Rossby Wave Example
!                       ========== ====== ==== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_COR
# define UV_ADV
# define ANA_GRID
# define ANA_INITIAL
# define AVERAGES
# define EW_PERIODIC
# define ANA_SMFLUX

# elif defined OVERFLOW
/*
!                       Gravitational/Overflow Example
!                       ====================== =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define UV_MIX_GEO
# define TS_DIF2
# define TS_MIX_GEO
# define SOLVE3D
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_BTFLUX

#elif defined UPWELLING
/*
!                       Upwelling Example
!                       ========= =======
*/
# define ETALON_CHECK
# undef OPENMP
# undef MPI
# define SOLVE3D
# define UV_COR
# define UV_ADV
# define ANA_GRID
# define ANA_INITIAL
# define AVERAGES
# define SALINITY
# define NONLIN_EOS
# define SPLIT_EOS
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_STFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_SMFLUX
# define LMD_MIXING
# define LMD_SKPP
# define LMD_BKPP
# define LMD_RIMIX
# define LMD_CONVEC
# define EW_PERIODIC

#elif defined VORTEX
/*
!                       Baroclinic Vortex Example (TEST AGRIF)
!                       ========== ====== ======= ===== ======
*/
# undef  ETALON_CHECK
# undef  OPENMP
# undef  MPI
# define AGRIF
# define AGRIF_2WAY
# define SOLVE3D
# define UV_COR
# define UV_ADV
# define ANA_STFLUX
# define ANA_SMFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_VMIX
# define SPONGE
# define TS_MIX_S
# define ZCLIMATOLOGY
# define M2CLIMATOLOGY
# define M3CLIMATOLOGY
# define TCLIMATOLOGY
# define ZNUDGING
# define M2NUDGING
# define M3NUDGING
# define TNUDGING
# define OBC_EAST
# define OBC_WEST
# define OBC_NORTH
# define OBC_SOUTH
# define OBC_M2FLATHER
# define OBC_TORLANSKI
# define OBC_M3ORLANSKI

#elif defined JET
/*
!                       Baroclinic JET Example
!                       ========== === =======
*/
# define ANA_JET
# undef  ETALON_CHECK
# undef  MPI
# undef  AGRIF
# undef  AGRIF_2WAY
# define SOLVE3D
# define UV_COR
# define UV_ADV
# define UV_VIS2
# define UV_MIX_S
# ifdef ANA_JET
#  define ANA_GRID
#  define ANA_INITIAL
# endif
# define ANA_STFLUX
# define ANA_SMFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_VMIX
# define EW_PERIODIC
# define CLIMATOLOGY
# ifdef CLIMATOLOGY
#  define ZCLIMATOLOGY
#  define M2CLIMATOLOGY
#  define M3CLIMATOLOGY
#  define TCLIMATOLOGY
#  define ZNUDGING
#  define M2NUDGING
#  define M3NUDGING
#  define TNUDGING
#  define ROBUST_DIAG
#  define ZONAL_NUDGING
#  ifdef ANA_JET
#   define ANA_SSH
#   define ANA_M2CLIMA
#   define ANA_M3CLIMA
#   define ANA_TCLIMA
#  endif
# endif
# define LMD_MIXING 
# ifdef  LMD_MIXING
#  undef  ANA_VMIX
#  define ANA_SRFLUX
#  undef  LMD_KPP
#  define LMD_RIMIX
#  define LMD_CONVEC
# endif 

#elif defined SHOREFACE
/*
!                       PLANAR BEACH Example
!                       ====== ===== =======
!
!   Uchiyama, Y., McWilliams, J.C. and Shchepetkin, A.F. (2010): 
!      Wave-current interaction in an oceanic circulation model with a 
!      vortex force formalism: Application to the surf zone.
!      Ocean Modelling Vol. 34:1-2, pp.16-35.
*/
# undef  ETALON_CHECK
# undef  OPENMP
# undef  MPI
# define SOLVE3D
# define UV_ADV
# undef  MASKING
# define NEW_S_COORD
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_SST
# define ANA_BTFLUX
# define NS_PERIODIC
# define OBC_WEST
# ifdef OBC_WEST
#  define OBC_M2CHARACT
#  define OBC_M3ORLANSKI
#  define OBC_TORLANSKI
#  define SPONGE
#  define TS_MIX_S
# endif
# define WET_DRY
# define MRL_WCI
# ifdef MRL_WCI
#  undef  WKB_WWAVE
#  undef  WAVE_ROLLER
#  undef  WAVE_STREAMING
#  undef  MRL_CEW
#  define WAVE_RAMP
# endif
# define LMD_MIXING
# define LMD_SKPP
# define LMD_BKPP
# undef  BBL
# undef  SEDIMENT

#elif defined RIP
/*
!                       Rip Current Example
!                       === ======= =======
!
!   Weir, B., Uchiyama, Y.. (2010): 
!      A vortex force analysis of the interaction of rip 
!      currents and surface gravity wave
!      JGR Vol. 116
*/
# undef  BISCA
# undef  DIAGNOSTICS_UV
# undef  OPENMP
# undef  MPI
# define SOLVE3D
# define UV_ADV
# define NEW_S_COORD
# ifndef BISCA
#  define ANA_GRID
# endif
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_SRFLUX
# define ANA_SST
# define ANA_BTFLUX
# ifndef BISCA
#  define NS_PERIODIC
# endif
# define OBC_WEST
# ifdef BISCA
#  define OBC_NORTH
#  define OBC_SOUTH
# endif
# ifdef OBC_WEST
#  define OBC_M2CHARACT
#  define OBC_M3ORLANSKI
#  define OBC_TORLANSKI
#  define SPONGE
# endif
# define WET_DRY
# define MRL_WCI
# ifdef MRL_WCI
#  define WKB_WWAVE
#  define WKB_OBC_WEST
#  define WAVE_ROLLER
#  define WAVE_STREAMING
#  define WAVE_FRICTION
#  define MRL_CEW
#  undef  WKB_KZ_FILTER
#  undef  WKB_TIME_FILTER
#  undef  WAVE_RAMP
# endif
# define LMD_MIXING
# define LMD_SKPP
# define LMD_BKPP
# define UV_VIS2
# ifdef UV_VIS2
#  define UV_MIX_S
#  define UV_VIS_SMAGO
# endif
# define BBL
# undef  SEDIMENT
# ifdef SEDIMENT
#  define ANA_SEDIMENT
#  undef  BED_ARMOR
#  undef  ANA_SPFLUX
#  undef  ANA_BPFLUX
# endif

#elif defined THACKER
/*
!                       Thacker Example
!                       ======= =======
!
! Thacker, W., (1981), Some exact solutions to the nonlinear shallow-water wave equations. 
! J. Fluid Mech., 107, 499–508.
*/
# undef  ETALON_CHECK
# undef  OPENMP
# undef  MPI
# define SOLVE3D
# define UV_COR
# define UV_ADV
# define WET_DRY
# define NEW_S_COORD
# define ANA_GRID
# define ANA_INITIAL
# define ANA_BTFLUX
# define ANA_SMFLUX
# define ANA_SRFLUX
# define ANA_STFLUX


#endif /* END OF CONFIGURATION CHOICE */

#include "set_global_definitions.h"
Does anyone encounter similar problem or have sulotion? Many thanks! :)

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blow up problems for an Internal Tide case

#2 Unread post by kate »

Many of us have blowing up problems from time to time. However, most of us on this forum are not familiar with that dialect of ROMS cppdefs. Is this the CROCO code you are using? I suggest you ask on the CROCO forums.

jinguangzhen
Posts: 6
Joined: Tue Aug 05, 2014 5:31 pm
Location: Ocean University of China

Re: Blow up problems for an Internal Tide case

#3 Unread post by jinguangzhen »

kate wrote:Many of us have blowing up problems from time to time. However, most of us on this forum are not familiar with that dialect of ROMS cppdefs. Is this the CROCO code you are using? I suggest you ask on the CROCO forums.
Hi, Kate. Thank you so much for your kind reply. The answer is yes, I am using the CROCO packages. Actually I think it might be an easy way to prepare ROMS inputs and many of ROMS users are using it. Any way, do you have any other suggestions to prepare ROMS input files for a new user?

Best regards,
Guangzhen

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blow up problems for an Internal Tide case

#4 Unread post by kate »

I was referring to your ROMS code, not the preprocessing code. I at least don't know what cppdefs to use with CROCO - some are familiar, some not. I just can't help you because it's not the code this forum was built to support.

geo_ytwin
Posts: 6
Joined: Tue Apr 13, 2010 4:16 pm
Location: South China Sea Institute of Oceanology, CAS

Re: Blow up problems for an Internal Tide case

#5 Unread post by geo_ytwin »

Set your ndtfast to 6 or 10, or even larger. As the 1 you set in the post, it means that the barotropic time step was 36 seconds. It is obviously too large.

jinguangzhen
Posts: 6
Joined: Tue Aug 05, 2014 5:31 pm
Location: Ocean University of China

Re: Blow up problems for an Internal Tide case

#6 Unread post by jinguangzhen »

kate wrote:I was referring to your ROMS code, not the preprocessing code. I at least don't know what cppdefs to use with CROCO - some are familiar, some not. I just can't help you because it's not the code this forum was built to support.
Thank you anyway,Kate. BTW, I have found out the gridpak toolbox for ROMS in the forum but haven't found any specific manuals online. I think it is used to prepare input files for ROMS? Would you please give me a method the get the manual? Thank you very much! :D

Best,
Guangzhen

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Blow up problems for an Internal Tide case

#7 Unread post by kate »

Well, gridpak is ancient and the manual is even older. There's a link to the manual from the github README. It's for making grids only. For all the boundary condition remapping, I use pyroms (Python 2), trying to transition to pycnal (Python 3). But really, there's no reason for you to give up your Matlab tools. As for CROCO, why not use it too? It's just that you have to look elsewhere for help with some of its features. The ROMS versions diverged maybe 15 years ago with rather different priorities for development since then.

jinguangzhen
Posts: 6
Joined: Tue Aug 05, 2014 5:31 pm
Location: Ocean University of China

Re: Blow up problems for an Internal Tide case

#8 Unread post by jinguangzhen »

geo_ytwin wrote:Set your ndtfast to 6 or 10, or even larger. As the 1 you set in the post, it means that the barotropic time step was 36 seconds. It is obviously too large.
Thank you so much geo_ytwin! This is the problem of BLOW UP in my case! Now the model can be run successfully!

Best,
Guangzhen

jinguangzhen
Posts: 6
Joined: Tue Aug 05, 2014 5:31 pm
Location: Ocean University of China

Re: Blow up problems for an Internal Tide case

#9 Unread post by jinguangzhen »

geo_ytwin wrote:Set your ndtfast to 6 or 10, or even larger. As the 1 you set in the post, it means that the barotropic time step was 36 seconds. It is obviously too large.
And wish you a happy new Year!

jinguangzhen
Posts: 6
Joined: Tue Aug 05, 2014 5:31 pm
Location: Ocean University of China

Re: Blow up problems for an Internal Tide case

#10 Unread post by jinguangzhen »

kate wrote:Well, gridpak is ancient and the manual is even older. There's a link to the manual from the github README. It's for making grids only. For all the boundary condition remapping, I use pyroms (Python 2), trying to transition to pycnal (Python 3). But really, there's no reason for you to give up your Matlab tools. As for CROCO, why not use it too? It's just that you have to look elsewhere for help with some of its features. The ROMS versions diverged maybe 15 years ago with rather different priorities for development since then.
Thank you, Kate! I am still trying to figure out how the Roms_tools work. I just think it maybe helpful to combine multiple tools to make the preparation easier and more convenient for me. I will learn more about your PYROMS. That must be interesting! Wish you a Happy New Year!
Best,
Guangzhen

Post Reply