i can't run model

[kobl1201]

i set my model

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

VARNAME = ../varinfo_cfsr.dat
! VARNAME = 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 == 197 ! Number of I-direction INTERIOR RHO-points
Mm == 183 ! Number of J-direction INTERIOR RHO-points
N == 50 ! 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 == 4 ! I-direction partition
NtileJ == 4 ! 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) == Cha Cha Cha Cha ! free-surface
LBC(isUbar) == Fla Fla Fla Fla ! 2D U-momentum
LBC(isVbar) == Fla Fla Fla Fla ! 2D V-momentum
LBC(isUvel) == RadNud RadNud RadNud RadNud ! 3D U-momentum
LBC(isVvel) == RadNud RadNud RadNud RadNud ! 3D V-momentum
LBC(isMtke) == Gra Gra Gra Gra ! mixing TKE

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

! Wec boundary conditions

LBC(isU2Sd) == Gra Gra Gra Gra ! 2D U-stokes
LBC(isV2Sd) == Gra Gra Gra Gra ! 2D V-stokes
LBC(isU3Sd) == Gra Gra Gra Gra ! 3D U-stokes
LBC(isV3Sd) == Gra Gra Gra Gra ! 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) == Cha Per Clo Per ! free-surface
ad_LBC(isUbar) == Fla Per Clo Per ! 2D U-momentum
ad_LBC(isVbar) == Fla Per Clo Per ! 2D U-momentum
ad_LBC(isUvel) == Gra Per Clo Per ! 3D U-momentum
ad_LBC(isVvel) == Gra Per Clo Per ! 3D V-momentum
ad_LBC(isMtke) == Gra Per Clo Per ! mixing TKE

ad_LBC(isTvar) == Gra Per Clo Per \ ! temperature
Gra Per Clo Per ! 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 == 172800 ! 365 days 1 day = 480 1440
DT == 180.0d0
NDTFAST == 30

! 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 == 172800
NSTA == 1
NFLT == 1
NINFO == 1

! Output history, average, diagnostic files parameters.

LDEFOUT == T
NHIS == 172800
NDEFHIS == 0 !175680
NTSAVG == 1
NAVG == 480 ! 30 days !480 ! 1 day
NDEFAVG == 0 !175680
p
NTSDIA == 1
NDIA == 1
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 == 20.0d0 20.0d0 ! m2/s
TNU4 == 0.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 == 100.0d0 ! 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 == 1.0d-6 1.0d-6 1.0d-6 1.0d-6 ! m2/s

ad_AKT_fac == 1.0d0 ! nondimensional

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

AKV_BAK == 1.0d-5 1.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.02d0 ! 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 == 1.0d0 ! surface stretching parameter
THETA_B == 4.0d0 ! bottom stretching parameter
TCLINE == 5.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 = 55927.0d0 ! days
TIDE_START = 0.0d0 ! days
TIME_REF = 18581117.00 ! yyyymmdd.dd

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

TNUDG == 10.0d0 10.0d0 ! days
ZNUDG == 10.0d0 ! days
M2NUDG == 10.0d0 ! days
M3NUDG == 10.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 == 10.0d0 ! nondimensional

! Linear equation of State parameters:

R0 == 1027.0d0 ! kg/m3
T0 == 10.0d0 ! Celsius
S0 == 35.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 == T ! sea-surface height
Lm2CLM == T ! 2D momentum
Lm3CLM == T ! 3D momentum

LtracerCLM == T T ! 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 == T ! 2D momentum
LnudgeM3CLM == T ! 3D momentum

LnudgeTCLM == T T ! 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) == F ! u 3D U-velocity
Hout(idVvel) == F ! 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) == F ! w 3D W-velocity
Hout(idOvel) == F ! 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) == F ! zeta free-surface
Hout(idBath) == F ! bath time-dependent bathymetry

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

Hout(idUair) == F ! Uwind surface U-wind
Hout(idVair) == F ! Vwind surface V-wind
Hout(idUairE) == F ! Uwind_eastward surface U-wind
Hout(idVairN) == F ! Vwind_northward surface V-wind
Hout(idUsms) == F ! sustr surface U-stress
Hout(idVsms) == F ! svstr surface V-stress
Hout(idUbms) == F ! bustr bottom U-stress
Hout(idVbms) == F ! 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) == T ! ubar_stokes 2D U-Stokes velocity
Hout(idV2Sd) == T ! vbar_stokes 2D V-Stokes velocity
Hout(idU3Sd) == T ! u_stokes 3D U-Stokes velocity
Hout(idV3Sd) == T ! v_stokes 3D V-Stokes velocity
Hout(idW3Sd) == T ! omega_stokes 3D Omega-Stokes velocity
Hout(idW3St) == T ! w_stokes 3D W-Stokes velocity
Hout(idWamp) == T ! Hwave wave height
Hout(idWlen) == T ! Lwave wave length-mean
Hout(idWlep) == T ! Lwavep wave length-peak
Hout(idWdir) == T ! Dwave wave direction
Hout(idWptp) == T ! Pwave_top wave surface period
Hout(idWpbt) == T ! Pwave_bot wave bottom period
Hout(idWorb) == T ! Uwave_rms wave bottom orbital velocity
Hout(idWbrk) == F ! Wave_break wave breaking (percent)
Hout(idUwav) == T ! uWave wave-depth avgeraged U-velocity
Hout(idVwav) == T ! 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) == F ! 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) == T ! dye_01, ... inert passive tracers

! Logical switches (TRUE/FALSE) to activate writing of exposed sediment
! layer properties into HISTORY output file. Currently, MBOTP properties
! are expected for the bottom boundary layer and/or sediment models:
!
! idBott( 1=isd50) grain_diameter mean grain diameter
! idBott( 2=idens) grain_density mean grain density
! idBott( 3=iwsed) settling_vel mean settling velocity
! idBott( 4=itauc) erosion_stress critical erosion stress
! idBott( 5=irlen) ripple_length ripple length
! idBott( 6=irhgt) ripple_height ripple height
! idBott( 7=ibwav) bed_wave_amp wave excursion amplitude
! idBott( 8=izdef) Zo_def default bottom roughness
! idBott( 9=izapp) Zo_app apparent bottom roughness
! idBott(10=izNik) Zo_Nik Nikuradse bottom roughness
! idBott(11=izbio) Zo_bio biological bottom roughness
! idBott(12=izbfm) Zo_bedform bed form bottom roughness
! idBott(13=izbld) Zo_bedload bed load bottom roughness
! idBott(14=izwbl) Zo_wbl wave bottom roughness
! idBott(15=iactv) active_layer_thickness active layer thickness
! idBott(16=ishgt) saltation saltation height
! idBott(17=imaxD) dep_net maximum inundation depth
! idBott(18=idnet) net erosion + dep Erosion or deposition
! idBott(19=idoff) tau critical offset dmix erodibility profile offset
! idBott(20=idslp) tau critical slope dmix or erodibility slope
! idBott(21=idtim) erodibility time scale erodibility profile restore time
! idBott(22=idbmx) diffusivity db_max Bed biodifusivity maximum
! idBott(23=idbmm) diffusivity db_m Bed biodifusivity minimum
! idBott(24=idbzs) diffusivity db_zs Bed biodifusivity zs
! idBott(25=idbzm) diffusivity db_zm Bed biodifusivity zm
! idBott(26=idbzp) diffusivity db_zphi Bed biodifusivity phi
! idBott(27=idprp) cohesive behavior cohesive behavior
!
! 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2
! 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7
!
Hout(idBott) == T T T T T T T T T F F F F F T F T T F F F F F F F F F

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

Aout(idUvel) == T ! u 3D U-velocity
Aout(idVvel) == T ! 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) == T ! w 3D W-velocity
Aout(idOvel) == F ! omega omega vertical velocity
Aout(idUbar) == T ! ubar 2D U-velocity
Aout(idVbar) == T ! 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) == T ! zeta free-surface
Aout(idBath) == F ! bath time-dependent bathymetry

Aout(idTvar) == T T ! 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) == F ! AKs vertical Salinity diffusion
Aout(idHsbl) == T ! Hsbl depth of surface boundary layer
Aout(idHbbl) == T ! 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/dhkim/NEW_MODEL/ROMS/Projects/TOPO/KODOS_roms_grid.nc
ININAME == ./INPUT/KODOS_ini_20120101.nc
ITLNAME == ocean_itl.nc
IRPNAME == ocean_irp.nc
IADNAME == ocean_iad.nc
FWDNAME == ocean_fwd.nc
ADSNAME == ocean_ads.nc
stop
! 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 == ./INPUT/KODOS_clm_2012_01.nc |
./INPUT/KODOS_clm_2012_02.nc |
./INPUT/KODOS_clm_2012_03.nc |
./INPUT/KODOS_clm_2012_04.nc |
./INPUT/KODOS_clm_2012_05.nc |
./INPUT/KODOS_clm_2012_06.nc

NBCFILES == 1 ! number of boundary files
BRYNAME == /home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/kodos_bdy_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/kodos_bdy_2012_OI.nc

! ./INPUT/KODOS_bdy_2012.nc
! ../2013/INPUT/KODOS_bdy_2013.nc

! BRYNAME == /home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/time_set/KODOS_bdy_2012.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 == 8 ! number of unique forcing files

FRCNAME == /home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_rain_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_rain_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_tair_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_tair_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_qair_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_qair_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_pair_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_pair_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_swrad_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_swrad_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2011/INPUT/cfsv_lwrad_2011_OI.nc |
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/cfsv_lwrad_2012_OI.nc \
/home/dhkim/NEW_MODEL/ROMS/Projects/2012/INPUT/coads_roms_frc.nc

! FRCNAME == ./Input_ROMS/KODOS_roms_frc.nc ! forcing file 1, grid 1

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

GSTNAME == ocean_gst_test1.nc
RSTNAME == ./Output_ROMS/KODOS_ocean_rst_2012_test1.nc
HISNAME == ./Output_ROMS/KODOS_ocean_his_2012_test1.nc
TLMNAME == ocean_tlm_test1.nc
TLFNAME == ocean_tlf_test1.nc
ADJNAME == ocean_adj_test1.nc
AVGNAME == ./Output_ROMS/KODOS_ocean_avg_2012_test1.nc
DIANAME == ocean_dia_test1.nc
STANAME == ocean_sta_test1.nc
FLTNAME == ocean_flt_test1.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


and i got erro message like this


GET_NGFLD - temperature southern boundary condition, t = 55927 00:00:00
(Rec=0000001, Index=1, File: kodos_bdy_2012_OI.nc)
(Tmin= 55927.0000 Tmax= 56292.0000)
(Min = 9.76048529E-01 Max = 2.75907154E+01)
GET_NGFLD - salinity southern boundary condition, t = 55927 00:00:00
(Rec=0000001, Index=1, File: kodos_bdy_2012_OI.nc)
(Tmin= 55927.0000 Tmax= 56292.0000)
(Min = 3.39031639E+01 Max = 3.51230469E+01)
GET_NGFLD - temperature northern boundary condition, t = 55927 00:00:00
(Rec=0000001, Index=1, File: kodos_bdy_2012_OI.nc)
(Tmin= 55927.0000 Tmax= 56292.0000)
(Min = 8.18639338E-01 Max = 2.61650333E+01)
GET_NGFLD - salinity northern boundary condition, t = 55927 00:00:00
(Rec=0000001, Index=1, File: kodos_bdy_2012_OI.nc)
(Tmin= 55927.0000 Tmax= 56292.0000)
(Min = 3.38812981E+01 Max = 3.47061501E+01)
GET_2DFLD - sea surface height climatology, t = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = -6.42439338E-02 Max = 2.62338230E-01)
GET_2DFLD - vertically integrated u-momentum climatologt = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = -8.00456479E-02 Max = 5.36643974E-02)
GET_2DFLD - vertically integrated v-momentum climatologt = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = -5.63130789E-02 Max = 5.74795417E-02)
GET_3DFLD - u-momentum component climatology, t = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = -1.07758546E+00 Max = 1.14394832E+00)
GET_3DFLD - v-momentum component climatology, t = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = -5.59971452E-01 Max = 9.15995300E-01)
GET_3DFLD - potential temperature climatology, t = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = 8.28390360E-01 Max = 2.71886845E+01)
GET_3DFLD - salinity climatology, t = 55928 00:00:00
(Rec=0000002, Index=2, File: KODOS_clm_2012_01.nc)
(Tmin= 55927.0000 Tmax= 55986.0000)
(Min = 3.36166954E+01 Max = 3.50767517E+01)

SET_NGFLD - current model time exceeds ending value for variable: zeta_west
TDAYS = 55927.0000
Data Tmin = 55927.0000 Data Tmax = 56292.0000
Data Tstr = 55926.0000 Data Tend = 55927.0000
TINTRP1 = 166886.0000 TINTRP2 = 166887.0000
FAC1 = 110960.0000 FAC2 = -110959.0000

Elapsed CPU time (seconds):

Node # 0 CPU: 1.141 Node # 1 CPU: 1.137
Node # 4 CPU: 1.128
Node # 8 CPU: 1.137
Node # 9 CPU: 1.133
Node # 10 CPU: 1.137
Node # 13 CPU: 1.137
Node # 15 CPU: 1.137
Node # 3 CPU: 1.141
Node # 5 CPU: 1.128
Node # 6 CPU: 1.116
Node # 7 CPU: 1.115
Node # 11 CPU: 1.137
Node # 12 CPU: 1.137
Node # 14 CPU: 1.137


Node # 2 CPU: 1.141
ROMS/TOMS - Output NetCDF summary for Grid 01:

Analytical header files used:

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

ROMS/TOMS - Input error ............. exit_flag: 2


ERROR: Abnormal termination: NetCDF INPUT.
REASON: No error

But I can't understand this erro message.

[rtopper]

The input file for zeta_west doesn't contain data for time = 55927.0000 (your start time) as indicated in your output:

SET_NGFLD - current model time exceeds ending value for variable: zeta_west
TDAYS = 55927.0000
Data Tmin = 55927.0000 Data Tmax = 56292.0000 <-- required input data time range (start time + 365 days)
Data Tstr = 55926.0000 Data Tend = 55927.0000 <-- time range of your input file

Check the input file that contains zeta_west. Make sure that it contains data for days 55927 - 56292.

[kobl1201]

when i checked zeta_west, this time range is 55927 to 56292.

But this model don't run.

How can i do?? please help me.

[salik]

It looks like you're using a revised version of varinfo.dat. Maybe check to see if the expected name for the zeta time variable is zeta_time?

In my file, the relevant piece of code is in line 1005 and looks like this:

Code: Select all

'zeta_west'                                        ! Input
  'free-surface western boundary condition'
  'meter'                                          ! [m]
  'zeta_west, scalar, series'
  'zeta_time'
  'idZbry(iwest)'
  'nulvar'
  1.0d0

If yours has a different value than "zeta_time", the code won't read your time variable from the input file. I had this exact problem when setting up on a new system recently