New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
thd_ice_2.F90 in branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2 – NEMO

source: branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/thd_ice_2.F90 @ 1858

Last change on this file since 1858 was 1858, checked in by gm, 14 years ago

ticket:#665 : step 1 - heat content of freezing-melting ice

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 4.8 KB
Line 
1MODULE thd_ice_2
2#if defined key_lim2
3   !!======================================================================
4   !!                       ***  MODULE thd_ice_2  ***
5   !! LIM 2.0 sea-ice :   Ice thermodynamics in 1D
6   !!=====================================================================
7   !! History : 2.0  ! 2002-11  (C. Ethe)  F90: Free form and module
8   !!           2.1  ! 2010-05  (Y. Aksenov, M. Vancoppenolle, G. Madec) add heat content exchanges
9   !!----------------------------------------------------------------------
10   !!   LIM 2.0, UCL-LOCEAN-IPSL (2005)
11   !! $Id$
12   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
13   !!----------------------------------------------------------------------
14   !! * Modules used
15   USE par_ice_2
16
17   IMPLICIT NONE
18   PRIVATE
19
20   !! * Share Module variables
21   REAL(wp) , PUBLIC ::   & !!! ** ice-thermo namelist (namicethd) **
22      hmelt   = -0.15  ,  &  !: maximum melting at the bottom
23      hicmin  = 0.2    ,  &  !: ice th. corr. to max. ener. in brine pocket
24      hiclim  = 0.05   ,  &  !: minimum ice thickness
25      amax    = 0.999  ,  &  !: maximum lead fraction
26      swiqst  = 1.0    ,  &  !: energy stored in brine pocket (1) or not (0)
27      sbeta   = 1.0    ,  &  !: numerical scheme for diffusion in ice
28      parlat  = 0.0    ,  &  !: percent. of energy used for lateral ablation
29      hakspl  = 0.5    ,  &  !: slope of distr. for Hakkinen-Mellro's lat. melt
30      hibspl  = 0.5    ,  &  !: slope of distribution for Hibler's lat. melt
31      exld    = 2.0    ,  &  !: exponent for leads-closure rate
32      hakdif  = 1.0    ,  &  !: coefficient for diffusions of ice and snow
33      thth    = 0.2    ,  &  !: thick. for comp. of eq. thermal conduct
34      hnzst   = 0.1    ,  &  !: thick. of the surf. layer in temp. comp.
35      parsub  = 1.0    ,  &  !: switch for snow sublimation or not
36      alphs   = 1.0          !: coef. for snow density when snow-ice formation
37
38   REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:   
39      hiccrit = (/0.3,0.3/)  !: ice th. for lateral accretion in the NH (SH) (m)
40
41   REAL(wp) , PUBLIC ::   &  !:
42      uscomi,             &  !: inverse of minimum lead fraction
43      cnscg                  !: ratio  rcpsn/rcpic
44
45   INTEGER , PUBLIC, DIMENSION(jpij) ::   &  !:
46      npb     ,   &   !: number of points where computations has to be done
47      npac            !: correspondance between the points
48
49   REAL(wp), PUBLIC, DIMENSION(jpij) ::   &  !:
50      qldif_1d    ,     &  !: corresponding to the 2D var  qldif
51      qcmif_1d    ,     &  !: corresponding to the 2D var  qcmif
52      thcm_1d     ,     &  !:    "                  "      thcm
53      fstbif_1d   ,     &  !:    "                  "      fstric
54      fltbif_1d   ,     &  !:    "                  "      ffltbif
55      fscbq_1d    ,     &  !:    "                  "      fscmcbq
56      qsr_ice_1d  ,     &  !:    "                  "      qsr_ice
57      fr1_i0_1d   ,     &  !:    "                  "      fr1_i0
58      fr2_i0_1d   ,     &  !:    "                  "      fr2_i0
59      qns_ice_1d  ,     &  !:    "                  "      qns_ice
60      qfvbq_1d    ,     &  !:    "                  "      qfvbq
61      sist_1d     ,     &  !:    "                  "      sist
62      tfu_1d      ,     &  !:    "                  "      tfu
63      sprecip_1d  ,     &  !:    "                  "      sprecip
64      h_snow_1d   ,     &  !:    "                  "      h_snow
65      h_ice_1d    ,     &  !:    "                  "      h_ice
66      frld_1d     ,     &  !:    "                  "      frld
67      qstbif_1d   ,     &  !:    "                  "      qstoif
68      fbif_1d     ,     &  !:    "                  "      fbif
69      rdm_ice_1d  ,     &  !:    "                  "      rdm_ice
70      rdq_ice_1d  ,     &  !:    "                  "      rdq_ice
71      rdm_snw_1d  ,     &  !:    "                  "      rdm_snw
72      rdq_snw_1d  ,     &  !:    "                  "      rdq_snw
73      qlbbq_1d    ,     &  !:    "                  "      qlbsbq
74      dmgwi_1d    ,     &  !:    "                  "      dmgwi
75      dvsbq_1d    ,     &  !:    "                  "      rdvosif
76      rdvomif_1d  ,     &  !:    "                  "      rdvomif
77      dvbbq_1d    ,     &  !:    "                  "      rdvobif
78      dvlbq_1d    ,     &  !:    "                  "      rdvolif
79      dvnbq_1d    ,     &  !:    "                  "      rdvolif
80      dqns_ice_1d ,     &  !:    "                  "      dqns_ice
81      qla_ice_1d  ,     &  !:    "                  "      qla_ice
82      dqla_ice_1d          !:    "                  "      dqla_ice
83
84   REAL(wp), PUBLIC, DIMENSION(jpij,jplayersp1) ::   &  !:
85      tbif_1d              !: corresponding to the 2D var  tbif
86
87   !!======================================================================
88#endif
89END MODULE thd_ice_2
Note: See TracBrowser for help on using the repository browser.