running error

[langlangxie]

when I are running ROMS,I run into a problem.follows:
Basin information for Grid 01:

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

Initial basin volumes: TotVolume = 6.3267250293E+15 m3
MinVolume = 8.3358286833E+08 m3
MaxVolume = 3.2792063556E+11 m3
Max/Min = 3.9338696609E+02

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

GET_2DFLD - surface u-momentum stress, t = 15 00:00:00
(Rec=0000001, Index=2, File: scs_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.41920103E-04 Max = -5.03152173E-06)
GET_2DFLD - surface v-momentum stress, t = 15 00:00:00
(Rec=0000001, Index=2, File: scs_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.37492202E-04 Max = -1.75192133E-08)
GET_2DFLD - solar shortwave radiation flux, t = 15 00:00:00
(Rec=0000001, Index=2, File: scs_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.31937231E-05 Max = 4.81916589E-05)
GET_2DFLD - surface net heat flux, t = 15 00:00:00
(Rec=0000001, Index=2, File: scs_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -9.77548439E-05 Max = 9.37606950E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 15 00:00:00
(Rec=0000001, Index=2, File: scs_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.18638232E-07 Max = 8.56997888E-08)
GET_NGFLD - free-surface western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - free-surface eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -7.97567557E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -7.97724727E+00 Max = 0.00000000E+00)
GET_NGFLD - free-surface northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.18108488E+00 Max = 0.00000000E+00)
GET_NGFLD - 2D u-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - 2D v-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.20794168E-01 Max = 2.39319701E-01)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 1.56021449E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.15372022E-01 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.45559910E-01 Max = 3.46445228E-01)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.43125938E-02 Max = 0.00000000E+00)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.75982703E-01 Max = 1.91234085E-01)
GET_NGFLD - 3D u-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - 3D v-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - 3D u-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.94580908E-01 Max = 3.88307495E-01)
GET_NGFLD - 3D v-momentum eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 6.09734925E-02)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.93077978E-01 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.93556123E-01 Max = 1.16396825E+00)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.73267105E-02 Max = 0.00000000E+00)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.13179043E-01 Max = 4.69278914E-01)
GET_NGFLD - temperature western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.58405853E+01 Max = 2.87204287E+01)
GET_NGFLD - salinity western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.99974794E+01 Max = 3.33290744E+01)
GET_NGFLD - temperature eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.18097712E+00 Max = 2.87141523E+01)
GET_NGFLD - salinity eastern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.37087321E+01 Max = 3.49626732E+01)
GET_NGFLD - temperature southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.65942121E+00 Max = 2.89971042E+01)
GET_NGFLD - salinity southern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.12381541E+01 Max = 3.47662936E+01)
GET_NGFLD - temperature northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 7.43037328E+00 Max = 2.24307879E+01)
GET_NGFLD - salinity northern boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.06074367E+01 Max = 3.47987568E+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 NaN 1.491859E+04 NaN 6.499725E+15
(89,001,25) 0.000000E+00 1.739733E-05 0.000000E+00 1.137579E-03

Elapsed CPU time (seconds):
I think my error is Cu and Cw is NaN,and I find some max to be NaN in my scs_bry.nc,is it making a error?
My borders are all open borders,and I always can't run next step.

[kate]

This is worrisome:

Code: Select all

GET_NGFLD - 2D u-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - 2D v-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)

The NaNs here will cause trouble, as you found. How did you create them?

[langlangxie]

This is worrisome:

Code: Select all

GET_NGFLD - 2D u-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)
GET_NGFLD - 2D v-momentum western boundary condition, t = 15 00:00:00
(Rec=0000001, Index=1, File: scs_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = NaN Max = NaN)

The NaNs here will cause trouble, as you found. How did you create them?

Thanks!
I use agrifroms to make Climate-state data files.I choose coads05 to make forcing filed,then choose woa2009 to make initial field and boundary field.what's more ,my area is 103°E-126°E,3°N-27°N,so I choose four directions to be open border.Then it create the NaNs.I think I can choose the way of interpolation to solve the problem,do you think ok? Can you tell me some other ways?

[kate]

For my newest domain, I used python scripts (which call scrip) to get initial and boundary conditions from a global HYCOM simulation. I think the model prefers boundary conditions from another model compared to something like a Levitus field.