DIAG speed trouble

[hpd14thu]

Dear all,

I'm trying to run ROMs by bulk forcing, river forcing and boundary values, but it blows up from the beginning. And ROMS doesn't show the details.

I wonder that which problems can result in this blow-up, then I can have a check.

Thank you!



----------------------------

mab

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -ftree-vectorize -ftree-loop-linear -funroll-loops -w -ffree-form -ffree-line-length-none -frecord-marker=4 -fconvert=big-endian

Input Script :

SVN Root URL : https:://myroms.org/svn/src
SVN Revision : 65M

Local Root : /home/ubuntu/app/COAWST_V3.2
Header Dir : /home/ubuntu/app/COAWST_V3.2/Projects/mab
Header file : mab.h
Analytical Dir: /home/ubuntu/app/COAWST_V3.2/Projects/mab

Resolution, Grid 01: 0166x0198x036, Parallel Nodes: 1, Tiling: 001x001


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

9216 ntimes Number of timesteps for 3-D equations.
150.000 dt Timestep size (s) for 3-D equations.
40 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.
120 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
1 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
360 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
2.0000E-01 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
2.0000E-01 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E-01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
0.000 wec_alpha WEC factor for roller/breaking energy distribution.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
2.5000E-02 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
2.0000E+00 blk_ZQ Height (m) of surface air humidity measurement.
2.0000E+00 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
1 lmd_Jwt Jerlov water type.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
5.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
50.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)
1.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
1.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Tnudg_SSS Nudging/relaxation time scale (days)
for sea surface salinity.
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
1.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.5000E+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.
10.000 T0 Background potential temperature (C) constant.
30.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T LuvSrc Turning ON momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
T LtracerSrc(01) Turning ON point Sources/Sinks on tracer 01: temp
T LtracerSrc(02) Turning ON 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(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(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt

T Aout(idSdif) Write out averaged vertical diffusion: AKt(isalt).

Output/Input Files:

Output Restart File: mab_rst.nc
Output History File: mab_his.nc
Output Averages File: mab_avg.nc
Input Grid File: /home/ubuntu/app/COAWST_V3.2/Projects/mab/mab_grid.nc
Input Nonlinear Initial File: /home/ubuntu/app/COAWST_V3.2/Projects/mab/mab_init.nc
Input Sources/Sinks File: /home/ubuntu/app/COAWST_V3.2/Projects/mab/mab_rivers.nc
Input Forcing File 01: /home/ubuntu/app/COAWST_V3.2/Projects/mab/mab_wind.nc
Input Boundary File 01: /home/ubuntu/app/COAWST_V3.2/Projects/mab/mab_boundary.nc

Tile partition information for Grid 01: 0166x0198x0036 tiling: 001x001

tile Istr Iend Jstr Jend Npts

Number of tracers: 2
0 1 166 1 198 1183248

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 166.50 0.50 198.50 RHO-points

0 1.00 166.00 0.50 198.50 U-points

0 0.50 166.50 1.00 198.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 354
HaloSizeJ(1) = 418
TileSide(1) = 203
TileSize(1) = 34713


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

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

zeta 1 Chapman Imp Chapman Imp Chapman Imp Chapman Imp

ubar 1 Flather Flather Flather Flather

vbar 1 Flather Flather Flather Flather

u 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud

v 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud

temp 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud

salt 1 Rad + Nud Rad + Nud Rad + Nud Rad + Nud

tke 1 Gradient Gradient Gradient Gradient

Activated C-preprocessing Options:

MAB mab
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
ATM_PRESS Impose atmospheric pressure onto sea surface.
AVERAGES Writing out time-averaged nonlinear model fields.
BULK_FLUXES Surface bulk fluxes parameterization.
CHARNOK Charnok surface roughness from wind stress.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EMINUSP Compute Salt Flux using E-P.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
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_LOGDRAG Logarithmic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.

Process Information:

Node # 0 (pid= 26220) 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

36 0.0000000 0.0000000 0.000 0.000 0.000 0.000
35 -0.0277778 -0.0001601 -0.053 -0.698 -1.798 -2.251
34 -0.0555556 -0.0006434 -0.107 -1.405 -4.471 -6.268
33 -0.0833333 -0.0014592 -0.160 -2.120 -8.045 -12.102
32 -0.1111111 -0.0026229 -0.214 -2.843 -12.561 -19.837
31 -0.1388889 -0.0041565 -0.267 -3.576 -18.079 -29.593
30 -0.1666667 -0.0060892 -0.321 -4.319 -24.678 -41.529
29 -0.1944444 -0.0084574 -0.375 -5.073 -32.457 -55.845
28 -0.2222222 -0.0113058 -0.428 -5.838 -41.536 -72.785
27 -0.2500000 -0.0146880 -0.482 -6.617 -52.062 -92.642
26 -0.2777778 -0.0186675 -0.536 -7.411 -64.205 -115.761
25 -0.3055556 -0.0233186 -0.589 -8.222 -78.167 -142.551
24 -0.3333333 -0.0287278 -0.643 -9.052 -94.182 -173.482
23 -0.3611111 -0.0349951 -0.697 -9.903 -112.522 -209.103
22 -0.3888889 -0.0422359 -0.751 -10.778 -133.499 -250.042
21 -0.4166667 -0.0505823 -0.805 -11.681 -157.470 -297.023
20 -0.4444444 -0.0601856 -0.859 -12.616 -184.846 -350.870
19 -0.4722222 -0.0712180 -0.914 -13.586 -216.093 -412.527
18 -0.5000000 -0.0838752 -0.968 -14.597 -251.740 -483.061
17 -0.5277778 -0.0983788 -1.023 -15.654 -292.389 -563.683
16 -0.5555556 -0.1149791 -1.077 -16.763 -338.717 -655.762
15 -0.5833333 -0.1339573 -1.132 -17.932 -391.486 -760.833
14 -0.6111111 -0.1556284 -1.187 -19.168 -451.549 -880.618
13 -0.6388889 -0.1803439 -1.243 -20.481 -519.857 -1017.038
12 -0.6666667 -0.2084930 -1.298 -21.879 -597.467 -1172.218
11 -0.6944444 -0.2405049 -1.354 -23.374 -685.539 -1348.504
10 -0.7222222 -0.2768484 -1.410 -24.977 -785.343 -1548.457
9 -0.7500000 -0.3180318 -1.467 -26.701 -898.257 -1774.856
8 -0.7777778 -0.3645994 -1.524 -28.559 -1025.755 -2030.672
7 -0.8055556 -0.4171263 -1.581 -30.567 -1169.394 -2319.050
6 -0.8333333 -0.4762095 -1.639 -32.739 -1330.792 -2643.251
5 -0.8611111 -0.5424547 -1.698 -35.089 -1511.589 -3006.584
4 -0.8888889 -0.6164579 -1.757 -37.634 -1713.399 -3412.306
3 -0.9166667 -0.6987799 -1.817 -40.386 -1937.742 -3863.481
2 -0.9444444 -0.7899141 -1.877 -43.359 -2185.953 -4362.806
1 -0.9722222 -0.8902454 -1.938 -46.562 -2459.076 -4912.382
0 -1.0000000 -1.0000000 -2.000 -50.000 -2757.723 -5513.446

Time Splitting Weights for Grid 01: ndtfast = 40 nfast = 56
==================================

Primary Secondary Accumulated to Current Step

1-0.0006038999279304 0.0250000000000000-0.0006038999279304 0.0250000000000000
2-0.0011053727845635 0.0250150974981983-0.0017092727124939 0.0500150974981983
3-0.0015044206959555 0.0250427318178123-0.0032136934084494 0.0750578293160106
4-0.0018010519301026 0.0250803423352112-0.0050147453385520 0.1001381716512218
5-0.0019952879837948 0.0251253686334638-0.0070100333223467 0.1252635402846856
6-0.0020871735042116 0.0251752508330587-0.0090972068265583 0.1504387911177443
7-0.0020767890452586 0.0252274301706640-0.0111739958718169 0.1756662212884082
8-0.0019642666586465 0.0252793498967954-0.0131382625304634 0.2009455711852036
9-0.0017498083197106 0.0253284565632616-0.0148880708501740 0.2262740277484652
10-0.0014337071879729 0.0253722017712543-0.0163217780381469 0.2516462295197196
11-0.0010163717024450 0.0254080444509537-0.0173381497405919 0.2770542739706733
12-0.0004983525116731 0.0254334537435148-0.0178365022522649 0.3024877277141881
13 0.0001196277614762 0.0254459125563066-0.0177168744907887 0.3279336402704947
14 0.0008366419202876 0.0254429218622697-0.0168802325705011 0.3533765621327644
15 0.0016515227599278 0.0254220058142625-0.0152287098105733 0.3787985679470269
16 0.0025628247984349 0.0253807177452643-0.0126658850121384 0.4041792856922913
17 0.0035687831729657 0.0253166471253035-0.0090971018391727 0.4294959328175947
18 0.0046672697013015 0.0252274275459793-0.0044298321378712 0.4547233603635741
19 0.0058557461086128 0.0251107458034468 0.0014259139707416 0.4798341061670209
20 0.0071312144194823 0.0249643521507315 0.0085571283902239 0.5047984583177524
21 0.0084901645151857 0.0247860717902444 0.0170472929054096 0.5295845301079968
22 0.0099285188562319 0.0245738176773648 0.0269758117616416 0.5541583477853615
23 0.0114415743701613 0.0243256047059590 0.0384173861318028 0.5784839524913206
24 0.0130239415046020 0.0240395653467049 0.0514413276364049 0.6025235178380255
25 0.0146694804455855 0.0237139668090899 0.0661108080819904 0.6262374846471154
26 0.0163712345011199 0.0233472297979502 0.0824820425831103 0.6495847144450656
27 0.0181213606500220 0.0229379489354222 0.1006034032331323 0.6725226633804878
28 0.0199110572560077 0.0224849149191717 0.1205144604891400 0.6950075782996595
29 0.0217304889470413 0.0219871384877715 0.1422449494361813 0.7169947167874310
30 0.0235687086599425 0.0214438762640955 0.1658136580961238 0.7384385930515265
31 0.0254135768502523 0.0208546585475969 0.1912272349463761 0.7592932515991234
32 0.0272516778673573 0.0202193191263406 0.2184789128137333 0.7795125707254640
33 0.0290682334948725 0.0195380271796567 0.2475471463086058 0.7990505979051207
34 0.0308470136562825 0.0188113213422849 0.2783941599648883 0.8178619192474056
35 0.0325702442858409 0.0180401460008778 0.3109644042507291 0.8359020652482834
36 0.0342185123647286 0.0172258898937318 0.3451829166154578 0.8531279551420151
37 0.0357706681224704 0.0163704270846136 0.3809535847379282 0.8694983822266287
38 0.0372037244036097 0.0154761603815518 0.4181573091415379 0.8849745426081804
39 0.0384927531996420 0.0145460672714616 0.4566500623411799 0.8995206098796420
40 0.0396107793462069 0.0135837484414705 0.4962608416873868 0.9131043583211125
41 0.0405286713855386 0.0125934789578153 0.5367895130729254 0.9256978372789278
42 0.0412150295941741 0.0115802621731769 0.5780045426670996 0.9372780994521047
43 0.0416360711759210 0.0105498864333225 0.6196406138430206 0.9478279858854273
44 0.0417555126200830 0.0095089846539245 0.6613961264631036 0.9573369705393517
45 0.0415344492249439 0.0084650968384224 0.7029305756880475 0.9658020673777742
46 0.0409312317865102 0.0074267356077988 0.7438618074745578 0.9732288029855730
47 0.0399013404525124 0.0064034548131361 0.7837631479270702 0.9796322577987091
48 0.0383972557416642 0.0054059213018233 0.8221604036687344 0.9850381791005324
49 0.0363683267281804 0.0044459899082817 0.8585287303969148 0.9894841690088140
50 0.0337606363915535 0.0035367817400771 0.8922893667884684 0.9930209507488911
51 0.0305168641315885 0.0026927658302883 0.9228062309200569 0.9957137165791794
52 0.0265761454486959 0.0019298442269986 0.9493823763687528 0.9976435608061780
53 0.0218739287894439 0.0012654405907812 0.9712563051581966 0.9989090013969592
54 0.0163418295573680 0.0007185923710451 0.9875981347155646 0.9996275937680043
55 0.0099074812890402 0.0003100466321109 0.9975056160046047 0.9999376404001151
56 0.0024943839953957 0.0000623595998849 1.0000000000000004 1.0000000000000000

ndtfast, nfast = 40 56 nfast/ndtfast = 1.40000

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

1.000000000000 1.039265571007 0.519632785503 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.21300

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

Minimum X-grid spacing, DXmin = 4.68165556E+00 km
Maximum X-grid spacing, DXmax = 5.28701590E+00 km
Minimum Y-grid spacing, DYmin = 5.03543630E+00 km
Maximum Y-grid spacing, DYmax = 5.03543631E+00 km
Minimum Z-grid spacing, DZmin = 5.34311187E-02 m
Maximum Z-grid spacing, DZmax = 6.01063990E+02 m

Minimum barotropic Courant Number = 4.59023827E-03
Maximum barotropic Courant Number = 2.39919205E-01
Maximum Coriolis Courant Number = 1.46732585E-02


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: mab_init.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = -5.15861928E-01 Max = 1.01138091E+00)
- vertically integrated v-momentum component
(Min = -4.19422179E-01 Max = 1.48187801E-01)
- u-momentum component
(Min = -8.79346848E-01 Max = 1.79157114E+00)
- v-momentum component
(Min = -1.20992875E+00 Max = 1.26725924E+00)
- potential temperature
(Min = -1.39353603E-01 Max = 3.04190636E+01)
- salinity
(Min = 0.00000000E+00 Max = 3.68937263E+01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 1.60000000E+01 Max = 6.40000000E+01)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 4.40000000E+01 Max = 1.66000000E+02)
GET_NGFLD - river runoff direction
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 2.77777778E-02 Max = 2.77777778E-02)
GET_NGFLD - river runoff mass transport, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_rivers.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.24755000E+02 Max = 6.96598200E+01)
GET_NGFLD - river runoff potential temperature, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_rivers.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 2.28000000E+01 Max = 2.74000000E+01)
GET_NGFLD - river runoff salinity, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_rivers.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.48879018E+00 Max = 9.48846082E+00)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -3.85694262E+00 Max = 8.27674899E+00)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 9.41213884E+02 Max = 1.02045113E+03)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 1.99932619E+01 Max = 2.83775444E+01)
GET_2DFLD - surface air relative humidity, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 6.49802349E-01 Max = 9.73909050E-01)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 5.84608824E-05 Max = 1.85498452E-04)
GET_2DFLD - net longwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -3.63744317E-05 Max = 8.11468262E-07)
GET_2DFLD - rain fall rate, t = 0 00:00:00
(Rec=0000001, Index=1, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -7.73038670E-05 Max = 6.48162364E-04)
GET_NGFLD - free-surface western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.87438655E-01)
GET_NGFLD - 2D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -3.00609842E-02 Max = 2.07886964E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.73493677E-01 Max = 3.81799936E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.91124165E-01 Max = 7.67858148E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -5.40703654E-01 Max = 6.88199878E-01)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.17271268E-01 Max = 5.99525034E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.79301214E-02 Max = 1.39712080E-01)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.22154427E-02 Max = 1.26435488E-01)
GET_NGFLD - 3D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.54327916E-02 Max = 4.86692816E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -6.77571893E-02 Max = 2.64625221E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -7.07323015E-01 Max = 1.51773870E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -5.39381504E-01 Max = 1.03494978E+00)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -7.26773083E-01 Max = 1.13102090E+00)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.80627805E-01 Max = 8.61527681E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.25189328E-01 Max = 1.75515458E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -7.26434067E-02 Max = 2.17351064E-01)
GET_NGFLD - temperature western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.93972244E+01)
GET_NGFLD - salinity western boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.63715286E+01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.84820652E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.66891594E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.95831852E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.67100487E+01)
GET_NGFLD - temperature northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.05079613E+01)
GET_NGFLD - salinity northern boundary condition, t = 0 00:00:00
(Rec=0000001, Index=2, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.45250435E+01)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 1.000000E-01 (Beckmann and Haidvogel)
rx1 = 6.279593E+00 (Haney)

Initial basin volumes: TotVolume = 1.5937623726E+15 m3
MinVolume = 1.2613937078E+06 m3
MaxVolume = 1.5952114150E+10 m3
Max/Min = 1.2646419632E+04

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

GET_NGFLD - river runoff mass transport, t = 1 00:00:00
(Rec=0000002, Index=1, File: mab_rivers.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.07764800E+02 Max = 6.99429900E+01)
GET_NGFLD - river runoff potential temperature, t = 1 00:00:00
(Rec=0000002, Index=1, File: mab_rivers.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 2.28000000E+01 Max = 2.86000000E+01)
GET_NGFLD - river runoff salinity, t = 1 00:00:00
(Rec=0000002, Index=1, File: mab_rivers.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.67755614E+00 Max = 9.68387638E+00)
GET_2DFLD - surface v-wind component, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -3.95093540E+00 Max = 7.99857274E+00)
GET_2DFLD - surface air pressure, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 9.41149331E+02 Max = 1.02057863E+03)
GET_2DFLD - surface air temperature, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 2.01950278E+01 Max = 2.84300235E+01)
GET_2DFLD - surface air relative humidity, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 5.93761673E-01 Max = 9.72107867E-01)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 6.46495803E-05 Max = 1.92980171E-04)
GET_2DFLD - net longwave radiation flux, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -3.78018656E-05 Max = 3.42369378E-07)
GET_2DFLD - rain fall rate, t = 0 00:30:00
(Rec=0000002, Index=2, File: mab_wind.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -6.36967807E-05 Max = 6.65809200E-04)
GET_NGFLD - free-surface western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.93966287E-01)
GET_NGFLD - 2D v-momentum western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.13958170E-02 Max = 2.13748202E-01)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.73708850E-01 Max = 3.80629003E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.49451026E-01 Max = 7.85718977E-01)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -5.48984230E-01 Max = 6.93641663E-01)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.17264503E-01 Max = 5.96684694E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.58427779E-02 Max = 1.48350596E-01)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.49377657E-02 Max = 1.22255221E-01)
GET_NGFLD - 3D u-momentum western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.09561586E-02 Max = 4.87650812E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -8.43785405E-02 Max = 2.72036523E-01)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -6.84598923E-01 Max = 1.49319148E+00)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -5.40083051E-01 Max = 1.02259660E+00)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -7.36855447E-01 Max = 1.12686682E+00)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -4.84302938E-01 Max = 8.56351137E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -2.02859089E-01 Max = 1.92000151E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = -8.37866589E-02 Max = 2.36102208E-01)
GET_NGFLD - temperature western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.94202271E+01)
GET_NGFLD - salinity western boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.63712006E+01)
GET_NGFLD - temperature eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.84867821E+01)
GET_NGFLD - salinity eastern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.66894455E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.95820427E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.67096710E+01)
GET_NGFLD - temperature northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 2.05952473E+01)
GET_NGFLD - salinity northern boundary condition, t = 0 00:30:00
(Rec=0000002, Index=1, File: mab_boundary.nc)
(Tmin= 0.0000 Tmax= 16.0000)
(Min = 0.00000000E+00 Max = 3.45262375E+01)

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

0 0 00:00:00 3.853334E-02 2.104765E+04 2.104769E+04 1.594152E+15
(093,078,28) 4.597547E-02 3.116953E-02 0.000000E+00 1.872812E+00
DEF_HIS - creating history file, Grid 01: mab_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: mab_avg.nc
1 0 00:02:30 NaN NaN NaN NaN
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 NaN
DIAG speed trouble NaN

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001

Elapsed CPU time (seconds):

Node # 0 CPU: 1.870
Total: 1.870

Nonlinear ocean model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.124 ( 6.6262 %)
Ocean state initialization ....................... 0.020 ( 1.0781 %)
Reading of input data ............................ 0.040 ( 2.1340 %)
Processing of input data ......................... 0.001 ( 0.0588 %)
Processing of output time averaged data .......... 0.000 ( 0.0034 %)
Computation of vertical boundary conditions ...... 0.001 ( 0.0573 %)
Computation of global information integrals ...... 0.014 ( 0.7221 %)
Writing of output data ........................... 0.233 (12.4585 %)
Model 2D kernel .................................. 0.211 (11.2898 %)
2D/3D coupling, vertical metrics ................. 0.022 ( 1.1747 %)
Omega vertical velocity .......................... 0.009 ( 0.5021 %)
Equation of state for seawater ................... 0.059 ( 3.1786 %)
Atmosphere-Ocean bulk flux parameterization ...... 0.064 ( 3.4340 %)
GLS vertical mixing parameterization ............. 0.468 (25.0481 %)
3D equations right-side terms .................... 0.027 ( 1.4612 %)
3D equations predictor step ...................... 0.097 ( 5.1795 %)
Pressure gradient ................................ 0.016 ( 0.8332 %)
Harmonic mixing of tracers, geopotentials ........ 0.025 ( 1.3593 %)
Harmonic stress tensor, S-surfaces ............... 0.013 ( 0.6728 %)
Corrector time-step for 3D momentum .............. 0.037 ( 1.9942 %)
Corrector time-step for tracers .................. 0.029 ( 1.5679 %)
Total: 1.511 80.8339

Nonlinear model message Passage profile, Grid: 01

Message Passage: 2D halo exchanges ............... 0.001 ( 0.0438 %)
Message Passage: 3D halo exchanges ............... 0.000 ( 0.0032 %)
Message Passage: 4D halo exchanges ............... 0.000 ( 0.0006 %)
Message Passage: data broadcast .................. 0.002 ( 0.1176 %)
Message Passage: data reduction .................. 0.000 ( 0.0020 %)
Message Passage: data gathering .................. 0.053 ( 2.8338 %)
Message Passage: data scattering.................. 0.021 ( 1.1433 %)
Total: 0.077 4.1443

All percentages are with respect to total time = 1.870


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/TOMS: DONE... Tuesday - September 10, 2019 - 6:37:41 PM

[jcwarner]

you have several options activated. suggest you try a more simple first approach and then add complexity. Try it without the river forcings and maybe without the volcons stuff and see if that works to start.
-j

[lizhicheng]

I also find the nan; I change the initial state some times, but still in the step 362 blow up.

360 2017-01-02 00:00:00.00 4.248238E-02 7.014599E+02 7.015024E+02 4.461919E+13
(264,057,01) 1.373269E-02 4.212748E-02 1.065118E+00 2.235862E+00
WRT_HIS - wrote history fields (Index=1,1) in record = 0000025
WRT_AVG - wrote averaged fields in record = 0000001
GET_NGFLD - river runoff mass transport, 2000-01-08 00:00:00.00
(Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc)
(Tmin= 0.0000 Tmax= 364.0000) t = 7.0000
(Min = -1.27772940E+23 Max = 2.81650723E-37)
GET_NGFLD - river runoff potential temperature, 2000-01-08 00:00:00.00
(Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc)
(Tmin= 0.0000 Tmax= 364.0000) t = 7.0000
(Min = -3.60683066E+28 Max = 1.24658598E+11)
GET_NGFLD - river runoff salinity, 2000-01-08 00:00:00.00
(Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc)
(Tmin= 0.0000 Tmax= 364.0000) t = 7.0000
(Min = -2.35343828-185 Max = -2.35343828-185)
GET_2DFLD - surface u-momentum stress, 2000-01-07 06:00:00.00
(Grid=01, Rec=0000026, Index=1, File: byecs_smsUstr_from_windu_ana6hr_1993_cycle.nc)
(Tmin= 0.0000 Tmax= 364.7500) t = 6.2500
(Min = -1.53729570E-06 Max = 5.40840912E-05) regrid = T
GET_2DFLD - surface v-momentum stress, 2000-01-07 06:00:00.00
(Grid=01, Rec=0000026, Index=1, File: byecs_smsVstr_from_windv_ana6hr_1993_cycle.nc)
(Tmin= 0.0000 Tmax= 364.7500) t = 6.2500
(Min = -9.69953848E-05 Max = -2.60725450E-06) regrid = T
361 2017-01-02 00:04:00.00 4.294329E-02 7.014896E+02 7.015326E+02 4.462479E+13
(264,057,01) 1.344975E-02 4.154651E-02 1.049825E+00 2.224724E+00
362 2017-01-02 00:08:00.00 NaN NaN NaN NaN
(264,057,01) 1.314301E-02 4.086271E-02 1.032439E+00 0.000000E+00
Found Error: 01 Line: 298 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 298 Source: ROMS/Drivers/nl_ocean.h

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):

[kate]

GET_NGFLD - river runoff mass transport, 2000-01-08 00:00:00.00
(Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc)
(Tmin= 0.0000 Tmax= 364.0000) t = 7.0000
(Min = -1.27772940E+23 Max = 2.81650723E-37)
GET_NGFLD - river runoff potential temperature, 2000-01-08 00:00:00.00
(Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc)
(Tmin= 0.0000 Tmax= 364.0000) t = 7.0000
(Min = -3.60683066E+28 Max = 1.24658598E+11)
GET_NGFLD - river runoff salinity, 2000-01-08 00:00:00.00
(Grid= 01, Rec=0000008, Index=1, File: create-river_daily2017_cycle_add9rivers_30layers-15-change-huanghe-1996coastl_dye.nc)
(Tmin= 0.0000 Tmax= 364.0000) t = 7.0000
(Min = -2.35343828-185 Max = -2.35343828-185)

How did you make your rivers file? I don't like your min/max values at all.

[lizhicheng]

THANKS for your reply;
I run one case sucessfully which uses current coastline before. However, I change the grid using old coastline, initial file and river file because of river location.