source: NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/cfgs/SHARED/namelist_ice_ref @ 12734

Last change on this file since 12734 was 12734, checked in by clem, 19 months ago

add the possibility to not use the heat contained in the leads to melt sea ice but wait for the ocean to warm up first (ln_leadhfx)

  • Property svn:mime-type set to text/x-fortran
File size: 19.0 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! SI3 namelist: 
3!!              1 - Generic parameters                 (nampar)
4!!              2 - Ice thickness discretization       (namitd)
5!!              3 - Ice dynamics                       (namdyn)
6!!              4 - Ice ridging/rafting                (namdyn_rdgrft)
7!!              5 - Ice rheology                       (namdyn_rhg)
8!!              6 - Ice advection                      (namdyn_adv)
9!!              7 - Ice surface boundary conditions    (namsbc)
10!!              8 - Ice thermodynamics                 (namthd)
11!!              9 - Ice heat diffusion                 (namthd_zdf)
12!!             10 - Ice lateral melting                (namthd_da)
13!!             11 - Ice growth in open water           (namthd_do)
14!!             12 - Ice salinity                       (namthd_sal)
15!!             13 - Ice melt ponds                     (namthd_pnd)
16!!             14 - Ice initialization                 (namini)
17!!             15 - Ice/snow albedos                   (namalb)
18!!             16 - Ice diagnostics                    (namdia)
19!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
20!
21!------------------------------------------------------------------------------
22&nampar         !   Generic parameters
23!------------------------------------------------------------------------------
24   jpl              =   5             !  number of ice  categories
25   nlay_i           =   2             !  number of ice  layers
26   nlay_s           =   1             !  number of snow layers (only 1 is working)
27   ln_virtual_itd   =   .false.       !  virtual ITD mono-category parameterization (jpl=1 only)
28                                      !     i.e. enhanced thermal conductivity & virtual thin ice melting
29   ln_icedyn        = .true.          !  ice dynamics (T) or not (F)
30   ln_icethd        = .true.          !  ice thermo   (T) or not (F)
31   rn_amax_n        =   0.997         !  maximum tolerated ice concentration NH
32   rn_amax_s        =   0.997         !  maximum tolerated ice concentration SH
33   cn_icerst_in     = "restart_ice"   !  suffix of ice restart name (input)
34   cn_icerst_out    = "restart_ice"   !  suffix of ice restart name (output)
35   cn_icerst_indir  = "."             !  directory to read   input ice restarts
36   cn_icerst_outdir = "."             !  directory to write output ice restarts
37/
38!------------------------------------------------------------------------------
39&namitd         !   Ice discretization
40!------------------------------------------------------------------------------
41   ln_cat_hfn       = .true.          !  ice categories are defined by a function following rn_himean**(-0.05)
42      rn_himean     =   2.0           !  expected domain-average ice thickness (m)
43   ln_cat_usr       = .false.         !  ice categories are defined by rn_catbnd below (m)
44      rn_catbnd     =   0.,0.45,1.1,2.1,3.7,6.0 
45   rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping
46/
47!------------------------------------------------------------------------------
48&namdyn         !   Ice dynamics
49!------------------------------------------------------------------------------
50   ln_dynALL        = .true.          !  dyn.: full ice dynamics                  (rheology + advection + ridging/rafting + correction)
51   ln_dynRHGADV     = .false.         !  dyn.: no ridge/raft & no corrections     (rheology + advection)
52   ln_dynADV1D      = .false.         !  dyn.: only advection 1D                  (Schar & Smolarkiewicz 1996 test case)
53   ln_dynADV2D      = .false.         !  dyn.: only advection 2D w prescribed vel.(rn_uvice + advection)
54      rn_uice       =   0.5           !        prescribed ice u-velocity
55      rn_vice       =   0.5           !        prescribed ice v-velocity
56   rn_ishlat        =   2.            !  lbc : free slip (0) ; partial slip (0-2) ; no slip (2) ; strong slip (>2)
57   ln_landfast_L16  = .false.         !  landfast: parameterization from Lemieux 2016
58      rn_depfra     =   0.125         !        fraction of ocean depth that ice must reach to initiate landfast
59                                      !          recommended range: [0.1 ; 0.25]
60      rn_icebfr     =  15.            !        maximum bottom stress per unit volume [N/m3]
61      rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1]
62      rn_tensile    =   0.05          !        isotropic tensile strength [0-0.5??]
63/
64!------------------------------------------------------------------------------
65&namdyn_rdgrft  !   Ice ridging/rafting
66!------------------------------------------------------------------------------
67          ! -- ice_rdgrft_strength -- !
68   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                     
69      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2]
70      rn_crhg       =   20.0          !     ice strength conc. parameter (-)
71                   ! -- ice_rdgrft -- !
72   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging
73              ! -- ice_rdgrft_prep -- !
74   ln_partf_lin     = .false.         !  Linear ridging participation function (Thorndike et al, 1975)
75      rn_gstar      =   0.15          !     fractional area of thin ice being ridged
76   ln_partf_exp     = .true.          !  Exponential ridging participation function (Lipscomb, 2007)
77      rn_astar      =   0.03          !     exponential measure of ridging ice fraction [set to 0.05 if hstar=100]
78   ln_ridging       = .true.          !  ridging activated (T) or not (F)
79      rn_hstar      =  25.0           !     determines the maximum thickness of ridged ice [m] (Hibler, 1980)
80      rn_porordg    =   0.3           !     porosity of newly ridged ice (Lepparanta et al., 1995)
81      rn_fsnwrdg    =   0.5           !     snow volume fraction that survives in ridging
82      rn_fpndrdg    =   1.0           !     pond fraction that survives in ridging (small a priori)
83   ln_rafting       = .true.          !  rafting activated (T) or not (F)
84      rn_hraft      =   0.75          !     threshold thickness for rafting [m]
85      rn_craft      =   5.0           !     squeezing coefficient used in the rafting function
86      rn_fsnwrft    =   0.5           !     snow volume fraction that survives in rafting
87      rn_fpndrft    =   1.0           !     pond fraction that survives in rafting (0.5 a priori)
88/
89!------------------------------------------------------------------------------
90&namdyn_rhg     !   Ice rheology
91!------------------------------------------------------------------------------
92   ln_rhg_EVP       = .true.          !  EVP rheology
93      ln_aEVP       = .true.          !     adaptive rheology (Kimmritz et al. 2016 & 2017)
94      rn_creepl     =   2.0e-9        !     creep limit [1/s]
95      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve         
96      nn_nevp       = 120             !     number of EVP subcycles                             
97      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast
98                                      !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300)
99/
100!------------------------------------------------------------------------------
101&namdyn_adv     !   Ice advection
102!------------------------------------------------------------------------------
103   ln_adv_Pra       = .true.         !  Advection scheme (Prather)
104   ln_adv_UMx       = .false.          !  Advection scheme (Ultimate-Macho)
105      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order)
106/
107!------------------------------------------------------------------------------
108&namsbc         !   Ice surface boundary conditions
109!------------------------------------------------------------------------------
110   rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-)
111   rn_blow_s        =   0.66          !  mesure of snow blowing into the leads
112                                      !     = 1 => no snow blowing, < 1 => some snow blowing
113   nn_flxdist       =  -1             !  Redistribute heat flux over ice categories
114                                      !     =-1  Do nothing (needs N(cat) fluxes)
115                                      !     = 0  Average N(cat) fluxes then apply the average over the N(cat) ice
116                                      !     = 1  Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity
117                                      !     = 2  Redistribute a single flux over categories
118   ln_cndflx        = .false.         !  Use conduction flux as surface boundary conditions (i.e. for Jules coupling)
119      ln_cndemulate = .false.         !     emulate conduction flux (if not provided in the inputs)
120/
121!------------------------------------------------------------------------------
122&namthd         !   Ice thermodynamics
123!------------------------------------------------------------------------------
124   ln_icedH         = .true.          !  activate ice thickness change from growing/melting (T) or not (F)
125   ln_icedA         = .true.          !  activate lateral melting param. (T) or not (F)
126   ln_icedO         = .true.          !  activate ice growth in open-water (T) or not (F)
127   ln_icedS         = .true.          !  activate brine drainage (T) or not (F)
128   !
129   ln_leadhfx       = .true.          !  heat in the leads is used to melt sea-ice before warming the ocean
130/
131!------------------------------------------------------------------------------
132&namthd_zdf     !   Ice heat diffusion
133!------------------------------------------------------------------------------
134   ln_zdf_BL99      = .true.          !  Heat diffusion follows Bitz and Lipscomb 1999
135   ln_cndi_U64      = .false.         !  sea ice thermal conductivity: k = k0 + beta.S/T            (Untersteiner, 1964)
136   ln_cndi_P07      = .true.          !  sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007)
137   rn_cnd_s         =   0.31          !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971)
138                                      !     Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)
139   rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice [1/m]
140/
141!------------------------------------------------------------------------------
142&namthd_da      !   Ice lateral melting
143!------------------------------------------------------------------------------
144   rn_beta          =   1.0           !  coef. beta for lateral melting param. Recommended range=[0.8-1.2]
145                                      !   => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2)
146                                      !      0.3 = best fit for western Fram Strait and Antarctica
147                                      !      1.4 = best fit for eastern Fram Strait     
148   rn_dmin          =   8.            !  minimum floe diameter for lateral melting param. Recommended range=[6-10]
149                                      !   => 6  vs 8m = +40% melting at the peak (A~0.5)
150                                      !      10 vs 8m = -20% melting
151/
152!------------------------------------------------------------------------------
153&namthd_do      !   Ice growth in open water
154!------------------------------------------------------------------------------
155   rn_hinew         =   0.1           !  thickness for new ice formation in open water (m), must be larger than rn_himin
156   ln_frazil        = .false.         !  Frazil ice parameterization (ice collection as a function of wind)
157      rn_maxfraz    =   1.0           !     maximum fraction of frazil ice collecting at the ice base
158      rn_vfraz      =   0.417         !     thresold drift speed for frazil ice collecting at the ice bottom (m/s)
159      rn_Cfraz      =   5.0           !     squeezing coefficient for frazil ice collecting at the ice bottom
160/
161!------------------------------------------------------------------------------
162&namthd_sal     !   Ice salinity
163!------------------------------------------------------------------------------
164   nn_icesal        =   2             !  ice salinity option
165                                      !     1: constant ice salinity (S=rn_icesal)
166                                      !     2: varying salinity parameterization S(z,t)
167                                      !     3: prescribed salinity profile S(z) (Schwarzacher 1959)
168   rn_icesal        =   4.            !      (nn_icesal=1) ice salinity (g/kg)
169   rn_sal_gd        =   5.            !      (nn_icesal=2) restoring ice salinity, gravity drainage (g/kg)
170   rn_time_gd       =   1.73e+6       !      (nn_icesal=2) restoring time scale,   gravity drainage  (s)
171   rn_sal_fl        =   2.            !      (nn_icesal=2) restoring ice salinity, flushing (g/kg)
172   rn_time_fl       =   8.64e+5       !      (nn_icesal=2) restoring time scale,   flushing (s)
173   rn_simax         =  20.            !  maximum tolerated ice salinity (g/kg)
174   rn_simin         =   0.1           !  minimum tolerated ice salinity (g/kg)
175/
176!------------------------------------------------------------------------------
177&namthd_pnd     !   Melt ponds
178!------------------------------------------------------------------------------
179   ln_pnd           = .false.         !  activate melt ponds or not
180     ln_pnd_H12     = .false.         !  activate evolutive melt ponds (from Flocco et al 2007,2010 & Holland et al 2012)
181       ln_pnd_lids  = .true.          !  ponds with frozen lids
182       ln_pnd_flush = .true.          !  ponds flushing trhu the ice 
183       rn_apnd_min  =   0.15          !  minimum ice fraction that contributes to melt pond. range: 0.0 -- 0.15 ??
184       rn_apnd_max  =   0.85          !  maximum ice fraction that contributes to melt pond. range: 0.7 -- 0.85 ??
185     ln_pnd_CST     = .false.         !  activate constant  melt ponds
186       rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC
187       rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC
188     ln_pnd_alb     = .false.         !  melt ponds affect albedo or not
189/
190!------------------------------------------------------------------------------
191&namini         !   Ice initialization
192!------------------------------------------------------------------------------
193   ln_iceini        = .true.          !  activate ice initialization (T) or not (F)
194   nn_iceini_file   =   0             !     0 = Initialise sea ice based on SSTs
195                                      !     1 = Initialise sea ice from single category netcdf file
196                                      !     2 = Initialise sea ice from multi category restart file
197   rn_thres_sst     =   2.0           !  max temp. above Tfreeze with initial ice = (sst - tfreeze)
198   rn_hti_ini_n     =   3.0           !  initial ice thickness       (m), North
199   rn_hti_ini_s     =   1.0           !        "            "             South
200   rn_hts_ini_n     =   0.3           !  initial snow thickness      (m), North
201   rn_hts_ini_s     =   0.3           !        "            "             South
202   rn_ati_ini_n     =   0.9           !  initial ice concentration   (-), North
203   rn_ati_ini_s     =   0.9           !        "            "             South
204   rn_smi_ini_n     =   6.3           !  initial ice salinity     (g/kg), North
205   rn_smi_ini_s     =   6.3           !        "            "             South
206   rn_tmi_ini_n     = 270.            !  initial ice temperature     (K), North
207   rn_tmi_ini_s     = 270.            !        "            "             South
208   rn_tsu_ini_n     = 270.            !  initial surface temperature (K), North
209   rn_tsu_ini_s     = 270.            !        "            "             South
210   rn_tms_ini_n     = 270.            !  initial snw temperature     (K), North
211   rn_tms_ini_s     = 270.            !        "            "             South
212   rn_apd_ini_n     =   0.2           !  initial pond fraction       (-), North
213   rn_apd_ini_s     =   0.2           !        "            "             South
214   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North
215   rn_hpd_ini_s     =   0.05          !        "            "             South
216   rn_hld_ini_n     =   0.0           !  initial pond lid depth      (m), North
217   rn_hld_ini_s     =   0.0           !        "            "             South
218   ! -- for nn_iceini_file = 1
219   sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', ''
220   sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', ''
221   sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', ''
222   sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', ''
223   sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', ''
224   sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', ''
225   sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', ''
226   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd)
227   sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
228   sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
229   sn_hld = 'NOT USED'              , -12 ,'hld'   ,  .false.  , .true., 'yearly'  , '' , '', ''
230   cn_dir='./'
231/
232!------------------------------------------------------------------------------
233&namalb         !   albedo parameters
234!------------------------------------------------------------------------------
235   !                                  !                          !  obs range (cloud-sky)
236   rn_alb_sdry      =   0.85          !  dry snow albedo         :  0.85 -- 0.87
237   rn_alb_smlt      =   0.75          !  melting snow albedo     :  0.72 -- 0.82
238   rn_alb_idry      =   0.60          !  dry ice albedo          :  0.54 -- 0.65
239   rn_alb_imlt      =   0.50          !  bare puddled ice albedo :  0.49 -- 0.58
240   rn_alb_dpnd      =   0.27          !  ponded ice albedo       :  0.10 -- 0.30
241/
242!------------------------------------------------------------------------------
243&namdia         !   Diagnostics
244!------------------------------------------------------------------------------
245   ln_icediachk     = .false.         !  check online heat, mass & salt budgets
246      !                               !   rate of ice spuriously gained/lost at each time step => rn_icechk=1 <=> 1.e-6 m/hour
247      rn_icechk_cel =  100.           !     check at each gridcell          (1.e-4m/h)=> stops the code if violated (and writes a file)
248      rn_icechk_glo =  1.             !     check over the entire ice cover (1.e-6m/h)=> only prints warnings
249   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F)
250   ln_icectl        = .false.         !  ice points output for debug (T or F)
251   iiceprt          =  10             !  i-index for debug
252   jiceprt          =  10             !  j-index for debug
253/
Note: See TracBrowser for help on using the repository browser.