big barotropic Courant Number

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
pavel_fayman
Posts: 7
Joined: Thu Feb 24, 2011 2:40 pm
Location: FERHRI

big barotropic Courant Number

#1 Unread post by pavel_fayman »

log file shows

Minimum X-grid spacing, DXmin = 1.73851754E+00 km
Maximum X-grid spacing, DXmax = 2.22556255E+00 km
Minimum Y-grid spacing, DYmin = 2.00000000E+00 km
Maximum Y-grid spacing, DYmax = 2.00000000E+00 km
Minimum Z-grid spacing, DZmin = 3.75000000E-01 m
Maximum Z-grid spacing, DZmax = 3.91957187E+02 m

Minimum barotropic Courant Number = 1.00000000E+20
Maximum barotropic Courant Number = -1.00000000E+20
Maximum Coriolis Courant Number = -1.00000000E+20



1.00000000E+20 - there is abnormally large number

output data of temperature, salinity and u,v-components are 1e+37

what is the reason?

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: big barotropic Courant Number

#2 Unread post by kate »

pavel_fayman wrote: output data of temperature, salinity and u,v-components are 1e+37

what is the reason?
Is this for the initial conditions? There's something wrong and you will have to investigate it further. Where do the initial conditions come from? File or ana? If a file, have you looked at it with ncview? Are the fields in bounds there?

pavel_fayman
Posts: 7
Joined: Thu Feb 24, 2011 2:40 pm
Location: FERHRI

Re: big barotropic Courant Number

#3 Unread post by pavel_fayman »

no problem with initial and boundary conditions - all ini and bound data take from netcdf files
next massage after Courant numbers in log files:


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(File: initial_2008_10_2.nc, Rec=0001, Index=1)
- free-surface
(Min = -1.63539832E-02 Max = 3.99238263E-01)
- vertically integrated u-momentum component
(Min = -1.37729413E-01 Max = 1.66987238E-01)
- vertically integrated v-momentum component
(Min = -1.52762200E-01 Max = 2.24287793E-01)
- u-momentum component
(Min = -4.30719233E-01 Max = 4.05836069E-01)
- v-momentum component
(Min = -5.65662029E-01 Max = 4.30758912E-01)
- potential temperature
(Min = 0.00000000E+00 Max = 2.06500482E+01)
- salinity
(Min = 0.00000000E+00 Max = 4.80799970E+01)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0001, Index=1, File: force_2008_10.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -5.52289494E+00 Max = 5.35037593E+00)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0001, Index=1, File: force_2008_10.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -6.43010344E+00 Max = 7.72222479E+00)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0001, Index=1, File: force_2008_10.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 8.15493645E+02 Max = 1.01727652E+03)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:00:00
(Rec=0001, Index=1, File: force_2008_10.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 6.34736878E-06 Max = 7.61811690E-05)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0001, Index=1, File: force_2008_10.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 1.05328449E+01 Max = 2.03014051E+01)
GET_2DFLD - surface air relative humidity, t = 0 00:00:00
(Rec=0001, Index=1, File: force_2008_10.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 3.82153687E-01 Max = 9.95000000E-01)
GET_NGFLD - free-surface western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.63539832E-02 Max = 3.84787902E-01)
GET_NGFLD - free-surface southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 0.00000000E+00 Max = 3.80767905E-01)
GET_NGFLD - free-surface northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.30527306E-02 Max = 1.69657139E-01)
GET_NGFLD - 2D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.88508456E-02 Max = 2.35897268E-02)
GET_NGFLD - 2D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.14395532E-02 Max = 1.67580893E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -2.94021266E-02 Max = 1.34843144E-02)
GET_NGFLD - 2D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -2.56481688E-02 Max = 6.53302891E-02)
GET_NGFLD - 2D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.37728168E-01 Max = 1.11989719E-02)
GET_NGFLD - 2D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -6.47755283E-03 Max = 6.94770140E-02)
GET_NGFLD - 3D u-momentum western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.49299064E-01 Max = 2.32547100E-01)
GET_NGFLD - 3D v-momentum western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -2.06744885E-01 Max = 2.79297984E-01)
GET_NGFLD - 3D u-momentum southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -1.71734042E-01 Max = 1.80543733E-01)
GET_NGFLD - 3D v-momentum southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -3.34541061E-01 Max = 1.40604646E-01)
GET_NGFLD - 3D u-momentum northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -2.87807229E-01 Max = 1.24616294E-01)
GET_NGFLD - 3D v-momentum northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = -5.46407429E-01 Max = 1.85562150E-01)
GET_NGFLD - temperature western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 8.70118895E-01 Max = 2.02672830E+01)
GET_NGFLD - salinity western boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 3.22908820E+01 Max = 3.48023484E+01)
GET_NGFLD - temperature southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 0.00000000E+00 Max = 2.02672830E+01)
GET_NGFLD - salinity southern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 0.00000000E+00 Max = 3.47697329E+01)
GET_NGFLD - temperature northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 0.00000000E+00 Max = 1.46583184E+01)
GET_NGFLD - salinity northern boundary condition, t = 0 00:00:00
(Rec=0001, Index=2, File: bound_2008_10_2.nc)
(Tmin= 0.0000 Tmax= 42.0000)
(Min = 0.00000000E+00 Max = 3.47718879E+01)

Post Reply