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

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

Initial revision

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