source: trunk/NEMO/LIM_SRC/ice.F90 @ 12

Last change on this file since 12 was 12, checked in by opalod, 17 years ago

CT : UPDATE001 : First major NEMO update

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.3 KB
Line 
1MODULE ice
2   !!======================================================================
3   !!                        ***  MODULE ice  ***
4   !! Sea Ice physics:  diagnostics variables of ice defined in memory
5   !!=====================================================================
6#if defined key_ice_lim
7   !!----------------------------------------------------------------------
8   !!   'key_ice_lim' :                                   LIM sea-ice model
9   !!----------------------------------------------------------------------
10   !! History :
11   !!   2.0  !  03-08  (C. Ethe)  F90: Free form and module
12   !!----------------------------------------------------------------------
13   !!  LIM 2.0, UCL-LODYC-IPSL (2002)
14   !!----------------------------------------------------------------------
15   !! * Modules used
16   USE par_ice          ! LIM sea-ice parameters
17
18   IMPLICIT NONE
19   PRIVATE
20
21   !! * Share Module variables
22   INTEGER , PUBLIC ::   & !!: ** ice-dynamic namelist (namicedyn) **
23      nbiter = 1      ,  &  !: number of sub-time steps for relaxation
24      nbitdr = 250          !: maximum number of iterations for relaxation
25
26   REAL(wp), PUBLIC ::   & !!: ** ice-dynamic namelist (namicedyn) **
27      epsd   = 1.0e-20,  &  !: tolerance parameter for dynamic
28      alpha  = 0.5    ,  &  !: coefficient for semi-implicit coriolis
29      dm     = 0.6e+03,  &  !: diffusion constant for dynamics
30      om     = 0.5    ,  &  !: relaxation constant
31      resl   = 5.0e-05,  &  !: maximum value for the residual of relaxation
32      cw     = 5.0e-03,  &  !: drag coefficient for oceanic stress
33      angvg  = 0.0    ,  &  !: turning angle for oceanic stress
34      pstar  = 1.0e+04,  &  !: first bulk-rheology parameter
35      c_rhg  = 20.0   ,  &  !: second bulk-rhelogy parameter
36      etamn  = 0.0e+07,  &  !: minimun value for viscosity
37      creepl = 2.0e-08,  &  !: creep limit
38      ecc    = 2.0    ,  &  !: eccentricity of the elliptical yield curve
39      ahi0   = 350.e0       !: sea-ice hor. eddy diffusivity coeff. (m2/s)
40
41   REAL(wp), PUBLIC ::   &  !:
42      usecc2          ,  &  !:  = 1.0 / ( ecc * ecc )
43      rhoco           ,  &  !: = rau0 * cw
44      sangvg, cangvg  ,  &  !: sin and cos of the turning angle for ocean stress
45      pstarh                !: pstar / 2.0
46
47   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::  &  !:
48      u_oce, v_oce,      &  !: surface ocean velocity used in ice dynamics
49      ahiu , ahiv ,      &  !: hor. diffusivity coeff. at ocean U- and V-points (m2/s)
50      pahu , pahv ,      &  !: ice hor. eddy diffusivity coef. at ocean U- and V-points
51      hsnm , hicm ,      &  !: mean snow and ice thicknesses
52      ust2s                 !: friction velocity
53
54
55   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
56      firic  ,   &  !: IR flux over the ice (only used for outputs)
57      fcsic  ,   &  !: Sensible heat flux over the ice (only used for outputs)
58      fleic  ,   &  !: Latent heat flux over the ice (only used for outputs)
59      qlatic ,   &  !: latent flux
60      rdvosif,   &  !: Variation of volume at surface (only used for outputs)
61      rdvobif,   &  !: Variation of ice volume at the bottom ice (only used for outputs)
62      fdvolif,   &  !: Total variation of ice volume (only used for outputs)
63      rdvonif,   &  !: Lateral Variation of ice volume (only used for outputs)
64      sist   ,   &  !: Sea-Ice Surface Temperature (Kelvin ??? degree ??? I don't know)
65      tfu    ,   &  !: Melting point temperature of sea water
66      hsnif  ,   &  !: Snow thickness
67      hicif  ,   &  !: Ice thickness
68      hicifp ,   &  !: Ice production/melting
69      frld   ,   &  !: Leads fraction = 1-a/totalarea
70      phicif ,   &  !: ice thickness  at previous time
71      pfrld  ,   &  !: Leads fraction at previous time 
72      qstoif ,   &  !: Energy stored in the brine pockets
73      fbif   ,   &  !: Heat flux at the ice base
74      rdmsnif,   &  !: Variation of snow mass
75      rdmicif,   &  !: Variation of ice mass
76      qldif  ,   &  !: heat balance of the lead (or of the open ocean)
77      qcmif  ,   &  !: Energy needed to bring the ocean surface layer until its freezing
78      fdtcn  ,   &  !: net downward heat flux from the ice to the ocean
79      qdtcn  ,   &  !: energy from the ice to the ocean
80      !             !  point (at a factor 2)
81      thcm   ,   &  !: part of the solar energy used in the lead heat budget
82      fstric ,   &  !: Solar flux transmitted trough the ice
83      ffltbif,   &  !: Array linked with the max heat contained in brine pockets (?)
84      fscmbq ,   &  !: Linked with the solar flux below the ice (?)
85      fsbbq  ,   &  !: Also linked with the solar flux below the ice (?)
86      qfvbq  ,   &  !: Array used to store energy in case of toral lateral ablation (?)
87      dmgwi         !: Variation of the mass of snow ice
88
89   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
90      albege ,   &  !: Albedo of the snow or ice (only for outputs)
91      albecn ,   &  !: Albedo of the ocean (only for outputs)
92      tauc   ,   &  !: Cloud optical depth
93      sdvt          !: u*^2/(Stress/density)
94
95
96   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
97      u_ice, v_ice,   &  !: two components of the ice velocity (m/s)
98      tio_u, tio_v       !: two components of the ice-ocean stress (N/m2)
99
100   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,nsmax) ::   &  !:
101      scal0              !: ???
102
103   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,nlayersp1) ::   &  !:
104      tbif          !: Temperature inside the ice/snow layer
105
106   REAL(wp), DIMENSION(jpi,jpj,0:kmax+1) ::    &  !:
107      reslum        !: Relative absorption of solar radiation in each ocean level
108
109   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
110         sxice, syice, sxxice, syyice, sxyice,      &  !: moments for advection
111         sxsn,  sysn,  sxxsn,  syysn,  sxysn,       &  !:
112         sxa,   sya,   sxxa,   syya,   sxya,        &  !:
113         sxc0,  syc0,  sxxc0,  syyc0,  sxyc0,       &  !:
114         sxc1,  syc1,  sxxc1,  syyc1,  sxyc1,       &  !:
115         sxc2,  syc2,  sxxc2,  syyc2,  sxyc2,       &  !:
116         sxst,  syst,  sxxst,  syyst,  sxyst           !:
117
118#else
119   !!----------------------------------------------------------------------
120   !!   Default option         Empty module            NO LIM sea-ice model
121   !!----------------------------------------------------------------------
122#endif
123
124   !!======================================================================
125END MODULE ice
Note: See TracBrowser for help on using the repository browser.