!
!  ROMS/TOMS Standard Input parameters.
!
!==============================================================================
!  Copyright (c) 2002 ROMS/TOMS Group, version 2.1                          ===
!========================================================= Hernan G. Arango ===
!                                                                             !
! Input parameters can be entered in ANY order, provided that the parameter   !
! KEYWORD (usually, upper case) is typed correctly followed by "="  or "=="   !
! symbols. Any comment lines are allowed and must begin with an exclamation   !
! mark (!) in column one.  Comments may  appear to the right of a parameter   !
! specification to improve documentation.  All comments will ignored during   !
! reading.  Blank lines are also allowed and ignored. Continuation lines in   !
! a parameter specification are allowed and must be preceded by a backslash   !
! (\).  In some instances, more than one value is required for a parameter.   !
! If fewer values are provided, the  last value  is assigned for the entire   !
! parameter array.  The multiplication symbol (*),  without blank spaces in   !
! between, is allowed for a parameter specification.  For example, in a two   !
! grids nested application:                                                   !
!                                                                             !
!    AKT_BAK == 2*1.0d-6  2*5.0d-6              ! m2/s                        !
!                                                                             !
! indicates that the first two entries of array AKT_BAK,  in fortran column-  !
! major order, will have the same value of "1.0d-6" for grid 1,  whereas the  !
! next two entries will have the same value of "5.0d-6" for grid 2.           !
!                                                                             !
! In multiple levels of nesting and/or multiple connected domains  step-ups,  !
! "Ngrids" entries are expected for some of these parameters.  In such case,  !
! the order of the entries for a parameter is extremely important.  It  must  !
! follow the same order (1:Ngrids) as in the state variable declaration. The  !
! USER may follow the above guidelines for specifying his/her values.  These  !
! parameters are marked by "==" plural symbol after the KEYWORD.              !
!                                                                             !
!==============================================================================
!

! Application title.

       TITLE = ROMS/TOMS 2.2 - Chesapeake Bay Coarse grid application

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

     VARNAME = /user1/jiangtao/ChesROMS/roms-2.2/External/varinfo.dat

! 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 == 8                              ! J-direction partition

! Time-Stepping parameters.

     NTIMES = 525600       !86400    !216000    ! 150 days 525600
          DT == 60.0d0
     NDTFAST == 20
       ERstr =  1
       ERend =  1

! Input/Output parameters.

       NRREC =  0
      LCYCLE == T
        NRST == 1440  !daily
        NSTA == 60    !hourly
        NFLT == 12
       NINFO == 1

! Output history and average files parameters.

     LDEFOUT == T
        NHIS == 1440
     NDEFHIS == 43200
      NTSAVG == 1
        NAVG == 180   !3 hourly
     NDEFAVG == 21600
      NTSDIA == 1
        NDIA == 1440
     NDEFDIA == 0

! Harmonic/biharmonic horizontal diffusion of tracer: [1:NAT,Ngrids].

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

! Harmononic/biharmonic, horizontal viscosity coefficient: [Ngrids].

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

! Vertical mixing coefficients for active tracers: [1:NAT,Ngrids]

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

! Vertical mixing coefficient for momentum: [Ngrids].

     AKV_BAK == 1.0d-5                          ! m2/s

! 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 == 0.d0                           !k-kl 
       GLS_M == 1.0d0
       GLS_N == 1.0d0
    GLS_Kmin == 5.0d-6
    GLS_Pmin == 1.0d-8

    GLS_CMU0 == 0.5544d0
      GLS_C1 == 0.9d0
      GLS_C2 == 0.5d0
     GLS_C3M == 2.53d0
     GLS_C3P == 1.0d0
    GLS_SIGK == 2.44d0
    GLS_SIGP == 2.44d0

! Constants used in momentum stress computation.
                                     
        RDRG == 0.0d0                      ! m/s
       RDRG2 == 3.0d-3                     ! nondimensional
         Zob == 0.0005d0                      ! m
         Zos == 0.02d0                     ! m

! Various parameters.

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

! Vertical S-coordinates parameters, [1:Ngrids].

     THETA_S == 6.0d0                      ! 
     THETA_B == 0.95d0                     ! 
      TCLINE == 10.0d0                     ! m
 
! Mean Density and time stamp.

        RHO0 =  1025.0d0                   ! kg/m3
      DSTART =  6940.0d0                    ! days                 
    TIME_REF =  19800101.0d0               ! yyyymmdd.dd

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

       TNUDG == 1.0d0  1.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 == 24.0d0                     ! nondimensional

! Linear equation of State parameters:

          R0 == 0.0d0                      ! kg/m3
          T0 == 15.0d0                     ! deg C
          S0 == 35.0d0                     ! PSU
       TCOEF == -0.28d0
       SCOEF ==  0.78d0

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

      GAMMA2 =  1.0d0

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

Hout(idUvel) == T                          ! 3D U-velocity
Hout(idVvel) == T                          ! 3D V-velocity
Hout(idWvel) == T                          ! 3D W-velocity
Hout(idOvel) == F                          ! omega vertical velocity
Hout(idUbar) == T                          ! 2D U-velocity
Hout(idVbar) == T                          ! 2D V-velocity
Hout(idFsur) == T                          ! free-surface

Hout(idTvar) == T T T                      ! temperature and salinity

Hout(idUsms) == T                          ! surface U-stress
Hout(idVsms) == T                          ! surface V-stress
Hout(idUbms) == T                          ! bottom U-stress
Hout(idVbms) == T                          ! bottom V-stress
Hout(idUbws) == F                          ! bottom U-wave stress
Hout(idVbws) == F                          ! bottom V-wave stress
Hout(idUbcs) == F                          ! bottom max wave-current U-stress
Hout(idVbcs) == F                          ! bottom max wave-current V-stress

Hout(idUbed) == F                          ! bed wave orbital U-velocity
Hout(idVbed) == F                          ! bed wave orbital V-velocity
Hout(idUbot) == F                          ! bottom U-momentum above bed
Hout(idVbot) == F                          ! bottom V-momentum above bed

Hout(idTsur) == T                          ! surface net heat flux
Hout(idLhea) == F                          ! latent heat flux
Hout(idShea) == F                          ! sensible heat flux
Hout(idLrad) == F                          ! longwave radiation flux
Hout(idSrad) == F                          ! shortwave radiation flux

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

! 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:
!
!   Hout(idBott(isd50)),  isd50 = 1        ! mean grain diameter
!   Hout(idBott(idens)),  idens = 2        ! mean grain density
!   Hout(idBott(iwsed)),  iwsed = 3        ! mean settling velocity
!   Hout(idBott(itauc)),  itauc = 4        ! critical erosion stress
!   Hout(idBott(irlen)),  irlen = 5        ! ripple length
!   Hout(idBott(irhgt)),  irhgt = 6        ! ripple height
!   Hout(idBott(ibwav)),  ibwav = 7        ! wave excursion amplitude
!   Hout(idBott(izdef)),  izdef = 8        ! default bottom roughness
!   Hout(idBott(izapp)),  izapp = 9        ! apparent bottom roughness
!   Hout(idBott(izNik)),  izNik = 10       ! Nikuradse bottom roughness
!   Hout(idBott(izbio)),  izbio = 11       ! biological bottom roughness
!   Hout(idBott(izbfm)),  izbfm = 12       ! bed form bottom roughness
!   Hout(idBott(izbld)),  izbld = 13       ! bed load bottom roughness
!   Hout(idBott(izwbl)),  izwbl = 14       ! wave bottom roughness
!   Hout(idBott(iactv)),  iactv = 15       ! active layer thickness
!   Hout(idBott(ishgt)),  ishgt = 16       ! saltation height

Hout(idBott) == T T T T T T T T T F F F F F F F

! Generic User parameters, [1:NUSER].

       NUSER =  0
        USER =  0.d0

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

     GRDNAME == chesroms_grd_smoother_msl.nc
     ININAME == chesroms_ini_1999.nc
     CLMNAME == ocean_clm.nc
     BRYNAME == chesroms_forcing_obc_1999m_ndw.nc
     FWDNAME == ocean_fwd.nc

! Input forcing NetCDF file name(s).  The USER has the option to enter
! several files names per each nested grid.  For example, the USER may
! have a different files for wind products, heat fluxes, rivers, 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 multiple forcing
! files per grid, enter first all the file names for grid 1, then grid 2,
! and so on.  Use a single line per entry with a continuation (\) symbol
! at the each entry, except the last one.

     NFFILES == 3                          ! number of forcing files

     FRCNAME == chesroms_forcing_rivers_wrivtemp_1999m_CnD.nc \
                chesroms_forcing_tides9_1999_S2N2K1adjust.nc  \
                chesroms_forcing_surface_1999.nc

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

     RSTNAME == chesroms_rst_Nz30.nc
     HISNAME == chesroms_his_Nz30.nc
     AVGNAME == chesroms_avg_Nz30.nc
     DIANAME == chesroms_dia_Nz30.nc
     STANAME == chesroms_sta_Nz30.nc
     FLTNAME == chesroms_flt_Nz30.nc

! Input ASCII parameter filenames.

     APARNAM =  assimilation.in
     SPOSNAM =  stations.in
     FPOSNAM =  floats.in
     BPARNAM =  bioFasham.in
     SPARNAM =  sediment.in
     USRNAME =  MyFile.dat

!
!  GLOSSARY:
!  =========
!
!------------------------------------------------------------------------------
! Application tile (string with a maximum of eighty characters).
!------------------------------------------------------------------------------
!
!  TITLE     Title of model Run.
!
!------------------------------------------------------------------------------
! Variable information file name (string with a maximum of eighty characters).
!------------------------------------------------------------------------------
!
!  VARNAME   Input/Output variable information file name.  This file needd to
!              be processed first so all information arrays and indices can be
!              initialized properly in "mod_ncparam.F".
!
!------------------------------------------------------------------------------
! Domain tile partition parameters.
!------------------------------------------------------------------------------
!
! Model tile decomposition parameters for serial and parallel configurations
! which are used to determine tile horizontal range indices (Istr,Iend) and
! (Jstr,Jend). In some computers, it is advantageous to have tile partitions
! in serial applications.
!
!  NtileI    Number of domain partitions in the I-direction (XI-coordinate).
!              It must be equal or greater than one.
!
!  NtileJ    Number of domain partitions in the J-direction (ETA-coordinate).
!              It must be equal or greater than one.
!
!  WARNING:  In shared-memory (OpenMP), the product of NtileI and NtileJ must
!            be a MULTIPLE of the number of parallel threads specified with
!            the OpenMP environmental variable OMP_NUM_THREADS.
!
!            In distributed-memory (MPI), the product of NtileI and NtileJ
!            must be EQUAL to the number of parallel nodes specified during
!            execution with the "mprun" or "mpirun" command.
!
!------------------------------------------------------------------------------
! Time-Stepping parameters.
!------------------------------------------------------------------------------
!
!  NTIMES    Total number time-steps in current run.  If 3D configuration,
!            NTIMES is the total of baroclinic time-steps.  If only 2D
!            configuration, NTIMES is the total of barotropic time-steps.
!
!  DT        Time-Step size in seconds.  If 3D configuration, DT is the
!            size of baroclinic time-step.  If only 2D configuration, DT
!            is the size of the barotropic time-step.
!
!  NDTFAST   Number of barotropic time-steps between each baroclinic time-step.
!            If only 2D configuration, NDTFAST should be unity since there is
!            not need to splitting time-stepping.
!
!  ERstr     Starting ensemble run (perturbation or iteration) number.
!
!  ERend     Ending   ensemble run (perturbation or iteration) number.
!
!------------------------------------------------------------------------------
! Input/Output parameters.
!------------------------------------------------------------------------------
!
!  NRREC     Switch to indicate re-start from a previous solution.  Use
!            NRREC=0 for new solutions. In a re-start solution, NRREC
!            is the time index of the re-start NetCDF file assigned for
!            initialization.  If NRREC is negative (said NRREC=-1), the
!            model will re-start from the most recent time record. That
!            is, the initialization record is assigned internally.
!            Notice that it is also possible to re-start from a history
!            or time-averaged NetCDF files.  If a history file is used
!            for re-start, it must contains all the necessary primitive
!            variables at all levels.
!
!  LCYCLE    Logical switch (T/F) used to recycle time records in output
!            re-start file.  If TRUE,  only the latest two re-start time
!            records are maintained.  If FALSE, all re-start fields are
!            saved every NRST time-steps without recycling.  The re-start
!            fields are written at all levels in double precision.
!
!  NRST      Number of time-steps between writing of re-start fields.
!
!  NSTA      Number of time-steps between writing data into stations file.
!            Station data is written at all levels.
!
!  NFLT      Number of time-steps between writing data into floats file.
!
!  NINFO     Number of time-steps between print of single line information
!            to standard output.  If also determines the interval between
!            computation of global energy diagnostics.
!
!------------------------------------------------------------------------------
!  Output history and average files parameters.
!------------------------------------------------------------------------------
!
!  LDEFOUT   Logical switch (T/F) used to create new output files when
!            initializing from a re-start file, abs(NRREC) > 0.  If TRUE
!            and applicable, a new history, average, diagnostic and
!            station files are created during the initialization stage.
!            If FALSE and applicable, data is appended to an existing
!            history, average, diagnostic and station files.  See also
!            parameters NDEFHIS, NDEFAVG and NDEFDIA below.
!
!  NHIS      Number of time-steps between writing fields into history file.
!
!  NDEFHIS   Number of time-steps between the creation of new history file.
!            If NDEFHIS=0, the model will only process one history file.
!            This feature is useful for long simulations when history files
!            get too large; it creates a new file every NDEFHIS time-steps.
!
!  NTSAVG    Starting time-step for the accumulation of output time-averaged
!            data.
!
!  NAVG      Number of time-steps between writing time-averaged data
!            into averages file.  Averaged date is written for all fields.
!
!  NDEFAVG   Number of time-steps between the creation of new average
!            file.  If NDEFAVG=0, the model will only process one average
!            file.  This feature is useful for long simulations when
!            average files get too large; it creates a new file every
!            NDEFAVG time-steps.
!
!  NTSDIA    Starting time-step for the accumulation of output time-averaged
!            diagnostics data.
!
!  NDIA      Number of time-steps between writing time-averaged diagnostics
!            data into diagnotics file.  Averaged date is written for all
!            fields.
!
!  NDEFDIA   Number of time-steps between the creation of new time-averaged
!            diagnostics file.  If NDEFDIA=0, the model will only process one
!            diagnostics file.  This feature is useful for long simulations
!            when diagnostics files get too large; it creates a new file
!            every NDEFDIA time-steps.
!
!------------------------------------------------------------------------------
! Harmonic/Biharmonic horizontal diffusion for active tracers.
!------------------------------------------------------------------------------
!
!  TNU2      Lateral, Laplacian, constant, mixing coefficient (m2/s) for
!            active tracer variables.  If variable horizontal diffusion is
!            activated, TNU2 is the mixing coefficient for the largest
!            grid-cell in the domain.
!
!  TNU4      Lateral, biharmonic, constant, mixing coefficient (m4/s) for
!            active tracer variables.  If variable horizontal diffusion is
!            activated, TNU4 is the mixing coefficient for the largest
!            grid-cell in the domain.
!
!------------------------------------------------------------------------------
! Harmonic/biharmonic horizontal viscosity coefficients.
!------------------------------------------------------------------------------
!
!  VISC2     Lateral, Laplacian, constant, mixing coefficient (m2/s) for
!            momentum.  If variable horizontal viscosity is activated, UVNU2
!            is the mixing coefficient for the largest grid-cell in the
!            domain.
!
!  VISC4     Lateral, biharmonic, constant mixing coefficient (m4/s) for
!            momentum. If variable horizontal viscosity is activated, UVNU4
!            is the mixing coefficient for the largest grid-cell in the
!            domain.
!
!------------------------------------------------------------------------------
! Vertical mixing coefficients for active tracers.
!------------------------------------------------------------------------------
!
!  AKT_BAK   Background vertical mixing coefficient (m2/s) for active
!            tracer variables.
!
!------------------------------------------------------------------------------
! Vertical mixing coefficient for momentum.
!------------------------------------------------------------------------------
!
!  AKV_BAK   Background vertical mixing coefficient (m2/s) for momentum.
!
!------------------------------------------------------------------------------
! Turbulent closure parameters.
!------------------------------------------------------------------------------
!
!  AKK_BAK   Background vertical mixing coefficient (m2/s) for turbulent
!            kinetic energy.
!
!  AKP_BAK   Background vertical mixing coefficient (m2/s) for turbulent
!            generic statistical field, "psi".
!
!  TKENU2    Lateral, harmonic, constant, mixing coefficient (m2/s) for
!            turbulent closure variables.
!
!  TKENU4    Lateral, biharmonic, constant mixing coefficient (m4/s) for
!            turbulent closure variables.
!
!------------------------------------------------------------------------------
! Generic length-scale turbulence closure parameters.
!------------------------------------------------------------------------------
!
!  GLS_P     Stability exponent (non-dimensional).
!
!  GLS_M     Turbulent kinetic energy exponent (non-dimensional).
!
!  GLS_N     Turbulent length scale exponent (non-dimensional).
!
!  GLS_Kmin  Minimum value of specific turbulent kinetic energy
!
!  GLS_Pmin  Minimum Value of dissipation.
!
! Closure independent constraint parameters (non-dimensional):
!
!  GLS_CMU0  Stability coefficient.
!
!  GLS_C1    Shear production coefficient.
!
!  GLS_C2    Dissipation coefficient.
!
!  GLS_C3M   Buoyancy production coefficient (minus).
!
!  GLS_C3P   Buoyancy production coefficient (plus).
!
!  GLS_SIGK  Constant Schmidt number (non-dimensional) for turbulent
!            kinetic energy diffusivity.
!
!  GLS_SIGP  Constant Schmidt number (non-dimensional) for turbulent
!            generic statistical field, "psi".
!
! Suggested values for various parameterizations:
!
!              MY2.5         K-epsilon    K-omega      K-omega      K-tao
!
!      GLS_P = 0.d0          3.0d0       -1.0d0       -1.0d0       -3.0d0
!      GLS_M = 1.d0          1.5d0        0.5d0        0.5d0        0.5d0
!      GLS_N = 1.d0         -1.0d0       -1.0d0       -1.0d0        1.0d0
!   GLS_Kmin = 5.0d-6        7.6d-6       7.6d-6       7.6d-6       7.6d-6
!   GLS_Pmin = 5.0d-6        1.0d-12      1.0d-12      1.0d-12      1.0d-12
!
!   GLS_CMU0 = 0.5544d0      0.5477d0     0.5477d0     0.5477d0     0.5477d0
!     GLS_C1 = 0.9d0         1.44d0       0.555d0      0.52d0       0.173d0
!     GLS_C2 = 0.5d0         1.92d0       0.833d0      0.8d0        0.225d0
!    GLS_C3M = 0.9d0        -0.4d0       -0.6d0       -0.6d0        0.0d0
!    GLS_C3P = 0.9d0         1.0d0        1.0d0        1.0d0        0.0d0
!   GLS_SIGK = 1.96d0        1.0d0        2.0d0        2.0d0        1.46d0
!   GLS_SIGP = 1.96d0        1.30d0       2.0d0        2.0d0       10.8d0
!
!
!------------------------------------------------------------------------------
! Constants used in the computation of momentum stress.
!------------------------------------------------------------------------------
!
! Bottom stress is computed using either the linear bottom drag coefficient,
! quadratic bottom drag coefficient, or logarithmic formulation via bottom
! roughness.  Choose only one formulation and set the other constants to zero.
! The model selects the formulation by checking which constant is not zero.
! The order of selection is: Zob, RDRG2, and RDRG. See "set_vbc" routine for
! details.
!
!  RDRG      Linear bottom drag coefficient (m/s).
!
!  RDRG2     Quadratic bottom drag coefficient.
!
!  Zob       Bottom roughness (m).
!
!  Zos       Surface roughness (m).
!
!------------------------------------------------------------------------------
! Jerlow Water type.
!------------------------------------------------------------------------------
!
!  WTYPE     Jerlov water type: an integer value from 1 to 5.
!
!------------------------------------------------------------------------------
! Body-force parameters. Used when CPP option BODYFORCE is activated.
!------------------------------------------------------------------------------
!
!  LEVSFRC   Deepest level to apply surface momentum stress as a body-force.
!
!  LEVBFRC   Shallowest level to apply bottom momentum stress as a body-force.
!
!------------------------------------------------------------------------------
! Vertical S-coordinates parameters.
!------------------------------------------------------------------------------
!
!  THETA_S   S-coordinate surface control parameter, [0 < theta_s < 20].
!
!  THETA_B   S-coordinate bottom  control parameter, [0 < theta_b < 1].
!
!  TCLINE    Width (m) of surface or bottom boundary layer in which
!            higher vertical resolution is required during stretching.
!
!            WARNING:  Users need to experiment with these parameters. We
!                      have found out that the model goes unstable with
!                      high values of THETA_S.  In steep and very tall
!                      topography, it is recommended to use THETA_S < 3.0.
!
!------------------------------------------------------------------------------
! Mean Density.
!------------------------------------------------------------------------------
!
!  RHO0      Mean density (Kg/m3) used when the Boussinesq approximation
!            is inferred.
!
!------------------------------------------------------------------------------
! Time Stamps.
!------------------------------------------------------------------------------
!
!  DSTART    Time stamp assigned to model initialization (days).  Usually
!            a Calendar linear coordinate, like modified Julian Day.  For
!            Example:
!
!                     Julian Day = 1  for  Nov 25, 0:0:0 4713 BCE
!            modified Julian Day = 1  for  May 24, 0:0:0 1968  CE GMT
!
!            It is called truncated or modified Julian day because an offset
!            of 2440000 needs to be added.
!
!  TIME_REF  Reference time (yyyymmdd.f) used to compute relative time:
!            elapsed time interval since reference-time.  The "units"
!            attribute takes the form "time-unit since reference-time".
!            This parameter also provides information about the calendar
!            used:
!
!            if TIME_REF = -2, model time and DSTART are in modified Julian
!            days units.  The "units" attribute is:
!
!                    'time-units since 1968-05-23 00:00:00 GMT'
!
!            if TIME_REF = -1, model time and DSTART are in a calendar
!            with 360 days in every year (30 days each month).  The "units"
!            attribute is:
!
!                    'time-units since 0000-01-01 00:00:00'
!
!            if TIME_REF = 0, model time and DSTART are in a common year
!            calendar with 365.25 days.  The "units" attribute is:
!
!                    'time-units since 0000-01-01 00:00:00'
!
!            if TIME_REF > 0, model time and DSTART are the elapsed time
!            units since specified reference time.  For example,
!            TIME_REF=20020115.5 will yield the following attribute:
!
!                    'time-units since 2002-01-15 12:00:00'
!
!------------------------------------------------------------------------------
! Nudging/relaxation time scales, inverse scales will be computed internally.
!------------------------------------------------------------------------------
!
! When passive/active open boundary conditions are activated, these nudging
! values correspond to the passive (outflow) nudging time scales.
!
!  TNUDG     Nudging time scale (days) for active tracer variables.
!            (1:NAT,1:Ngrids) values are expected.
!
!  ZNUDG     Nudging time scale (days) for free-surface.
!
!  M2NUDG    Nudging time scale (days) for 2D momentum.
!
!  M3NUDG    Nudging time scale (days) for 3D momentum.
!
!  OBCFAC    Factor between passive (outflow) and active (inflow) open
!              boundary conditions.  The nudging time scales for the
!              active (inflow) conditions are obtained by multiplying
!              the passive values by OBCFAC. If OBCFAC > 1, nudging on
!              inflow is stronger than on outflow (recommended).
!
!------------------------------------------------------------------------------
! Linear equation of State parameters.
!------------------------------------------------------------------------------
!
!  R0        Background density value (Kg/m3) used in Linear Equation of
!            State.
!
!  T0        Background potential temperature (Celsius) constant.
!
!  S0        Background salinity (PSU) constant.
!
!  TCOEF     Thermal expansion coefficient in Linear Equation of State.
!
!  SCOEF     Saline contraction coefficient in Linear Equation of State.
!
!------------------------------------------------------------------------------
! Slipperiness parameter.
!------------------------------------------------------------------------------
!
!  GAMMA2    Slipperiness variable, either 1.0 (free slip) or -1.0 (no slip).
!
!------------------------------------------------------------------------------
! Logical switches (T/F) to activate writing of fields into HISTORY file.
!------------------------------------------------------------------------------
!
!  Hout(idUvel)  Write out 3D U-velocity component.
!  Hout(idVvel)  Write out 3D V-velocity component.
!  Hout(idWvel)  Write out 3D W-velocity component.
!  Hout(idOvel)  Write out 3D omega vertical velocity.
!  Hout(idUbar)  Write out 2D U-velocity component.
!  Hout(idVbar)  Write out 2D V-velocity component.
!  Hout(idFsur)  Write out free-surface.
!
!  Hout(idTvar)  Write out active tracers: temperature and salinity
!
!  Hout(idUsms)  Write out surface U-momentum stress.
!  Hout(idVsms)  Write out surface V-momentum stress.
!  Hout(idUbms)  Write out bottom  U-momentum stress.
!  Hout(idVbms)  Write out bottom  V-momentum stress.
!  Hout(idUbws)  Write out wind-induced, bottom U-wave stress.
!  Hout(idVbws)  Write out wind-induced, bottom V-wave stress.
!
!  Hout(idUbed)  Write out wind-induced, bed wave orbital U-velocity.
!  Hout(idVbed)  Write out wind-induced, bed wave orbital V-velocity.
!  Hout(idUbot)  Write out bottom U-momentum above bed.
!  Hout(idVbot)  Write out bottom V-momentum above bed.
!  Hout(idUbcs)  Write out bottom maximum wave and current U-stress.
!  Hout(idVbcs)  Write out bottom maximum wave and current U-stress.
!
!  Hout(idTsur)  Write out surface net heat flux.
!  Hout(idLhea)  Write out latent heat flux.
!  Hout(idShea)  Write out sensible heat flux.
!  Hout(idLrad)  Write out long-wave radiation flux.
!  Hout(idSrad)  Write out short-wave radiation flux.
!
!  Hout(idDano)  Write out density anomaly.
!  Hout(idVvis)  Write out vertical viscosity coefficient.
!  Hout(idTdif)  Write out vertical diffusion coefficient of temperature.
!  Hout(idSdif)  Write out vertical diffusion coefficient of salinity.
!  Hout(idHsbl)  Write out depth of oceanic surface boundary layer.
!  Hout(idHbbl)  Write out depth of oceanic bottom boundary layer.
!  Hout(idMtke)  Write out turbulent kinetic energy.
!  Hout(idMtls)  Write out turbulent kinetic energy times length scale.
!
!  Hout(idBott)  Write out exposed sediment layer properties, 1:MBOTP.
!
!------------------------------------------------------------------------------
! Generic User parameters.
!------------------------------------------------------------------------------
!
!  NUSER     Number of User parameters to consider (integer).
!  USER      Vector containing user parameters (real array).
!
!------------------------------------------------------------------------------
! Input/output NetCDF file names (string with a maximum of eighty characters).
!------------------------------------------------------------------------------
!
!  RSTNAME   Output re-start file name.
!  HISNAME   Output history file name.
!  AVGNAME   Output averages file name.
!  DIANAME   Output diagnostics file name.
!  STANAME   Output stations file name.
!  FLTNAME   Output floats file name.
!
!  GRDNAME   Input grid file name.
!  ININAME   Input initial conditions file name. It can be a re-start file.
!  FRCNAME   Input forcing fields file name.
!  CLMNAME   Input climatology fields file name.
!  BRYNAME   Input climatology data file name.
!  FWDNAME   Input forward solution fields file name.
!
!------------------------------------------------------------------------------
! Input ASCII parameters file names.
!------------------------------------------------------------------------------
!
!  APARNAM   Input assimilation parameters file name.
!  SPOSNAM   Input stations positions file name.
!  FPOSNAM   Input initial drifters positions file name.
!  BPARNAM   Input biological parameters file name.
!  SPARNAM   Input sediment transport parameters file name.
!  USRNAME   USER's input generic file name.
!
