source: NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/cfgs/SHARED/namelist_ice_ref @ 11518

Last change on this file since 11518 was 11518, checked in by clem, 14 months ago

add the final touch to the famous gaston's branch. More precisely, add the possibility to have melt ponds as input file when using bdy

  • Property svn:mime-type set to text/x-fortran
File size: 18.4 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] - L16=0.125 - home=0.15
60      rn_icebfr     =  15.            !        ln_landfast_L16:  maximum bottom stress per unit volume [N/m3]
61                                      !        ln_landfast_home: maximum bottom stress per unit area of contact [N/m2]
62                                      !          recommended range: ?? L16=15 - home=10
63      rn_lfrelax    =   1.e-5         !        relaxation time scale to reach static friction [s-1]
64      rn_tensile    =   0.2           !        ln_landfast_L16: isotropic tensile strength
65/
66!------------------------------------------------------------------------------
67&namdyn_rdgrft  !   Ice ridging/rafting
68!------------------------------------------------------------------------------
69          ! -- ice_rdgrft_strength -- !
70   ln_str_H79       = .true.          !  ice strength param.: Hibler_79   => P = pstar*<h>*exp(-c_rhg*A)                     
71      rn_pstar      =   2.0e+04       !     ice strength thickness parameter [N/m2]
72      rn_crhg       =   20.0          !     ice strength conc. parameter (-)
73                   ! -- ice_rdgrft -- !
74   rn_csrdg         =   0.5           !  fraction of shearing energy contributing to ridging
75              ! -- ice_rdgrft_prep -- !
76   ln_partf_lin     = .false.         !  Linear ridging participation function (Thorndike et al, 1975)
77      rn_gstar      =   0.15          !     fractional area of thin ice being ridged
78   ln_partf_exp     = .true.          !  Exponential ridging participation function (Lipscomb, 2007)
79      rn_astar      =   0.03          !     exponential measure of ridging ice fraction [set to 0.05 if hstar=100]
80   ln_ridging       = .true.          !  ridging activated (T) or not (F)
81      rn_hstar      =  25.0           !     determines the maximum thickness of ridged ice [m] (Hibler, 1980)
82      rn_porordg    =   0.3           !     porosity of newly ridged ice (Lepparanta et al., 1995)
83      rn_fsnwrdg    =   0.5           !     snow volume fraction that survives in ridging
84      rn_fpndrdg    =   1.0           !     pond fraction that survives in ridging (small a priori)
85   ln_rafting       = .true.          !  rafting activated (T) or not (F)
86      rn_hraft      =   0.75          !     threshold thickness for rafting [m]
87      rn_craft      =   5.0           !     squeezing coefficient used in the rafting function
88      rn_fsnwrft    =   0.5           !     snow volume fraction that survives in rafting
89      rn_fpndrft    =   1.0           !     pond fraction that survives in rafting (0.5 a priori)
90/
91!------------------------------------------------------------------------------
92&namdyn_rhg     !   Ice rheology
93!------------------------------------------------------------------------------
94   ln_rhg_EVP       = .true.          !  EVP rheology
95      ln_aEVP       = .false.         !     adaptive rheology (Kimmritz et al. 2016 & 2017)
96      rn_creepl     =   2.0e-9        !     creep limit [1/s]
97      rn_ecc        =   2.0           !     eccentricity of the elliptical yield curve         
98      nn_nevp       = 120             !     number of EVP subcycles                             
99      rn_relast     =   0.333         !     ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast
100                                      !        advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300)
101/
102!------------------------------------------------------------------------------
103&namdyn_adv     !   Ice advection
104!------------------------------------------------------------------------------
105   ln_adv_Pra       = .false.         !  Advection scheme (Prather)
106   ln_adv_UMx       = .true.          !  Advection scheme (Ultimate-Macho)
107      nn_UMx        =   5             !     order of the scheme for UMx (1-5 ; 20=centered 2nd order)
108/
109!------------------------------------------------------------------------------
110&namsbc         !   Ice surface boundary conditions
111!------------------------------------------------------------------------------
112   rn_cio           =   5.0e-03       !  ice-ocean drag coefficient (-)
113   rn_blow_s        =   0.66          !  mesure of snow blowing into the leads
114                                      !     = 1 => no snow blowing, < 1 => some snow blowing
115   nn_flxdist       =  -1             !  Redistribute heat flux over ice categories
116                                      !     =-1  Do nothing (needs N(cat) fluxes)
117                                      !     = 0  Average N(cat) fluxes then apply the average over the N(cat) ice
118                                      !     = 1  Average N(cat) fluxes then redistribute over the N(cat) ice using T-ice and albedo sensitivity
119                                      !     = 2  Redistribute a single flux over categories
120   ln_cndflx        = .false.         !  Use conduction flux as surface boundary conditions (i.e. for Jules coupling)
121      ln_cndemulate = .false.         !     emulate conduction flux (if not provided in the inputs)
122/
123!------------------------------------------------------------------------------
124&namthd         !   Ice thermodynamics
125!------------------------------------------------------------------------------
126   ln_icedH         = .true.          !  activate ice thickness change from growing/melting (T) or not (F)
127   ln_icedA         = .true.          !  activate lateral melting param. (T) or not (F)
128   ln_icedO         = .true.          !  activate ice growth in open-water (T) or not (F)
129   ln_icedS         = .true.          !  activate brine drainage (T) or not (F)
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 Holland et al 2012)
181     ln_pnd_CST     = .false.         !  activate constant  melt ponds
182       rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC
183       rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC
184     ln_pnd_alb     = .false.         !  melt ponds affect albedo or not
185/
186!------------------------------------------------------------------------------
187&namini         !   Ice initialization
188!------------------------------------------------------------------------------
189   ln_iceini        = .true.          !  activate ice initialization (T) or not (F)
190   ln_iceini_file   = .false.         !  netcdf file provided for initialization (T) or not (F)
191   rn_thres_sst     =   2.0           !  max temp. above Tfreeze with initial ice = (sst - tfreeze)
192   rn_hti_ini_n     =   3.0           !  initial ice thickness       (m), North
193   rn_hti_ini_s     =   1.0           !        "            "             South
194   rn_hts_ini_n     =   0.3           !  initial snow thickness      (m), North
195   rn_hts_ini_s     =   0.3           !        "            "             South
196   rn_ati_ini_n     =   0.9           !  initial ice concentration   (-), North
197   rn_ati_ini_s     =   0.9           !        "            "             South
198   rn_smi_ini_n     =   6.3           !  initial ice salinity     (g/kg), North
199   rn_smi_ini_s     =   6.3           !        "            "             South
200   rn_tmi_ini_n     = 270.            !  initial ice temperature     (K), North
201   rn_tmi_ini_s     = 270.            !        "            "             South
202   rn_tsu_ini_n     = 270.            !  initial surface temperature (K), North
203   rn_tsu_ini_s     = 270.            !        "            "             South
204   rn_tms_ini_n     = 270.            !  initial snw temperature     (K), North
205   rn_tms_ini_s     = 270.            !        "            "             South
206   rn_apd_ini_n     =   0.2           !  initial pond fraction       (-), North
207   rn_apd_ini_s     =   0.2           !        "            "             South
208   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North
209   rn_hpd_ini_s     =   0.05          !        "            "             South
210   ! -- for ln_iceini_file = T
211   sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', ''
212   sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', ''
213   sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', ''
214   sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', ''
215   sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', ''
216   sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', ''
217   sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', ''
218   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd)
219   sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
220   sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
221   cn_dir='./'
222/
223!------------------------------------------------------------------------------
224&namalb         !   albedo parameters
225!------------------------------------------------------------------------------
226   !                                  !                          !  obs range (cloud-sky)
227   rn_alb_sdry      =   0.85          !  dry snow albedo         :  0.85 -- 0.87
228   rn_alb_smlt      =   0.75          !  melting snow albedo     :  0.72 -- 0.82
229   rn_alb_idry      =   0.60          !  dry ice albedo          :  0.54 -- 0.65
230   rn_alb_imlt      =   0.50          !  bare puddled ice albedo :  0.49 -- 0.58
231   rn_alb_dpnd      =   0.27          !  ponded ice albedo       :  0.10 -- 0.30
232/
233!------------------------------------------------------------------------------
234&namdia         !   Diagnostics
235!------------------------------------------------------------------------------
236   ln_icediachk     = .false.         !  check online the heat, mass & salt budgets at each time step
237      !                               !     rate of ice spuriously gained/lost. For ex., rn_icechk=1. <=> 1mm/year, rn_icechk=0.1 <=> 1mm/10years                                   
238      rn_icechk_cel =  1.             !     check at any gridcell           => stops the code if violated (and writes a file)
239      rn_icechk_glo =  0.1            !     check over the entire ice cover => only prints warnings
240   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F)
241   ln_icectl        = .false.         !  ice points output for debug (T or F)
242   iiceprt          =  10             !  i-index for debug
243   jiceprt          =  10             !  j-index for debug
244/
Note: See TracBrowser for help on using the repository browser.