ROMS BLOWUP MaxDensity

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
markus.qs
Posts: 9
Joined: Mon Jun 21, 2021 12:40 pm
Location: Ranold AS

ROMS BLOWUP MaxDensity

#1 Unread post by 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
Posts: 71
Joined: Fri Nov 19, 2010 2:33 pm
Location: University of Aegean

Re: ROMS BLOWUP MaxDensity

#2 Unread post by 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
Posts: 9
Joined: Mon Jun 21, 2021 12:40 pm
Location: Ranold AS

Re: ROMS BLOWUP MaxDensity

#3 Unread post by 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
Posts: 1
Joined: Fri Apr 29, 2022 1:39 am
Location: IIT DELHI

Re: ROMS BLOWUP MaxDensity

#4 Unread post by 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

Post Reply