1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
2 | !! NEMO/LIM3 : 1 - dynamics/advection/thermo (namicerun) |
---|
3 | !! namelists 2 - ice intialisation (namiceini) |
---|
4 | !! 3 - ice dynamics (namicedyn) |
---|
5 | !! 4 - ice advection (namicetrp) |
---|
6 | !! 5 - thermodynamics (namicethd) |
---|
7 | !! 6 - ice salinity (namicesal) |
---|
8 | !! 7 - mechanical redistribution of ice (namiceitdme) |
---|
9 | !! 8 - ice diagnostics (namicedia) |
---|
10 | !! 9 - ice outputs (namiceout) |
---|
11 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
12 | |
---|
13 | !----------------------------------------------------------------------- |
---|
14 | &namicerun ! Share parameters for dynamics/advection/thermo |
---|
15 | !----------------------------------------------------------------------- |
---|
16 | jpl = 5 ! number of ice categories |
---|
17 | nlay_i = 2 ! number of ice layers |
---|
18 | nlay_s = 1 ! number of snow layers |
---|
19 | cn_icerst_in = "restart_ice" ! suffix of ice restart name (input) |
---|
20 | cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) |
---|
21 | ln_limdyn = .true. ! ice dynamics (T) or thermodynamics only (F) |
---|
22 | amax = 0.999 ! maximum tolerated ice concentration |
---|
23 | ln_nicep = .false. ! ice points output for debug (yes or no) |
---|
24 | ln_limdiahsb = .false. ! check the heat and salt budgets (T) or not (F) |
---|
25 | ln_limdiaout = .true. ! output the heat and salt budgets (T) or not (F) |
---|
26 | / |
---|
27 | !----------------------------------------------------------------------- |
---|
28 | &namiceitd ! ice thickness distribution |
---|
29 | !----------------------------------------------------------------------- |
---|
30 | nn_itdshp = 1 ! shape of ice categories distribution following: |
---|
31 | ! 1 = tanh function |
---|
32 | ! 2 = h^(-alpha) |
---|
33 | rn_itmean = 2.5 ! mean thickness of the domain (used to compute the distribution, nn_itdshp = 2 only) |
---|
34 | / |
---|
35 | !----------------------------------------------------------------------- |
---|
36 | &namiceini ! ice initialization |
---|
37 | !----------------------------------------------------------------------- |
---|
38 | ln_limini = .true. ! activate ice initialization (T) or not (F) |
---|
39 | thres_sst = 2.0 ! threshold water temperature for initial sea ice |
---|
40 | hts_ini_n = 0.3 ! initial snow thickness (North) |
---|
41 | hts_ini_s = 0.3 ! " " (South) |
---|
42 | hti_ini_n = 3.0 ! initial ice thickness (North) |
---|
43 | hti_ini_s = 1.0 ! " " (South) |
---|
44 | ati_ini_n = 0.9 ! initial ice concentration (North) |
---|
45 | ati_ini_s = 0.9 ! " " (South) |
---|
46 | smi_ini_n = 6.3 ! initial ice salinity (North) |
---|
47 | smi_ini_s = 6.3 ! " " (South) |
---|
48 | tmi_ini_n = 270. ! initial ice/snw temp (North) |
---|
49 | tmi_ini_s = 270. ! initial ice/snw temp (South) |
---|
50 | / |
---|
51 | !----------------------------------------------------------------------- |
---|
52 | &namicedyn ! ice dynamics |
---|
53 | !----------------------------------------------------------------------- |
---|
54 | cw = 5.0e-03 ! ice-ocean drag coefficient |
---|
55 | pstar = 2.0e+04 ! 1st bulk-rheology parameter (N/m) |
---|
56 | c_rhg = 20.0 ! 2nd bulk-rhelogy parameter |
---|
57 | creepl = 1.0e-12 ! creep limit (s-1) |
---|
58 | ecc = 2.0 ! eccentricity of the elliptical yield curve |
---|
59 | ahi0 = 350.e0 ! horizontal eddy diffusivity coefficient for sea-ice (m2/s) - DO NOT CHANGE IT |
---|
60 | ! the value is set up for orca2 and is adapted automatically in the code |
---|
61 | ! depending on the resolution of your experiment |
---|
62 | nevp = 120 ! number of EVP subcycles |
---|
63 | relast = 0.333 ! ratio of elastic timescale to ice time step (1/3 if nevp=120, 1/9 if nevp=300) |
---|
64 | / |
---|
65 | !----------------------------------------------------------------------- |
---|
66 | &namicethd ! ice thermodynamics |
---|
67 | !----------------------------------------------------------------------- |
---|
68 | hiccrit = 0.1 ! thickness for new ice formation in open water |
---|
69 | fraz_swi = 0 ! use frazil ice collection thickness as a function of wind (1.0) or not (0.0) |
---|
70 | maxfrazb = 0.0 ! maximum portion of frazil ice collecting at the ice bottom |
---|
71 | vfrazb = 0.417 ! thresold drift speed for frazil ice collecting at the ice bottom |
---|
72 | Cfrazb = 5.0 ! squeezing coefficient for frazil ice collecting at the ice bottom |
---|
73 | hiclim = 0.10 ! minimum thickness of the first thickness category, must be smaller than hiccrit |
---|
74 | betas = 0.6 ! exponent in lead-ice fractionation of snow precipitation 0.66 |
---|
75 | ! betas = 1 -> equipartition, betas < 1 -> more on leads |
---|
76 | kappa_i = 1.0 ! extinction radiation parameter in sea ice (m-1) |
---|
77 | nconv_i_thd = 50 ! maximal number of iterations for heat diffusion computation |
---|
78 | maxer_i_thd = 0.0001 ! temperature error tolerance after heat diffusion |
---|
79 | thcon_i_swi = 1 ! sea ice thermal conductivity, Untersteiner_64 (0) or Pringle_et_al_07 (1) |
---|
80 | nn_monocat = 0 ! virtual ITD mono-category parameterizations (1) or not (0); 1 ok for jpl=1 only |
---|
81 | ! 2=replace ridging by piling, 3=activate G(he) only, 4=activate lateral melting only) --- those are temporary test options |
---|
82 | / |
---|
83 | !----------------------------------------------------------------------- |
---|
84 | &namicesal ! ice salinity |
---|
85 | !----------------------------------------------------------------------- |
---|
86 | num_sal = 2 ! S=bulk_sal (1), S(z,t) - (2), S(z) Schwarzacher_59 (3) |
---|
87 | bulk_sal = 4.0 ! if 1 is used, it represents the ice salinity |
---|
88 | sal_G = 5.00 ! restoring salinity, gravity drainage (g/kg) |
---|
89 | time_G = 1.728e+6 ! restoring time scale, gravity drainage (s) |
---|
90 | sal_F = 2.00 ! restoring salinity, flushing (g/kg) |
---|
91 | time_F = 8.640e+5 ! restoring time, flushing (s) |
---|
92 | s_i_max = 20.0 ! maximum tolerated ice salinity (g/kg) |
---|
93 | s_i_min = 0.1 ! minimum tolerated ice salinity (g/kg) |
---|
94 | / |
---|
95 | !----------------------------------------------------------------------- |
---|
96 | &namiceitdme ! parameters for mechanical redistribution of ice |
---|
97 | !----------------------------------------------------------------------- |
---|
98 | ridge_scheme_swi = 0 ! ice strength parameteriztaion, Hibler_79 (0), Rothrock_75 (1) |
---|
99 | Cs = 0.50 ! ratio of shearing energy contributing to ridging |
---|
100 | Cf = 17.0 ! ratio of ridging work to potential energy change in ridging |
---|
101 | fsnowrdg = 0.5 ! snow volume fraction that survives in ridging |
---|
102 | fsnowrft = 0.5 ! snow volume fraction that survives in rafting |
---|
103 | Gstar = 0.15 ! fractional area of thin ice being ridged (partfun_swi = 0) |
---|
104 | astar = 0.05 ! measures fractional area of thin ice being ridged (partfun_swi = 1) |
---|
105 | Hstar = 100.0 ! determines the maximum thickness of ridged ice (m) |
---|
106 | raft_swi = 1 ! rafting activated (1) or not (0) |
---|
107 | hparmeter = 0.75 ! threshold thickness for rafting (m) |
---|
108 | Craft = 5.0 ! coefficient used in the rafting function |
---|
109 | ridge_por = 0.3 ! porosity of newly ridged ice |
---|
110 | partfun_swi = 1 ! participation function: linear, Thorndike et al 75 (0) or exponential, Lipscomb 07 (1) |
---|
111 | brinstren_swi = 0 ! ice strength is a function of brine volume fraction (1) or not (0) |
---|
112 | / |
---|