ROMS BLOWUP MaxDensity

[markus.qs]

Hello,
I'm having a difficult time trying to run the ROMS model on a local lake. The error states that the blow-up is because of the maxdensity, here 7.3*10⁵. How do I correct it?
Below is some of the output from ROMS:



Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 1.21456751E-01 km Water points = 1.21619998E-01 km
Maximum X-grid spacing, DXmax = 1.21961210E-01 km Water points = 1.21961210E-01 km
Minimum Y-grid spacing, DYmin = 5.96607629E-02 km Water points = 5.97458276E-02 km
Maximum Y-grid spacing, DYmax = 5.99181318E-02 km Water points = 5.99181318E-02 km
Minimum Z-grid spacing, DZmin = 3.33486014E-02 m Water points = 5.43715807E-02 m
Maximum Z-grid spacing, DZmax = 3.33111113E+01 m Water points = 3.33111113E+01 m

Minimum barotropic Courant Number = 1.32849609E-02
Maximum barotropic Courant Number = 1.26711469E-01
Maximum Coriolis Courant Number = 1.25369904E-03

Horizontal mixing scaled by grid area squared root, MAXVAL(grdscl) = 8.54850153E-02 km (Water points)


AND THE ERROR:

Found Error: 01 Line: 314 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 301 Source: ROMS/Drivers/nl_roms.h, ROMS_run

Blowing-up: Saving latest model state into RESTART file
REASON: MaxDensity = 7.3E+25

WRT_RST_NF90 - wrote re-start fields (Index=2,2) in record = 1



Fyi, the bathymetry is smoothed and ROMS manage to run the first timestep before blown-up.

Best regard
Markus QS

[ymamoutos]

Greetings

Is it possible to post your log file?
What are the options you are using
at the *.h? What is your dt and dtfast?

Giannis

[markus.qs]

Thanks for the reply
I figured it out eventually. The shortwave radiation values generated from the forcingfile were way to high, and therefore the total energy became extremely negative after just one time step

Sorry for the disturbance

[ast212183]

hello markus.qs ,

can you please tell what changes have you made to the header file/input file in order to solve the issue. I am very new to ROMS and I am also getting similar error.


ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 1
number of time records written in RESTART file = 1

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_nudgcoef.h

MAIN: Abnormal termination: BLOWUP.
REASON: MaxDensity = 3.7E+02