Blowing-up: Saving latest model state into RESTART file

[Scarlett]

Hi modelers:
-
I'm working with roms-3.7 and when I run my model,it blow up at seven time steps. If I change the passage of time, I only get a blow up delayed.

-
I will add my out.txt, my ocean.in and my .h file.
-
If somebody can help me I will be so gratefull for ever
-
Thanks in advance.
Mar.Mo.

[kate]

I've written the generic response here. Perhaps you should tell us what you are trying to achieve with a periodic channel and wetting and drying - but look to see how it's going bad first.

[lalvarez]

Hello Scarlett

You should also take a look at this

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

It is greater than 0.4, perhaps you need to smooth your grid

Can you share a map or figure of your area?

Saluditos

[Scarlett]

Hi!
-
Many thanks to both for your help.
-
I had smoothing the grid ... step by step to the extent that I smoothing at 2000m deep, because I only have interesting in the first 200 m of the surface, but I still have the same problem ...
-
My Beckmann value still been 9.943090e-1 although the grid is so smoothing



------------------------------------------------------------------------------------------------------------------------------------------------------------------
...
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt

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 Dout(M2rate) Write out 2D momentum acceleration.
T Dout(M2pgrd) Write out 2D momentum pressure gradient.
T Dout(M2fcor) Write out 2D momentum Coriolis force.
T Dout(M2hadv) Write out 2D momentum horizontal advection.
T Dout(M2xadv) Write out 2D momentum horizontal X-advection.
T Dout(M2yadv) Write out 2D momentum horizontal Y-advection.
T Dout(M2hvis) Write out 2D momentum horizontal viscosity.
T Dout(M2xvis) Write out 2D momentum horizontal X-viscosity.
T Dout(M2yvis) Write out 2D momentum horizontal Y-viscosity.
T Dout(M2sstr) Write out 2D momentum surface stress.
T Dout(M2bstr) Write out 2D momentum bottom stress.

T Dout(M3rate) Write out 3D momentum acceleration.
T Dout(M3pgrd) Write out 3D momentum pressure gradient.
T Dout(M3fcor) Write out 3D momentum Coriolis force.
T Dout(M3hadv) Write out 3D momentum horizontal advection.
T Dout(M3xadv) Write out 3D momentum horizontal X-advection.
T Dout(M3yadv) Write out 3D momentum horizontal Y-advection.
T Dout(M3vadv) Write out 3D momentum vertical advection.
T Dout(M3hvis) Write out 3D momentum horizontal viscosity.
T Dout(M3xvis) Write out 3D momentum horizontal X-viscosity.
T Dout(M3yvis) Write out 3D momentum horizontal Y-viscosity.
T Dout(M3vvis) Write out 3D momentum vertical viscosity.

T Dout(iTrate) Write out rate of change of tracer 01: temp
T Dout(iTrate) Write out rate of change of tracer 02: salt
T Dout(iThadv) Write out horizontal advection, tracer 01: temp
T Dout(iThadv) Write out horizontal advection, tracer 02: salt
T Dout(iTxadv) Write out horizontal X-advection, tracer 01: temp
T Dout(iTxadv) Write out horizontal X-advection, tracer 02: salt
T Dout(iTyadv) Write out horizontal Y-advection, tracer 01: temp
T Dout(iTyadv) Write out horizontal Y-advection, tracer 02: salt
T Dout(iTvadv) Write out vertical advection, tracer 01: temp
T Dout(iTvadv) Write out vertical advection, tracer 02: salt
T Dout(iThdif) Write out horizontal diffusion, tracer 01: temp
T Dout(iThdif) Write out horizontal diffusion, tracer 02: salt
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 01: temp
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 02: salt
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 01: temp
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 02: salt
T Dout(iTvdif) Write out vertical diffusion, tracer 01: temp
T Dout(iTvdif) Write out vertical diffusion, tracer 02: salt

Output/Input Files:

Output Restart File: /home/ahumada/Roms_tools/Run/SCRATCH/ocean_rst.nc
Output History File: /home/ahumada/Roms_tools/Run/SCRATCH/ocean_his.nc
Output Diagnostics File: /home/ahumada/Roms_tools/Run/SCRATCH/ocean_dia.nc
Input Grid File: /home/ahumada/roms-3.7/ROMS/External/tehua4srp5_roms_grd.nc

Tile partition information for Grid 01: 0088x0088x0032 tiling: 002x012

tile Istr Iend Jstr Jend Npts

0 1 44 1 4 5632
1 45 88 1 4 5632
2 1 44 5 12 11264
3 45 88 5 12 11264
4 1 44 13 20 11264
5 45 88 13 20 11264
6 1 44 21 28 11264
7 45 88 21 28 11264
8 1 44 29 36 11264
9 45 88 29 36 11264
10 1 44 37 44 11264
11 45 88 37 44 11264
12 1 44 45 52 11264
13 45 88 45 52 11264
14 1 44 53 60 11264
15 45 88 53 60 11264
16 1 44 61 68 11264
17 45 88 61 68 11264
18 1 44 69 76 11264
19 45 88 69 76 11264
20 1 44 77 84 11264
21 45 88 77 84 11264
22 1 44 85 88 5632
23 45 88 85 88 5632

Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 -1.50 91.50 0.50 90.50 RHO-points
1 -2.50 90.50 0.50 90.50 RHO-points
2 -1.50 91.50 -0.50 90.50 RHO-points
3 -2.50 90.50 -0.50 90.50 RHO-points
4 -1.50 91.50 -0.50 90.50 RHO-points
5 -2.50 90.50 -0.50 90.50 RHO-points
6 -1.50 91.50 -0.50 90.50 RHO-points
7 -2.50 90.50 -0.50 90.50 RHO-points
8 -1.50 91.50 -0.50 90.50 RHO-points
9 -2.50 90.50 -0.50 90.50 RHO-points
10 -1.50 91.50 -0.50 90.50 RHO-points
11 -2.50 90.50 -0.50 90.50 RHO-points
12 -1.50 91.50 -0.50 90.50 RHO-points
13 -2.50 90.50 -0.50 90.50 RHO-points
14 -1.50 91.50 -0.50 90.50 RHO-points
15 -2.50 90.50 -0.50 90.50 RHO-points
16 -1.50 91.50 -0.50 90.50 RHO-points
17 -2.50 90.50 -0.50 90.50 RHO-points
18 -1.50 91.50 -0.50 90.50 RHO-points
19 -2.50 90.50 -0.50 90.50 RHO-points
20 -1.50 91.50 -0.50 90.50 RHO-points
21 -2.50 90.50 -0.50 90.50 RHO-points
22 -1.50 91.50 -0.50 89.50 RHO-points
23 -2.50 90.50 -0.50 89.50 RHO-points

0 -2.00 91.50 0.50 90.50 U-points
1 -2.50 90.00 0.50 90.50 U-points
2 -2.00 91.50 -0.50 90.50 U-points
3 -2.50 90.00 -0.50 90.50 U-points
4 -2.00 91.50 -0.50 90.50 U-points
5 -2.50 90.00 -0.50 90.50 U-points
6 -2.00 91.50 -0.50 90.50 U-points
7 -2.50 90.00 -0.50 90.50 U-points
8 -2.00 91.50 -0.50 90.50 U-points
9 -2.50 90.00 -0.50 90.50 U-points
10 -2.00 91.50 -0.50 90.50 U-points
11 -2.50 90.00 -0.50 90.50 U-points
12 -2.00 91.50 -0.50 90.50 U-points
13 -2.50 90.00 -0.50 90.50 U-points
14 -2.00 91.50 -0.50 90.50 U-points
15 -2.50 90.00 -0.50 90.50 U-points
16 -2.00 91.50 -0.50 90.50 U-points
17 -2.50 90.00 -0.50 90.50 U-points
18 -2.00 91.50 -0.50 90.50 U-points
19 -2.50 90.00 -0.50 90.50 U-points
20 -2.00 91.50 -0.50 90.50 U-points
21 -2.50 90.00 -0.50 90.50 U-points
22 -2.00 91.50 -0.50 89.50 U-points
23 -2.50 90.00 -0.50 89.50 U-points

0 -1.50 91.50 0.00 90.50 V-points
1 -2.50 90.50 0.00 90.50 V-points
2 -1.50 91.50 -0.50 90.50 V-points
3 -2.50 90.50 -0.50 90.50 V-points
4 -1.50 91.50 -0.50 90.50 V-points
5 -2.50 90.50 -0.50 90.50 V-points
6 -1.50 91.50 -0.50 90.50 V-points
7 -2.50 90.50 -0.50 90.50 V-points
8 -1.50 91.50 -0.50 90.50 V-points
9 -2.50 90.50 -0.50 90.50 V-points
10 -1.50 91.50 -0.50 90.50 V-points
11 -2.50 90.50 -0.50 90.50 V-points
12 -1.50 91.50 -0.50 90.50 V-points
13 -2.50 90.50 -0.50 90.50 V-points
14 -1.50 91.50 -0.50 90.50 V-points
15 -2.50 90.50 -0.50 90.50 V-points
16 -1.50 91.50 -0.50 90.50 V-points
17 -2.50 90.50 -0.50 90.50 V-points
18 -1.50 91.50 -0.50 90.50 V-points
19 -2.50 90.50 -0.50 90.50 V-points
20 -1.50 91.50 -0.50 90.50 V-points
21 -2.50 90.50 -0.50 90.50 V-points
22 -1.50 91.50 -0.50 89.00 V-points
23 -2.50 90.50 -0.50 89.00 V-points

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:

tehuantepec
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_INITIAL Analytical initial conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_SMFLUX Analytical kinematic surface momentum flux.
ANA_SRFLUX Analytical kinematic shortwave radiation flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
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.
FSOBC_REDUCED Using free-surface data in reduced physics conditions
GLS_MIXING Generic Length-Scale turbulence closure.
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction.
MASKING Land/Sea masking.
MIX_S_TS Mixing of tracers along constant S-surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
OUT_DOUBLE Double precision output fields in NetCDF files.
_OPENMP OpenMP parallel shared-memory directives.
PERFECT_RESTART Processing perfect restart variables.
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.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion.
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity.
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.
WET_DRY Wetting and drying activated.

Process Information:

Thread # 11 (pid= 5394) is active.
Thread # 18 (pid= 5394) is active.
Thread # 7 (pid= 5394) is active.
Thread # 0 (pid= 5394) is active.
Thread # 20 (pid= 5394) is active.
Thread # 15 (pid= 5394) is active.
Thread # 5 (pid= 5394) is active.
Thread # 17 (pid= 5394) is active.
Thread # 9 (pid= 5394) is active.
Thread # 22 (pid= 5394) is active.
Thread # 13 (pid= 5394) is active.
Thread # 1 (pid= 5394) is active.
Thread # 3 (pid= 5394) is active.
Thread # 23 (pid= 5394) is active.
Thread # 16 (pid= 5394) is active.
Thread # 6 (pid= 5394) is active.
Thread # 4 (pid= 5394) is active.
Thread # 12 (pid= 5394) is active.
Thread # 14 (pid= 5394) is active.
Thread # 8 (pid= 5394) is active.
Thread # 10 (pid= 5394) is active.
Thread # 21 (pid= 5394) is active.
Thread # 19 (pid= 5394) is active.
Thread # 2 (pid= 5394) is active.

INITIAL: Configuring and initializing forward nonlinear model ...
*******

Vertical S-coordinate System, Grid 01:

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

32 0.0000000 0.0000000 0.000 0.000 0.000 0.000
31 -0.0312500 -0.0002027 -0.157 -0.157 -0.511 -0.714
30 -0.0625000 -0.0008157 -0.317 -0.317 -1.431 -2.245
29 -0.0937500 -0.0018538 -0.478 -0.478 -2.773 -4.622
28 -0.1250000 -0.0033418 -0.642 -0.642 -4.563 -7.894
27 -0.1562500 -0.0053155 -0.808 -0.808 -6.837 -12.133
26 -0.1875000 -0.0078222 -0.977 -0.977 -9.640 -17.432
25 -0.2187500 -0.0109217 -1.148 -1.148 -13.034 -23.911
24 -0.2500000 -0.0146880 -1.323 -1.323 -17.091 -31.717
23 -0.2812500 -0.0192105 -1.502 -1.502 -21.901 -41.028
22 -0.3125000 -0.0245961 -1.685 -1.685 -27.570 -52.057
21 -0.3437500 -0.0309716 -1.874 -1.874 -34.223 -65.055
20 -0.3750000 -0.0384860 -2.067 -2.067 -42.010 -80.320
19 -0.4062500 -0.0473137 -2.268 -2.268 -51.104 -98.199
18 -0.4375000 -0.0576579 -2.476 -2.476 -61.707 -119.095
17 -0.4687500 -0.0697544 -2.693 -2.693 -74.054 -143.479
16 -0.5000000 -0.0838752 -2.919 -2.919 -88.415 -171.891
15 -0.5312500 -0.1003338 -3.158 -3.158 -105.102 -204.955
14 -0.5625000 -0.1194890 -3.410 -3.410 -124.473 -243.386
13 -0.5937500 -0.1417503 -3.678 -3.678 -146.935 -287.998
12 -0.6250000 -0.1675822 -3.963 -3.963 -172.949 -339.716
11 -0.6562500 -0.1975084 -4.269 -4.269 -203.038 -399.581
10 -0.6875000 -0.2321148 -4.598 -4.598 -237.784 -468.761
9 -0.7187500 -0.2720516 -4.954 -4.954 -277.835 -548.548
8 -0.7500000 -0.3180318 -5.340 -5.340 -323.899 -640.362
7 -0.7812500 -0.3708273 -5.760 -5.760 -376.745 -745.738
6 -0.8125000 -0.4312588 -6.219 -6.219 -437.190 -866.311
5 -0.8437500 -0.5001797 -6.720 -6.720 -506.082 -1003.778
4 -0.8750000 -0.5784503 -7.267 -7.267 -584.279 -1159.851
3 -0.9062500 -0.6669010 -7.866 -7.866 -672.605 -1336.184
2 -0.9375000 -0.7662811 -8.519 -8.519 -771.808 -1534.266
1 -0.9687500 -0.8771915 -9.230 -9.230 -882.484 -1755.294
0 -1.0000000 -1.0000000 -10.000 -10.000 -1005.000 -2000.000

Time Splitting Weights for Grid 01: 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.0022367614257357 0.4026337758325831
13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801076
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029188
16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
19 0.0199444086685725 0.0311389577709445 0.0857756755402367 0.6298591744489550
20 0.0227631639997064 0.0304741441486588 0.1085388395399431 0.6603333185976138
21 0.0256737146312910 0.0297153720153352 0.1342125541712341 0.6900486906129490
22 0.0286498597812016 0.0288595815276255 0.1628624139524357 0.7189082721405746
23 0.0316613792205220 0.0279045862015855 0.1945237931729577 0.7468128583421600
24 0.0346736416507075 0.0268492068942347 0.2291974348236651 0.7736620652363948
25 0.0376471948657328 0.0256934188392112 0.2668446296893979 0.7993554840756060
26 0.0405373376992232 0.0244385123436867 0.3073819673886211 0.8237939964192927
27 0.0432936737565711 0.0230872677537126 0.3506756411451922 0.8468812641730054
28 0.0458596469320356 0.0216441452951603 0.3965352880772278 0.8685254094681656
29 0.0481720587108284 0.0201154903974257 0.4447073467880562 0.8886408998655914
30 0.0501605672561820 0.0185097551070648 0.4948679140442381 0.9071506549726561
31 0.0517471682814030 0.0168377361985254 0.5466150823256412 0.9239883911711815
32 0.0528456577069106 0.0151128305891453 0.5994607400325517 0.9391012217603267
33 0.0533610761022577 0.0133513086655816 0.6528218161348094 0.9524525304259084
34 0.0531891349131379 0.0115726061288397 0.7060109510479473 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213234 0.9738247715198165
36 0.0503158038019031 0.0080591141492892 0.8085423793232265 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079418 0.9882658063583315
38 0.0431818225418188 0.0048034616164019 0.8990779740497606 0.9930692679747335
39 0.0376397765791564 0.0033640675316746 0.9367177506289170 0.9964333355064081
40 0.0305543017255206 0.0021094083123694 0.9672720523544377 0.9985427438187775
41 0.0217382098544505 0.0010909315881854 0.9890102622088881 0.9996336754069629
42 0.0109897377911119 0.0003663245930371 1.0000000000000000 1.0000000000000000

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

Metrics information for Grid 01:
===============================

Minimum X-grid spacing, DXmin = 0.00000000E+00 km
Maximum X-grid spacing, DXmax = 1.24434963E+01 km
Minimum Y-grid spacing, DYmin = 0.00000000E+00 km
Maximum Y-grid spacing, DYmax = 1.10642938E+01 km
Minimum Z-grid spacing, DZmin = 0.00000000E+00 m
Maximum Z-grid spacing, DZmax = 2.44706031E+02 m

Minimum barotropic Courant Number = 1.41636580E-02
Maximum barotropic Courant Number = 1.70393192E-01
Maximum Coriolis Courant Number = 1.27838438E-02


Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 1.6884675908E+15 m3
MinVolume = 2.5425025520E+07 m3
MaxVolume = 3.5401977069E+10 m3
Max/Min = 1.3924067467E+03

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


TIME-STEP YYYY-MM-DD hh:mm:ss.ss KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed

0 0001-01-01 00:00:00.00 0.000000E+00 9.721168E+03 9.721168E+03 1.690551E+15
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: /home/ahumada/Roms_tools/Run/SCRATCH/ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) in record = 0000001
DEF_DIAGS - creating diagnostics file, Grid 01: /home/ahumada/Roms_tools/Run/SCRATCH/ocean_dia.nc
1 0001-01-01 00:05:00.00 1.134897E-03 9.732794E+03 9.732795E+03 1.688434E+15
(01,56,32) 2.204008E-02 6.220827E-03 0.000000E+00 1.124059E+00
2 0001-01-01 00:10:00.00 2.759096E-02 9.731637E+03 9.731665E+03 1.688210E+15
(88,56,32) 1.494777E-02 0.000000E+00 2.927881E+04 6.782240E+00
3 0001-01-01 00:15:00.00 5.606112E-02 9.730155E+03 9.730211E+03 1.687893E+15
(88,56,32) 2.743207E-02 0.000000E+00 1.056139E+05 8.911713E+00
4 0001-01-01 00:20:00.00 9.784686E-02 9.728696E+03 9.728794E+03 1.687556E+15
(88,56,32) 6.252911E-02 0.000000E+00 1.437877E+05 1.133421E+01
5 0001-01-01 00:25:00.00 1.488974E-01 9.727229E+03 9.727378E+03 1.687204E+15
(88,56,32) 7.666016E-02 0.000000E+00 1.657057E+05 1.452295E+01
6 0001-01-01 00:30:00.00 2.098312E-01 9.725725E+03 9.725934E+03 1.686839E+15
(88,56,32) 1.273641E-02 0.000000E+00 1.805207E+05 2.504551E+01

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,1) in record = 0000001

Elapsed CPU time (seconds):

Thread # 19 CPU: 1.634
Thread # 20 CPU: 1.635
Thread # 16 CPU: 1.635
Thread # 5 CPU: 1.635
Thread # 7 CPU: 1.635
Thread # 9 CPU: 1.635
Thread # 13 CPU: 1.635
Thread # 1 CPU: 1.635
Thread # 10 CPU: 1.634
Thread # 22 CPU: 1.635
Thread # 21 CPU: 1.634
Thread # 15 CPU: 1.635
Thread # 3 CPU: 1.635
Thread # 23 CPU: 1.635
Thread # 0 CPU: 1.635
Thread # 18 CPU: 1.635
Thread # 6 CPU: 1.635
Thread # 4 CPU: 1.635
Thread # 11 CPU: 1.635
Thread # 8 CPU: 1.635
Thread # 14 CPU: 1.635
Thread # 2 CPU: 1.633
Thread # 17 CPU: 1.635
Thread # 12 CPU: 1.635
Total: 39.229

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 26.134 (66.6193 %)
Ocean state initialization ....................... 2.067 ( 5.2697 %)
Reading of input data ............................ 0.000 ( 0.0000 %)
Processing of input data ......................... 0.005 ( 0.0122 %)
Processing of output time averaged data .......... 0.599 ( 1.5270 %)
Computation of vertical boundary conditions ...... 0.003 ( 0.0064 %)
Computation of global information integrals ...... 0.042 ( 0.1072 %)
Writing of output data ........................... 0.085 ( 0.2171 %)
Model 2D kernel .................................. 3.614 ( 9.2137 %)
2D/3D coupling, vertical metrics ................. 0.122 ( 0.3108 %)
Omega vertical velocity .......................... 0.064 ( 0.1640 %)
Equation of state for seawater ................... 0.130 ( 0.3311 %)
GLS vertical mixing parameterization ............. 1.224 ( 3.1214 %)
3D equations right-side terms .................... 0.287 ( 0.7310 %)
3D equations predictor step ...................... 0.742 ( 1.8919 %)
Pressure gradient ................................ 0.123 ( 0.3133 %)
Harmonic mixing of tracers, S-surfaces ........... 0.088 ( 0.2240 %)
Harmonic stress tensor, S-surfaces ............... 0.099 ( 0.2525 %)
Corrector time-step for 3D momentum .............. 0.499 ( 1.2710 %)
Corrector time-step for tracers .................. 0.379 ( 0.9672 %)
Total: 36.307 92.5510

All percentages are with respect to total time = 39.229


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_initial.h
ROMS/Functionals/ana_m2obc.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_srflux.h
ROMS/Functionals/ana_stflux.h

ROMS/TOMS: DONE... Wednesday - June 14, 2017 - 5:17:35 PM
------------------------------------------------------------------------------------------------------------------------------------
I have no idea why the problem continue... Do you think can help me with this problem?
-
Thanks in advance =)
-
Mar.Mo.

[kate]

Your bathymetry is still very cliff-like and also inconsistent with this:

Code: Select all

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

I would start with a finite minimum depth (20m??) and turn off wetting and drying to see if you can get that running first.

[Scarlett]

Hi Kate:
So many thanks for your help, it was so useful, it run at 20 m deep!
-
I've been doing some sensitive test, and It's run at 180m maximum deep. But, What would happen if I want to considered a realistic bathymetry?How can I avoid this inconsistency, without smoothing so much the bathymetry? What do I need to considered?
I want to use 500 m deep
-
I'm working with:
-
r=0.5
n_filter_deep_topo=2
n_filter_final=2

theta_s=5
theta_b=.4
Tcline=10
N=32

Vtransform == 2
Vstretching == 4
-
Could you help me? Annex figure with real bathymetry
Thanks in advance.
-
Mar.Mo.

[kate]

Did you change the boundary conditions?

The way to model a steeper bathymetry is more horizontal resolution. With bathymetry smoothing there are trade-offs between realistic geometry and good model behavior. You just have to think about what you are resolving and what you need to resolve to get the answers you want. Plus try it and see.