source: branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_2/thd_ice_2.F90 @ 2590

MODULE thd_ice_2
#if defined key_lim2
   !!======================================================================
   !!                       ***  MODULE thd_ice_2  ***
   !! LIM 2.0 sea-ice :   Ice thermodynamics in 1D
   !!=====================================================================
   !! History :
   !!   2.0  !  02-11  (C. Ethe)  F90: Free form and module
   !!----------------------------------------------------------------------
   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
   !! $Id$
   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
   !!----------------------------------------------------------------------
   !! * Modules used
   USE par_ice_2

   IMPLICIT NONE
   PRIVATE

   PUBLIC thd_ice_alloc_2 ! Routine called by nemogcm.F90

   !! * Share Module variables
   REAL(wp) , PUBLIC ::   & !!! ** ice-thermo namelist (namicethd) **
      hmelt   = -0.15  ,  &  !: maximum melting at the bottom
      hicmin  = 0.2    ,  &  !: ice th. corr. to max. ener. in brine pocket
      hiclim  = 0.05   ,  &  !: minimum ice thickness
      amax    = 0.999  ,  &  !: maximum lead fraction
      swiqst  = 1.0    ,  &  !: energy stored in brine pocket (1) or not (0)
      sbeta   = 1.0    ,  &  !: numerical scheme for diffusion in ice 
      parlat  = 0.0    ,  &  !: percent. of energy used for lateral ablation
      hakspl  = 0.5    ,  &  !: slope of distr. for Hakkinen-Mellro's lat. melt
      hibspl  = 0.5    ,  &  !: slope of distribution for Hibler's lat. melt
      exld    = 2.0    ,  &  !: exponent for leads-closure rate
      hakdif  = 1.0    ,  &  !: coefficient for diffusions of ice and snow
      thth    = 0.2    ,  &  !: thick. for comp. of eq. thermal conduct
      hnzst   = 0.1    ,  &  !: thick. of the surf. layer in temp. comp.
      parsub  = 1.0    ,  &  !: switch for snow sublimation or not
      alphs   = 1.0          !: coef. for snow density when snow-ice formation

   REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:   
      hiccrit = (/0.3,0.3/)  !: ice th. for lateral accretion in the NH (SH) (m)

   REAL(wp) , PUBLIC ::   &  !:
      uscomi,             &  !: inverse of minimum lead fraction
      cnscg                  !: ratio  rcpsn/rcpic

   INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   &  !:
      npb     ,   &   !: number of points where computations has to be done
      npac            !: correspondance between the points

   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   &  !: 
      qldif_1d    ,     &  !: corresponding to the 2D var  qldif
      qcmif_1d    ,     &  !: corresponding to the 2D var  qcmif
      thcm_1d     ,     &  !:    "                  "      thcm
      fstbif_1d   ,     &  !:    "                  "      fstric
      fltbif_1d   ,     &  !:    "                  "      ffltbif
      fscbq_1d    ,     &  !:    "                  "      fscmcbq
      qsr_ice_1d  ,     &  !:    "                  "      qsr_ice
      fr1_i0_1d   ,     &  !:    "                  "      fr1_i0
      fr2_i0_1d   ,     &  !:    "                  "      fr2_i0
      qns_ice_1d  ,     &  !:    "                  "      qns_ice
      qfvbq_1d    ,     &  !:    "                  "      qfvbq
      sist_1d     ,     &  !:    "                  "      sist
      tfu_1d      ,     &  !:    "                  "      tfu
      sprecip_1d  ,     &  !:    "                  "      sprecip
      h_snow_1d   ,     &  !:    "                  "      h_snow
      h_ice_1d    ,     &  !:    "                  "      h_ice
      frld_1d     ,     &  !:    "                  "      frld
      qstbif_1d   ,     &  !:    "                  "      qstoif
      fbif_1d     ,     &  !:    "                  "      fbif
      rdmicif_1d  ,     &  !:    "                  "      rdmicif
      rdmsnif_1d  ,     &  !:    "                  "      rdmsnif
      qlbbq_1d    ,     &  !:    "                  "      qlbsbq
      dmgwi_1d    ,     &  !:    "                  "      dmgwi
      dvsbq_1d    ,     &  !:    "                  "      rdvosif
      rdvomif_1d  ,     &  !:    "                  "      rdvomif
      dvbbq_1d    ,     &  !:    "                  "      rdvobif
      dvlbq_1d    ,     &  !:    "                  "      rdvolif
      dvnbq_1d    ,     &  !:    "                  "      rdvolif
      dqns_ice_1d ,     &  !:    "                  "      dqns_ice
      qla_ice_1d  ,     &  !:    "                  "      qla_ice
      dqla_ice_1d          !:    "                  "      dqla_ice

   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   &  !:
      tbif_1d              !: corresponding to the 2D var  tbif

   !!======================================================================
 CONTAINS

   FUNCTION thd_ice_alloc_2()
      USE in_out_manager, ONLY: ctl_warn
      IMPLICIT none
      INTEGER :: thd_ice_alloc_2
      ! Local vars
      INTEGER :: ierr(4)

      ierr(:) = 0

      ALLOCATE(npb(jpij), npac(jpij),                             &
               qldif_1d(jpij), qcmif_1d(jpij), thcm_1d(jpij),     &
               fstbif_1d(jpij), fltbif_1d(jpij), fscbq_1d(jpij),  &
               qsr_ice_1d(jpij),fr1_i0_1d(jpij), fr2_i0_1d(jpij), Stat=ierr(1))

      ALLOCATE(qns_ice_1d(jpij), qfvbq_1d(jpij), sist_1d(jpij), tfu_1d(jpij), &
               sprecip_1d(jpij), h_snow_1d(jpij),h_ice_1d(jpij),frld_1d(jpij),&
               qstbif_1d(jpij),  fbif_1d(jpij),  Stat=ierr(2))

      ALLOCATE(rdmicif_1d(jpij), rdmsnif_1d(jpij), qlbbq_1d(jpij),   &
               dmgwi_1d(jpij)  , dvsbq_1d(jpij)  , rdvomif_1d(jpij), &
               dvbbq_1d(jpij)  , dvlbq_1d(jpij)  , dvnbq_1d(jpij)  , &
               Stat=ierr(3))

      ALLOCATE(dqns_ice_1d(jpij) ,qla_ice_1d(jpij), dqla_ice_1d(jpij), &
               tbif_1d(jpij, jplayersp1), Stat=ierr(4))

      thd_ice_alloc_2 = MAXVAL(ierr)

      IF(thd_ice_alloc_2 /= 0)THEN
         CALL ctl_warn('thd_ice_alloc_2: failed to allocate arrays.')
      END IF

   END FUNCTION thd_ice_alloc_2

#endif
END MODULE thd_ice_2