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_lim3 in branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00 – NEMO

source: branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_ice_lim3 @ 4345

Last change on this file since 4345 was 4345, checked in by clem, 10 years ago

clean the glob_sum diags and make scalar.nc work

File size: 14.5 KB
Line 
1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2!! NEMO/LIM3 :  1 - dynamics/advection/thermo          (namicerun)
3!! namelists    2 - ice intialisation                  (namiceini)
4!!              3 - ice dynamic                        (namicedyn)
5!!              4 - ice advection                      (namicetrp)
6!!              5 - thermodynamic                      (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   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input)
17   cn_icerst_out = "restart_ice"      !  suffix of ice restart name (output)
18   ln_limdyn   = .true.    !  ice dynamics (T) or thermodynamics only (F)
19   amax        = 0.999      !  maximum ice concentration
20   cai         =  1.40e-3  !  atmospheric drag over sea ice
21   cao         =  1.00e-3  !  atmospheric drag over ocean
22   ln_nicep    = .false.   !  Ice points output for debug (yes or no)
23   ln_limdiahsb  = .false.    !  check the heat and salt budgets (T) or not (F)
24   ln_limdiaout  = .true.    !  output the heat and salt budgets (T) or not (F)
25/
26!-----------------------------------------------------------------------
27&namiceini     !   ice initialisation
28!-----------------------------------------------------------------------
29   ttest       =  2.0      !  threshold water temperature for initial sea ice
30   hninn       =  0.3      !  initial snow thickness in the north
31   hnins       =  0.1      !        "            "          south
32   hginn       =  3.50     !  initial undeformed ice thickness in the north
33   hgins       =  1.0      !        "            "              "     south
34   aginn       =  0.95     !  initial undeformed ice concentration in the north
35   agins       =  0.9      !        "            "              "         south
36   sinn        =  6.301    !  initial salinity in the north
37   sins        =  6.301    !        "            "    south
38/
39!-----------------------------------------------------------------------
40&namicedyn     !   ice dynamic
41!-----------------------------------------------------------------------
42   epsd        =   1.0e-20 !  tolerance parameter
43   alpha       =   0.5     !  coefficient for semi-implicit coriolis
44   dm          =   0.6e+03 !  diffusion constant for dynamics
45   nbiter      =   1       !  number of sub-time steps for relaxation
46   nbitdr      = 500       !  maximum number of iterations for relaxation
47   om          =   0.5     !  relaxation constant
48   resl        =   5.0e-05 !  maximum value for the residual of relaxation
49   cw          =   5.0e-03 !  drag coefficient for oceanic stress
50   angvg       =   0.0     !  turning angle for oceanic stress
51   pstar       =   2.0e+04 !  1st bulk-rheology parameter
52   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter
53   etamn       =   0.0e+07 !  minimun value for viscosity
54   creepl      =   1.0e-12 !  creep limit
55   ecc         =   2.0     !  eccentricity of the elliptical yield curve
56   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s]
57   nevp        = 300       !  number of iterations for subcycling in EVP, SB, 400
58   telast      =9600.0     !  timescale for elastic waves, SB, 720.0
59   alphaevp    =   1.0     !  coefficient for the solution of internal ice stresses
60   hminrhg     =   0.001   !  ice volume (a*h in m) below which ice velocity equal ocean velocity
61/
62!-----------------------------------------------------------------------
63&namicetrp     !   ice transport
64!-----------------------------------------------------------------------
65   bound       =   0.      !  boundary conditions (=0.0 no-slip, =1.0 free-slip)
66/
67!-----------------------------------------------------------------------
68&namicethd     !   ice thermodynamic
69!-----------------------------------------------------------------------
70   hmelt       = -0.15     !  maximum melting at the bottom
71   hiccrit     = 0.1 , 0.1 !  ice thickness for lateral accretion in the Northern (Southern) Hemisphere
72                           !         caution 1.0, 1.0 best value to be used!!! (gilles G.)  ????
73   fraz_swi    = 0.0       !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0)
74   maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom
75   vfrazb      = 0.4166667 !  thresold drift speed for frazil ice collecting at the ice bottom
76   Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom
77   hicmin      = 0.2       !  ice thickness corr. to max. energy stored in brine pocket
78   hiclim      = 0.10      !  minimum ice thickness
79   sbeta       = 1.        !  numerical caracteritic of the scheme for diffusion in ice
80                           !        Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)
81   parlat      = 0.0       !  percentage of energy used for lateral ablation
82   hakspl      = 0.25      !  slope of distr. for Hakkinen-Mellor's lateral melting
83   hibspl      = 0.5       !  slope of distribution for Hibler's lateral melting
84   exld        = 2.0       !  exponent for leads-closure rate
85   hakdif      = 1.0       !  coefficient for diffusions of ice and snow
86   thth        = 0.2       !  threshold thickness for comp. of eq. thermal conductivity
87   hnzst       = 0.1       !  thickness of the surf. layer in temp. computation
88   parsub      = 1.0       !  switch for snow sublimation or not
89   alphs       = 1.0       !  coefficient for snow density when snow ice formation
90   betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66
91                           !        betas = 1 -> equipartition, betas < 1 -> more on leads
92   kappa_i     = 1.0       !  extinction radiation parameter in sea ice (1.0)
93   nconv_i_thd = 50        !  maximal number of iterations for heat diffusion computation
94   maxer_i_thd = 0.0001    !  maximal error in temperature for heat diffusion computation
95   thcon_i_swi = 1         !  switch for computation of thermal conductivity in the ice
96                           !        (0) Untersteiner (1964), (1) Pringle et al. (2007)
97/
98!-----------------------------------------------------------------------
99&namicesal     !   ice salinity
100!-----------------------------------------------------------------------
101   num_sal     =  2        !  salinity option: 1 -> S = bulk_sal
102                           !                   2 -> S = S(z,t) with a simple parameterization
103                           !                   3 -> S = S(z) profile of Scwharzacher [1959]
104                           !                   4 -> S = S(h) Cox and Weeks [1974]
105   bulk_sal    =  4.0      !  if 1 is used, it represents the ice salinity
106   sal_G       =  5.00     !  restoring salinity for GD
107   time_G      =  1.728e+6 !  restoring time for GD
108   sal_F       =  2.00     !  restoring salinity for flushing
109   time_F      =  8.640e+5 !  restoring time for flushing
110   s_i_max     = 20.0      !  Maximum salinity
111   s_i_min     =  0.1      !  Minimum tolerated ice salinity
112   s_i_0       =  3.5      !  1st salinity for salinity profile
113   s_i_1       =  4.5      !  2nd salinity for salinity profile
114/
115!-----------------------------------------------------------------------
116&namiceitdme   !   parameters for mechanical redistribution of ice
117!-----------------------------------------------------------------------
118   ridge_scheme_swi = 0      !  which ridging scheme using (1=Rothrock,else=Hibler79)
119   Cs               =   0.50 !  shearing energy contribution to ridging
120   Cf               =  17.0  !  ratio of ridging work to PE change in ridging
121   fsnowrdg         =   0.5  !  snow fraction that survives in ridging
122   fsnowrft         =   0.5  !  snow fraction that survives in rafting
123   Gstar            =   0.15 !  fractional area of thin ice being ridged
124   astar            =   0.05 !  equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice)
125   Hstar            = 100.0  !  parameter determining the maximum thickness of ridged ice
126   raftswi          =   1    !  rafting or not
127   hparmeter        =   0.75 !  threshold thickness for rafting or not
128   Craft            =   5.0  !  coefficient used in the rafting function
129   ridge_por        =   0.3  !  initial porosity of the ridged ice (typically 0.30)
130   sal_max_ridge    =  15.0  !  maximum ridged ice salinity
131   partfun_swi      =   1    !  participation function linear, TH75 (0) or exponential Letal07 (1)
132   transfun_swi     =   0    !  transfer function uniform of H80 (0) or exponential Letal07 (1)
133   brinstren_swi    =   0    !  (1) use brine volume to diminish ice strength
134/
135!-----------------------------------------------------------------------
136&namicedia     !   ice diagnostics
137!-----------------------------------------------------------------------
138   fmtinf      ='1PE13.5 ' !  format of the output values
139   nfrinf      = 4         !  number of variables written in one line
140   ntmoy       = 1         !  instantaneous values of ice evolution or averaging
141   ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging
142/
143!!-----------------------------------------------------------------------
144!&namicehsb       !  Heat and salt budgets
145!!-----------------------------------------------------------------------
146!
147!/
148!-----------------------------------------------------------------------
149&namiceout     !   parameters for outputs
150!-----------------------------------------------------------------------
151   noumef      =   43      !  number of fields
152   add_diag_swi=    1      !  1 -> diagnose distribution in thickness space
153                           !  0 -> only simple diagnostics
154!
155!           !         title of the field           !  name     !   units   !  save  ! multipl. ! additive !
156!           !                                      !           !           ! or not !  factor  !  factor  !
157   field_1  = 'Ice concentration                  ', 'iiceconc', '%       ',    1   ,  1.0     ,    0.0
158   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0
159   field_3  = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0
160   field_4  = 'Daily bottom thermo ice production ', 'iicebopr', 'km3/day ',    1   , 1.0e-9   ,    0.0
161   field_5  = 'Daily dynamic ice production       ', 'iicedypr', 'km3/day ',    1   , 1.0e-9   ,    0.0
162   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0
163   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0
164   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0
165   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15
166   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0
167   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0
168   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0
169   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0
170   field_14 = 'Heat flux due to brine release     ', 'iocehebr', 'w/m2    ',    1   ,  1.0     ,    0.0
171   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0
172   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0
173   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0
174   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0
175   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/m2/d ',    1   ,  1.0     ,    0.0
176   field_20 = 'Mean ice salinity                  ', 'iicesali', 'psu     ',    1   ,  1.0     ,    0.0
177   field_21 = 'Mean ice age                       ', 'iiceages', 'years   ',    1   ,  0.002739,    0.0
178   field_22 = 'Daily lateral thermo ice prod.     ', 'iicelapr', 'km3/day ',    1   ,1.0e-9    ,    0.0
179   field_23 = 'Daily snowice ice production       ', 'iicesipr', 'km3/day ',    1   ,1.0e-9    ,    0.0
180   field_24 = 'Mean ice temperature               ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15
181   field_25 = 'Ice total heat content             ', 'iiceheco', '10^9 J  ',    1   ,  1.0     ,    0.0
182   field_26 = 'Ice surface temperature            ', 'iicesurt', 'C       ',    1   ,  1.0     , -273.15
183   field_27 = 'Snow temperature                   ', 'isnotem2', 'C       ',    1   ,  1.0     , -273.15
184   field_28 = 'Fsbri - brine salt flux            ', 'iicefsbr', 'kg/m2/d ',    1   ,  1.0     ,    0.0
185   field_29 = 'Fseqv - equivalent FW salt flux    ', 'iicefseq', 'kg/m2/d ',    1   ,  1.0     ,    0.0
186   field_30 = 'Brine volume                       ', 'ibrinvol', '%       ',    1   ,  100.0   ,    0.0
187   field_31 = 'Frazil ice collection thickness    ', 'iicecolf', 'm       ',    1   ,  1.0     ,    0.0
188   field_32 = 'Ice strength                       ', 'iicestre', 'N/m     ',    1   ,  0.001   ,    0.0
189   field_33 = 'Ice velocity                       ', 'iicevelo', 'm/s     ',    1   ,  1.0     ,    0.0
190   field_34 = 'Surface melt                       ', 'iicesume', 'km3/day ',    1   ,1.0e-9    ,    0.0
191   field_35 = 'Bottom melt                        ', 'iicebome', 'km3/day ',    1   ,1.0e-9    ,    0.0
192   field_36 = 'Divergence                         ', 'iicedive', '10-8s-1 ',    1   ,  1.0e8   ,    0.0
193   field_37 = 'Shear                              ', 'iiceshea', '10-8s-1 ',    1   ,  1.0e8   ,    0.0
194   field_38 = 'Daily resultant ice prod/melt      ', 'iicerepr', 'km3/day ',    1   ,  1.0e-9  ,    0.0
195   field_39 = 'Ice volume                         ', 'iicevolu', 'km3     ',    1   ,  1.0e-9  ,    0.0
196   field_40 = 'Snow volume                        ', 'isnowvol', 'km3     ',    1   ,  1.0e-9  ,    0.0
197   field_41 = 'Fsrpo - salt flux from ridg/raft   ', 'iicefsrp', 'kg/m2/d ',    1   ,  1.0     ,    0.0
198   field_42 = 'Fsres - salt flux from limupdate   ', 'iicefsre', 'kg/m2/d ',    1   ,  1.0     ,    0.0
199   field_43 = 'Ice volume transport               ', 'iicevtrp', 'km3/day ',    1   ,1.0e-9    ,    0.0
200/
201
Note: See TracBrowser for help on using the repository browser.