pkg/thsice/THSICE_VARS.h

87ea84cac6 Jean*0001 #ifdef ALLOW_THSICE
                0002
                0003 C     !ROUTINE: THSICE_VARS.h
                0004 C -------------------------------
                0005 C   THSICE_VARS.h
                0006 C  variable for thermodynamics - Sea-Ice model
                0007 C -------------------------------
                0008
                0009 C-- COMMON /THSICE_VARS/  state variables of sea-ice model :
                0010 C   iceMask   :: sea-ice fraction: no ice=0, grid all ice 1  []
                0011 C   iceHeight :: depth of ice layer        [m]
                0012 C   snowHeight:: depth of snow layer       [m]
                0013 C   Tsrf      :: temperature at surface     [oC]
                0014 C   Tice1     :: temperature of ice layer 1 [oC]
                0015 C   Tice2     :: temperature of ice layer 2 [oC]
bd8a00c1fc Jean*0016 C   Qice1     :: enthalpy of ice layer 1  [J/kg]
                0017 C   Qice2     :: enthalpy of ice layer 2  [J/kg]
87ea84cac6 Jean*0018 C   snowAge   :: snow age                   [s]
                0019       COMMON /THSICE_VARS/
                0020      &       iceMask, iceHeight, snowHeight,
                0021      &       Tsrf, Tice1, Tice2,
                0022      &       Qice1, Qice2, snowAge
                0023
                0024       _RL iceMask(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0025       _RL iceHeight (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0026       _RL snowHeight(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0027       _RL Tsrf   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0028       _RL Tice1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0029       _RL Tice2  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0030       _RL Qice1  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0031       _RL Qice2  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0032       _RL snowAge(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0033
                0034 C-- COMMON /THSICE_FLUX/ hold fluxes needed for thsice model
                0035 C   sHeating :: surf heating left to melt snow or ice (= Atmos-Conduction)
                0036 C   flxCndBt :: heat flux conducted through the ice to bottom surface
a4e6fa7055 Jean*0037 C   siceAlb  :: area weighted sea-ice albedo           [0-1]
ce354ad541 Jean*0038 C   icAlbNIR :: seaice(+snow) surface Near Infra-Red Albedo
7269783f6f Jean*0039 C atmospheric fluxes (change along the time-stepping):
49c48d1136 Jean*0040 C   icFlxSW  :: short-wave heat flux (+=down) over sea-ice
7269783f6f Jean*0041 C               (downward SW / net SW @ surface / net SW below sea-ice)
                0042 C   icFlxAtm :: Atmospheric surf. heat flux over sea-ice [W/m2] (+=down)
                0043 C               (over sea-ice only / weighted by ice-fraction)
                0044 C   icFrwAtm :: fresh-water flux (E-P) from the atmosphere [kg/m2/s] (+=up)
                0045 C               ( ice Evap only / ice E-P / ice - ocean weighted E-P )
de836be2bc Jean*0046 C   adjustFrW :: global adjustment to surface fresh-water flux [kg/m2/s]
87ea84cac6 Jean*0047       COMMON / THSICE_FLUX /
                0048      &       sHeating, flxCndBt,
c8bafbe9ad Jean*0049      &       siceAlb, icAlbNIR,
de836be2bc Jean*0050      &       icFlxSW, icFlxAtm, icFrwAtm,
                0051      &       adjustFrW
87ea84cac6 Jean*0052
                0053       _RL sHeating(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0054       _RL flxCndBt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
a4e6fa7055 Jean*0055       _RL siceAlb (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
ce354ad541 Jean*0056       _RL icAlbNIR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
7269783f6f Jean*0057       _RL icFlxSW (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0058       _RL icFlxAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0059       _RL icFrwAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
de836be2bc Jean*0060       _RL adjustFrW
87ea84cac6 Jean*0061
f6439efa9f Jean*0062 C-- COMMON /THSICE_DYN_R/  variables used with sea-ice advection/diffusion
                0063 C   oceFWfx   :: fresh water flux to the ocean  [kg/m^2/s]
ba0b047096 Mart*0064 C   oceSflx   :: salt flux to the ocean         [g/m^2/s]
f6439efa9f Jean*0065 C   oceQnet   :: heat flux to the ocean         [W/m^2]
                0066 C---
                0067 C    Note :: when ice volume is too small to be kept, ice & snow is melt
                0068 C        and fresh water, salt and heat are returned to the ocean.
                0069 C        For now, those fluxes are stored separately and will try to find
                0070 C        out how to incorporate them more naturally in the usual forcing.
                0071 C---
                0072       COMMON /THSICE_DYN_R/
                0073      &       oceFWfx, oceSflx, oceQnet
                0074
                0075       _RL oceFWfx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0076       _RL oceSflx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0077       _RL oceQnet(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0078
87ea84cac6 Jean*0079 C-- COMMON / THSICE_OCEMXLAYER / oceanic mixed layer state
                0080 C   hOceMxL :: thickness   of the ocean mixed layer [m]
                0081 C   tOceMxL :: temperature in the ocean mixed layer [oC]
ba0b047096 Mart*0082 C   sOceMxL :: salinity    in the ocean mixed layer [g/kg]
87ea84cac6 Jean*0083 C   v2ocMxL :: velocity (square) in the mixed layer [m2/s2]
                0084       COMMON / THSICE_OCEMXLAYER /
                0085      &       hOceMxL, tOceMxL, sOceMxL, v2ocMxL
                0086
                0087       _RL hOceMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0088       _RL tOceMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0089       _RL sOceMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0090       _RL v2ocMxL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0091
                0092 #endif /* ALLOW_THSICE */
df4e8f7bcf Ed H*0093
                0094 CEH3 ;;; Local Variables: ***
                0095 CEH3 ;;; mode:fortran ***
                0096 CEH3 ;;; End: ***