Hi
I have changed 'theta S', 'theta b' and 'hc' for my domain few time and time step but every time I run model face to following error (BLOWUP:
South China Sea
1080 ntimes Total number of timesteps for 3D equations.
2400.00 dt Timestep [sec] for 3D equations
45 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
Bulk Data File: ROMS_FILES/roms_blk.nc
Climatology File: ROMS_FILES/roms_clm.nc
Boundary File: ROMS_FILES/roms_bry.nc
Initial State File: ROMS_FILES/roms_ini.nc Record: 1
Restart File: ROMS_FILES/roms_rst.nc nrst = 1080 rec/file: -1
History File: ROMS_FILES/roms_his.nc Create new: T nwrt = 180 rec/file = 0
1 ntsavg Starting timestep for the accumulation of output
time-averaged data.
180 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.
F write UBAR 2D U-momentum component.
F write VBAR 2D V-momentum component.
F write U 3D U-momentum component.
F write V 3D V-momentum component.
F write T(1) Tracer of index 1.
F write T(2) Tracer of index 2.
F write RHO Density anomaly.
F write Omega Omega vertical velocity.
F write W True vertical velocity.
F write Akv Vertical viscosity.
F write Akt Vertical diffusivity for temperature.
F write Aks Vertical diffusivity for salinity.
F write Hbl Depth of KPP-model boundary layer.
F write Bostr Bottom Stress.
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 Bostr Bottom Stress.
1025.0000 rho0 Boussinesq approximation mean density, kg/m3.
0.000E+00 visc2 Horizontal Laplacian mixing coefficient [m2/s]
for momentum.
0.000E+00 tnu2(1) Horizontal Laplacian mixing coefficient (m2/s)
for tracer 1.
0.000E+00 tnu2(2) Horizontal Laplacian mixing coefficient (m2/s)
for tracer 2.
0.000E+00 rdrg Linear bottom drag coefficient (m/si).
0.000E+00 rdrg2 Quadratic bottom drag coefficient.
1.000E-02 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.
1.00E+05 x_sponge Thickness of sponge and/or nudging layer (m)
800.00 v_sponge Viscosity in sponge layer (m2/s)
1.157E-05 tauT_in Nudging coefficients [sec^-1]
3.215E-08 tauT_out Nudging coefficients [sec^-1]
1.157E-06 tauM_in Nudging coefficients [sec^-1]
3.215E-08 tauM_out Nudging coefficients [sec^-1]
WARNING: Unrecognized keyword: diagnostics --> DISREGARDED.
WARNING: Unrecognized keyword: diag_avg --> DISREGARDED.
WARNING: Unrecognized keyword: diagnostics_bio --> DISREGARDED.
WARNING: Unrecognized keyword: diagbio_avg --> 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.
Activated C-preprocessing Options:
REGIONAL
SCS
OPENMP
AGRIF_2WAY
OBC_EAST
OBC_WEST
OBC_SOUTH
SOLVE3D
UV_COR
UV_ADV
CURVGRID
SPHERICAL
MASKING
AVERAGES
AVERAGES_K
SALINITY
NONLIN_EOS
SPLIT_EOS
QCORRECTION
SFLX_CORR
DIURNAL_SRFLUX
BULK_FLUX
BULK_EP
BULK_SMFLUX
BULK_WVEC
BULK_WSTR
SPONGE
CLIMATOLOGY
ZCLIMATOLOGY
M2CLIMATOLOGY
M3CLIMATOLOGY
TCLIMATOLOGY
ZNUDGING
M2NUDGING
M3NUDGING
TNUDGING
FRC_BRY
Z_FRC_BRY
M2_FRC_BRY
M3_FRC_BRY
T_FRC_BRY
ANA_BSFLUX
ANA_BTFLUX
UV_VIS2
MIX_GP_UV
TS_DIF2
MIX_GP_TS
CLIMAT_TS_MIXH
LMD_MIXING
LMD_SKPP
LMD_BKPP
LMD_RIMIX
LMD_CONVEC
OBC_M2FLATHER
OBC_VOLCONS
OBC_M2ORLANSKI
OBC_M3ORLANSKI
OBC_TORLANSKI
M2FILTER_COSINE
Linux 2.6.23.17-88.fc7 i686
NUMBER OF THREADS: 2 BLOCKING: 2 x 4.
-11
Spherical grid detected.
-10
hmin hmax grdmin grdmax Cu_min Cu_max
75.000000 4321.321078 .181203344E+05 .185221205E+05 0.11045571 0.85701886
volume=9.914351461535635000000E+14 open_cross=1.412581681829876899719E+09
-9
Vertical S-coordinate System:
level S-coord Cs-curve at_hmin over_slope at_hmax
35 0.0000000 0.0000000 0.000 0.000 0.000
34 -0.0285714 -0.0008540 -0.341 -2.154 -3.968
33 -0.0571429 -0.0017332 -0.684 -4.364 -8.044
32 -0.0857143 -0.0026635 -1.030 -6.685 -12.340
31 -0.1142857 -0.0036722 -1.382 -9.178 -16.975
30 -0.1428571 -0.0047891 -1.740 -11.908 -22.076
29 -0.1714286 -0.0060471 -2.107 -14.946 -27.785
28 -0.2000000 -0.0074832 -2.486 -18.374 -34.263
27 -0.2285714 -0.0091398 -2.880 -22.285 -41.690
26 -0.2571429 -0.0110656 -3.291 -26.785 -50.279
25 -0.2857143 -0.0133175 -3.723 -31.998 -60.273
24 -0.3142857 -0.0159616 -4.180 -38.069 -71.959
23 -0.3428571 -0.0190760 -4.669 -45.170 -85.671
22 -0.3714286 -0.0227524 -5.193 -53.500 -101.807
21 -0.4000000 -0.0270990 -5.761 -63.297 -120.833
20 -0.4285714 -0.0322440 -6.382 -74.841 -143.300
19 -0.4571429 -0.0383389 -7.063 -88.463 -169.863
18 -0.4857143 -0.0455632 -7.819 -104.557 -201.295
17 -0.5142857 -0.0541298 -8.661 -123.588 -238.514
16 -0.5428571 -0.0642911 -9.607 -146.108 -282.608
15 -0.5714286 -0.0763464 -10.677 -172.772 -334.868
14 -0.6000000 -0.0906508 -11.892 -204.358 -396.825
13 -0.6285714 -0.1076258 -13.281 -241.788 -470.295
12 -0.6571429 -0.1277714 -14.877 -286.156 -557.435
11 -0.6857143 -0.1516811 -16.716 -338.760 -660.803
10 -0.7142857 -0.1800593 -18.847 -401.142 -783.436
9 -0.7428571 -0.2137420 -21.322 -475.130 -928.939
8 -0.7714286 -0.2537216 -24.206 -562.898 -1101.589
7 -0.8000000 -0.3011757 -27.576 -667.021 -1306.465
6 -0.8285714 -0.3575023 -31.523 -790.558 -1549.593
5 -0.8571429 -0.4243610 -36.155 -937.141 -1838.128
4 -0.8857143 -0.5037212 -41.599 -1111.080 -2180.561
3 -0.9142857 -0.5979209 -48.008 -1317.490 -2586.972
2 -0.9428571 -0.7097353 -55.561 -1562.443 -3069.325
1 -0.9714286 -0.8424583 -64.474 -1853.148 -3641.823
0 -1.0000000 -1.0000000 -75.000 -2198.161 -4321.321
Time splitting: ndtfast = 45 nfast = 67
-8
-7
Maximum grid stiffness ratios: rx0 =0.2301486227 rx1 = 2.7446044817
GET_INITIAL -- Processing data for time = 0.000 record = 1
-6
-5
-4
-3
-2
GET_TCLIMA -- Read climatology of tracer 1 for time = 345.0
GET_TCLIMA -- Read climatology of tracer 1 for time = 15.00
GET_TCLIMA -- Read climatology of tracer 2 for time = 345.0
GET_TCLIMA -- Read climatology of tracer 2 for time = 15.00
GET_UCLIMA -- Read momentum climatology for time = 345.0
GET_UCLIMA -- Read momentum climatology for time = 15.00
GET_SSH - Read SSH climatology for time = 345.0
GET_SSH - Read SSH climatology for time = 15.00
GET_BRY -- Read all boundary data for time = 345.0
GET_BRY -- Read all boundary data for time = 15.00
GET_SMFLUX -- Read surface momentum stresses for time = 345.0
GET_SMFLUX -- Read surface momentum stresses for time = 15.00
GET_BULK -- Read fields for bulk formula for time = 345.0
GET_BULK -- Read fields for bulk formula for time = 15.00
-1
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
0
MAIN: started time-steping.
STEP time[DAYS] KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME trd
0 0.00000 0.000000000E+00 7.9516822E+00 7.9516822E+00 9.9143515E+14 1
1 0.02778 NaN NaN NaN NaN 1
MAIN - number of records written into history file(s): 1
number of records written into restart file(s): 0
number of records written into averages file(s): 0
MAIN: Abnormal termination: BLOWUP
I have confused completely and tired. Everybody that tell me cause of problem and can help me for resolve problem I thanks so much.
Best Regards
F.Daryabor
RUN
Re: RUN
The warnings you get about stations, floats, diagnostics, etc. could be due to confusion about ROMS versions. The Rutgers repository ROMS has STATIONS, FLOATS, SEDIMENT, etc, but you seem to be using the AGRIF ROMS.
What is your horizontal grid resolution? Does that ROMS report anything about Courant numbers? Is there any way you could start simpler and add some things back in later, for instance the climatology nudging?
What is your horizontal grid resolution? Does that ROMS report anything about Courant numbers? Is there any way you could start simpler and add some things back in later, for instance the climatology nudging?
Re: RUN
Hi
Thank So Much Kate
I use Agrif ROMS and in this ROMS version I have STATION, FLOATS and SEDIMENT options, but I don't selection them. The model resolution is:
Resolution: 1/6 deg
Create the grid file...
LLm = 95
MMm = 84
Fill the grid file...
Compute the metrics...
Min dx=18.1122 km - Max dx=18.5221 km
Min dy=18.1285 km - Max dy=18.5221 km
Fill the grid file...
Add topography...
ROMS resolution : 18.4 km
Topography data resolution : 3.69 km
Topography resolution halved 3 times
New topography resolution : 29.5 km
Theta S =6 , Theta b =0 and hc=10
hmin=75 and hmax_coast=500
----------------------------------------------------------------------------
Kate you sent me 'ifdefs.F' script for fixed '#if' error, but I don't know how employ it for removing '#if' error. please you guide me how I employ it.
thanks so much for any things.
Best Regards
F.Daryabor
Thank So Much Kate
I use Agrif ROMS and in this ROMS version I have STATION, FLOATS and SEDIMENT options, but I don't selection them. The model resolution is:
Resolution: 1/6 deg
Create the grid file...
LLm = 95
MMm = 84
Fill the grid file...
Compute the metrics...
Min dx=18.1122 km - Max dx=18.5221 km
Min dy=18.1285 km - Max dy=18.5221 km
Fill the grid file...
Add topography...
ROMS resolution : 18.4 km
Topography data resolution : 3.69 km
Topography resolution halved 3 times
New topography resolution : 29.5 km
Theta S =6 , Theta b =0 and hc=10
hmin=75 and hmax_coast=500
----------------------------------------------------------------------------
Kate you sent me 'ifdefs.F' script for fixed '#if' error, but I don't know how employ it for removing '#if' error. please you guide me how I employ it.
thanks so much for any things.
Best Regards
F.Daryabor
Re: RUN
For your blow up error, your timestep seems long compared to your grid spacing. I would try a shorter timestep, maybe half what you have. If that works, explore options between the two values (maybe 1800 s).
It will give output something like:
allowing you to see if all the #ifs have matching #endifs.
The thing should be renamed simply 'ifdefs'. To use it:Kate you sent me 'ifdefs.F' script for fixed '#if' error, but I don't know how employ it for removing '#if' error. please you guide me how I employ it.
Code: Select all
ifdefs main.FCode: Select all
cygnus.arsc.edu 195% ifdefs ana_biology_goanpz.h
# ifdef DEPAVG
# ifdef IRON_LIMIT
# endif
# else
# ifdef IRON_LIMIT
# endif
# endif /* DEPAVG */
#ifdef GAK1D
#endif