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.
namelist_ice_ref in NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/cfgs/SHARED – NEMO

source: NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/cfgs/SHARED/namelist_ice_ref @ 11777

Last change on this file since 11777 was 11777, checked in by acc, 4 years ago

Branch 2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps. Begin merge prepartions by merging changes on the trunk since this branch creation. Do this in staged commits in case of issues. First stage is to merge and commit changes to cfgs directory (cd cfgs; svn merge -r 10721:11740 svn+ssh://acc@forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/NEMO/trunk/cfgs ./)

  • Property svn:mime-type set to text/x-fortran
File size: 18.1 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!------------------------------------------------------------------------------
130&namthd_zdf     !   Ice heat diffusion
131!------------------------------------------------------------------------------
132   ln_zdf_BL99      = .true.          !  Heat diffusion follows Bitz and Lipscomb 1999
133   ln_cndi_U64      = .false.         !  sea ice thermal conductivity: k = k0 + beta.S/T            (Untersteiner, 1964)
134   ln_cndi_P07      = .true.          !  sea ice thermal conductivity: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007)
135   rn_cnd_s         =   0.31          !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971)
136                                      !     Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)
137   rn_kappa_i       =   1.0           !  radiation attenuation coefficient in sea ice [1/m]
138/
139!------------------------------------------------------------------------------
140&namthd_da      !   Ice lateral melting
141!------------------------------------------------------------------------------
142   rn_beta          =   1.0           !  coef. beta for lateral melting param. Recommended range=[0.8-1.2]
143                                      !   => decrease = more melt and melt peaks toward higher concentration (A~0.5 for beta=1 ; A~0.8 for beta=0.2)
144                                      !      0.3 = best fit for western Fram Strait and Antarctica
145                                      !      1.4 = best fit for eastern Fram Strait     
146   rn_dmin          =   8.            !  minimum floe diameter for lateral melting param. Recommended range=[6-10]
147                                      !   => 6  vs 8m = +40% melting at the peak (A~0.5)
148                                      !      10 vs 8m = -20% melting
149/
150!------------------------------------------------------------------------------
151&namthd_do      !   Ice growth in open water
152!------------------------------------------------------------------------------
153   rn_hinew         =   0.1           !  thickness for new ice formation in open water (m), must be larger than rn_himin
154   ln_frazil        = .false.         !  Frazil ice parameterization (ice collection as a function of wind)
155      rn_maxfraz    =   1.0           !     maximum fraction of frazil ice collecting at the ice base
156      rn_vfraz      =   0.417         !     thresold drift speed for frazil ice collecting at the ice bottom (m/s)
157      rn_Cfraz      =   5.0           !     squeezing coefficient for frazil ice collecting at the ice bottom
158/
159!------------------------------------------------------------------------------
160&namthd_sal     !   Ice salinity
161!------------------------------------------------------------------------------
162   nn_icesal        =   2             !  ice salinity option
163                                      !     1: constant ice salinity (S=rn_icesal)
164                                      !     2: varying salinity parameterization S(z,t)
165                                      !     3: prescribed salinity profile S(z) (Schwarzacher 1959)
166   rn_icesal        =   4.            !      (nn_icesal=1) ice salinity (g/kg)
167   rn_sal_gd        =   5.            !      (nn_icesal=2) restoring ice salinity, gravity drainage (g/kg)
168   rn_time_gd       =   1.73e+6       !      (nn_icesal=2) restoring time scale,   gravity drainage  (s)
169   rn_sal_fl        =   2.            !      (nn_icesal=2) restoring ice salinity, flushing (g/kg)
170   rn_time_fl       =   8.64e+5       !      (nn_icesal=2) restoring time scale,   flushing (s)
171   rn_simax         =  20.            !  maximum tolerated ice salinity (g/kg)
172   rn_simin         =   0.1           !  minimum tolerated ice salinity (g/kg)
173/
174!------------------------------------------------------------------------------
175&namthd_pnd     !   Melt ponds
176!------------------------------------------------------------------------------
177   ln_pnd           = .false.         !  activate melt ponds or not
178     ln_pnd_H12     = .false.         !  activate evolutive melt ponds (from Holland et al 2012)
179     ln_pnd_CST     = .false.         !  activate constant  melt ponds
180       rn_apnd      =   0.2           !     prescribed pond fraction, at Tsu=0 degC
181       rn_hpnd      =   0.05          !     prescribed pond depth,    at Tsu=0 degC
182     ln_pnd_alb     = .false.         !  melt ponds affect albedo or not
183/
184!------------------------------------------------------------------------------
185&namini         !   Ice initialization
186!------------------------------------------------------------------------------
187   ln_iceini        = .true.          !  activate ice initialization (T) or not (F)
188   ln_iceini_file   = .false.         !  netcdf file provided for initialization (T) or not (F)
189   rn_thres_sst     =   2.0           !  max temp. above Tfreeze with initial ice = (sst - tfreeze)
190   rn_hti_ini_n     =   3.0           !  initial ice thickness       (m), North
191   rn_hti_ini_s     =   1.0           !        "            "             South
192   rn_hts_ini_n     =   0.3           !  initial snow thickness      (m), North
193   rn_hts_ini_s     =   0.3           !        "            "             South
194   rn_ati_ini_n     =   0.9           !  initial ice concentration   (-), North
195   rn_ati_ini_s     =   0.9           !        "            "             South
196   rn_smi_ini_n     =   6.3           !  initial ice salinity     (g/kg), North
197   rn_smi_ini_s     =   6.3           !        "            "             South
198   rn_tmi_ini_n     = 270.            !  initial ice temperature     (K), North
199   rn_tmi_ini_s     = 270.            !        "            "             South
200   rn_tsu_ini_n     = 270.            !  initial surface temperature (K), North
201   rn_tsu_ini_s     = 270.            !        "            "             South
202   rn_tms_ini_n     = 270.            !  initial snw temperature     (K), North
203   rn_tms_ini_s     = 270.            !        "            "             South
204   rn_apd_ini_n     =   0.2           !  initial pond fraction       (-), North
205   rn_apd_ini_s     =   0.2           !        "            "             South
206   rn_hpd_ini_n     =   0.05          !  initial pond depth          (m), North
207   rn_hpd_ini_s     =   0.05          !        "            "             South
208   ! -- for ln_iceini_file = T
209   sn_hti = 'Ice_initialization'    , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', ''
210   sn_hts = 'Ice_initialization'    , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', ''
211   sn_ati = 'Ice_initialization'    , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', ''
212   sn_smi = 'Ice_initialization'    , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', ''
213   sn_tmi = 'Ice_initialization'    , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', ''
214   sn_tsu = 'Ice_initialization'    , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', ''
215   sn_tms = 'NOT USED'              , -12 ,'tms'   ,  .false.  , .true., 'yearly'  , '' , '', ''
216   !      melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd)
217   sn_apd = 'NOT USED'              , -12 ,'apd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
218   sn_hpd = 'NOT USED'              , -12 ,'hpd'   ,  .false.  , .true., 'yearly'  , '' , '', ''
219   cn_dir='./'
220/
221!------------------------------------------------------------------------------
222&namalb         !   albedo parameters
223!------------------------------------------------------------------------------
224   !                                  !                          !  obs range (cloud-sky)
225   rn_alb_sdry      =   0.85          !  dry snow albedo         :  0.85 -- 0.87
226   rn_alb_smlt      =   0.75          !  melting snow albedo     :  0.72 -- 0.82
227   rn_alb_idry      =   0.60          !  dry ice albedo          :  0.54 -- 0.65
228   rn_alb_imlt      =   0.50          !  bare puddled ice albedo :  0.49 -- 0.58
229   rn_alb_dpnd      =   0.27          !  ponded ice albedo       :  0.10 -- 0.30
230/
231!------------------------------------------------------------------------------
232&namdia         !   Diagnostics
233!------------------------------------------------------------------------------
234   ln_icediachk     = .false.         !  check online heat, mass & salt budgets
235      !                               !   rate of ice spuriously gained/lost at each time step => rn_icechk=1 <=> 1.e-6 m/hour
236      rn_icechk_cel =  100.           !     check at each gridcell          (1.e-4m/h)=> stops the code if violated (and writes a file)
237      rn_icechk_glo =  1.             !     check over the entire ice cover (1.e-6m/h)=> only prints warnings
238   ln_icediahsb     = .false.         !  output the heat, mass & salt budgets (T) or not (F)
239   ln_icectl        = .false.         !  ice points output for debug (T or F)
240   iiceprt          =  10             !  i-index for debug
241   jiceprt          =  10             !  j-index for debug
242/
Note: See TracBrowser for help on using the repository browser.