ROMS/TOMS - Blows up ................ exit_flag: 1

[JIA]

I got the problem of blowing up.
I also check that there are any nan-data in the initial and grid file.

I'm testing WRF-ROMS, but I checked that the WRF only case works well. so I think this error probably came from the ROMS.


the below is .out file

=================================================================================
Coupled Input File name = ./Projects/Sandy/coupling_sandy.in
Coupled Input File name = ./Projects/Sandy/coupling_sandy.in

Model Coupling:


Ocean Model MPI nodes: 000 - 000

Atmos Model MPI nodes: 001 - 001

Model Input Parameters: ROMS/TOMS version 3.7
Thursday - July 16, 2015 - 11:23:14 PM
-----------------------------------------------------------------------------
Quilting with 1 groups of 0 I/O tasks.
Ntasks in X 1, ntasks in Y 1
--- WARNING: traj_opt is zero, but num_traj is not zero; setting num_traj to zero.
--- NOTE: grid_fdda is 0 for domain 1, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain 1, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain 1, setting obs nudging interval and ending time to 0 for that domain.
bl_pbl_physics /= 4, implies mfshconv must be 0, resetting
--- NOTE: num_soil_layers has been set to 4
WRF V3.6.1 MODEL
*************************************
Parent domain
ids,ide,jds,jde 1 130 1 133
ims,ime,jms,jme -4 135 -4 138
ips,ipe,jps,jpe 1 130 1 133
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate

Idealized Coastal shelf

Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /usr/local/mpich2/intel14_15/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -g -check uninit -ftrapuv -traceback -I/home/JIA/MCT/include -I../WRF/main -I../WRF/external/esmf_time_f90 -I../WRF/frame -I../WR

Input Script :

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

Local Root : /home/JIA/COAWST_nonest
Header Dir : /home/JIA/COAWST_nonest/Projects/Sandy
Header file : sandy.h
Analytical Dir: /home/JIA/COAWST_nonest/Projects/Sandy

Resolution, Grid 01: 0127x0130x016, Parallel Nodes: 1, Tiling: 001x001


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

5760 ntimes Number of timesteps for 3-D equations.
10.000 dt Timestep size (s) for 3-D equations.
28 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.
432 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).
864 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
432 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
0.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
0.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
5.0000E-02 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).
1 lmd_Jwt Jerlov water type.
1 Vtransform S-coordinate transformation equation.
1 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)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
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).
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 Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idUbms) Write out bottom U-momentum stress.
T Hout(idVbms) Write out bottom V-momentum stress.
T Hout(idVvis) Write out vertical viscosity: AKv.
T Hout(idMtke) Write out turbulent kinetic energy.
T Hout(idMtls) Write out turbulent generic length-scale.

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

Output/Input Files:

Output Restart File: ocean_rst.nc
Output History File: ocean_joe_tc_coarse_his.nc
Output Averages File: ocean_avg.nc
Physical parameters File:
Input Grid File: Projects/Sandy/Sandy_roms_grid.nc
Input Nonlinear Initial File: Projects/Sandy/Sandy_ini.nc

Tile partition information for Grid 01: 0127x0130x0016 tiling: 001x001

tile Istr Iend Jstr Jend Npts

Number of tracers: 2
0 1 127 1 130 264160

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

tile Xmin Xmax Ymin Ymax grid

0 0.50 127.50 0.50 130.50 RHO-points

0 1.00 127.00 0.50 130.50 U-points

0 0.50 127.50 1.00 130.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 274
HaloSizeJ(1) = 282
TileSide(1) = 135
TileSize(1) = 17685


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

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

zeta 1 Radiation Radiation Radiation Radiation

ubar 1 Radiation Radiation Radiation Radiation

vbar 1 Radiation Radiation Radiation Radiation

u 1 Radiation Radiation Radiation Radiation

v 1 Radiation Radiation Radiation Radiation

temp 1 Radiation Radiation Radiation Radiation

salt 1 Radiation Radiation Radiation Radiation

tke 1 Radiation Radiation Radiation Radiation

Activated C-preprocessing Options:

SANDY Idealized Coastal shelf
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
ATM_PRESS Impose atmospheric pressure onto sea surface.
ATM2OCN_FLUXES Surface fluxes from ATM Model.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MCT_LIB Using Model Coupling Toolkit library.
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.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
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.
WRF_COUPLING WRF model coupling.

Process Information:

Node # 0 (pid= 25876) is active.
alloc_space_field: domain 1, 496007556 bytes allocated
med_initialdata_input: calling input_input
Timing for processing wrfinput file (stream 0) for domain 1: 0.05231 elapsed seconds
INPUT LandUse = "USGS"
LANDUSE TYPE = "USGS" FOUND 33 CATEGORIES 2 SEASONS WATER CATEGORY = 16 SNOW CATEGORY = 24
INITIALIZE THREE Noah LSM RELATED TABLES
LANDUSE TYPE = USGS FOUND 27 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND 19 CATEGORIES

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

16 0.0000000 -0.0000000 0.000 0.000 -0.000 -0.000
15 -0.0625000 -0.0048948 213.291 213.291 200.048 186.805
14 -0.1250000 -0.0119914 426.582 426.582 394.139 361.696
13 -0.1875000 -0.0230908 639.873 639.873 577.400 514.928
12 -0.2500000 -0.0409944 853.164 853.164 742.253 631.343
11 -0.3125000 -0.0696342 1066.454 1066.454 878.059 689.664
10 -0.3750000 -0.1133185 1279.745 1279.745 973.162 666.579
9 -0.4375000 -0.1742175 1493.036 1493.036 1021.691 550.346
8 -0.5000000 -0.2489214 1706.327 1706.327 1032.871 359.414
7 -0.5625000 -0.3284417 1919.618 1919.618 1031.019 142.421
6 -0.6250000 -0.4042644 2132.909 2132.909 1039.172 -54.565
5 -0.6875000 -0.4744488 2346.200 2346.200 1062.579 -221.041
4 -0.7500000 -0.5437741 2559.491 2559.491 1088.310 -382.870
3 -0.8125000 -0.6205544 2772.781 2772.781 1093.872 -585.037
2 -0.8750000 -0.7145184 2986.072 2986.072 1052.944 -880.185
1 -0.9375000 -0.8366254 3199.363 3199.363 935.874 -1327.615
0 -1.0000000 -1.0000000 3412.654 3412.654 707.155 -1998.344

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

Primary Secondary Accumulated to Current Step

1-0.0008549076326633 0.0357142857142857-0.0008549076326633 0.0357142857142857
2-0.0014558918889834 0.0357448181297380-0.0023107995216466 0.0714591038440237
3-0.0018029740646893 0.0357968142686302-0.0041137735863359 0.1072559181126539
4-0.0018962369765050 0.0358612061995120-0.0060100105628409 0.1431171243121659
5-0.0017358959479127 0.0359289289486729-0.0077459065107536 0.1790460532608388
6-0.0013223981892219 0.0359909252325269-0.0090683046999755 0.2150369784933658
7-0.0006565505719430 0.0360381537392848-0.0097248552719185 0.2510751322326506
8 0.0002603242025330 0.0360616019739971-0.0094645310693855 0.2871367342066477
9 0.0014262030399494 0.0360523046810495-0.0080383280294361 0.3231890388876972
10 0.0028381494958941 0.0360013688581941-0.0052001785335421 0.3591904077458913
11 0.0044920724242039 0.0359000063761979-0.0007081061093381 0.3950904141220893
12 0.0063824562310637 0.0357395752181906 0.0056743501217256 0.4308299893402799
13 0.0085020627348001 0.0355116303527955 0.0141764128565257 0.4663416196930755
14 0.0108416046313699 0.0352079852551241 0.0250180174878956 0.5015496049481996
15 0.0133893905655432 0.0348207850897180 0.0384074080534388 0.5363703900379175
16 0.0161309418077813 0.0343425925695200 0.0545383498612201 0.5707129826074375
17 0.0190485805368085 0.0337664875049564 0.0735869303980287 0.6044794701123939
18 0.0221209897278798 0.0330861810572133 0.0957079201259085 0.6375656511696072
19 0.0253227446467419 0.0322961457097890 0.1210306647726504 0.6698617968793962
20 0.0286238159492900 0.0313917619724053 0.1496544807219404 0.7012535588518015
21 0.0319890443869185 0.0303694828313593 0.1816435251088589 0.7316230416831608
22 0.0353775871175665 0.0292270169603979 0.2170211122264254 0.7608500586435587
23 0.0387423356224582 0.0279635317061991 0.2557634478488836 0.7888135903497578
24 0.0420293052285376 0.0265798768625399 0.2977927530774212 0.8153934672122977
25 0.0451769962365983 0.0250788302472350 0.3429697493140195 0.8404722974595327
26 0.0481157266551073 0.0234653660959279 0.3910854759691268 0.8639376635554605
27 0.0507669365397241 0.0217469472868169 0.4418524125088509 0.8856846108422775
28 0.0530424639385140 0.0199338424103982 0.4948948764473649 0.9056184532526756
29 0.0548437924428562 0.0180394686983084 0.5497386688902210 0.9236579219509841
30 0.0560612703440466 0.0160807618253492 0.6057999392342676 0.9397386837763333
31 0.0565733013955950 0.0140785735987762 0.6623732406298626 0.9538172573751095
32 0.0562455071812177 0.0120580985489335 0.7186187478110803 0.9658753559240429
33 0.0549298610885233 0.0100493304353186 0.7735486088996035 0.9759246863593615
34 0.0524637938883947 0.0080875496821570 0.8260124027879983 0.9840122360415184
35 0.0486692709200650 0.0062138427575715 0.8746816737080633 0.9902260787990900
36 0.0433518408818877 0.0044756545104263 0.9180335145899510 0.9947017333095163
37 0.0362996562278024 0.0029273744789303 0.9543331708177534 0.9976291077884466
38 0.0272824651694943 0.0016309581850802 0.9816156359872477 0.9992600659735268
39 0.0160505752842493 0.0006565844290269 0.9976662112714970 0.9999166504025536
40 0.0023337887285028 0.0000833495974465 0.9999999999999998 1.0000000000000002

ndtfast, nfast = 28 40 nfast/ndtfast = 1.42857

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

1.000000000000 1.049755056052 0.524877528026 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.18257

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

Minimum X-grid spacing, DXmin = 8.94129798E+00 km
Maximum X-grid spacing, DXmax = 8.99699403E+00 km
Minimum Y-grid spacing, DYmin = 8.94321390E+00 km
Maximum Y-grid spacing, DYmax = 8.99643161E+00 km
Minimum Z-grid spacing, DZmin = -2.13290881E+02 m
Maximum Z-grid spacing, DZmax = 6.70728644E+02 m

Minimum barotropic Courant Number = 1.44052346E-05
Maximum barotropic Courant Number = 1.03136365E-02
Maximum Coriolis Courant Number = 9.75576986E-04


NLM: GET_STATE - Read state initial conditions, t = 56493 00:00:00
(Grid 01, File: Sandy_ini.nc, Rec=0001, Index=1)
- free-surface
(Min = -6.77090883E-02 Max = 6.85521662E-01)
- vertically integrated u-momentum component
(Min = -6.95050062E+05 Max = 3.75780688E+05)
- vertically integrated v-momentum component
(Min = -8.86146375E+05 Max = 6.06092438E+05)
- u-momentum component
(Min = -2.70264500E+06 Max = 1.20221945E+05)
- v-momentum component
(Min = -3.31302975E+06 Max = 2.43865234E+05)
- potential temperature
(Min = -4.22937750E+06 Max = 6.49815160E+07)
- salinity
(Min = -2.62840350E+06 Max = 2.62479775E+06)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 1.157955E+02 (Beckmann and Haidvogel)
rx1 = 8.751518E+04 (Haney)

Initial basin volumes: TotVolume = -2.5813502425E+14 m3
MinVolume = -1.7124735277E+10 m3
MaxVolume = 3.5765252533E+10 m3
Max/Min = -2.0885141846E+00
## WRF grid 1 sent data to ROMS grid 1
## ROMS grid 1 recv data from WRF grid 1


## ROMS grid 1 sent data to WRF grid 1



NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: ******** - 00005760)

## WRF grid 1 recvd data from ROMS grid 1

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

******* 56493 00:00:00 2.907356E+11 -3.271950E+04 2.907356E+11 -9.934884E+13
(121,117,01) 2.836295E+03 3.544466E+03 0.000000E+00 4.072956E+06
DIAG speed trouble 4072955.95644003
DIAG speed ijk 243865.234375000 119 62 2

Elapsed CPU time (seconds):

Node # 0 CPU: 1.658
Total: 1.658

Nonlinear ocean model elapsed time profile:

Allocation and array initialization .............. 0.050 ( 3.0157 %)
Ocean state initialization ....................... 0.011 ( 0.6634 %)
Computation of global information integrals ...... 0.003 ( 0.1809 %)
2D/3D coupling, vertical metrics ................. 0.007 ( 0.4222 %)
Omega vertical velocity .......................... 0.002 ( 0.1206 %)
Equation of state for seawater ................... 0.011 ( 0.6635 %)
Total: 0.084 5.0663

Nonlinear model message Passage profile:

Message Passage: data scattering.................. 0.005 ( 0.3015 %)
Message Passage: boundary data gathering ......... 0.003 ( 0.1810 %)
Total: 0.008 0.4825

All percentages are with respect to total time = 1.658

ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_stflux.h

ROMS/TOMS - Blows up ................ exit_flag: 1


MAIN: Abnormal termination: BLOWUP.
Timing for Writing wrfout_d01_2013-07-20_00:00:00 for domain 1: 0.16819 elapsed seconds
d01 2013-07-20_00:00:00 Input data processed for aux input 4 for domain 1
Timing for processing lateral boundary for domain 1: 0.00987 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 130 JS 1 JE 133
WRF NUMBER OF TILES = 1

===================================================================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= EXIT CODE: 134
= CLEANING UP REMAINING PROCESSES
= YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Aborted (signal 6)
This typically refers to a problem with your application.
Please see the FAQ page for debugging suggestions

================================================================================



Thus, my cpp file

#define ROMS_MODEL
#undef NESTING
#define WRF_MODEL
#undef SWAN_MODEL
#define MCT_LIB
#define MCT_INTERP_OC2AT
#undef MCT_INTERP_WV2AT
#undef MCT_INTERP_OC2WV

#if defined WRF_MODEL && defined SWAN_MODEL
# define DRAGLIM_DAVIS
# define COARE_TAYLOR_YELLAND
#endif

#ifdef ROMS_MODEL
/* Physics + numerics */
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define MIX_S_UV
# undef TS_FIXED
# define TS_U3HADVECTION
# define TS_C4VADVECTION
# undef TS_MPDATA

# undef SSW_BBL
# ifdef SSW_BBL
# define SSW_CALC_ZNOT
# /*define ANA_SEDIMENT*/
# else
# define UV_LOGDRAG
# endif
# if !defined SWAN_MODEL && defined SSW_BBL
# define ANA_WWAVE
# endif

# define DJ_GRADPS
# define TS_DIF2
# define MIX_GEO_TS
# define CURVGRID

# define SALINITY
# define SOLVE3D
# define SPLINES
# define AVERAGES
# define NONLIN_EOS

/* Grid and Initial */
# define MASKING

/* Forcing */
# ifdef WRF_MODEL
# undef BULK_FLUXES
# define ATM2OCN_FLUXES
# define ANA_SSFLUX
# undef LONGWAVE_OUT
# define ATM_PRESS
# else
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX
# endif
# define ANA_BTFLUX
# define ANA_BSFLUX
# define EMINUSP
# define SOLAR_SOURCE

/* Turbulence closure */
# define GLS_MIXING
# undef MY25_MIXING
# define AKLIMIT

# if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif

/* Output */
# undef DIAGNOSTICS_UV
# undef DIAGNOSTICS_TS
#endif

[rtopper]

- potential temperature
(Min = -4.22937750E+06 Max = 6.49815160E+07)
- salinity
(Min = -2.62840350E+06 Max = 2.62479775E+06)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 1.157955E+02 (Beckmann and Haidvogel)
rx1 = 8.751518E+04 (Haney)

Initial basin volumes: TotVolume = -2.5813502425E+14 m3
MinVolume = -1.7124735277E+10 m3
MaxVolume = 3.5765252533E+10 m3
Max/Min = -2.0885141846E+00

Based on your output I would say that your ROMS input files are causing the problem. Negative basin volumes and temperature/salinity in the +/- 1e6 - 1e7 range are not reasonable.

EDIT: Just noticed your post is already 2 weeks old, but hope the reply still helps

[lalvarez]

Hello

I guess you need to smooth your grid rx0 is too big