issue to handle masked surface forcing files, model is blown

[tara]

Hello Roms community,

I am using the tropflux daily climatology data forcings (Tair,swrad, lwrad, Qair, Pair from coreV2) and GPCP data for rainfall rate and Uwind & Vwind from quikscat data. The current setup is a trial to run a 10 years climatological run for the domain [35E to 90E & 25S to 20N] and resolution is 1/4. But model is blowing in the second time step itself. description of log file is under mentioned.


ndtfast, nfast = 50 69 nfast/ndtfast = 1.38000

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

1.000000000000 1.033396903681 0.516698451840 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.22720

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

Minimum X-grid spacing, DXmin = 2.51801363E+01 km
Maximum X-grid spacing, DXmax = 2.77831716E+01 km
Minimum Y-grid spacing, DYmin = 2.52032577E+01 km
Maximum Y-grid spacing, DYmax = 2.77831405E+01 km
Minimum Z-grid spacing, DZmin = 1.25073687E+00 m
Maximum Z-grid spacing, DZmax = 7.52012906E+02 m

Minimum barotropic Courant Number = 3.40694121E-02
Maximum barotropic Courant Number = 2.98301057E-01
Maximum Coriolis Courant Number = 7.37608015E-02

Horizontal mixing scaled by grid size, GRDMAX = 2.77831249E+01 km

Minimum horizontal diffusion coefficient = 2.09930786E+02 m2/s
Maximum horizontal diffusion coefficient = 2.31526041E+02 m2/s

Minimum horizontal viscosity coefficient = 2.09930786E+02 m2/s
Maximum horizontal viscosity coefficient = 2.31526041E+02 m2/s

NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: roms_ini.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 = 6.12046405E-01 Max = 2.96272088E+01)
- salinity
(Min = 3.07662600E+01 Max = 4.06677452E+01)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: uwind.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -8.86058779E+00 Max = 1.00000000E+20)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: vwind.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -8.67955616E+00 Max = 1.00000000E+20)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0000001, Index=1, File: surface_pres.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 1.00951002E+03 Max = 1.02011134E+03)
GET_2DFLD - solar shortwave radiation flux, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.97214003E-05 Max = 7.47180066E-05)
GET_2DFLD - net longwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: longwave.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -2.81343400E-05 Max = 2.44820516E+13)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0000001, Index=1, File: airtemp.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 2.06291232E+01 Max = 1.00000000E+20)
GET_2DFLD - surface air specific humidity, t = 0 00:00:00
(Rec=0000001, Index=1, File: specific_hum.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 9.47512199E+00 Max = 1.00000000E+20)
GET_2DFLD - rain fall rate, t = 0 00:00:00
(Rec=0000001, Index=1, File: precip.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 0.00000000E+00 Max = 3.03253871E-04)
GET_2DFLD - sea surface salinity climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.75694853E+01 Max = 3.96281540E+01)
GET_NGFLD - free-surface eastern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.71301593E-01 Max = 1.30383335E-01)
GET_NGFLD - free-surface southern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.71301593E-01 Max = 8.02520623E+00)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.30252428E-01 Max = 1.18584624E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.45121988E-02 Max = 1.28864014E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.11475038E-02 Max = 7.97645937E-03)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.93501570E-01 Max = 7.76984871E-02)
GET_NGFLD - temperature eastern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 7.27592240E-01 Max = 2.86222639E+01)
GET_NGFLD - salinity eastern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.02278060E+01 Max = 3.58199364E+01)
GET_NGFLD - temperature southern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 6.13947362E-01 Max = 2.68585280E+01)
GET_NGFLD - salinity southern boundary condition, t = 345 00:00:00
(Rec=0000012, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.44397525E+01 Max = 3.58216616E+01)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 2.154182E-01 (Beckmann and Haidvogel)
rx1 = 4.478907E+00 (Haney)

Initial basin volumes: TotVolume = 8.4646060247E+16 m3
MinVolume = 8.5321855646E+08 m3
MaxVolume = 5.8048028930E+11 m3
Max/Min = 6.8034184782E+02

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 08064720)

GET_2DFLD - surface u-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: uwind.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -9.42937975E+00 Max = 1.00000000E+20)
GET_2DFLD - surface v-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: vwind.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -8.49320788E+00 Max = 1.00000000E+20)
GET_2DFLD - surface air pressure, t = 1 00:00:00
(Rec=0000002, Index=2, File: surface_pres.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 1.00934552E+03 Max = 1.01951950E+03)
GET_2DFLD - solar shortwave radiation flux, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 4.14254984E-05 Max = 7.63300206E-05)
GET_2DFLD - net longwave radiation flux, t = 1 00:00:00
(Rec=0000002, Index=2, File: longwave.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -2.80249039E-05 Max = 2.44820516E+13)
GET_2DFLD - surface air temperature, t = 1 00:00:00
(Rec=0000002, Index=2, File: airtemp.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 2.04292793E+01 Max = 1.00000000E+20)
GET_2DFLD - surface air specific humidity, t = 1 00:00:00
(Rec=0000002, Index=2, File: specific_hum.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 9.55675877E+00 Max = 1.00000000E+20)
GET_2DFLD - rain fall rate, t = 1 00:00:00
(Rec=0000002, Index=2, File: precip.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = 0.00000000E+00 Max = 2.52505480E-04)
GET_2DFLD - sea surface salinity climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.87745855E+01 Max = 3.86804531E+01)
GET_NGFLD - free-surface eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.49242404E-01 Max = 1.21437297E-01)
GET_NGFLD - free-surface southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.35165429E-01 Max = 8.02420879E+00)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.33005595E-01 Max = 1.00078927E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.78489296E-02 Max = 1.75680212E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.46073857E-02 Max = 7.39300480E-03)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.06468452E-01 Max = 6.61020167E-02)
GET_NGFLD - temperature eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 7.27592240E-01 Max = 2.88058121E+01)
GET_NGFLD - salinity eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.07662600E+01 Max = 3.58525282E+01)
GET_NGFLD - temperature southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 6.13947362E-01 Max = 2.75359423E+01)
GET_NGFLD - salinity southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.44510161E+01 Max = 3.58558299E+01)

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

0 0 00:00:00 0.000000E+00 1.993934E+04 1.993934E+04 8.519198E+16
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: output/ocean_his_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
1 0 00:20:00 NaN NaN NaN NaN
(219,185,01) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=2,2) into time record = 0000001

Elapsed CPU time (seconds):

Node # 18 CPU: 3.657
Node # 22 CPU: 3.643
Node # 26 CPU: 3.648
Node # 30 CPU: 3.675
Node # 34 CPU: 3.645
Node # 0 CPU: 2.426
Node # 6 CPU: 3.672
Node # 24 CPU: 3.667
Node # 2 CPU: 3.656
Node # 4 CPU: 3.676
Node # 10 CPU: 3.670
Node # 14 CPU: 3.671
Node # 16 CPU: 3.668
Node # 20 CPU: 3.632
Node # 28 CPU: 3.637

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

Analytical header files used:

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


ROMS/TOMS: DONE... Friday - July 24, 2015 - 9:04:10 PM
Node # 11 CPU: 3.665
Node # 15 CPU: 3.676
Node # 19 CPU: 3.644
Node # 23 CPU: 3.647
Node # 27 CPU: 3.606
Node # 31 CPU: 3.622
Node # 7 CPU: 3.651
Node # 17 CPU: 3.668
Node # 21 CPU: 3.656
Node # 25 CPU: 3.644
Node # 29 CPU: 3.678
Node # 33 CPU: 3.614
Node # 35 CPU: 3.623
Node # 1 CPU: 3.661
Node # 3 CPU: 3.664
Node # 5 CPU: 3.646
Node # 9 CPU: 3.669
Node # 13 CPU: 3.672
Node # 8 CPU: 3.652
Node # 12 CPU: 3.679
Node # 32 CPU: 3.637

please suggest some probable solution.


Best Regards-
Tara

[kate]

There are many ways to get ROMS to blow up that quickly. I would check the forcing files first, though:

GET_2DFLD - surface u-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: uwind.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -9.42937975E+00 Max = 1.00000000E+20)
GET_2DFLD - surface v-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: vwind.nc)
(Tmin= 0.0000 Tmax= 364.0000)
(Min = -8.49320788E+00 Max = 1.00000000E+20)

What is that 10E+20? A special land value?

[tara]

Hello Kate,

I am very grateful for your instant reply, yes you are correct, my tropflux data that I am using for Tair, Qair and quikscat(Uwind,Vwind) data both are masked over land , and my doubt was also on these masked values present in land points. But I am surprise that if model is not dealing with land points then how these values are creating problems and the question comes that ...

Q1> If I need to use this kind of masked land points data for forcing, how to process the data or provide to the model in a suitable manner, so that model will not blow?

Actually, I tried to solve this issue by providing some different domains so that these points will not come in the calculation (Avoiding NaN i.e. in the first time step). I tried this different domain setups & all of them was over Indian Ocean and less than the spatial range what my forcing files have.

Q2> can you kindly suggest me whether this issue can cause for blown up,if yes then how to make this data compatible with model?


With Best Regards--
Tara

[kate]

I typically use MERRA, which has values everywhere. However, the land values are appropriate to say winds over land, not winds over water. My colleague Raphael has introduced us to a tool called sosie which we use to process the MERRA files, extrapolating water values into the land. This allows ROMS to do its interpolation to the ROMS ocean points using valid "wet" values even when the nearest wind point is actually a land point.

Another option would be to look into the ROMS interpolation code and perhaps replace it with something smarter about masked values. I've been hearing of groups switching over to ESMF now that it provides interpolation routines. I'd be interested in your mods if you take this route.

[tara]

Dear Kate,

Thank you very much for your kind help.

[tara]

We resolved this issue by adding attribute "_FillValue = 1.e+20;" to the uwind, vwind variables in forcings files. By default it was having missing_value attribute.
I just followed the same procedure for all my required forcings files.
Now model is working fine.

I got this _FillValue attribute information from here viewtopic.php?f=14&t=1078

Thanks & Regards,
Tara