!********+*********+*********+*********+*********+*********+*********+*
!*   This file determines the adjustable run-time parameters for the  *
!* Sea Ice Simulator (SIS), versions 2 and later.  Where              *
!* appropriate, MKS units are used.                                   *
!********+*********+*********+*********+*********+*********+*********+*
! >>> Global coupled ice-ocean SIS2 5 layer ice, 63 layer ocean test case.

!  Specify properties of the physical domain.
OMEGA = 7.2921e-5               !   The rotation rate of the earth in s-1.
ROTATION = "2omegasinlat"       !   Rotation approximation (2omegasinlat,
                                !  betaplane or USER).
G_EARTH = 9.80                  !    G_EARTH is the Earth's gravitational
                                !  acceleration, in m s-2.
GRID_CONFIG = "mosaic"          !    Method for defining horizontal grid
                                !  = file|cartesian|spherical|mercator
                                !  file - read grid from file "GRID_FILE"
                                !  cartesian - a Cartesian grid
                                !  spherical - a spherical grid
                                !  mercator  - a Mercator grid
GRID_FILE = "ocean_hgrid.nc"
                                !    Name of file to read horizontal grid data
INPUTDIR = "INPUT"
                                !    INPUTDIR is a directory in which NetCDF
                                !  input files might be found.

!  Specify the numerical domain.
NIGLOBAL = 120                  !    NIGLOBAL and NJGLOBAL are the number of thickness
NJGLOBAL = 84                   !  grid points in the zonal and meridional
                                !  directions of the physical domain.
NCAT_ICE = 5                    !    The number of sea ice thickness categories.
NIHALO = 1                      !    NIHALO and NJHALO are the number of halo
NJHALO = 1                      !  points on each side in the x- and y-directions.
                                !  If static memory allocation is used, these
                                !  must match the values of NIHALO_ and NJHALO_
                                !  in MOM_memory.h.
IO_LAYOUT = 1, 1                !    The processor layout to be used, or 0,0 to automatically
                                !  set the io_layout to be the same as the layout.

REENTRANT_X = True              !    If defined, the domain is zonally reentrant.
REENTRANT_Y = False             !    If defined, the domain is meridionally
                                !  reentrant.
TRIPOLAR_N = True               !    Use tripolar connectivity at the northern
                                !  edge of the domain.  With TRIPOLAR_N, NIGLOBAL
                                !  must be even.

DEBUG = False                   !    If true, write out verbose debugging data.
SEND_LOG_TO_STDOUT = False      !    If true write out log information to stdout.
ICE_STATS_INTERVAL = 7          !    The interval in units of TIMEUNIT between
                                ! writes of the globally summed ice statistics
                                ! and conservation checks.
EULERIAN_TSURF = False          ! If true, use previous calculations of the ice-top surface
                                ! skin temperature for tsurf at the start of atmospheric
                                ! time stepping, including interpolating between tsurf
                                ! values from other categories in the same location.
RECATEGORIZE_ICE = False        ! If true, readjust the distribution into ice thickness
                                ! categories after advection.


!  Specify the time integration scheme.
DT_ICE_DYNAMICS = 7200.0        ! The time step used for the slow ice dynamics, 
                                ! including stepping the continuity equation and
                                ! interactions between the ice mass field and 
                                ! velocities, in s.
DT_RHEOLOGY = 100.0             ! The sub-cycling time step, in s, for iterating
                                ! the rheology and ice momentum equations. If
                                ! DT_RHEOLOGY is negative, the time step is set
                                ! via NSTEPS_DYN.

!  Specify the ice properties.
ICE_STRENGTH_PSTAR = 2.75e4     !    ICE_STRENGTH_PSTAR is a constant in the
                                !  expression for the ice strength, P* in
                                !  Hunke & Dukowics '97, in Pa.
RHO_OCEAN = 1030.0              !    RHO_OCEAN is used in the Boussinesq
                                !  approximation to calculations of pressure and
                                !  pressure gradients, in units of kg m-3.
RHO_ICE = 905.0                 !    RHO_ICE is the nominal density of sea ice
                                !  used in SIS, in units of kg m-3.
RHO_SNOW = 330.0                !    RHO_SNOW is the nominal density of snow
                                !  used in SIS, in units of kg m-3.
CP_ICE = 2100.0                 !   [J kg-1 K-1] default = 2100.0
                                ! The heat capacity of fresh ice, approximated as a
                                ! constant.
CP_SEAWATER = 4200.0            !   [J kg-1 K-1] default = 4200.0
                                ! The heat capacity of sea water, approximated as a
                                ! constant.
CP_BRINE = 2100.0               !   [J kg-1 K-1] default = 2100.0
                                ! The heat capacity of water in brine pockets within the
                                ! sea-ice, approximated as a constant.  CP_BRINE and
                                ! CP_WATER should be equal, but for computational
                                ! convenience CP_BRINE has often been set equal to CP_ICE.
! SNOW_ALBEDO = 0.88              !    The albedo of dry snow atop sea ice, ND.
! ICE_ALBEDO  = 0.68              !    The albedo of dry bare sea ice, ND.
ICE_BULK_SALINITY = 5.0         !    The fixed bulk salinity of sea ice, in g/kg.
SIS2_FILLING_FRAZIL = False     !   If true, apply frazil to fill as many categories as
                                ! possible to fill in a uniform (minimum) amount of ice
                                ! in all the thinnest categories. Otherwise the frazil is
                                ! always assigned to a single category.

CONSTANT_COSZEN_IC = 0.0        ! A constant value to use to initialize the
                                ! cosine of the solar zenith angle for the first
                                ! radiation step, or a negative number to use
                                ! the current time and astronomy.

!  Specify the properties of the surface forcing.
ADD_DIURNAL_SW = True           !    If true, add a synthetic diurnal cycle to
                                !  the shortwave radiation.

!  These are used for debugging the SIS_input files.
REPORT_UNUSED_PARAMS = True     !   [Boolean] default = False
                                ! If true, report any parameter lines that are not used
                                ! in the run.
FATAL_UNUSED_PARAMS = True      !   [Boolean] default = False
                                ! If true, kill the run if there are any unused
                                ! parameters.