Model blows up at 1st time step

[lyzheng]

Hi, All,

I am using latest version of ROMS to run a case with one boundary (south). The initial conditions are 0 for seal level and velocities and climatological data for TS. The external forces include winds, heat fluxes, and rivers. Along the open boundary, the boundary file provides the sea level and transports as 0 (since there is not data) and TS as climatological data (constant during simulation period). The south edge open boundary conditions are SOUTH_FSCHAPMAN for sea level, SOUTH_M2FLATHER for transports, SOUTH_M3RADIATION for 3D velocities, and SOUTH_TRADIATION and SOUTH_TNUDGING for temperature and salinity.

When I run the model, it bows up at the 1st time step. The output file shows following message:

I STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME

0 41225 00:00:00 0.000000E+00 1.260187E+04 1.260187E+04 4.315980E+14
DEF_HIS - creating history file: WFS_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file: WFS_avr.nc
DEF_DIAGS - creating diagnostics file: WFS_dia.nc
DEF_STATION - creating stations file: WFS_sta.nc
DEF_FLOATS - creating floats file: WFS_flt.nc
1 41225 00:02:00 6.129442E-08 NaN NaN 4.315980E+14
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):


I check the input data from the output file and they are correct as input:

NLM: GET_STATE - Read state initial conditions, t = 41225 00:00:00
(File: ic72.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 5.09499979E+00 Max = 3.02528725E+01)
- salinity
(Min = 3.48903580E+01 Max = 3.65991669E+01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 4.00000000E+00 Max = 1.16000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 1.20000000E+01 Max = 7.50000000E+01)
GET_NGFLD - river runoff direction
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 0.00000000E+00 Max = 2.00000000E-01)
GET_NGFLD - river runoff mass transport, t = 289 03:45:00
(Rec=0010, Index=2, File: river_clim2.nc)
(Tmin= 15.2188 Tmax= 350.0312)
(Min = 5.29000000E+01 Max = 8.80840000E+03)
GET_2DFLD - surface u-wind component, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -3.98408008E+00 Max = 1.93082631E+00)
GET_2DFLD - surface v-wind component, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -3.10085535E+00 Max = 1.77354598E+00)
GET_2DFLD - surface air pressure, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 1.00793756E+03 Max = 1.01813959E+03)
GET_2DFLD - cloud fraction, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 7.00000003E-02 Max = 7.00000003E-02)
GET_2DFLD - solar shortwave radiation flux, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 6.84567836E-12 Max = 6.84567836E-12)
GET_2DFLD - surface air temperature, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 2.06800423E+01 Max = 2.76649666E+01)
GET_2DFLD - surface air relative humidity, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 5.63671708E-01 Max = 9.72829878E-01)
GET_2DFLD - rain fall rate, t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 2.65815161E-05)
GET_2DFLD - sea surface temperature climatology, t = 41225 00:00:00
(Rec=0001, Index=1, File: sst72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 2.38468552E+01 Max = 2.87542210E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 41225 00:00:00
(Rec=0001, Index=1, File: sst72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -1.15740752E-05 Max = -1.15740752E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 41225 00:00:00
(Rec=0001, Index=1, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -8.67638432E-08 Max = -1.15116658E-10)
GET_NGFLD - free-surface southern boundary condition, t = 41225 00:00:00
(Rec=0001, Index=2, File: bc_nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 41225 00:00:00
(Rec=0001, Index=2, File: bc_nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 41225 00:00:00
(Rec=0001, Index=2, File: bc_nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_3DFLD - potential temperature climatology, t = 41225 00:00:00
(Rec=0001, Index=1, File: nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 5.09499979E+00 Max = 3.02528725E+01)
GET_3DFLD - salinity climatology, t = 41225 00:00:00
(Rec=0001, Index=1, File: nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 3.48903580E+01 Max = 3.65991669E+01)

Maximum grid stiffness ratios: rx0 = 0.000000E+00 (Beckmann and Haidvogel)
rx1 = 0.000000E+00 (Haney)


Initial basin volumes: TotVolume = 0.0000000000E+00 m3
MinVolume = 1.0000000000E+20 m3
MaxVolume = 0.0000000000E+00 m3
Max/Min = 0.0000000000E+00

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00263520)
GET_NGFLD - river runoff mass transport, t = 319 14:15:00
(Rec=0011, Index=1, File: river_clim2.nc)
(Tmin= 15.2188 Tmax= 350.0312)
(Min = 4.30000000E+00 Max = 1.06833000E+04)
GET_2DFLD - surface u-wind component, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -5.01726055E+00 Max = 2.08626461E+00)
GET_2DFLD - surface v-wind component, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -1.03191841E+00 Max = 2.52696967E+00)
GET_2DFLD - surface air pressure, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 1.00754950E+03 Max = 1.01946692E+03)
GET_2DFLD - cloud fraction, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 7.00000003E-02 Max = 8.00000012E-01)
GET_2DFLD - solar shortwave radiation flux, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 6.84567836E-12 Max = 7.64232168E-07)
GET_2DFLD - surface air temperature, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 1.90295238E+01 Max = 2.71164627E+01)
GET_2DFLD - surface air relative humidity, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 6.26445651E-01 Max = 9.87701476E-01)
GET_2DFLD - rain fall rate, t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 6.98466320E-05)
GET_2DFLD - sea surface temperature climatology, t = 41225 06:00:00
(Rec=0002, Index=2, File: sst72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 2.20039940E+01 Max = 2.87542210E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 41225 06:00:00
(Rec=0002, Index=2, File: sst72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -1.15740752E-05 Max = -1.15740752E-05)
GET_2DFLD - surface net freswater flux, (E-P), t = 41225 06:00:00
(Rec=0002, Index=2, File: atm72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = -1.21943334E-07 Max = 2.19937052E-09)
GET_NGFLD - free-surface southern boundary condition, t = 41591 00:00:00
(Rec=0002, Index=1, File: bc_nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 41591 00:00:00
(Rec=0002, Index=1, File: bc_nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 41591 00:00:00
(Rec=0002, Index=1, File: bc_nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_3DFLD - potential temperature climatology, t = 41591 00:00:00
(Rec=0002, Index=2, File: nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 5.09499979E+00 Max = 3.02528725E+01)
GET_3DFLD - salinity climatology, t = 41591 00:00:00
(Rec=0002, Index=2, File: nest72.nc)
(Tmin= 41225.0000 Tmax= 41591.0000)
(Min = 3.48903580E+01 Max = 3.65991669E+01)

It is very likely caused by incorrect model setup. Does someone has the same thing happen before? Any comments are greatly appreciated.

Zheng

[alazhar]

I see something weird in your basin volumes report where your max volume is zero.

Initial basin volumes: TotVolume = 0.0000000000E+00 m3
MinVolume = 1.0000000000E+20 m3
MaxVolume = 0.0000000000E+00 m3
Max/Min = 0.0000000000E+00

maybe you need to recheck your grid setup.

[lyzheng]

Thanks, Alazhar,

The grid input file can be correctly applied to ROMS version 3.0. I run version 3.0 using the same grid input file and it gets:

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


Initial basin volumes: TotVolume = 4.30769841105073E+14 m3
MinVolume = 1.08564579831466E+06 m3
MaxVolume = 2.94063866737793E+10 m3
Max/Min = 2.70865384634926E+04

but if I run version 3.8, it gets:

Maximum grid stiffness ratios: rx0 = 0.000000E+00 (Beckmann and Haidvogel)
rx1 = 0.000000E+00 (Haney)


Initial basin volumes: TotVolume = 0.0000000000E+00 m3
MinVolume = 1.0000000000E+20 m3
MaxVolume = 0.0000000000E+00 m3
Max/Min = 0.0000000000E+00

Does this mean that for versions 3.0 and 3.8, the grid input files are different?

Thanks.

Zheng