all result = 0?

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

all result = 0?

#1 Unread post by hpd14thu »

Hello ,every one,I just use a grid file, a initial file and a easy wind file, in order to calculate the flow forced by the wind. But the result( u,v ,salt...) is zero. I tried a lot ,but it didn't work. I really need your help, thank you!

Here is my in.file and headfile

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#2 Unread post by hpd14thu »

! Application title.

TITLE = ZFIRST

! C-preprocessing Flag.

MyAppCPP = ZFIRST

! Input variable information file name. This file needs to be processed
! first so all information arrays can be initialized properly.

VARNAME = /home/hpd14/Documents/COAWST_V3.2/ROMS/External/varinfo.dat

! Number of nested grids.

Ngrids = 1

! Number of grid nesting layers. This parameter is used to allow refinement
! and composite grid combinations.

NestLayers = 1

! Number of grids in each nesting layer [1:NestLayers].

GridsInLayer = 1

! Grid dimension parameters. See notes below in the Glossary for how to set
! these parameters correctly.

Lm == 48 ! Number of I-direction INTERIOR RHO-points
Mm == 48 ! Number of J-direction INTERIOR RHO-points
N == 20 ! Number of vertical levels
ND == 0 ! Number of wave directional bins

Nbed = 0 ! Number of sediment bed layers

NAT = 2 ! Number of active tracers (usually, 2)
NPT = 0 ! Number of inactive passive tracers
NCS = 0 ! Number of cohesive (mud) sediment tracers
NNS = 0 ! Number of non-cohesive (sand) sediment tracers

! Domain decomposition parameters for serial, distributed-memory or
! shared-memory configurations used to determine tile horizontal range
! indices (Istr,Iend) and (Jstr,Jend), [1:Ngrids].

NtileI == 1 ! I-direction partition
NtileJ == 1 ! J-direction partition

! Set lateral boundary conditions keyword. Notice that a value is expected
! for each boundary segment per nested grid for each state variable.
!
! Each tracer variable requires [1:4,1:NAT+NPT,Ngrids] values. Otherwise,
! [1:4,1:Ngrids] values are expected for other variables. The boundary
! order is: 1=west, 2=south, 3=east, and 4=north. That is, anticlockwise
! starting at the western boundary.
!
! The keyword is case insensitive and usually has three characters. However,
! it is possible to have compound keywords, if applicable. For example, the
! keyword "RadNud" implies radiation boundary condition with nudging. This
! combination is usually used in active/passive radiation conditions.
!
! Keyword Lateral Boundary Condition Type
!
! Cha Chapman_implicit (free-surface)
! Che Chapman_explicit (free-surface)
! Cla Clamped
! Clo Closed
! Fla Flather (2D momentum) _____N_____ j=Mm
! Gra Gradient | 4 |
! Nes Nested (refinement) | |
! Nud Nudging 1 W E 3
! Per Periodic | |
! Rad Radiation |_____S_____|
! Red Reduced Physics (2D momentum) 2 j=1
! Shc Shchepetkin (2D momentum) i=1 i=Lm
!
! W S E N
! e o a o
! s u s r
! t t t t
! h h
!
! 1 2 3 4

LBC(isFsur) == Clo Clo Clo Clo ! free-surface
LBC(isUbar) == Clo Clo Clo Clo ! 2D U-momentum
LBC(isVbar) == Clo Clo Clo Clo ! 2D V-momentum
LBC(isUvel) == Clo Clo Clo Clo ! 3D U-momentum
LBC(isVvel) == Clo Clo Clo Clo ! 3D V-momentum
LBC(isMtke) == Clo Clo Clo Clo ! mixing TKE

LBC(isTvar) == Clo Clo Clo Clo \ ! temperature
Clo Clo Clo Clo ! salinity

! Wec boundary conditions

LBC(isU2Sd) == Clo Clo Clo Clo ! 2D U-stokes
LBC(isV2Sd) == Clo Clo Clo Clo ! 2D V-stokes
LBC(isU3Sd) == Clo Clo Clo Clo ! 3D U-stokes
LBC(isV3Sd) == Clo Clo Clo Clo ! 3D V-stokes

! Ice boundary conditions

LBC(isAice) == Clo Clo Clo Clo ! ice concentration
LBC(isHice) == Clo Clo Clo Clo ! ice thickness
LBC(isHsno) == Clo Clo Clo Clo ! snow thickness
LBC(isTice) == Clo Clo Clo Clo ! ice temperature
LBC(isSfwat)== Clo Clo Clo Clo ! surface water
LBC(isSig11)== Clo Clo Clo Clo ! sigma-11
LBC(isSig12)== Clo Clo Clo Clo ! sigma-12
LBC(isSig22)== Clo Clo Clo Clo ! sigma-22
LBC(isUice) == Clo Clo Clo Clo ! ice U-momentum
LBC(isVice) == Clo Clo Clo Clo ! ice V-momentum

! Adjoint-based algorithms can have different lateral boundary
! conditions keywords.

ad_LBC(isFsur) == Clo Clo Clo Clo ! free-surface
ad_LBC(isUbar) == Clo Clo Clo Clo ! 2D U-momentum
ad_LBC(isVbar) == Clo Clo Clo Clo ! 2D U-momentum
ad_LBC(isUvel) == Clo Clo Clo Clo ! 3D U-momentum
ad_LBC(isVvel) == Clo Clo Clo Clo ! 3D V-momentum
ad_LBC(isMtke) == Clo Clo Clo Clo ! mixing TKE

ad_LBC(isTvar) == Clo Clo Clo Clo \ ! temperature
Clo Clo Clo Clo ! salinity

! Set lateral open boundary edge volume conservation switch for
! nonlinear model and adjoint-based algorithms. Usually activated
! with radiation boundary conditions to enforce global mass
! conservation, except if tidal forcing is enabled. [1:Ngrids].

VolCons(west) == F ! western boundary
VolCons(east) == F ! eastern boundary
VolCons(south) == F ! southern boundary
VolCons(north) == F ! northern boundary

ad_VolCons(west) == F ! western boundary
ad_VolCons(east) == F ! eastern boundary
ad_VolCons(south) == F ! southern boundary
ad_VolCons(north) == F ! northern boundary

! Time-Stepping parameters.

NTIMES == 8640
DT == 0.02d0
NDTFAST == 20

! Model iteration loops parameters.

ERstr = 1
ERend = 1
Nouter = 1
Ninner = 1
Nintervals = 1

! Number of eigenvalues (NEV) and eigenvectors (NCV) to compute for the
! Lanczos/Arnoldi problem in the Generalized Stability Theory (GST)
! analysis. NCV must be greater than NEV (see documentation below).

NEV = 2 ! Number of eigenvalues
NCV = 10 ! Number of eigenvectors

! Input/Output parameters.

NRREC == 0
LcycleRST == T
NRST == 120
NSTA == 1
NFLT == 1
NINFO == 1

! Output history, average, diagnostic files parameters.

LDEFOUT == T
NHIS == 1
NDEFHIS == 0
NTSAVG == 1
NAVG == 360
NDEFAVG == 0
NTSDIA == 1
NDIA == 360
NDEFDIA == 0

! Output tangent linear and adjoint models parameters.

LcycleTLM == F
NTLM == 720
NDEFTLM == 0
LcycleADJ == F
NADJ == 720
NDEFADJ == 0
NSFF == 720
NOBC == 720

! GST output and check pointing restart parameters.

LmultiGST = F ! one eigenvector per file
LrstGST = F ! GST restart switch
MaxIterGST = 500 ! maximum number of iterations
NGST = 10 ! check pointing interval

! Relative accuracy of the Ritz values computed in the GST analysis.

Ritz_tol = 1.0d-15

! Harmonic/biharmonic horizontal diffusion of tracer for nonlinear model
! and adjoint-based algorithms: [1:NAT+NPT,Ngrids].

TNU2 == 0.0d0 0.0d0 ! m2/s
TNU4 == 2.0d0 0.0d0 ! m4/s

ad_TNU2 == 0.0d0 0.0d0 ! m2/s
ad_TNU4 == 0.0d0 0.0d0 ! m4/s

! Harmonic/biharmonic, horizontal viscosity coefficient for nonlinear model
! and adjoint-based algorithms: [Ngrids].

VISC2 == 0.05d0 ! m2/s
VISC4 == 0.0d0 ! m4/s

ad_VISC2 == 0.0d0 ! m2/s
ad_VISC4 == 0.0d0 ! m4/s

! Logical switches (TRUE/FALSE) to increase/decrease horizontal viscosity
! and/or diffusivity in specific areas of the application domain (like
! sponge areas) for the desired application grid.

LuvSponge == F ! horizontal momentum
LtracerSponge == F F ! temperature, salinity, inert

! Vertical mixing coefficients for tracers in nonlinear model and
! basic state scale factor in adjoint-based algorithms: [1:NAT+NPT,Ngrids]

AKT_BAK == 5.0d-6 5.0d-6 5.0d-6 5.0d-6 ! m2/s

ad_AKT_fac == 1.0d0 1.0d0 ! nondimensional

! Vertical mixing coefficient for momentum for nonlinear model and
! basic state scale factor in adjoint-based algorithms: [Ngrids].

AKV_BAK == 5.0d-5 5.0d-5 ! m2/s

ad_AKV_fac == 1.0d0 ! nondimensional

! Turbulent closure parameters.

AKK_BAK == 5.0d-6 ! m2/s
AKP_BAK == 5.0d-6 ! m2/s
TKENU2 == 0.0d0 ! m2/s
TKENU4 == 0.0d0 ! m4/s

! Generic length-scale turbulence closure parameters.

GLS_P == 3.0d0 ! K-epsilon
GLS_M == 1.5d0
GLS_N == -1.0d0
GLS_Kmin == 7.6d-6
GLS_Pmin == 1.0d-12

GLS_CMU0 == 0.5477d0
GLS_C1 == 1.44d0
GLS_C2 == 1.92d0
GLS_C3M == -0.4d0
GLS_C3P == 1.0d0
GLS_SIGK == 1.0d0
GLS_SIGP == 1.30d0

! Constants used in surface turbulent kinetic energy flux computation.

CHARNOK_ALPHA == 1400.0d0 ! Charnok surface roughness
ZOS_HSIG_ALPHA == 0.5d0 ! roughness from wave amplitude
SZ_ALPHA == 0.25d0 ! roughness from wave dissipation
CRGBAN_CW == 100.0d0 ! Craig and Banner wave breaking
WEC_ALPHA == 0.0d0 ! 0: all wave dissip goes to break and none to roller.
! 1: all wave dissip goes to roller and none to breaking.

! Constants used in momentum stress computation.

RDRG == 3.0d-04 ! m/s
RDRG2 == 0.025d0 ! nondimensional
Zob == 0.02d0 ! m
Zos == 0.5d0 ! m

! Height (m) of atmospheric measurements for Bulk fluxes parameterization.

BLK_ZQ == 2.0d0 ! air humidity
BLK_ZT == 2.0d0 ! air temperature
BLK_ZW == 10.0d0 ! winds

! Minimum depth for wetting and drying.

DCRIT == 0.10d0 ! m

! Various parameters.

WTYPE == 1
LEVSFRC == 15
LEVBFRC == 1

! Set vertical, terrain-following coordinates transformation equation and
! stretching function (see below for details), [1:Ngrids].

Vtransform == 2 ! transformation equation
Vstretching == 4 ! stretching function

! Vertical S-coordinates parameters (see below for details), [1:Ngrids].

THETA_S == 8.0d0 ! surface stretching parameter
THETA_B == 0.5d0 ! bottom stretching parameter
TCLINE == 20.0d0 ! critical depth (m)

! Mean Density and Brunt-Vaisala frequency.

RHO0 = 1025.0d0 ! kg/m3
BVF_BAK = 1.0d-5 ! 1/s2

! Time-stamp assigned for model initialization, reference time
! origin for tidal forcing, and model reference time for output
! NetCDF units attribute.

DSTART = 0.0d0 ! days
TIDE_START = 0.0d0 ! days
TIME_REF = 0.0d0 ! yyyymmdd.dd

! Nudging/relaxation time scales, inverse scales will be computed
! internally, [1:Ngrids].

TNUDG == 0.0d0 0.0d0 ! days
ZNUDG == 0.0d0 ! days
M2NUDG == 0.0d0 ! days
M3NUDG == 0.0d0 ! days

! Factor between passive (outflow) and active (inflow) open boundary
! conditions, [1:Ngrids]. If OBCFAC > 1, nudging on inflow is stronger
! than on outflow (recommended).

OBCFAC == 0.0d0 ! nondimensional

! Linear equation of State parameters:

R0 == 1027.0d0 ! kg/m3
T0 == 10.0d0 ! Celsius
S0 == 30.0d0 ! nondimensional
TCOEF == 1.7d-4 ! 1/Celsius
SCOEF == 7.6d-4 ! nondimensional

! Slipperiness parameter: 1.0 (free slip) or -1.0 (no slip)

GAMMA2 == 1.0d0

! Logical switches (TRUE/FALSE) to activate horizontal momentum transport
! point Sources/Sinks (like river runoff transport) and mass point
! Sources/Sinks (like volume vertical influx), [1:Ngrids].

LuvSrc == F ! horizontal momentum transport
LwSrc == F ! volume vertical influx

! Logical switches (TRUE/FALSE) to activate tracers point Sources/Sinks
! (like river runoff) and to specify which tracer variables to consider:
! [1:NAT+NPT,Ngrids]. See glossary below for details.

LtracerSrc == F F ! temperature, salinity, inert

! Logical switches (TRUE/FALSE) to read and process climatology fields.
! See glossary below for details.

LsshCLM == F ! sea-surface height
Lm2CLM == F ! 2D momentum
Lm3CLM == F ! 3D momentum

LtracerCLM == F F ! temperature, salinity, inert

! Logical switches (TRUE/FALSE) to nudge the desired climatology field(s).
! If not analytical climatology fields, users need to turn ON the logical
! switches above to process the fields from the climatology NetCDF file
! that are needed for nudging. See glossary below for details.

LnudgeM2CLM == F ! 2D momentum
LnudgeM3CLM == F ! 3D momentum

LnudgeTCLM == F F ! temperature, salinity, inert

! Starting (DstrS) and ending (DendS) day for adjoint sensitivity forcing.
! DstrS must be less or equal to DendS. If both values are zero, their
! values are reset internally to the full range of the adjoint integration.

DstrS == 0.0d0 ! starting day
DendS == 0.0d0 ! ending day

! Starting and ending vertical levels of the 3D adjoint state variables
! whose sensitivity is required.

KstrS == 1 ! starting level
KendS == 1 ! ending level

! Logical switches (TRUE/FALSE) to specify the adjoint state variables
! whose sensitivity is required.

Lstate(isFsur) == F ! free-surface
Lstate(isUbar) == F ! 2D U-momentum
Lstate(isVbar) == F ! 2D V-momentum
Lstate(isUvel) == F ! 3D U-momentum
Lstate(isVvel) == F ! 3D V-momentum

Lstate(isTvar) == F F ! NT tracers

! Logical switches (TRUE/FALSE) to specify the state variables for
! which Forcing Singular Vectors or Stochastic Optimals is required.

Fstate(isFsur) == F ! free-surface
Fstate(isUbar) == F ! 2D U-momentum
Fstate(isVbar) == F ! 2D V-momentum
Fstate(isUvel) == F ! 3D U-momentum
Fstate(isVvel) == F ! 3D V-momentum
Fstate(isTvar) == F F ! NT tracers

Fstate(isUstr) == T ! surface U-stress
Fstate(isVstr) == T ! surface V-stress
Fstate(isTsur) == F F ! NT surface tracers flux

! Stochastic Optimals time decorrelation scale (days) assumed for
! red noise processes.

SO_decay == 2.0d0 ! days

! Stochastic Optimals surface forcing standard deviation for
! dimensionalization.

SO_sdev(isFsur) == 1.0d0 ! free-surface
SO_sdev(isUbar) == 1.0d0 ! 2D U-momentum
SO_sdev(isVbar) == 1.0d0 ! 2D V-momentum
SO_sdev(isUvel) == 1.0d0 ! 3D U-momentum
SO_sdev(isVvel) == 1.0d0 ! 3D V-momentum
SO_sdev(isTvar) == 1.0d0 1.0d0 ! NT tracers

SO_sdev(isUstr) == 1.0d0 ! surface U-stress
SO_sdev(isVstr) == 1.0d0 ! surface V-stress
SO_sdev(isTsur) == 1.0d0 1.0d0 ! NT surface tracers flux

! Logical switches (TRUE/FALSE) to activate writing of fields into
! HISTORY output file.

Hout(idUvel) == T ! u 3D U-velocity
Hout(idVvel) == T ! v 3D V-velocity
Hout(idu3dE) == F ! u_eastward 3D U-eastward at RHO-points
Hout(idv3dN) == F ! v_northward 3D V-northward at RHO-points
Hout(idWvel) == T ! w 3D W-velocity
Hout(idOvel) == T ! omega omega vertical velocity
Hout(idUbar) == T ! ubar 2D U-velocity
Hout(idVbar) == T ! vbar 2D V-velocity
Hout(idu2dE) == F ! ubar_eastward 2D U-eastward at RHO-points
Hout(idv2dN) == F ! vbar_northward 2D V-northward at RHO-points
Hout(idFsur) == T ! zeta free-surface
Hout(idBath) == T ! bath time-dependent bathymetry

Hout(idTvar) == T T ! temp, salt temperature and salinity

Hout(idUair) == T ! Uwind surface U-wind
Hout(idVair) == T ! Vwind surface V-wind
Hout(idUairE) == F ! Uwind_eastward surface U-wind
Hout(idVairN) == F ! Vwind_northward surface V-wind
Hout(idUsms) == T ! sustr surface U-stress
Hout(idVsms) == T ! svstr surface V-stress
Hout(idUbms) == T ! bustr bottom U-stress
Hout(idVbms) == T ! bvstr bottom V-stress

Hout(idUbrs) == F ! bustrc bottom U-current stress
Hout(idVbrs) == F ! bvstrc bottom V-current stress
Hout(idUbws) == F ! bustrw bottom U-wave stress
Hout(idVbws) == F ! bvstrw bottom V-wave stress
Hout(idUbcs) == F ! bustrcwmax bottom max wave-current U-stress
Hout(idVbcs) == F ! bvstrcwmax bottom max wave-current V-stress
Hout(idUVwc) == F ! bstrcwmax bottom max wave-current stress magnitude

Hout(idUbot) == F ! Ubot bed wave orbital U-velocity
Hout(idVbot) == F ! Vbot bed wave orbital V-velocity
Hout(idUbur) == F ! Ur bottom U-velocity above bed
Hout(idVbvr) == F ! Vr bottom V-velocity above bed

Hout(idW2xx) == F ! Sxx_bar WEC_Mellor 2D Sxx radiation stress
Hout(idW2xy) == F ! Sxy_bar WEC_Mellor 2D Sxy radiation stress
Hout(idW2yy) == F ! Syy_bar WEC_Mellor 2D Syy radiation stress
Hout(idW3xx) == F ! Sxx WEC_Mellor 3D Sxx radiation stress
Hout(idW3xy) == F ! Sxy WEC_Mellor 3D Sxy radiation stress
Hout(idW3yy) == F ! Syy WEC_Mellor 3D Syy radiation stress
Hout(idW3zx) == F ! Szx WEC_Mellor 3D Szx radiation stress
Hout(idW3zy) == F ! Szy WEC_Mellor 3D Szy radiation stress

Hout(idWztw) == F ! zetaw WEC_VF quasi-static sea level adjustment
Hout(idWqsp) == F ! qsp WEC_VF quasi-static pressure
Hout(idWbeh) == F ! bh WEC_VF Bernoulli head

Hout(idU2rs) == F ! ubar_Wecstress WEC 2D U-stress
Hout(idV2rs) == F ! vbar_Wecstress WEC 2D V-stress
Hout(idU3rs) == F ! u_Rstress WEC 3D U-stress
Hout(idV3rs) == F ! v_Rstress WEC 3D V-stress

Hout(idU2Sd) == F ! ubar_stokes 2D U-Stokes velocity
Hout(idV2Sd) == F ! vbar_stokes 2D V-Stokes velocity
Hout(idU3Sd) == F ! u_stokes 3D U-Stokes velocity
Hout(idV3Sd) == F ! v_stokes 3D V-Stokes velocity
Hout(idW3Sd) == F ! omega_stokes 3D Omega-Stokes velocity
Hout(idW3St) == F ! w_stokes 3D W-Stokes velocity
Hout(idWamp) == F ! Hwave wave height
Hout(idWlen) == F ! Lwave wave length-mean
Hout(idWlep) == F ! Lwavep wave length-peak
Hout(idWdir) == F ! Dwave wave direction
Hout(idWptp) == F ! Pwave_top wave surface period
Hout(idWpbt) == F ! Pwave_bot wave bottom period
Hout(idWorb) == F ! Uwave_rms wave bottom orbital velocity
Hout(idWbrk) == F ! Wave_break wave breaking (percent)
Hout(idUwav) == F ! uWave wave-depth avgeraged U-velocity
Hout(idVwav) == F ! vWave wave-depth avgeraged V-velocity
Hout(idWdif) == F ! Dissip_fric wave dissipation due to bottom friction
Hout(idWdib) == F ! Dissip_break wave dissipation due to breaking
Hout(idWdiw) == F ! Dissip_wcap wave dissipation due to white capping
Hout(idWdis) == F ! Dissip_roller wave roller dissipation
Hout(idWrol) == F ! rollA wave roller action density

Hout(idRunoff) == F ! Runoff surface runoff from land
Hout(idPair) == T ! Pair surface air pressure
Hout(idUair) == F ! Uair surface U-wind component
Hout(idVair) == F ! Vair surface V-wind component

Hout(idTsur) == F F ! shflux, ssflux surface net heat and salt flux
Hout(idLhea) == F ! latent latent heat flux
Hout(idShea) == F ! sensible sensible heat flux
Hout(idLrad) == F ! lwrad longwave radiation flux
Hout(idSrad) == F ! swrad shortwave radiation flux
Hout(idEmPf) == F ! EminusP E-P flux
Hout(idevap) == F ! evaporation evaporation rate
Hout(idrain) == F ! rain precipitation rate

Hout(idDano) == F ! rho density anomaly
Hout(idVvis) == F ! AKv vertical viscosity
Hout(idTdif) == F ! AKt vertical T-diffusion
Hout(idSdif) == F ! AKs vertical Salinity diffusion
Hout(idHsbl) == F ! Hsbl depth of surface boundary layer
Hout(idHbbl) == F ! Hbbl depth of bottom boundary layer
Hout(idMtke) == F ! tke turbulent kinetic energy
Hout(idMtls) == F ! gls turbulent length scale

! Logical switches (TRUE/FALSE) to activate writing of extra inert passive
! tracers other than biological and sediment tracers. An inert passive tracer
! is one that it is only advected and diffused. Other processes are ignored.
! These tracers include, for example, dyes, pollutants, oil spills, etc.
! NPT values are expected. However, these switches can be activated using
! compact parameter specification.

Hout(inert) == F ! dye_01, ... inert passive tracers

! Logical switches (TRUE/FALSE) to activate writing of time-averaged
! fields into AVERAGE output file.

Aout(idUvel) == F ! u 3D U-velocity
Aout(idVvel) == F ! v 3D V-velocity
Aout(idu3dE) == F ! u_eastward 3D U-eastward at RHO-points
Aout(idv3dN) == F ! v_northward 3D V-northward at RHO-points
Aout(idWvel) == F ! w 3D W-velocity
Aout(idOvel) == F ! omega omega vertical velocity
Aout(idUbar) == F ! ubar 2D U-velocity
Aout(idVbar) == F ! vbar 2D V-velocity
Aout(idu2dE) == F ! ubar_eastward 2D U-eastward at RHO-points
Aout(idv2dN) == F ! vbar_northward 2D V-northward at RHO-points
Aout(idFsur) == F ! zeta free-surface
Aout(idBath) == F ! bath time-dependent bathymetry

Aout(idTvar) == F F ! temp, salt temperature and salinity

Aout(idUair) == F ! Uwind surface U-wind
Aout(idVair) == F ! Vwind surface V-wind
Aout(idUairE) == F ! Uwind_eastward surface U-wind
Aout(idVairN) == F ! Vwind_northward surface V-wind
Aout(idUsms) == F ! sustr surface U-stress
Aout(idVsms) == F ! svstr surface V-stress
Aout(idUbms) == F ! bustr bottom U-stress
Aout(idVbms) == F ! bvstr bottom V-stress

Aout(idUbrs) == F ! bustrc bottom U-current stress
Aout(idVbrs) == F ! bvstrc bottom V-current stress
Aout(idUbws) == F ! bustrw bottom U-wave stress
Aout(idVbws) == F ! bvstrw bottom V-wave stress
Aout(idUbcs) == F ! bustrcwmax bottom max wave-current U-stress
Aout(idVbcs) == F ! bvstrcwmax bottom max wave-current V-stress
Aout(idUVwc) == F ! bstrcwmax bottom max wave-current stress magnitude

Aout(idUbot) == F ! Ubot bed wave orbital U-velocity
Aout(idVbot) == F ! Vbot bed wave orbital V-velocity
Aout(idUbur) == F ! Ur bottom U-velocity above bed
Aout(idVbvr) == F ! Vr bottom V-velocity above bed

Aout(idW2xx) == F ! Sxx_bar WEC_Mellor 2D Sxx radiation stress
Aout(idW2xy) == F ! Sxy_bar WEC_Mellor 2D Sxy radiation stress
Aout(idW2yy) == F ! Syy_bar WEC_Mellor 2D Syy radiation stress
Aout(idW3xx) == F ! Sxx WEC_Mellor 3D Sxx radiation stress
Aout(idW3xy) == F ! Sxy WEC_Mellor 3D Sxy radiation stress
Aout(idW3yy) == F ! Syy WEC_Mellor 3D Syy radiation stress
Aout(idW3zx) == F ! Szx WEC_Mellor 3D Szx radiation stress
Aout(idW3zy) == F ! Szy WEC_Mellor 3D Szy radiation stress

Aout(idWztw) == F ! zetaw WEC_VF quasi-static sea level adjustment
Aout(idWqsp) == F ! qsp WEC_VF quasi-static pressure
Aout(idWbeh) == F ! bh WEC_VF Bernoulli head

Aout(idU2rs) == F ! ubar_Wecstress WEC 2D U-stress
Aout(idV2rs) == F ! vbar_Wecstress WEC 2D V-stress
Aout(idU3rs) == F ! u_Rstress WEC 3D U-stress
Aout(idV3rs) == F ! v_Rstress WEC 3D V-stress

Aout(idU2Sd) == F ! ubar_stokes 2D U-Stokes velocity
Aout(idV2Sd) == F ! vbar_stokes 2D V-Stokes velocity
Aout(idU3Sd) == F ! u_stokes 3D U-Stokes velocity
Aout(idV3Sd) == F ! v_stokes 3D V-Stokes velocity
Aout(idW3Sd) == F ! omega_stokes 3D Omega-Stokes velocity
Aout(idW3St) == F ! w_stokes 3D W-Stokes velocity
Aout(idWamp) == F ! Hwave wave height
Aout(idWlen) == F ! Lwave wave length-mean
Aout(idWlep) == F ! Lwavep wave length-peak
Aout(idWdir) == F ! Dwave wave direction
Aout(idWptp) == F ! Pwave_top wave surface period
Aout(idWpbt) == F ! Pwave_bot wave bottom period
Aout(idWorb) == F ! Uwave_rms wave bottom orbital velocity
Aout(idWbrk) == F ! Wave_break wave breaking (percent)
Aout(idUwav) == F ! uWave wave-depth avgeraged U-velocity
Aout(idVwav) == F ! vWave wave-depth avgeraged V-velocity
Aout(idWdif) == F ! Dissip_fric wave dissipation due to bottom friction
Aout(idWdib) == F ! Dissip_break wave dissipation due to breaking
Aout(idWdiw) == F ! Dissip_wcap wave dissipation due to white capping
Aout(idWdis) == F ! Dissip_roller wave roller dissipation
Aout(idWrol) == F ! rollA wave roller action density

Aout(idRunoff) == F ! Runoff surface runoff from land
Aout(idPair) == F ! Pair surface air pressure
Aout(idUair) == F ! Uair surface U-wind component
Aout(idVair) == F ! Vair surface V-wind component

Aout(idTsur) == F F ! shflux, ssflux surface net heat and salt flux
Aout(idLhea) == F ! latent latent heat flux
Aout(idShea) == F ! sensible sensible heat flux
Aout(idLrad) == F ! lwrad longwave radiation flux
Aout(idSrad) == F ! swrad shortwave radiation flux
Aout(idevap) == F ! evaporation evaporation rate
Aout(idrain) == F ! rain precipitation rate

Aout(idDano) == F ! rho density anomaly
Aout(idVvis) == F ! AKv vertical viscosity
Aout(idTdif) == F ! AKt vertical T-diffusion
Aout(idSdif) == T ! AKs vertical Salinity diffusion
Aout(idHsbl) == F ! Hsbl depth of surface boundary layer
Aout(idHbbl) == F ! Hbbl depth of bottom boundary layer

Aout(id2dRV) == F ! pvorticity_bar 2D relative vorticity
Aout(id3dRV) == F ! pvorticity 3D relative vorticity
Aout(id2dPV) == F ! rvorticity_bar 2D potential vorticity
Aout(id3dPV) == F ! rvorticity 3D potential vorticity

Aout(idu3dD) == F ! u_detided detided 3D U-velocity
Aout(idv3dD) == F ! v_detided detided 3D V-velocity
Aout(idu2dD) == F ! ubar_detided detided 2D U-velocity
Aout(idv2dD) == F ! vbar_detided detided 2D V-velocity
Aout(idFsuD) == F ! zeta_detided detided free-surface

Aout(idTrcD) == F F ! temp_detided, ... detided temperature and salinity

Aout(idHUav) == F ! Huon u-volume flux, Huon
Aout(idHVav) == F ! Hvom v-volume flux, Hvom
Aout(idUUav) == F ! uu quadratic <u*u> term
Aout(idUVav) == F ! uv quadratic <u*v> term
Aout(idVVav) == F ! vv quadratic <v*v> term
Aout(idU2av) == F ! ubar2 quadratic <ubar*ubar> term
Aout(idV2av) == F ! vbar2 quadratic <vbar*vbar> term
Aout(idZZav) == F ! zeta2 quadratic <zeta*zeta> term

Aout(idTTav) == F F ! temp_2, ... quadratic <t*t> tracer terms
Aout(idUTav) == F F ! u_temp, ... quadratic <u*t> tracer terms
Aout(idVTav) == F F ! v_temp, ... quadratic <v*t> tracer terms
Aout(iHUTav) == F F ! Huon_temp, ... tracer volume flux, <Huon*t>
Aout(iHVTav) == F F ! Hvom_temp, ... tracer volume flux, <Hvom*t>

! Logical switches (TRUE/FALSE) to activate writing of extra inert passive
! tracers other than biological and sediment tracers into the AVERAGE file.

Aout(inert) == F ! dye_01, ... inert passive tracers

! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
! 2D momentum (ubar,vbar) diagnostic terms into DIAGNOSTIC output file.

Dout(M2rate) == F ! ubar_accel, ... acceleration
Dout(M2pgrd) == F ! ubar_prsgrd, ... pressure gradient
Dout(M2fcor) == F ! ubar_cor, ... Coriolis force
Dout(M2hadv) == F ! ubar_hadv, ... horizontal total advection
Dout(M2xadv) == F ! ubar_xadv, ... horizontal XI-advection
Dout(M2yadv) == F ! ubar_yadv, ... horizontal ETA-advection
Dout(M2hrad) == F ! ubar_hrad, ... horizontal total wec_mellor radiation stress
Dout(M2hvis) == F ! ubar_hvisc, ... horizontal total viscosity
Dout(M2xvis) == F ! ubar_xvisc, ... horizontal XI-viscosity
Dout(M2yvis) == F ! ubar_yvisc, ... horizontal ETA-viscosity
Dout(M2sstr) == F ! ubar_sstr, ... surface stress
Dout(M2bstr) == F ! ubar_bstr, ... bottom stress
Dout(M2hjvf) == F ! 2D wec_vf horizontal J vortex force
Dout(M2kvrf) == F ! 2D wec_vf K vortex force
Dout(M2fsco) == F ! 2D wec_vf coriolis-stokes
Dout(M2bstm) == F ! 2D wec_vf bottom streaming
Dout(M2sstm) == F ! 2D wec_vf surface streaming
Dout(M2wrol) == F ! 2D wec_vf wave roller accel
Dout(M2wbrk) == F ! 2D wec_vf wave breaking
Dout(M2zeta) == F ! 2D wec_vf Eulerian sea level adjustment
Dout(M2zetw) == F ! 2D wec_vf quasi-static sea level adjustment
Dout(M2zqsp) == F ! 2D wec_vf quasi-static pressure
Dout(M2zbeh) == F ! 2D wec_vf Bernoulli head

! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
! 3D momentum (u,v) diagnostic terms into DIAGNOSTIC output file.

Dout(M3rate) == F ! u_accel, ... acceleration
Dout(M3pgrd) == F ! u_prsgrd, ... pressure gradient
Dout(M3fcor) == F ! u_cor, ... Coriolis force
Dout(M3hadv) == F ! u_hadv, ... horizontal total advection
Dout(M3xadv) == F ! u_xadv, ... horizontal XI-advection
Dout(M3yadv) == F ! u_yadv, ... horizontal ETA-advection
Dout(M3vadv) == F ! u_vadv, ... vertical advection
Dout(M3hrad) == F ! u_hrad, ... horizontal total wec_mellor radiation stress
Dout(M3vrad) == F ! v_hrad, ... vertical total wec_mellor radiation stress
Dout(M3hvis) == F ! u_hvisc, ... horizontal total viscosity
Dout(M3xvis) == F ! u_xvisc, ... horizontal XI-viscosity
Dout(M3yvis) == F ! u_yvisc, ... horizontal ETA-viscosity
Dout(M3vvis) == F ! u_vvisc, ... vertical viscosity
Dout(M3vjvf) == F ! 3D wec_vf vertical J vortex force
Dout(M3hjvf) == F ! 3D wec_vf horizontal J vortex force
Dout(M3kvrf) == F ! 3D wec_vf K vortex force
Dout(M3fsco) == F ! 3D wec_vf coriolis-stokes
Dout(M3bstm) == F ! 3D wec_vf bottom streaming
Dout(M3sstm) == F ! 3D wec_vf surface streaming
Dout(M3wrol) == F ! 3D wec_vf wave roller accel
Dout(M3wbrk) == F ! 3D wec_vf wave breaking

! Logical switches (TRUE/FALSE) to activate writing of time-averaged,
! active (temperature and salinity) and passive (inert) tracer diagnostic
! terms into DIAGNOSTIC output file: [1:NAT+NPT,Ngrids].

Dout(iTrate) == F F ! temp_rate, ... time rate of change
Dout(iThadv) == F F ! temp_hadv, ... horizontal total advection
Dout(iTxadv) == F F ! temp_xadv, ... horizontal XI-advection
Dout(iTyadv) == F F ! temp_yadv, ... horizontal ETA-advection
Dout(iTvadv) == F F ! temp_vadv, ... vertical advection
Dout(iThdif) == F F ! temp_hdiff, ... horizontal total diffusion
Dout(iTxdif) == F F ! temp_xdiff, ... horizontal XI-diffusion
Dout(iTydif) == F F ! temp_ydiff, ... horizontal ETA-diffusion
Dout(iTsdif) == F F ! temp_sdiff, ... horizontal S-diffusion
Dout(iTvdif) == F F ! temp_vdiff, ... vertical diffusion

! Generic User parameters, [1:NUSER].

NUSER = 0
USER = 0.d0

! NetCDF-4/HDF5 compression parameters for output files.

NC_SHUFFLE = 1 ! if non-zero, turn on shuffle filter
NC_DEFLATE = 1 ! if non-zero, turn on deflate filter
NC_DLEVEL = 1 ! deflate level [0-9]

! Input NetCDF file names, [1:Ngrids].

GRDNAME == /home/hpd14/Documents/COAWST_V3.2/Projects/zfirst/zfirst_grid_cco.nc
ININAME == /home/hpd14/Documents/COAWST_V3.2/Projects/zfirst/zfirst_init_cco.nc
! ITLNAME == ocean_itl.nc
! IRPNAME == ocean_irp.nc
! IADNAME == ocean_iad.nc
! FWDNAME == ocean_fwd.nc
! ADSNAME == ocean_ads.nc

! Nesting grids connectivity data: contact points information. This
! NetCDF file is special and complex. It is currently generated using
! the script "matlab/grid/contact.m" from the Matlab repository.

! NGCNAME = ocean_ngc.nc

! Input lateral boundary conditions and climatology file names. The
! USER has the option to split input data time records into several
! NetCDF files (see prologue instructions above). If so, use a single
! line per entry with a vertical bar (|) symbol after each entry,
! except the last one.

NCLMFILES == 1 ! number of climate files
! CLMNAME == ocean_clm.nc

NBCFILES == 1 ! number of boundary files
! BRYNAME == ocean_bry.nc

! Input climatology nudging coefficients file name.

! NUDNAME == ocean_nud.nc

! Input Sources/Sinks forcing (like river runoff) file name.

! SSFNAME == ocean_rivers.nc

! Input forcing NetCDF file name(s). The USER has the option to enter
! several file names for each nested grid. For example, the USER may
! have different files for wind products, heat fluxes, tides, etc.
! The model will scan the file list and will read the needed data from
! the first file in the list containing the forcing field. Therefore,
! the order of the file names is very important. If using multiple forcing
! files per grid, first enter all the file names for grid 1, then grid 2,
! and so on. It is also possible to split input data time records into
! several NetCDF files (see prologue instructions above). Use a single line
! per entry with a continuation (\) or vertical bar (|) symbol after each
! entry, except the last one.

NFFILES == 1 ! number of forcing files

FRCNAME == /home/hpd14/Documents/COAWST_V3.2/Projects/zfirst/zfirst_wind_cco.nc

! Output NetCDF file names, [1:Ngrids].

GSTNAME == ocean_gst.nc
RSTNAME == ocean_rst.nc
HISNAME == ocean_zfirst_his.nc
TLMNAME == ocean_tlm.nc
TLFNAME == ocean_tlf.nc
ADJNAME == ocean_adj.nc
AVGNAME == ocean_avg.nc
DIANAME == ocean_dia.nc
STANAME == ocean_sta.nc
FLTNAME == ocean_flt.nc

! Input ASCII parameter filenames.

APARNAM = ROMS/External/s4dvar.in
SPOSNAM = ROMS/External/stations.in
FPOSNAM = ROMS/External/floats.in
BPARNAM = ROMS/External/bioFasham.in
SPARNAM = ROMS/External/sediment.in
USRNAME = ROMS/External/MyFile.dat

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#3 Unread post by hpd14thu »

/*
** svn $Id: sandy.h 25 2007-04-09 23:43:58Z jcwarner $
*******************************************************************************
** Copyright (c) 2002-2007 The ROMS/TOMS Group **
** Licensed under a MIT/X style license **
** See License_ROMS.txt **
*******************************************************************************
**
** Options for Sandy Test.
**
** Application flag: SANDY
*/

#define ROMS_MODEL
#undef NESTING
/*#define WRF_MODEL
#define SWAN_MODEL
#define MCT_LIB
#define MCT_INTERP_OC2AT
#define MCT_INTERP_WV2AT
#define MCT_INTERP_OC2WV*/

#if defined WRF_MODEL && defined SWAN_MODEL
# define DRAGLIM_DAVIS
# define COARE_TAYLOR_YELLAND
#endif

#ifdef ROMS_MODEL
/* Physics + numerics */
# define UV_ADV
# define UV_COR
# define UV_VIS2
# define MIX_S_UV
# undef TS_FIXED
# define TS_U3HADVECTION
# define TS_C4VADVECTION
# undef TS_MPDATA

# undef SSW_BBL
# ifdef SSW_BBL
# define SSW_CALC_ZNOT
/* # define ANA_SEDIMENT */
# else
# define UV_LOGDRAG
# endif
# if !defined SWAN_MODEL && defined SSW_BBL
# define ANA_WWAVE
# endif

# define DJ_GRADPS
# define TS_DIF2
# define MIX_GEO_TS
# define CURVGRID
# define SALINITY
# define SOLVE3D
# define AVERAGES
# define NONLIN_EOS


/* Grid and Initial */
# define MASKING

/* Forcing */

# ifdef WRF_MODEL
# undef BULK_FLUXES
# define ATM2OCN_FLUXES
# define ANA_SSFLUX
# undef LONGWAVE_OUT
# define ATM_PRESS
# else
# define BULK_FLUXES
/* # define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SSFLUX */
# define ATM_PRESS
/*# define LONGWAVE_OUT*/
# endif
# define ANA_BTFLUX
# define ANA_BSFLUX
# define EMINUSP
# define SOLAR_SOURCE

/* Turbulence closure */
# define GLS_MIXING
# undef MY25_MIXING
# define AKLIMIT

# if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# endif

/* Output */
# undef DIAGNOSTICS_UV
# undef DIAGNOSTICS_TS
#endif

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#4 Unread post by hpd14thu »

And here is my result:


Lateral Boundary Conditions: NLM
============================

Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------

zeta 1 Closed Closed Closed Closed

ubar 1 Closed Closed Closed Closed

vbar 1 Closed Closed Closed Closed

u 1 Closed Closed Closed Closed

v 1 Closed Closed Closed Closed

temp 1 Closed Closed Closed Closed

salt 1 Closed Closed Closed Closed

tke 1 Closed Closed Closed Closed

Activated C-preprocessing Options:

ZFIRST ZFIRST
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
ATM_PRESS Impose atmospheric pressure onto sea surface.
AVERAGES Writing out time-averaged nonlinear model fields.
BULK_FLUXES Surface bulk fluxes parameterization.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EMINUSP Compute Salt Flux using E-P.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLAR_SOURCE Solar Radiation Source Term.
SOLVE3D Solving 3D Primitive Equations.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.

Process Information:

Node # 0 (pid= 61073) is active.

INITIAL: Configuring and initializing forward nonlinear model ...
*******

Vertical S-coordinate System, Grid 01:

level S-coord Cs-curve Z at hmin at hc half way at hmax

20 0.0000000 0.0000000 0.000 0.000 0.000 0.000
19 -0.0500000 -0.0000692 -0.501 -0.501 -0.501 -0.501
18 -0.1000000 -0.0002879 -1.003 -1.003 -1.003 -1.003
17 -0.1500000 -0.0006915 -1.507 -1.507 -1.507 -1.507
16 -0.2000000 -0.0013455 -2.013 -2.013 -2.013 -2.013
15 -0.2500000 -0.0023555 -2.524 -2.524 -2.524 -2.524
14 -0.3000000 -0.0038848 -3.039 -3.039 -3.039 -3.039
13 -0.3500000 -0.0061801 -3.562 -3.562 -3.562 -3.562
12 -0.4000000 -0.0096110 -4.096 -4.096 -4.096 -4.096
11 -0.4500000 -0.0147275 -4.647 -4.647 -4.647 -4.647
10 -0.5000000 -0.0223460 -5.223 -5.223 -5.223 -5.223
9 -0.5500000 -0.0336719 -5.837 -5.837 -5.837 -5.837
8 -0.6000000 -0.0504756 -6.505 -6.505 -6.505 -6.505
7 -0.6500000 -0.0753361 -7.253 -7.253 -7.253 -7.253
6 -0.7000000 -0.1119642 -8.120 -8.120 -8.120 -8.120
5 -0.7500000 -0.1655988 -9.156 -9.156 -9.156 -9.156
4 -0.8000000 -0.2434198 -10.434 -10.434 -10.434 -10.434
3 -0.8500000 -0.3548031 -12.048 -12.048 -12.048 -12.048
2 -0.9000000 -0.5110137 -14.110 -14.110 -14.110 -14.110
1 -0.9500000 -0.7235564 -16.736 -16.736 -16.736 -16.736
0 -1.0000000 -1.0000000 -20.000 -20.000 -20.000 -20.000

Time Splitting Weights for Grid 01: ndtfast = 20 nfast = 29
==================================

Primary Secondary Accumulated to Current Step

1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651644 0.0501156998674241-0.0034654566136460 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000

ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000

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

1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.14200

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

Minimum X-grid spacing, DXmin = 1.01116846E-03 km
Maximum X-grid spacing, DXmax = 1.01116846E-03 km
Minimum Y-grid spacing, DYmin = 9.87495621E-04 km
Maximum Y-grid spacing, DYmax = 9.87495621E-04 km
Minimum Z-grid spacing, DZmin = 5.00691658E-01 m
Maximum Z-grid spacing, DZmax = 3.26443623E+00 m

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


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: zfirst_init_cco.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 = 1.80000000E+01 Max = 1.80000000E+01)
- salinity
(Min = 3.50000000E+01 Max = 3.50000000E+01)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E+00 Max = 8.00000000E+00)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.01300000E+03 Max = 1.01300000E+03)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.80000000E+01 Max = 1.80000000E+01)
GET_2DFLD - surface air relative humidity, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E-01 Max = 8.00000000E-01)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - net longwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - rain fall rate, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)

Basin information for Grid 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 - 00008640)

GET_2DFLD - surface u-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E+00 Max = 8.00000000E+00)
GET_2DFLD - surface v-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air pressure, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.01300000E+03 Max = 1.01300000E+03)
GET_2DFLD - surface air temperature, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.80000000E+01 Max = 1.80000000E+01)
GET_2DFLD - surface air relative humidity, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E-01 Max = 8.00000000E-01)
GET_2DFLD - solar shortwave radiation flux, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - net longwave radiation flux, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - rain fall rate, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)

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 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: ocean_zfirst_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: ocean_avg.nc
1 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000002
2 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000003
3 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000004
4 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000005
5 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000006
6 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000007
7 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000008
8 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000009
9 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000010
10 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000011
11 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000012
12 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000013
13 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000014
14 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000015
15 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000016
16 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000017
17 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000018
18 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000019
19 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000020
20 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000021
21 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000022
22 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000023
23 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000024
24 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000025
25 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000026
26 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000027
27 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000028
28 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000029
29 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000030
30 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000031
31 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000032
32 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000033
33 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000034
34 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000035
35 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000036
36 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000037
37 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000038
38 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000039
39 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000040
40 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000041
41 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000042
42 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000043
43 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000044
44 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000045
45 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000046
46 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000047
47 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000048
48 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000049
49 0 00:00:00 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000050
50 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000051
51 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000052
52 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000053
53 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000054
54 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000055
55 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000056
56 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000057
57 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000058
58 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000059
59 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000060
60 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000061
61 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000062
62 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000063
63 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000064
64 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000065
65 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000066
66 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000067
67 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000068
68 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000069
69 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000070
70 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000071
71 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000072
72 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000073
73 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000074
74 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000075
75 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000076
76 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000077
77 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000078
78 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000079
79 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000080
80 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000081
81 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E+04
(00,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000082
82 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E

......... all result is zero...

【and when I checked the calculation process in his file , I found that the sus ,the svs..in each step is 10^34, I don't know if it means that the calculation is blow out.】

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#5 Unread post by hpd14thu »

The forcing file just contains frc_bulk variables:
Uwind , Vwind , Pair..
But I don't use ustress adn vstress

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

Re: all result = 0?

#6 Unread post by kate »

【and when I checked the calculation process in his file , I found that the sus ,the svs..in each step is 10^34, I don't know if it means that the calculation is blow out.】
It means you have some serious problems. I'm not sure what, but I'd look carefully at all the inputs and outputs.

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#7 Unread post by hpd14thu »

kate wrote:
【and when I checked the calculation process in his file , I found that the sus ,the svs..in each step is 10^34, I don't know if it means that the calculation is blow out.】
It means you have some serious problems. I'm not sure what, but I'd look carefully at all the inputs and outputs.
Thank you for your reply, does it mean that there are some problems in my force nc file? I mean, I really don't know which file I should check or what I can do next.

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#8 Unread post by hpd14thu »

Here is the file which used to generae force.nc :

%1) Enter name of netcdf forcing file to be created.

% If it already exists it will be overwritten!!.

% fn='zfirst_wind_cco.nc';

fn='zfirst_wind_cco.nc';



%2) If you want to read the forcing data from another netcdf file, then

% enter read_forcing=1 and provide a name of history or rst file to read data from.

% Otherwise set read_forcing=0.

read_forcing=0;

data_file='coawst_us_20121024_13.nc';



%3) Enter start time of forcing file, in seconds.

% This time needs to be consistent with model time (ie dstart and time_ref).

% See *.in files for more detail.

if (read_forcing)

time=ncread(data_file,'wind_time');

t=length(wind_time);

wind_time=wind_time(t);

else

% wind_time=0;

% Pair_time=0:2 / 288:2;

% Tair_time=0:2 / 288:2;

% Qair_time=0:2 / 288:2;

% swrad_time=0:2 / 288:2;

% lwrad_time=0:2 / 288:2;

% rain_time=0:2 / 288:2;

wind_time = 0;

Pair_time = 0;

Tair_time = 0;

Qair_time = 0;

swrad_time = 0;

lwrad_time = 0;

rain_time = 0;

end



%4) Obtain grid information.

grid_file='/home/hpd14/Documents/COAWST_V3.2/Projects/zfirst/zfirst_grid_cco.nc' %<-enter name of grid here

x_rho=ncread(grid_file,'x_rho');

y_rho=ncread(grid_file,'y_rho');

[ix iy]=size(x_rho);

t=length(wind_time);

xrho = ix;

yrho = iy;



Gout=get_roms_grid(grid_file);



%5) Forcing Uwind, Vwind.



% Forcing values for Uwind, Vwind.

display('Forcing Uwind, Vwind')

%

if (read_forcing)

Uwind=ncread(data_file,'Uwind',[1 1 t],[Inf Inf 1]);

Uwind(isnan(Uwind))=0;

%

Vwind=ncread(data_file,'Vwind',[1 1 t],[Inf Inf 1]);

Vwind(isnan(Vwind))=0;

else



Uwind(1:xrho,1:yrho,1:length(wind_time)) = 8;

Vwind(1:xrho,1:yrho,1:length(wind_time)) = 0;

end



%6) Forcing Pair.



% Forcing values for Pair.

display('Forcing Pair')

%

if (read_forcing)

Pair=ncread(data_file,'Pair',[1 1 t],[Inf Inf 1]);

Pair(isnan(Pair))=0;
%

else
Pair(1:xrho,1:yrho,1:length(Pair_time)) = 1013;





%7) Forcing Tair.



% Forcing values for Tair.

display('Forcing Tair')

%

if (read_forcing)

Tair=ncread(data_file,'Tair',[1 1 t],[Inf Inf 1]);

Tair(isnan(Tair))=0;

%

else

Tair(1:xrho,1:yrho,1:length(Tair_time)) = 18;

end



%8) Forcing Qair.



% Forcing values for Qair.

display('Forcing Qair')

%

if (read_forcing)

Qair=ncread(data_file,'Qair',[1 1 t],[Inf Inf 1]);

Qair(isnan(Qair))=0;

%

else

Qair(1:xrho,1:yrho,1:length(Qair_time)) = 80;

end



%9) Forcing swrad.



% Forcing values for swrad.

display('Forcing swrad')

%

if (read_forcing)

swrad=ncread(data_file,'swrad',[1 1 t],[Inf Inf 1]);

swrad(isnan(swrad))=0;

%

else

swrad(1:xrho,1:yrho,1:length(swrad_time)) = 0;

end



%8) Forcing lwrad.



% Forcing values for lwrad.

display('Forcing lwrad')

%

if (read_forcing)

lwrad=ncread(data_file,'lwrad',[1 1 t],[Inf Inf 1]);

lwrad(isnan(lwrad))=0;

%

else

lwrad(1:xrho,1:yrho,1:length(lwrad_time)) = 0;

end



%9) Forcing rain.



% Forcing values for rain.

display('Forcing rain')

%

if (read_forcing)

rain=ncread(data_file,'rain',[1 1 t],[Inf Inf 1]);

rain(isnan(rain))=0;

%

else

rain(1:xrho,1:yrho,1:length(rain_time)) = 0;

end



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% END of USER INPUT %

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%create forcing file

time = 0;

create_roms_forcings(x_rho,y_rho,time,fn, 'Uwind', 'Vwind', 'Pair', 'Tair', 'Qair' ,'swrad', 'lwrad', 'rain')



%now write the data from the arrays to the netcdf file

disp(' ## Filling Variables in netcdf file with data...')



ncwrite(fn,'wind_time',wind_time);

ncwrite(fn,'Pair_time',Pair_time);

ncwrite(fn,'Tair_time',Tair_time);

ncwrite(fn,'Qair_time',Qair_time);

ncwrite(fn,'swrad_time',swrad_time);

ncwrite(fn,'lwrad_time',lwrad_time);

ncwrite(fn,'rain_time',rain_time);

ncwrite(fn,'x',x_rho);

ncwrite(fn,'y',y_rho);

ncwrite(fn,'Uwind',Uwind);

ncwrite(fn,'Vwind',Vwind);

ncwrite(fn,'Pair',Pair);

ncwrite(fn,'Qair',Qair);

ncwrite(fn,'Tair',Tair);

ncwrite(fn,'swrad',swrad);

ncwrite(fn,'lwrad',lwrad);

ncwrite(fn,'rain',rain);



%close file

disp(['created ', fn])

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#9 Unread post by hpd14thu »

kate wrote:
【and when I checked the calculation process in his file , I found that the sus ,the svs..in each step is 10^34, I don't know if it means that the calculation is blow out.】
It means you have some serious problems. I'm not sure what, but I'd look carefully at all the inputs and outputs.

And here is my force.cdl

netcdf zfirst_wind_cco {
dimensions:
xrho = 50 ;
yrho = 50 ;
time = 3 ;
wind_time = 3 ;
Pair_time = 3 ;
Tair_time = 3 ;
Qair_time = 3 ;
rain_time = 3 ;
swrad_time = 3 ;
lwrad_time = 3 ;
variables:
double time(time) ;
time:long_name = "atmospheric forcing time" ;
time:units = "days" ;
time:field = "time, scalar, series" ;
double x(yrho, xrho) ;
x:long_name = "x" ;
x:units = "meter" ;
x:field = "xp, scalar, series" ;
double y(yrho, xrho) ;
y:long_name = "y" ;
y:units = "meter" ;
y:field = "yp, scalar, series" ;
double wind_time(wind_time) ;
wind_time:long_name = "wind_time" ;
wind_time:units = "days" ;
wind_time:field = "Uwind_time, scalar, series" ;
double Uwind(wind_time, yrho, xrho) ;
Uwind:long_name = "surface u-wind component" ;
Uwind:units = "meter second-1" ;
Uwind:field = "Uwind, scalar, series" ;
Uwind:coordinates = "x y" ;
Uwind:time = "wind_time" ;
double Vwind(wind_time, yrho, xrho) ;
Vwind:long_name = "surface v-wind component" ;
Vwind:units = "meter second-1" ;
Vwind:field = "Vwind, scalar, series" ;
Vwind:coordinates = "x y" ;
Vwind:time = "wind_time" ;
double Pair_time(Pair_time) ;
Pair_time:long_name = "Pair_time" ;
Pair_time:units = "days" ;
Pair_time:field = "Pair_time, scalar, series" ;
double Pair(Pair_time, yrho, xrho) ;
Pair:long_name = "surface air pressure" ;
Pair:units = "millibar" ;
Pair:field = "Pair, scalar, series" ;
Pair:coordinates = "x y" ;
Pair:time = "Pair_time" ;
double Tair_time(Tair_time) ;
Tair_time:long_name = "Tair_time" ;
Tair_time:units = "days" ;
Tair_time:field = "Tair_time, scalar, series" ;
double Tair(Tair_time, yrho, xrho) ;
Tair:long_name = "surface air temperature" ;
Tair:units = "Celsius" ;
Tair:field = "Tair, scalar, series" ;
Tair:coordinates = "x y" ;
Tair:time = "Tair_time" ;
double Qair_time(Qair_time) ;
Qair_time:long_name = "Qair_time" ;
Qair_time:units = "days" ;
Qair_time:field = "Qair_time, scalar, series" ;
double Qair(Qair_time, yrho, xrho) ;
Qair:long_name = "surface air relative humidity" ;
Qair:units = "percentage" ;
Qair:field = "Qair, scalar, series" ;
Qair:coordinates = "x y" ;
Qair:time = "Qair_time" ;
double rain_time(rain_time) ;
rain_time:long_name = "rain_time" ;
rain_time:units = "days" ;
rain_time:field = "rain_time, scalar, series" ;
double rain(rain_time, yrho, xrho) ;
rain:long_name = "rain fall rate" ;
rain:units = "kilogram meter-2 second-1" ;
rain:field = "Qair, scalar, series" ;
rain:coordinates = "x y" ;
rain:time = "rain_time" ;
double swrad_time(swrad_time) ;
swrad_time:long_name = "swrad_time" ;
swrad_time:units = "days" ;
swrad_time:field = "swrad_time, scalar, series" ;
double swrad(swrad_time, yrho, xrho) ;
swrad:long_name = "solar shortwave radiation" ;
swrad:units = "Watts meter-2" ;
swrad:positive_value = "downward flux, heating" ;
swrad:negative_value = "upward flux, cooling" ;
swrad:field = "swrad, scalar, series" ;
swrad:coordinates = "x y" ;
swrad:time = "swrad_time" ;
double lwrad_time(lwrad_time) ;
lwrad_time:long_name = "lwrad_time" ;
lwrad_time:units = "days" ;
lwrad_time:field = "lwrad_time, scalar, series" ;
double lwrad(lwrad_time, yrho, xrho) ;
lwrad:long_name = "solar longwave radiation" ;
lwrad:units = "Watts meter-2" ;
lwrad:positive_value = "downward flux, heating" ;
lwrad:negative_value = "upward flux, cooling" ;
lwrad:field = "lwrad, scalar, series" ;
lwrad:coordinates = "x y" ;
lwrad:time = "lwrad_time" ;

// global attributes:
:history = "Created by create_roms_forcings on 23-Apr-2018 23:15:12" ;
:type = "bulk fluxes forcing file" ;
}

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

Re: all result = 0?

#10 Unread post by kate »

I don't know what to check either, but you need to be looking at the values in all the input and output files. Where is that weird 1e34 coming from? I like ncview for viewing these files.

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#11 Unread post by hpd14thu »

kate wrote:I don't know what to check either, but you need to be looking at the values in all the input and output files. Where is that weird 1e34 coming from? I like ncview for viewing these files.
Attachments
zfirst_wind_cco.nc
(511.73 KiB) Downloaded 244 times
zfirst_init_cco.nc
Thank you for your reply! Here is my input file.I have checked all of them . It seems that all values are normal. But the values in output files are all zero(not 10^34)


And I wonder that if there are somes problems in my head file, which leads to the problem : my wind values have been read but not used to calculated
(1.5 MiB) Downloaded 258 times

jklinck
Posts: 34
Joined: Mon Jun 30, 2003 2:29 pm
Location: CCPO/ODU, USA

Re: all result = 0?

#12 Unread post by jklinck »

Given the values in the wind forcing file (10^34), I think these are all "missing values" for netcdf. Are you sure that you assigned variables to the forcing file? I would plot the wind speeds in the forcing file to see that they are what you wanted. If the forcing file is ok, then make sure that the surface stress values in the model are correct. I did not look carefully at your first output file, but I assume that you don't have masking. So pick a grid index and print the values of the winds, other met parameters, and the wind stress at that point for every model time step to see that they are non-zero. And run the model for about 10 long time steps until you figure out why the forcing is not working.

Advice:
For a new simulation, run the model for 10 steps (saving every step) to see that something happens and that all the parts work without failure. Then run 100 steps (saving every tenth) to see that things are going well. Then run for a day or two (saving every hour or so). By this time, you should know if things are working, so you can start making longer simulations.

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#13 Unread post by hpd14thu »

jklinck wrote:Given the values in the wind forcing file (10^34), I think these are all "missing values" for netcdf. Are you sure that you assigned variables to the forcing file? I would plot the wind speeds in the forcing file to see that they are what you wanted. If the forcing file is ok, then make sure that the surface stress values in the model are correct. I did not look carefully at your first output file, but I assume that you don't have masking. So pick a grid index and print the values of the winds, other met parameters, and the wind stress at that point for every model time step to see that they are non-zero. And run the model for about 10 long time steps until you figure out why the forcing is not working.

Advice:
For a new simulation, run the model for 10 steps (saving every step) to see that something happens and that all the parts work without failure. Then run 100 steps (saving every tenth) to see that things are going well. Then run for a day or two (saving every hour or so). By this time, you should know if things are working, so you can start making longer simulations.

Hello,jklinck! I'm an undergraduate student. None of my classmates could use this model, so I can't find anybody in real life to help me. I really thank you for your help, it's sooo kind of you.

I just want to learn this model. So I try a easy case. In the case, I just set a closed region(50m*50m), all intitial conditions=0, and a surface uwind effect on the region: Uwind=8 m/s. But after I run it, all results=0. According to your suggestion,I check my forcing file,but I consider that all values are normal.

Attached pls find the forcing file for your necessary,Thank you!

And if possible, could you please provide me a little case about the Uwind, maybe I can compare our cases and find sth wrong in my case.
Attachments
zfirst_wind_cco.nc
(511.73 KiB) Downloaded 244 times

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

Re: all result = 0?

#14 Unread post by kate »

The winds file looks fine to me. What I would do next is add a bunch of print statements to the bulk flux routine, where winds are converted to stresses. Or watch it in a debugger - this should be small enough to run in serial mode so that gdb should work, right?

jklinck
Posts: 34
Joined: Mon Jun 30, 2003 2:29 pm
Location: CCPO/ODU, USA

Re: all result = 0?

#15 Unread post by jklinck »

To add to Kate's suggestion, I would run this model for about 10 time steps having compiled with debug on in the compile script (however you compile the model). Run in scalar mode (one processor). Make sure that you turn on array bounds checking. This makes sure that arrays are defined and indexes are properly specified. The fortran optimizer sometimes damages the code (with O3 and higher levels of optimization).

As Kate suggests, put a print statement in the place were wind in converted to stress.

Do you have BULK_FLUXES defined? the routine Nonlinear/bulk_flux.F is the place were winds are converted to stress. That would be the place to print values to make sure the stress is nonzero.

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#16 Unread post by hpd14thu »

jklinck wrote:To add to Kate's suggestion, I would run this model for about 10 time steps having compiled with debug on in the compile script (however you compile the model). Run in scalar mode (one processor). Make sure that you turn on array bounds checking. This makes sure that arrays are defined and indexes are properly specified. The fortran optimizer sometimes damages the code (with O3 and higher levels of optimization).

As Kate suggests, put a print statement in the place were wind in converted to stress.

Do you have BULK_FLUXES defined? the routine Nonlinear/bulk_flux.F is the place were winds are converted to stress. That would be the place to print values to make sure the stress is nonzero.
Thank you so much! I wonder that how can I define arrays?
And I have BULK_FLUXES.

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

Re: all result = 0?

#17 Unread post by kate »

John is saying to recompile with USE_DEBUG set to on. That should turn on array-bounds checking and the rest.

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: all result = 0?

#18 Unread post by wilkin »

Reading the log file output your quoted
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000082
82 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E
... that says to me the model has run for only 1 second of simulated time. If so, you're not going to see much action. What is your model DT?

You say you made a 50x50 m grid, which is very small so perhaps you have a tiny DT. Make your grid larger so you can run with a DT that lets something happen on geophysical timescales.

Your forcing files has times 0, 1, 2 in days, so you can run at most 2 days which should be enough to see something.

I see nothing obviously wrong with your forcing file. Your initial conditions have water temperature = 0 (freezing) but that of itself won't explain the stresses.

It appears you may have configured a complex system with sediments and boundary layers. I suggest you undefine those options for now and at least get a simple ocean circulation running with only the bulk fluxes computation setting the stresses.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#19 Unread post by hpd14thu »

wilkin wrote:Reading the log file output your quoted
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000082
82 0 00:00:01 0.000000E+00 9.608405E+01 9.608405E+01 4.601201E
... that says to me the model has run for only 1 second of simulated time. If so, you're not going to see much action. What is your model DT?

You say you made a 50x50 m grid, which is very small so perhaps you have a tiny DT. Make your grid larger so you can run with a DT that lets something happen on geophysical timescales.

Your forcing files has times 0, 1, 2 in days, so you can run at most 2 days which should be enough to see something.

I see nothing obviously wrong with your forcing file. Your initial conditions have water temperature = 0 (freezing) but that of itself won't explain the stresses.

It appears you may have configured a complex system with sediments and boundary layers. I suggest you undefine those options for now and at least get a simple ocean circulation running with only the bulk fluxes computation setting the stresses.
Dear Mr.wilkin,

Thank you for your help, it's really so kind of you! I have deleted the sediments and boundary layers in my init.file , but it did't work.

My dt = 0.02s ,and I think it's small enough. It has taken me 2 weeks to deal with this case, but I really didn't find anything abnormal.

Here is my project. I Look forward to your reply!

Paddy
Attachments
zfirst.zip
(79.28 KiB) Downloaded 221 times

bzhang
Posts: 11
Joined: Wed Mar 26, 2003 9:25 pm
Location: CICS/ESSIC University of Maryland

Re: all result = 0?

#20 Unread post by bzhang »

You probably set all sea grids as land pixels. So all the velocity/kinetic energy are zeros, but total potential energy and volume are not since the calculation of these diagnostics are not masked off. Can you make a plot your land/sea masks to check?
-Bin

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#21 Unread post by hpd14thu »

bzhang wrote:You probably set all sea grids as land pixels. So all the velocity/kinetic energy are zeros, but total potential energy and volume are not since the calculation of these diagnostics are not masked off. Can you make a plot your land/sea masks to check?
-Bin
Hello,Bin

Thank you for your help and suggestions! I see all of my mask in grid.nc , the values are all zero. So maybe the mask is normal? There should be some other problems?

Here is my grid.nc.

And I wonder that if possible, could you please provide a ease case? I don't find any wind-induced cases . I just want to try a easiest case which just a closed area and a constant wind , but I have dealed with it for 2 weeks and it did't work.

Thank you!
Attachments
zfirst_grid_cco.nc
(389.3 KiB) Downloaded 210 times

User avatar
hpd14thu
Posts: 68
Joined: Tue May 01, 2018 3:56 pm
Location: Tsinghua University

Re: all result = 0?

#22 Unread post by hpd14thu »

Hello,here is the part about reading the forcing file in result.out.

I wonder that if there is any problems in reading part.

I just set one variable : Uwind = 8, all of other values =0

LM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: zfirst_init_cco.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 = 0.00000000E+00 Max = 0.00000000E+00)
- salinity
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface u-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E+00 Max = 8.00000000E+00)
GET_2DFLD - surface v-wind component, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air pressure, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.01300000E+03 Max = 1.01300000E+03)
GET_2DFLD - surface air temperature, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.80000000E+01 Max = 1.80000000E+01)
GET_2DFLD - surface air relative humidity, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E-01 Max = 8.00000000E-01)
GET_2DFLD - solar shortwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - net longwave radiation flux, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - rain fall rate, t = 0 00:00:00
(Rec=0000001, Index=1, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)

Basin information for Grid 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 - 00008640)

GET_2DFLD - surface u-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E+00 Max = 8.00000000E+00)
GET_2DFLD - surface v-wind component, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface air pressure, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.01300000E+03 Max = 1.01300000E+03)
GET_2DFLD - surface air temperature, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 1.80000000E+01 Max = 1.80000000E+01)
GET_2DFLD - surface air relative humidity, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 8.00000000E-01 Max = 8.00000000E-01)
GET_2DFLD - solar shortwave radiation flux, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - net longwave radiation flux, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - rain fall rate, t = 1 00:00:00
(Rec=0000002, Index=2, File: zfirst_wind_cco.nc)
(Tmin= 0.0000 Tmax= 2.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)

bzhang
Posts: 11
Joined: Wed Mar 26, 2003 9:25 pm
Location: CICS/ESSIC University of Maryland

Re: all result = 0?

#23 Unread post by bzhang »

Make sure sea masks set to one,not zero. Change this first,then check others step by step.

Post Reply