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.
ice.F90 in trunk/NEMO/LIM_SRC_3 – NEMO

source: trunk/NEMO/LIM_SRC_3/ice.F90 @ 921

Last change on this file since 921 was 921, checked in by rblod, 16 years ago

Correct indentation and print for debug in LIM3, see ticket #134, step I

File size: 29.6 KB
Line 
1MODULE ice
2#if defined key_lim3
3   !!----------------------------------------------------------------------
4   !!   'key_lim3' :                                   LIM3 sea-ice model
5   !!----------------------------------------------------------------------
6   !! History :
7   !!   2.0  !  03-08  (C. Ethe)  F90: Free form and module
8   !!   3.0  !  08-03  (M. Vancoppenolle) : LIM3 !
9   !!----------------------------------------------------------------------
10   !!  LIM 3.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          ! LIM sea-ice parameters
16
17   IMPLICIT NONE
18   PRIVATE
19   !!
20   !!======================================================================
21   !!                        ***  MODULE ice  ***
22   !!
23   !!                           **************
24   !!                           * L I M  3.0 *
25   !!                           **************
26   !!
27   !!                         ''in ice we trust''
28   !!
29   !!                   This module contains the sea ice
30   !!                 diagnostics variables of ice defined
31   !!                             in memory
32   !!
33   !!======================================================================
34   !!
35   !! LIM3 by the use of sweat, agile fingers and sometimes brain juice,
36   !!  was developed in Louvain-la-Neuve by :
37   !!
38   !!    * Martin Vancoppenolle (UCL-ASTR, Belgium)
39   !!    * Sylvain Bouillon (UCL-ASTR, Belgium)
40   !!    * Miguel Angel Morales Maqueda (POL, UK)
41   !!
42   !! Based on extremely valuable earlier work by
43   !!
44   !!    * Thierry Fichefet
45   !!    * Hugues Goosse
46   !!
47   !! The following persons also contributed to the code in various ways
48   !!
49   !!    * Gurvan Madec, Claude Talandier, Christian Ethe
50   !!      and Rachid Benshila (LOCEAN-IPSL, France)
51   !!    * Xavier Fettweis (UCL-ASTR), Ralph Timmermann (AWI, Germany)
52   !!    * Bill Lipscomb (LANL), Cecilia Bitz (UWa)
53   !!      and Elisabeth Hunke (LANL), USA.
54   !!
55   !! (c) UCL-ASTR, 2005-2008
56   !!
57   !! For more info, the interested user is kindly invited to consult the
58   !! following references
59   !!    For model description and validation :
60   !!    * Vancoppenolle et al., Ocean Modelling, 2008a.
61   !!    * Vancoppenolle et al., Ocean Modelling, 2008b.
62   !!
63   !!    For a specific description of EVP :
64   !!    * Bouillon et al., in prep for 2008.
65   !!
66   !!    Or the reference manual, that should be available by 2009
67   !!
68   !!======================================================================
69   !!                                                                     |
70   !!            *****************************************                |
71   !!            *                                       *                |
72   !! ************ I C E   S T A T E   V A R I A B L E S **************** |
73   !!            *                                       *                |
74   !!            *****************************************                |
75   !!                                                                     |
76   !! Introduction :                                                      |
77   !! --------------                                                      |
78   !!                                                                     |
79   !! Every ice-covered grid cell is characterized by a series of state   |
80   !! variables. To account for unresolved spatial variability in ice     |
81   !! thickness, the ice cover in divided in ice thickness categories.    |
82   !!                                                                     |
83   !! Sea ice state variables depend on the ice thickness category        |
84   !!                                                                     |
85   !! Those variables are divided into two groups                         |
86   !! * Extensive (or global) variables.                                  |
87   !!   These are the variables that are transported by all means         |
88   !! * Intensive (or equivalent) variables.                              |
89   !!   These are the variables that are either physically more           |
90   !!   meaningful and/or used in ice thermodynamics                      |
91   !!                                                                     |
92   !! Routines in limvar.F90 perform conversions                          |
93   !!  - lim_var_glo2eqv  : from global to equivalent variables           |
94   !!  - lim_var_eqv2glo  : from equivalent to global variables           |
95   !!                                                                     |
96   !! For various purposes, the sea ice state variables have sometimes    |
97   !! to be aggregated over all ice thickness categories. This operation  |
98   !! is done in :                                                        |
99   !!  - lim_var_agg                                                      |
100   !!                                                                     |
101   !! in icestp.F90, the routines that compute the changes in the ice     |
102   !! state variables are called                                          |
103   !! - lim_dyn : ice dynamics                                            |
104   !! - lim_trp : ice transport                                           |
105   !! - lim_itd_me : mechanical redistribution (ridging and rafting)      |
106   !! - lim_thd : ice halo-thermodynamics                                 |
107   !! - lim_itd_th : thermodynamic changes in ice thickness distribution  |
108   !!                and creation of new ice                              |
109   !!                                                                     |
110   !! See the associated routines for more information                    |
111   !!                                                                     |
112   !! List of ice state variables :                                       |
113   !! -----------------------------                                       |
114   !!                                                                     |
115   !!-------------|-------------|---------------------------------|-------|
116   !!   name in   |   name in   |              meaning            | units |
117   !! 2D routines | 1D routines |                                 |       |
118   !!-------------|-------------|---------------------------------|-------|
119   !!                                                                     |
120   !! ******************************************************************* |
121   !! ***         Dynamical variables (prognostic)                    *** |
122   !! ******************************************************************* |
123   !!                                                                     |
124   !! u_ice       |      -      |    Comp. U of the ice velocity  | m/s   |
125   !! v_ice       |      -      |    Comp. V of the ice velocity  | m/s   |
126   !!                                                                     |
127   !! ******************************************************************* |
128   !! ***         Category dependent state variables (prognostic)     *** |
129   !! ******************************************************************* |
130   !!                                                                     |
131   !! ** Global variables                                                 |
132   !!                                                                     |
133   !!-------------|-------------|---------------------------------|-------|
134   !! a_i         | a_i_b       |    Ice concentration            |       |
135   !! v_i         |      -      |    Ice volume per unit area     | m     |
136   !! v_s         |      -      |    Snow volume per unit area    | m     |
137   !! smv_i       |      -      |    Sea ice salt content         | ppt.m |
138   !! oa_i        !      -      !    Sea ice areal age content    | day   |
139   !! e_i         !      -      !    Ice enthalpy                 | 10^9 J|
140   !!      -      ! q_i_b       !    Ice enthalpy per unit vol.   | J/m3  |
141   !! e_s         !      -      !    Snow enthalpy                | 10^9 J|
142   !!      -      ! q_s_b       !    Snow enthalpy per unit vol.  | J/m3  |
143   !!                                                                     |
144   !!-------------|-------------|---------------------------------|-------|
145   !!                                                                     |
146   !! ** Equivalent variables                                             |
147   !!                                                                     |
148   !!-------------|-------------|---------------------------------|-------|
149   !!                                                                     |
150   !! ht_i        | ht_i_b      |    Ice thickness                | m     |
151   !! ht_s        ! ht_s_b      |    Snow depth                   | m     |
152   !! sm_i        ! sm_i_b      |    Sea ice bulk salinity        ! ppt   |
153   !! s_i         ! s_i_b       |    Sea ice salinity profile     ! ppt   |
154   !! o_i         !      -      |    Sea ice Age                  ! days  |
155   !! t_i         ! t_i_b       |    Sea ice temperature          ! K     |
156   !! t_s         ! t_s_b       |    Snow temperature             ! K     |
157   !! t_su        ! t_su_b      |    Sea ice surface temperature  ! K     |
158   !!                                                                     |
159   !! notes: the ice model only sees a bulk (i.e., vertically averaged)   |
160   !!        salinity, except in thermodynamic computations, for which    |
161   !!        the salinity profile is computed as a function of bulk       |
162   !!        salinity                                                     |
163   !!                                                                     |
164   !!        the sea ice surface temperature is not associated to any     |
165   !!        heat content. Therefore, it is not a state variable and      |
166   !!        does not have to be advected. Nevertheless, it has to be     |
167   !!        computed to determine whether the ice is melting or not      |
168   !!                                                                     |
169   !! ******************************************************************* |
170   !! ***         Category-summed state variables (diagnostic)        *** |
171   !! ******************************************************************* |
172   !! at_i        | at_i_b      |    Total ice concentration      |       |
173   !! vt_i        |      -      |    Total ice vol. per unit area | m     |
174   !! vt_s        |      -      |    Total snow vol. per unit ar. | m     |
175   !! smt_i       |      -      |    Mean sea ice salinity        | ppt   |
176   !! tm_i        |      -      |    Mean sea ice temperature     | K     |
177   !! ot_i        !      -      !    Sea ice areal age content    | day   |
178   !! et_i        !      -      !    Total ice enthalpy           | 10^9 J|
179   !! et_s        !      -      !    Total snow enthalpy          | 10^9 J|
180   !! bv_i        !      -      !    Mean relative brine volume   | ???   |
181   !!                                                                     |
182   !!                                                                     |
183   !!=====================================================================
184
185   LOGICAL, PUBLIC ::    &
186      con_i = .false.           ! switch for conservation test
187
188   !!--------------------------------------------------------------------------
189   !! * Share Module variables
190   !!--------------------------------------------------------------------------
191   INTEGER , PUBLIC ::   &     !!: ** ice-dynamic namelist (namicedyn) **
192      nbiter = 1      ,  &     !: number of sub-time steps for relaxation
193      nbitdr = 250    ,  &     !: maximum number of iterations for relaxation
194      nevp   = 400    ,  &     !: number of iterations for subcycling
195      nlay_i = 5               !: number of layers in the ice
196
197   REAL(wp), PUBLIC ::   &     !!: ** ice-dynamic namelist (namicedyn) **
198      epsd   = 1.0e-20,  &     !: tolerance parameter for dynamic
199      alpha  = 0.5    ,  &     !: coefficient for semi-implicit coriolis
200      dm     = 0.6e+03,  &     !: diffusion constant for dynamics
201      om     = 0.5    ,  &     !: relaxation constant
202      resl   = 5.0e-05,  &     !: maximum value for the residual of relaxation
203      cw     = 5.0e-03,  &     !: drag coefficient for oceanic stress
204      angvg  = 0.0    ,  &     !: turning angle for oceanic stress
205      pstar  = 1.0e+04,  &     !: determines ice strength (N/M), Hibler JPO79
206      c_rhg  = 20.0   ,  &     !: determines changes in ice strength
207      etamn  = 0.0e+07,  &     !: minimun value for viscosity : has to be 0
208      creepl = 2.0e-08,  &     !: creep limit : has to be under 1.0e-9
209      ecc    = 2.0    ,  &  !: eccentricity of the elliptical yield curve
210      ahi0   = 350.e0 ,  &  !: sea-ice hor. eddy diffusivity coeff. (m2/s)
211      telast = 2880.0 ,  &  !: timescale for elastic waves (s) !SB
212      alphaevp = 1.0  ,  &  !: coeficient of the internal stresses !SB
213      unit_fac = 1.0e9      !: conversion factor for ice / snow enthalpy
214
215   REAL(wp), PUBLIC ::   & !!: ** ice-salinity namelist (namicesal) **
216      s_i_max  = 20.0 ,  &  !: maximum ice salinity (ppt)
217      s_i_min  =  0.1 ,  &  !: minimum ice salinity (ppt)
218      s_i_0    =  3.5 ,  &  !: 1st sal. value for the computation of sal .prof.
219                                !: (ppt)
220      s_i_1    =  4.5 ,  &  !: 2nd sal. value for the computation of sal .prof.
221                                !: (ppt)
222      sal_G    = 5.00 ,  &  !: restoring salinity for gravity drainage
223                                !: (ppt)
224      sal_F    = 2.50 ,  &  !: restoring salinity for flushing
225                                !: (ppt)
226      time_G   = 1.728e+06,&!: restoring time constant for gravity drainage
227                                !: (= 20 days, in s)
228      time_F   = 8.640e+05,&!: restoring time constant for gravity drainage
229                                !: (= 10 days, in s)
230      bulk_sal = 4.0        !: bulk salinity (ppt) in case of constant salinity
231
232   INTEGER , PUBLIC ::   & !!: ** ice-salinity namelist (namicesal) **
233      num_sal  = 1    ,  &  !: salinity configuration used in the model
234                                !: 1 - s constant in space and time
235                                !: 2 - prognostic salinity (s(z,t))
236                                !: 3 - salinity profile, constant in time
237                                !: 4 - salinity variations affect only ice
238                                !      thermodynamics
239      sal_prof = 1    ,  &  !: salinity profile or not
240      thcon_i_swi = 1       !: thermal conductivity of Untersteiner (1964) (1) or
241   !: Pringle et al (2007) (2)
242
243   REAL(wp), PUBLIC ::   & !!: ** ice-mechanical redistribution namelist (namiceitdme)
244      Cs = 0.25       ,  & !!: fraction of shearing energy contributing to ridging           
245      Cf = 17.0       ,  & !!: ratio of ridging work to PE loss
246      fsnowrdg = 0.5  ,  & !!: fractional snow loss to the ocean during ridging
247      fsnowrft = 0.5  ,  & !!: fractional snow loss to the ocean during ridging
248      Gstar = 0.15    ,  & !!: fractional area of young ice contributing to ridging
249      astar = 0.05    ,  & !!: equivalent of G* for an exponential participation function
250      Hstar = 100.0   ,  & !!: thickness that determines the maximal thickness of ridged
251                                !!: ice
252      hparmeter = 0.75,  & !!: threshold thickness (m) for rafting / ridging
253      Craft = 5.0     ,  & !!: coefficient for smoothness of the hyperbolic tangent in rafting
254      ridge_por = 0.0 ,  & !!: initial porosity of ridges (0.3 regular value)
255      sal_max_ridge = 15.0, & !!: maximum ridged ice salinity (ppt)
256      betas    = 1.0      , & !:: coef. for partitioning of snowfall between leads and sea ice
257      kappa_i  = 1.0      , & !!: coefficient for the extinction of radiation
258                                !!: Grenfell et al. (2006) (m-1)
259      nconv_i_thd = 50    , & !!: maximal number of iterations for heat diffusion
260      maxer_i_thd = 1.0e-4    !!: maximal tolerated error (C) for heat diffusion
261
262   INTEGER , PUBLIC ::   & !!: ** ice-mechanical redistribution namelist (namiceitdme)
263      ridge_scheme_swi = 0, & !!: scheme used for ice ridging
264      raftswi          = 1, & !!: rafting of ice or not                       
265      partfun_swi      = 1, & !!: participation function Thorndike et al. JGR75 (0)
266                                !!: or Lipscomb et al. JGR07 (1)
267      transfun_swi     = 0, & !!: transfer function of Hibler, MWR80 (0)
268                                !!: or Lipscomb et al., 2007 (1)
269      brinstren_swi    = 0    !!: use brine volume to diminish ice strength
270
271   REAL(wp), PUBLIC ::   &  !:
272      usecc2          ,  &  !:  = 1.0 / ( ecc * ecc )
273      rhoco           ,  &  !: = rau0 * cw
274      sangvg, cangvg  ,  &  !: sin and cos of the turning angle for ocean stress
275      pstarh                !: pstar / 2.0
276
277   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::  &  !:
278      u_oce, v_oce    ,  &  !: surface ocean velocity used in ice dynamics
279      ahiu , ahiv     ,  &  !: hor. diffusivity coeff. at ocean U- and V-points (m2/s)
280      pahu , pahv     ,  &  !: ice hor. eddy diffusivity coef. at ocean U- and V-points
281      ust2s, hicol    ,  &  !: friction velocity, ice collection thickness accreted in leads
282      strength        ,  &  !: ice strength
283      strp1, strp2    ,  &  !: strength at previous time steps
284      stress1_i       ,  &  !: first stress tensor element
285      stress2_i       ,  &  !: second stress tensor element
286      stress12_i      ,  &  !: diagonal stress tensor element
287      delta_i         ,  &  !: Delta factor for the ice rheology (see Flato and Hibler 95) [s-1] -> limrhg.F90
288      divu_i          ,  &  !: Divergence of the velocity field [s-1] -> limrhg.F90
289      shear_i               !: Shear of the velocity field [s-1] -> limrhg.F90
290
291   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
292      firic  ,   &  !: IR flux over the ice (only used for outputs)
293      fcsic  ,   &  !: Sensible heat flux over the ice (only used for outputs)
294      fleic  ,   &  !: Latent heat flux over the ice (only used for outputs)
295      qlatic ,   &  !: latent flux
296      rdvosif,   &  !: Variation of volume at surface (only used for outputs)
297      rdvobif,   &  !: Variation of ice volume at the bottom ice (only used for outputs)
298      fdvolif,   &  !: Total variation of ice volume (only used for outputs)
299      rdvonif,   &  !: Lateral Variation of ice volume (only used for outputs)
300      sist   ,   &  !: Average Sea-Ice Surface Temperature (Kelvin ??? degree ??? I don't know)
301      t_bo   ,   &  !: Sea-Ice bottom temperature (Kelvin)     
302      hicifp ,   &  !: Ice production/melting
303                                !obsolete... can be removed
304      frld   ,   &  !: Leads fraction = 1-a/totalarea REFERS TO LEAD FRACTION everywhere
305                                !: except in the OUTPUTS!!!!
306      pfrld  ,   &  !: Leads fraction at previous time 
307      phicif ,   &  !: Old ice thickness
308      fbif   ,   &  !: Heat flux at the ice base
309      rdmsnif,   &  !: Variation of snow mass
310      rdmicif,   &  !: Variation of ice mass
311      qldif  ,   &  !: heat balance of the lead (or of the open ocean)
312      qcmif  ,   &  !: Energy needed to bring the ocean surface layer until its freezing
313      fdtcn  ,   &  !: net downward heat flux from the ice to the ocean
314      qdtcn  ,   &  !: energy from the ice to the ocean
315      fstric ,   &  !: transmitted solar radiation under ice
316      fscmbq ,   &  !: associated with lead chipotage with solar flux
317      ffltbif,   &  !: Array linked with the max heat contained in brine pockets (?)
318      fsbbq  ,   &  !: Also linked with the solar flux below the ice (?)
319      qfvbq  ,   &  !: Array used to store energy in case of toral lateral ablation (?)
320      dmgwi  ,   &  !: Variation of the mass of snow ice
321      fsalt_res, &  !: Residual salt flux due to correction of ice thickness
322      fsbri  ,   &  !: Salt flux due to brine rejection
323      fsalt_rpo, &  !: Salt flux associated with porous ridged ice formation
324      fheat_rpo, &  !: Heat flux associated with porous ridged ice formation
325      fhbri  ,   &  !: heat flux due to brine rejection
326      fmmec  ,   &  !: Mass flux due to snow loss during compression
327      fseqv  ,   &  !: Equivalent salt flux due to ice growth/melt
328      fheat_res, &  !: Residual heat flux due to correction of ice thickness
329      fhmec         !: Heat flux due to snow loss during compression
330
331   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::       &  !:
332      albege ,   &  !: Albedo of the snow or ice (only for outputs)
333      albecn ,   &  !: Albedo of the ocean (only for outputs)
334      tauc          !: Cloud optical depth
335
336   ! temporary arrays for dummy version of the code
337   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
338      dh_i_surf2D, dh_i_bott2D, fstbif, fsup2D, focea2D, q_s
339
340   INTEGER, PUBLIC, DIMENSION(jpi, jpj, jpl) ::          &   !::
341      patho_case ! number of the pathological case (if any, of course)
342
343   !!--------------------------------------------------------------------------
344   !! * Ice global state variables
345   !!--------------------------------------------------------------------------
346   !! Variables defined for each ice category
347   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl)        ::   &  !:
348      ht_i   ,   &  !: Ice thickness (m)
349      a_i    ,   &  !: Ice fractional areas (concentration)
350      v_i    ,   &  !: Ice volume per unit area (m)
351      v_s    ,   &  !: Snow volume per unit area(m)
352      ht_s   ,   &  !: Snow thickness (m)
353      t_su   ,   &  !: Sea-Ice Surface Temperature (K)
354      sm_i   ,   &  !: Sea-Ice Bulk salinity (ppt)
355      smv_i  ,   &  !: Sea-Ice Bulk salinity times volume per area (ppt.m)
356                                !: this is an extensive variable that has to be transported
357      o_i    ,   &  !: Sea-Ice Age (days)
358      ov_i   ,   &  !: Sea-Ice Age times volume per area (days.m)
359      oa_i          !: Sea-Ice Age times ice area (days)
360
361   !! Variables summed over all categories, or associated to
362   !! all the ice in a single grid cell
363   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
364      u_ice, v_ice,   &  !: two components of the ice velocity (m/s)
365      tio_u, tio_v,   &  !: two components of the ice-ocean stress (N/m2)
366      vt_i        ,   &  !: ice total volume per unit area (m)
367      vt_s        ,   &  !: snow total volume per unit area (m)
368      at_i        ,   &  !: ice total fractional area (ice concentration)
369      ato_i       ,   &  !: total open water fractional area (1-at_i)
370      et_i        ,   &  !: total ice heat content
371      et_s        ,   &  !: total snow heat content
372      ot_i        ,   &  !: mean age over all categories
373      tm_i        ,   &  !: mean ice temperature over all categories
374      bv_i        ,   &  !: brine volume averaged over all categories
375      smt_i              !: mean sea ice salinity averaged over all categories
376
377   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpm) ::   &  !:
378      at_i_typ    ,   &  !: total area contained in each ice type
379      vt_i_typ           !: total volume contained in each ice type
380
381   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpsmax) ::   &  !:
382      scal0              !: ???
383
384   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,nlay_s,jpl) :: & !:
385      t_s,            &  !: Snow temperatures (K)
386      e_s
387
388   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) ::   &  !: ! go to trash
389      e_i_cat
390
391   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jkmax,jpl) :: & !:
392      t_i,            &  !: Ice temperatures     [ Kelvins     ]
393      e_i,            &  !: Ice thermal contents [ Joules*10^9 ]
394      s_i                !: Ice salinities
395
396   REAL(wp), DIMENSION(jpi,jpj,0:jpkmax+1) ::    &  !:
397      reslum        !: Relative absorption of solar radiation in each ocean level
398
399   !!--------------------------------------------------------------------------
400   !! * Moments for advection
401   !!--------------------------------------------------------------------------
402   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   &  !:
403      sxopw, syopw, sxxopw, syyopw, sxyopw          !: open water in sea ice
404
405   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) ::   &  !:
406      sxice, syice, sxxice, syyice, sxyice,      &  !: ice thickness moments for advection
407      sxsn,  sysn,  sxxsn,  syysn,  sxysn,       &  !: snow thickness
408      sxa,   sya,   sxxa,   syya,   sxya,        &  !: lead fraction
409      sxc0,  syc0,  sxxc0,  syyc0,  sxyc0,       &  !: snow thermal content
410      sxsal, sysal, sxxsal, syysal, sxysal,      &  !: ice salinity
411      sxage, syage, sxxage, syyage, sxyage          !: ice age
412
413   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jkmax,jpl) ::   &  !:
414      sxe ,  sye ,  sxxe ,  syye ,  sxye            !: ice layers heat content
415
416   !!--------------------------------------------------------------------------
417   !! * Old values of global variables
418   !!--------------------------------------------------------------------------
419   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) ::   &  !:
420      old_v_s, old_v_i,                          &  !: snow and ice volumes
421      old_a_i, old_smv_i, old_oa_i
422   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,nlay_s,jpl) ::   &  !:
423      old_e_s                                       !: snow heat content
424   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jkmax,jpl) ::   &  !:
425      old_e_i                          !: ice temperatures
426   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !: ice velocity (gv6 and gv7)
427      old_u_ice, old_v_ice
428
429   !!--------------------------------------------------------------------------
430   !! * Increment of global variables
431   !!--------------------------------------------------------------------------
432   ! thd refers to changes induced by thermodynamics
433   ! trp   ''         ''     ''       advection (transport of ice)
434   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) ::   &  !:
435      d_a_i_thd  , d_a_i_trp ,                   &  !: icefractions                 
436      d_v_s_thd  , d_v_s_trp,                    &  !: snow volume
437      d_v_i_thd  , d_v_i_trp,                    &  !: ice  volume
438      d_smv_i_thd, d_smv_i_trp,                  &     
439      d_sm_i_fl  , d_sm_i_gd ,                   &
440      d_sm_i_se  , d_sm_i_si , d_sm_i_la ,       &
441      d_oa_i_thd , d_oa_i_trp, s_i_newice
442
443   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,nlay_s,jpl) ::  &  !:
444      d_e_s_thd, d_e_s_trp
445
446   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jkmax,jpl) ::   &  !:
447      d_e_i_thd, d_e_i_trp
448
449   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::       &  !: ice velocity
450      d_u_ice_dyn, d_v_ice_dyn
451
452   !!--------------------------------------------------------------------------
453   !! * Ice thickness distribution variables
454   !!--------------------------------------------------------------------------
455   ! REMOVE
456   INTEGER, PUBLIC, DIMENSION(jpl)                ::   &  !:
457      ice_types      !: Vector making the connection between types and categories
458
459   INTEGER, PUBLIC, DIMENSION(jpm,2)              ::   &  !:
460      ice_cat_bounds !: Matrix containing the integer upper and
461   !: lower boundaries of ice thickness categories
462
463   ! REMOVE
464   INTEGER, PUBLIC, DIMENSION(jpm)                ::   &  !:
465      ice_ncat_types !: Vector containing the number of thickness categories in each ice type
466
467   REAL(wp), PUBLIC, DIMENSION(0:jpl)             ::   &  !:
468      hi_max          !: Boundary of ice thickness categories in thickness space
469
470   REAL(wp), PUBLIC, DIMENSION(jpl)               ::   &  !:
471      hi_mean         !: Mean ice thickness in catgories
472
473   ! REMOVE
474   REAL(wp), PUBLIC, DIMENSION(0:jpl,jpm)         ::   &  !:
475      hi_max_typ     !: Boundary of ice thickness categories
476   !:in thickness space (same but specific for each ice type)
477
478   !!--------------------------------------------------------------------------
479   !! * Ice Run
480   !!--------------------------------------------------------------------------
481   !! Namelist namicerun read in iceini
482   LOGICAL , PUBLIC  ::     & !!! ** init namelist (namicerun) **
483      ln_limdyn   = .TRUE., & !: flag for ice dynamics (T) or not (F)
484      ln_nicep    = .TRUE.    !: flag for sea-ice points output (T) or not (F)
485   REAL(wp), PUBLIC  ::   &  !:
486      hsndif = 0.e0    ,  &  !: computation of temp. in snow (0) or not (9999)
487      hicdif = 0.e0    ,  &  !: computation of temp. in ice (0) or not (9999)
488      cai    = 1.40e-3 ,  &  !: atmospheric drag over sea ice
489      cao    = 1.00e-3       !: atmospheric drag over ocean
490   REAL(wp), PUBLIC, DIMENSION(2)  ::  &  !:
491      acrit  = (/ 1.e-06 , 1.e-06 /)    !: minimum fraction for leads in
492   !                                   !  north and south hemisphere
493
494   !!--------------------------------------------------------------------------
495   !! * Ice diagnostics
496   !!--------------------------------------------------------------------------
497   !! Check if everything down here is necessary
498   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !: volume of ice formed in the leads
499      v_newice
500   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) ::   &  !: thermodynamic growth rates
501      dv_dt_thd,  &
502      izero, fstroc, fhbricat
503   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   & 
504      diag_sni_gr,                           & ! snow ice growth
505      diag_lat_gr,                           & ! lateral ice growth
506      diag_bot_gr,                           & ! bottom ice growth
507      diag_dyn_gr,                           & ! dynamical ice growth
508      diag_bot_me,                           & ! vertical bottom melt
509      diag_sur_me                              ! vertical surface melt
510   INTEGER , PUBLIC ::   jiindx, jjindx        !: indexes of the debugging point
511
512#else
513   !!----------------------------------------------------------------------
514   !!   Default option         Empty module            NO LIM sea-ice model
515   !!----------------------------------------------------------------------
516#endif
517
518   !!======================================================================
519END MODULE ice
Note: See TracBrowser for help on using the repository browser.