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/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/LIM_SRC_2 – NEMO

source: branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/LIM_SRC_2/thd_ice_2.F90 @ 2833

Last change on this file since 2833 was 2715, checked in by rblod, 13 years ago

First attempt to put dynamic allocation on the trunk

  • Property svn:keywords set to Id
File size: 6.3 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 :
8   !!   2.0  !  02-11  (C. Ethe)  F90: Free form and module
9   !!----------------------------------------------------------------------
10   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
11   !! $Id$
12   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
13   !!----------------------------------------------------------------------
14   USE par_ice_2
15
16   IMPLICIT NONE
17   PRIVATE
18
19   PUBLIC thd_ice_alloc_2 ! Routine called by nemogcm.F90
20
21   !! * Share Module variables
22   REAL(wp) , PUBLIC ::   & !!! ** ice-thermo namelist (namicethd) **
23      hmelt   = -0.15  ,  &  !: maximum melting at the bottom
24      hicmin  = 0.2    ,  &  !: ice th. corr. to max. ener. in brine pocket
25      hiclim  = 0.05   ,  &  !: minimum ice thickness
26      amax    = 0.999  ,  &  !: maximum lead fraction
27      swiqst  = 1.0    ,  &  !: energy stored in brine pocket (1) or not (0)
28      sbeta   = 1.0    ,  &  !: numerical scheme for diffusion in ice
29      parlat  = 0.0    ,  &  !: percent. of energy used for lateral ablation
30      hakspl  = 0.5    ,  &  !: slope of distr. for Hakkinen-Mellro's lat. melt
31      hibspl  = 0.5    ,  &  !: slope of distribution for Hibler's lat. melt
32      exld    = 2.0    ,  &  !: exponent for leads-closure rate
33      hakdif  = 1.0    ,  &  !: coefficient for diffusions of ice and snow
34      thth    = 0.2    ,  &  !: thick. for comp. of eq. thermal conduct
35      hnzst   = 0.1    ,  &  !: thick. of the surf. layer in temp. comp.
36      parsub  = 1.0    ,  &  !: switch for snow sublimation or not
37      alphs   = 1.0          !: coef. for snow density when snow-ice formation
38
39   REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:   
40      hiccrit = (/0.3,0.3/)  !: ice th. for lateral accretion in the NH (SH) (m)
41
42   REAL(wp) , PUBLIC ::   &  !:
43      uscomi,             &  !: inverse of minimum lead fraction
44      cnscg                  !: ratio  rcpsn/rcpic
45
46   INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   &  !:
47      npb     ,   &   !: number of points where computations has to be done
48      npac            !: correspondance between the points
49
50   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   &  !:
51      qldif_1d    ,     &  !: corresponding to the 2D var  qldif
52      qcmif_1d    ,     &  !: corresponding to the 2D var  qcmif
53      thcm_1d     ,     &  !:    "                  "      thcm
54      fstbif_1d   ,     &  !:    "                  "      fstric
55      fltbif_1d   ,     &  !:    "                  "      ffltbif
56      fscbq_1d    ,     &  !:    "                  "      fscmcbq
57      qsr_ice_1d  ,     &  !:    "                  "      qsr_ice
58      fr1_i0_1d   ,     &  !:    "                  "      fr1_i0
59      fr2_i0_1d   ,     &  !:    "                  "      fr2_i0
60      qns_ice_1d  ,     &  !:    "                  "      qns_ice
61      qfvbq_1d    ,     &  !:    "                  "      qfvbq
62      sist_1d     ,     &  !:    "                  "      sist
63      tfu_1d      ,     &  !:    "                  "      tfu
64      sprecip_1d  ,     &  !:    "                  "      sprecip
65      h_snow_1d   ,     &  !:    "                  "      h_snow
66      h_ice_1d    ,     &  !:    "                  "      h_ice
67      frld_1d     ,     &  !:    "                  "      frld
68      qstbif_1d   ,     &  !:    "                  "      qstoif
69      fbif_1d     ,     &  !:    "                  "      fbif
70      rdmicif_1d  ,     &  !:    "                  "      rdmicif
71      rdmsnif_1d  ,     &  !:    "                  "      rdmsnif
72      qlbbq_1d    ,     &  !:    "                  "      qlbsbq
73      dmgwi_1d    ,     &  !:    "                  "      dmgwi
74      dvsbq_1d    ,     &  !:    "                  "      rdvosif
75      rdvomif_1d  ,     &  !:    "                  "      rdvomif
76      dvbbq_1d    ,     &  !:    "                  "      rdvobif
77      dvlbq_1d    ,     &  !:    "                  "      rdvolif
78      dvnbq_1d    ,     &  !:    "                  "      rdvolif
79      dqns_ice_1d ,     &  !:    "                  "      dqns_ice
80      qla_ice_1d  ,     &  !:    "                  "      qla_ice
81      dqla_ice_1d          !:    "                  "      dqla_ice
82
83   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   &  !:
84      tbif_1d              !: corresponding to the 2D var  tbif
85
86   !!----------------------------------------------------------------------
87   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
88   !! $Id$
89   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
90   !!----------------------------------------------------------------------
91 CONTAINS
92
93   INTEGER FUNCTION thd_ice_alloc_2()
94      !!----------------------------------------------------------------------
95      USE lib_mpp        ! MPP library
96      INTEGER :: ierr(4)
97      !!----------------------------------------------------------------------
98      !
99      ierr(:) = 0
100      !
101      ALLOCATE( npb(jpij), npac(jpij),                             &
102         &      qldif_1d(jpij), qcmif_1d(jpij), thcm_1d(jpij),     &
103         &      fstbif_1d(jpij), fltbif_1d(jpij), fscbq_1d(jpij),  &
104         &      qsr_ice_1d(jpij),fr1_i0_1d(jpij), fr2_i0_1d(jpij), Stat=ierr(1))
105         !
106      ALLOCATE( qns_ice_1d(jpij), qfvbq_1d(jpij), sist_1d(jpij), tfu_1d(jpij), &
107         &      sprecip_1d(jpij), h_snow_1d(jpij),h_ice_1d(jpij),frld_1d(jpij),&
108         &      qstbif_1d(jpij),  fbif_1d(jpij),  Stat=ierr(2))
109         !
110      ALLOCATE( rdmicif_1d(jpij), rdmsnif_1d(jpij), qlbbq_1d(jpij),   &
111         &      dmgwi_1d(jpij)  , dvsbq_1d(jpij)  , rdvomif_1d(jpij), &
112         &      dvbbq_1d(jpij)  , dvlbq_1d(jpij)  , dvnbq_1d(jpij)  , &
113         &      Stat=ierr(3))
114         !
115      ALLOCATE( dqns_ice_1d(jpij) ,qla_ice_1d(jpij), dqla_ice_1d(jpij), &
116         &      tbif_1d(jpij, jplayersp1), Stat=ierr(4))
117         !
118      thd_ice_alloc_2 = MAXVAL(ierr)
119      IF( thd_ice_alloc_2 /= 0 )   CALL ctl_warn('thd_ice_alloc_2: failed to allocate arrays')
120      !
121   END FUNCTION thd_ice_alloc_2
122
123#endif
124   !!======================================================================
125END MODULE thd_ice_2
Note: See TracBrowser for help on using the repository browser.