Long time hindcast instability

Information about ongoing ROMS/TOMS applications

Moderators: arango, robertson

Post Reply
Message
Author
LuoZhengpu
Posts: 4
Joined: Tue Aug 02, 2022 2:23 pm
Location: Ocean University of China

Long time hindcast instability

#1 Unread post by LuoZhengpu »

Using Hourly CFSR atmospheric bulk forcing to drive ROMS hindcast from 2004 to 2020, a sudden occurrence of NAN occurred at 2007 during the calculation process. Through experiments, it was found tthat it was a problem with the wind file (Replacing with ERA5 wind can run normally), but there was no NAN in the CFSR wind file, indicating no outlier data. What could be the reason for the instability of the model? :cry:
Is it because the Vwind speed is too high?But in fact, there seemed to be a higher speed before, but it ran smoothly.
By the way, the time step is 60s. Model reselution is about 8km.

Here is the output before it blowed up:
1667880 2007-03-04 06:00:00.00 1.538431E-02 2.705588E+02 2.705742E+02 1.765453E+13
(104,095,32) 5.740918E-04 1.071267E-03 1.655787E-01 9.994749E-01
1667890 2007-03-04 06:10:00.00 1.546097E-02 2.704776E+02 2.704931E+02 1.765144E+13
(104,095,32) 4.842636E-04 1.076456E-03 1.607030E-01 1.002374E+00
1667900 2007-03-04 06:20:00.00 1.557903E-02 2.703980E+02 2.704136E+02 1.764826E+13
(063,117,30) 4.022391E-04 2.680786E-03 1.561503E-01 1.003446E+00
1667910 2007-03-04 06:30:00.00 1.572910E-02 2.703206E+02 2.703363E+02 1.764501E+13
(063,117,30) 3.782827E-04 2.698790E-03 1.563011E-01 1.003125E+00
1667920 2007-03-04 06:40:00.00 1.590136E-02 2.702458E+02 2.702617E+02 1.764172E+13
(063,117,30) 3.561241E-04 2.717476E-03 1.566878E-01 1.001529E+00
1667930 2007-03-04 06:50:00.00 1.608707E-02 2.701742E+02 2.701903E+02 1.763841E+13
(063,117,30) 3.375654E-04 2.735594E-03 1.573522E-01 9.987905E-01
GET_2DFLD_NF90 - surface u-wind component, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_wind_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = -1.37525229E+01 Max = 1.46642400E+01) regrid = T
GET_2DFLD_NF90 - surface v-wind component, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_wind_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = -2.21711187E+01 Max = 2.04555164E+01) regrid = T
GET_2DFLD_NF90 - surface air pressure, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_surface_air_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = 9.90562060E+02 Max = 1.01734676E+03) regrid = T
GET_2DFLD_NF90 - solar shortwave radiation flux, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_flux_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = 2.51117417E-06 Max = 6.22807026E-05) regrid = T
GET_2DFLD_NF90 - net longwave radiation flux, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_flux_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = 6.90479480E-05 Max = 9.74249733E-05) regrid = T
GET_2DFLD_NF90 - surface air temperature, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_surface_air_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = -5.02027887E+00 Max = 1.77817922E+01) regrid = T
GET_2DFLD_NF90 - surface air relative humidity, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_surface_air_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = 7.61973396E-01 Max = 9.86484414E-01) regrid = T
GET_2DFLD_NF90 - rain fall rate, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_bulk_flux_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = 1.33951269E-06 Max = 2.79875110E-03) regrid = T
GET_2DFLD_NF90 - sea surface temperature climatology, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_frc_sst_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = -4.48853350E+00 Max = 1.56885500E+01) regrid = F
GET_2DFLD_NF90 - surface net heat flux sensitivity to SST, 2007-03-04 08:00:00.00
(Grid=01, Rec=1497, Index=2, File: BYECS_HINDCAST_frc_sst_hourly_2007_CFSR.nc)
(Tmin= 1096.0000 Tmax= 1460.9583) t = 1158.3333
(Min = -1.82112137E-05 Max = -2.64627197E-06) regrid = F
1667940 2007-03-04 07:00:00.00 1.627521E-02 2.701063E+02 2.701226E+02 1.763511E+13
(063,117,30) 3.238365E-04 2.752176E-03 1.582398E-01 9.947457E-01
1667950 2007-03-04 07:10:00.00 1.646089E-02 2.700424E+02 2.700589E+02 1.763183E+13
(063,117,29) 3.044345E-04 2.635372E-03 1.611834E-01 9.909538E-01
1667960 2007-03-04 07:20:00.00 1.663722E-02 2.699829E+02 2.699995E+02 1.762860E+13
(063,117,29) 3.471054E-04 2.666293E-03 1.672286E-01 9.874487E-01
1667970 2007-03-04 07:30:00.00 NaN NaN NaN NaN
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 NaN
Found Error: 01 Line: 321 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 297 Source: ROMS/Drivers/nl_roms.h, ROMS_run

Blowing-up: Saving latest model state into RESTART file
REASON: KEchar = NaN, PEchar = NaN

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


Looking forward to your answer! :D

Post Reply