!=======================================================================
!BOP
!
! !MODULE: ice_constants - sets physical constants
!
! !DESCRIPTION:
!
! This module defines a variety of physical and numerical constants
! used throughout the ice model \\
!
! Code originally based on constants.F in POP
!
! !REVISION HISTORY:
!  SVN:$Id: ice_constants.F90 37 2006-11-29 18:06:44Z eclare $
!
! author Elizabeth C. Hunke, LANL
!
! !INTERFACE:


      module ice_constants 39,2
!
! !USES:
!
      use shr_const_mod
      use ice_kinds_mod
!
!EOP
!
      implicit none
      save

      !-----------------------------------------------------------------
      ! physical constants
      !-----------------------------------------------------------------

      real (kind=dbl_kind), parameter :: &
         pi        = SHR_CONST_PI    ,&! pi
         gravit    = SHR_CONST_G     ,&! gravitational acceleration (m/s^2)
         secday    = SHR_CONST_CDAY  ,&! seconds in calendar day
         omega     = SHR_CONST_OMEGA ,&! angular velocity of earth (rad/sec)
         radius    = SHR_CONST_REARTH,&! earth radius (m)
         rhos      = 330.0_dbl_kind  ,&! density of snow (kg/m^3)
         rhoi      = SHR_CONST_RHOICE,&! density of ice (kg/m^3)
         rhow      = SHR_CONST_RHOSW ,&! density of seawater (kg/m^3)
         rhofresh  = SHR_CONST_RHOFW ,&! density of fresh water (kg/m^3)
         zvir      = SHR_CONST_ZVIR  ,&! rh2o/rair - 1.0
         vonkar    = SHR_CONST_KARMAN,&! von Karman constant
         cp_air    = SHR_CONST_CPDAIR,&! specific heat of air (J/kg/K)
         cp_wv     = SHR_CONST_CPWV  ,&! specific heat of water vapor (J/kg/K)
         ! (Briegleb JGR 97 11475-11485  July 1992)
         emissivity = 0.95_dbl_kind  ,&! emissivity of snow and ice
         stefan_boltzmann = SHR_CONST_STEBOL,&!  W/m^2/K^4
         Tffresh   = SHR_CONST_TKFRZ ,&! freezing temp of fresh ice (K)
         cp_ice    = SHR_CONST_CPICE ,&! specific heat of fresh ice (J/kg/K)
         cp_ocn    = SHR_CONST_CPSW  ,&! specific heat of ocn    (J/kg/K)
         depressT  = 0.054_dbl_kind  ,&! Tf:brine salinity ratio (C/ppt)
         Lsub      = SHR_CONST_LATSUB,&! latent heat, sublimation freshwater (J/kg)
         Lvap      = SHR_CONST_LATVAP,&! latent heat, vaporization freshwater (J/kg)
         Lfresh    = SHR_CONST_LATICE,&! latent heat of melting of fresh ice (J/kg)
         Timelt    = SHR_CONST_TKFRZ-SHR_CONST_TKFRZ,&! melting temp. ice top surface  (C)
         Tsmelt    = SHR_CONST_TKFRZ-SHR_CONST_TKFRZ,&! melting temp. snow top surface (C)
         ice_ref_salinity = SHR_CONST_ICE_REF_SAL ,&! (psu)
!        ocn_ref_salinity = SHR_CONST_OCN_REF_SAL ,&! (psu)
         albocn = 0.06_dbl_kind                   ,&! ocean albedo
         dragio = 0.00536_dbl_kind                ,&! ice-ocn drag coefficient
!        rho_air   = SHR_CONST_RHODAIR,&! ambient air density (kg/m^3)
         spval_dbl = SHR_CONST_SPVAL  ,&! special value
         snowpatch = 0.005_dbl_kind     ! parameter for fractional snow area (m)

      real (kind=real_kind), parameter :: &
         spval     = 1.0e30_real_kind   ! special value for netCDF output

      real (kind=dbl_kind), parameter :: &
         iceruf   = 0.0005_dbl_kind   ,&! ice surface roughness (m)

         ! (Ebert, Schramm and Curry JGR 100 15965-15975 Aug 1995)
         kappav = 1.4_dbl_kind ,&! vis extnctn coef in ice, wvlngth<700nm (1/m)
         kappan = 17.6_dbl_kind,&! vis extnctn coef in ice, wvlngth<700nm (1/m)

         kice   = 2.03_dbl_kind  ,&! thermal conductivity of fresh ice(W/m/deg)
         kseaice= 2.00_dbl_kind  ,&! thermal conductivity of sea ice (W/m/deg)
                                   ! (used in zero layer thermodynamics option)
         ksno   = 0.30_dbl_kind  ,&! thermal conductivity of snow  (W/m/deg)
         zref   = 10._dbl_kind   ,&! reference height for stability (m)
         Tocnfrz= -34.0_dbl_kind*depressT  ! freezing temp of seawater (C),
                                           ! used as Tsfcn for open water

      ! weights for albedos 
! 4 Jan 2007 BPB  Following are appropriate for complete cloud
! in a summer polar atmosphere with 1.5m bare sea ice surface:
! .636/.364 vis/nir with only 0.5% direct for each band.
      real (kind=dbl_kind), parameter :: &           ! currently used only
         awtvdr = 0.00318_dbl_kind, &! visible, direct  ! for history and
         awtidr = 0.00182_dbl_kind, &! near IR, direct
         awtvdf = 0.63282_dbl_kind, &! visible, diffuse  ! diagnostics
         awtidf = 0.36218_dbl_kind   ! near IR, diffuse

      real (kind=dbl_kind), parameter :: &
         hs0   = 0.03_dbl_kind,     &! parameter for delta-Eddington snow frac
         hsmin = 0.0001_dbl_kind     ! minimum snow thickness for dEdd

      real (kind=dbl_kind), parameter :: &
         qqqice  = 11637800._dbl_kind   ,&! for qsat over ice
         TTTice  = 5897.8_dbl_kind      ,&! for qsat over ice
         qqqocn  = 627572.4_dbl_kind    ,&! for qsat over ocn
         TTTocn  = 5107.4_dbl_kind        ! for qsat over ocn

      ! these are currently set so as to have no effect on the decomposition
      real (kind=dbl_kind), parameter :: &
         shlat  =  30.0_dbl_kind   ,&! artificial masking edge (deg)
         nhlat  = -30.0_dbl_kind     ! artificial masking edge (deg)
   
      !-----------------------------------------------------------------
      ! numbers
      !-----------------------------------------------------------------

      real (kind=dbl_kind), parameter :: &
        c0   = 0.0_dbl_kind, &
        c1   = 1.0_dbl_kind, &
        c1p5 = 1.5_dbl_kind, &
        c2   = 2.0_dbl_kind, &
        c3   = 3.0_dbl_kind, &
        c4   = 4.0_dbl_kind, &
        c5   = 5.0_dbl_kind, &
        c6   = 6.0_dbl_kind, &
        c7   = 7.0_dbl_kind, &
        c8   = 8.0_dbl_kind, &
        c9   = 9.0_dbl_kind, &
        c10  = 10.0_dbl_kind, &
        c12  = 12.0_dbl_kind, &
        c15  = 15.0_dbl_kind, &
        c16  = 16.0_dbl_kind, &
        c20  = 20.0_dbl_kind, &
        c25  = 25.0_dbl_kind, &
        c90  = 90.0_dbl_kind, &
        c100 = 100.0_dbl_kind, &
        c180 = 180.0_dbl_kind, &
        c360 = 360.0_dbl_kind, &
        c365 = 365.0_dbl_kind, &
        c3600= 3600.0_dbl_kind, &
        c1000= 1000.0_dbl_kind, &
        p001 = 0.001_dbl_kind, &
        p01  = 0.01_dbl_kind, &
        p1   = 0.1_dbl_kind, &
        p2   = 0.2_dbl_kind, &
        p4   = 0.4_dbl_kind, &
        p5   = 0.5_dbl_kind, &
        p6   = 0.6_dbl_kind, &
        p05  = 0.05_dbl_kind, &
        p15  = 0.15_dbl_kind, &
        p25  = 0.25_dbl_kind, &
        p75  = 0.75_dbl_kind, &
        p166 = c1/c6, &
        p333 = c1/c3, &
        p666 = c2/c3, &
        p111 = c1/c9, &
        p055 = p111*p5, &
        p027 = p055*p5, &
        p222 = c2/c9, &
        eps04  = 1.0e-4_dbl_kind, &
        eps11  = 1.0e-11_dbl_kind, &
        eps12  = 1.0e-12_dbl_kind, &
        eps13  = 1.0e-13_dbl_kind, &
        eps15  = 1.0e-15_dbl_kind, &
        eps16  = 1.0e-16_dbl_kind, &
        puny   = eps11, &
        bignum = 1.0e+30_dbl_kind, &
        pih    = p5*pi, &
        pi2    = c2*pi

      !-----------------------------------------------------------------
      ! location of fields for staggered grids
      !-----------------------------------------------------------------

      integer (int_kind), parameter :: &   
        field_loc_unknown  =  0, & 
        field_loc_noupdate = -1, & 
        field_loc_center   =  1, & 
        field_loc_NEcorner =  2, & 
        field_loc_Nface    =  3, & 
        field_loc_Eface    =  4, &
        field_loc_Wface    =  5


      !-----------------------------------------------------------------
      ! field type attribute - necessary for handling
      ! changes of direction across tripole boundary
      !-----------------------------------------------------------------

      integer (int_kind), parameter :: &   
        field_type_unknown  =  0, & 
        field_type_noupdate = -1, & 
        field_type_scalar   =  1, & 
        field_type_vector   =  2, & 
        field_type_angle    =  3

      !-----------------------------------------------------------------
      ! conversion factors
      !-----------------------------------------------------------------

      real (kind=dbl_kind), parameter :: &
        cm_to_m       = 0.01_dbl_kind   ,&! cm to meters
        m_to_cm       = 100._dbl_kind   ,&! meters to cm
        m2_to_km2     = 1.e-6_dbl_kind  ,&! m^2 to km^2
        kg_to_g       = 1000._dbl_kind  ,&! kilograms to grams
        mps_to_cmpdy  = 8.64e6_dbl_kind ,&! m per s to cm per day
        rad_to_deg    = 180._dbl_kind/pi  ! degree-radian conversion

#ifndef USE_ESMF
      integer (kind=int_kind), parameter :: &
         ESMF_SUCCESS = 0   ! otherwise ESMF defines this parameter
#endif

      ! useful for debugging
      integer (kind=int_kind), parameter :: &
         mtest = -999, itest = 1, jtest = 1, ntest = 1, btest = 1
!         mtest = 2, itest = 50, jtest = 53, ntest = 1, btest = 1

!=======================================================================

      end module ice_constants

!=======================================================================