source: CMIP6/DAMIP/CM61-LR-hist-GHG-03.1950/Debug/CM61-LR-hist-GHG-03.1950_19830101_19831231_namelist_ice @ 4385

Last change on this file since 4385 was 4385, checked in by ggalod, 5 years ago

sauvegarde exp DAMIP
File size: 13.7 KB
Line 
1namelist_ice_cfg namelist_ice_ref
2
3_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
4|  0   namelist_ice_cfg
5- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
60!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
70!! NEMO/LIM-3 : Ice configuration namelist. Overwrites SHARED/namelist_ice_lim3_ref
80!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
90
100!-----------------------------------------------------------------------
110&namicerun     !   Share parameters for dynamics/advection/thermo
120!-----------------------------------------------------------------------
130   cn_icerst_in  = "restart_ice_in" !  suffix of ice restart name (input)
140   rn_amax_n     = 0.997            !  maximum tolerated ice concentration NH
150   rn_amax_s     = 0.950            !  maximum tolerated ice concentration SH
160   ln_limdiahsb  = .false.          !  check the heat and salt budgets (T) or not (F)
170   ln_limdiaout  = .true.            !  output the heat and salt budgets (T) or not (F)
180/
190!-----------------------------------------------------------------------
200&namiceini     !   ice initialisation
210!-----------------------------------------------------------------------
220/
230!-----------------------------------------------------------------------
240&namicedyn     !   ice dynamic
250!-----------------------------------------------------------------------
260/
270!------------------------------------------------------------------------------
280&namicehdf     !   Ice horizontal diffusion
290!------------------------------------------------------------------------------
300/
310!-----------------------------------------------------------------------
320&namicethd     !   ice thermodynamic
330!-----------------------------------------------------------------------
340   rn_cdsn     = 0.50              !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971)
350                                   !  Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)
360/
370!-----------------------------------------------------------------------
380&namicesal     !   ice salinity
390!-----------------------------------------------------------------------
400/
410!-----------------------------------------------------------------------
420&namiceitdme   !   parameters for mechanical redistribution of ice
430!-----------------------------------------------------------------------
440   rn_astar    =   0.03            !  exponential measure of ridging ice fraction (nn_partfun = 1)
450   rn_hstar    = 25.0             !  determines the maximum thickness of ridged ice (m) (Hibler, 1980)
460/
470!-----------------------------------------------------------------------
480&namicedia     !   ice diagnostics
490!-----------------------------------------------------------------------
500/
510!------------------------------------------------------------------------------
520&namiceitd     !   Ice discretization
530!------------------------------------------------------------------------------
540/
55
56_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _
57|  1   namelist_ice_ref
58- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
591!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
601!! LIM3 namelist : 
611!!              1 - Generic parameters                 (namicerun)
621!!              2 - Ice initialization                 (namiceini)
631!!              3 - Ice discretization                 (namiceitd)
641!!              4 - Ice dynamics and transport         (namicedyn)
651!!              5 - Ice thermodynamics                 (namicethd)
661!!              6 - Ice salinity                       (namicesal)
671!!              7 - Ice mechanical redistribution      (namiceitdme)
681!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
691!
701!------------------------------------------------------------------------------
711&namicerun     !   Generic parameters
721!------------------------------------------------------------------------------
731   jpl            =    5           !  number of ice  categories
741   nlay_i         =    2           !  number of ice  layers
751   nlay_s         =    1           !  number of snow layers (only 1 is working)
761   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input)
771   cn_icerst_indir = "."           !  directory from which to read input ice restarts
781   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output)
791   cn_icerst_outdir = "."          !  directory in which to write output ice restarts
801   ln_limdyn     = .true.          !  ice dynamics (T) or thermodynamics only (F)
811   rn_amax_n     = 0.999           !  maximum tolerated ice concentration NH
821   rn_amax_s     = 0.999           !  maximum tolerated ice concentration SH
831   ln_limdiahsb  = .false.         !  check the heat and salt budgets (T) or not (F)
841   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F)
851   ln_icectl     = .false.         !  ice points output for debug (T or F)
861   iiceprt       = 10              !  i-index for debug
871   jiceprt       = 10              !  j-index for debug
881/
891!------------------------------------------------------------------------------
901&namiceini     !   Ice initialization
911!------------------------------------------------------------------------------
921   ln_iceini      = .true.         !  activate ice initialization (T) or not (F)
931   rn_thres_sst   =  2.0           !  maximum water temperature with initial ice (degC)
941   rn_hts_ini_n   =  0.3           !  initial real snow thickness (m), North
951   rn_hts_ini_s   =  0.3           !        "            "             South
961   rn_hti_ini_n   =  3.0           !  initial real ice thickness  (m), North
971   rn_hti_ini_s   =  1.0           !        "            "             South
981   rn_ati_ini_n   =  0.9           !  initial ice concentration   (-), North
991   rn_ati_ini_s   =  0.9           !        "            "             South
1001   rn_smi_ini_n   =  6.3           !  initial ice salinity     (g/kg), North
1011   rn_smi_ini_s   =  6.3           !        "            "             South
1021   rn_tmi_ini_n   =  270.          !  initial ice/snw temperature (K), North
1031   rn_tmi_ini_s   =  270.          !        "            "             South
1041/
1051!------------------------------------------------------------------------------
1061&namiceitd     !   Ice discretization
1071!------------------------------------------------------------------------------
1081   nn_catbnd      =    2           !  computation of ice category boundaries based on
1091                                   !      1: tanh function
1101                                   !      2: h^(-alpha), function of rn_himean
1111   rn_himean      =    2.0         !  expected domain-average ice thickness (m), nn_catbnd = 2 only
1121/
1131!------------------------------------------------------------------------------
1141&namicedyn     !   Ice dynamics and transport
1151!------------------------------------------------------------------------------
1161   nn_icestr      =    0           !  ice strength parameteriztaion                     
1171                                   !     0: Hibler_79     P = pstar*<h>*exp(-c_rhg*A)
1181                                   !     1: Rothrock_75   P = Cf*coeff*integral(wr.h^2)   
1191   ln_icestr_bvf  =    .false.     !  ice strength function brine volume (T) or not (F)     
1201   rn_pe_rdg      =   17.0         !  ridging work divided by pot. energy change in ridging, if nn_icestr = 1
1211   rn_pstar       =    2.0e+04     !  ice strength thickness parameter (N/m2), nn_icestr = 0
1221   rn_crhg        =   20.0         !  ice strength conc. parameter (-), nn_icestr = 0       
1231   rn_cio         =    5.0e-03     !  ice-ocean drag coefficient           (-)             
1241   rn_creepl      =    1.0e-12     !  creep limit (s-1)                                   
1251   rn_ecc         =    2.0         !  eccentricity of the elliptical yield curve         
1261   nn_nevp        =  120           !  number of EVP subcycles                             
1271   rn_relast      =    0.333       !  ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast
1281                                   !     advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300)
1291/
1301!------------------------------------------------------------------------------
1311&namicehdf     !   Ice horizontal diffusion
1321!------------------------------------------------------------------------------
1331   nn_ahi0        =    -1          !  horizontal diffusivity computation
1341                                   !    -1: no diffusion (bypass limhdf)
1351                                   !     0: use rn_ahi0_ref
1361                                   !     1: use rn_ahi0_ref x mean grid cell length / ( 2deg mean grid cell length )
1371                                   !     2: use rn_ahi0_ref x grid cell length      / ( 2deg mean grid cell length )
1381   rn_ahi0_ref    = 350.0          !  horizontal sea ice diffusivity (m2/s)
1391                                   !     if nn_ahi0 > 0, rn_ahi0_ref is the reference value at a nominal 2 deg resolution
1401   nn_convfrq     = 5              !  convergence check frequency of the Crant-Nicholson scheme (perf. optimization)
1411/
1421!------------------------------------------------------------------------------
1431&namicethd     !   Ice thermodynamics
1441!------------------------------------------------------------------------------
1451   rn_hnewice  = 0.1               !  thickness for new ice formation in open water (m)
1461   ln_frazil   = .false.           !  use frazil ice collection thickness as a function of wind (T) or not (F)
1471   rn_maxfrazb = 1.0               !  maximum fraction of frazil ice collecting at the ice base
1481   rn_vfrazb   = 0.417             !  thresold drift speed for frazil ice collecting at the ice bottom (m/s)
1491   rn_Cfrazb   = 5.0               !  squeezing coefficient for frazil ice collecting at the ice bottom
1501   rn_himin    = 0.10              !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice
1511   rn_betas    = 0.66              !  exponent in lead-ice repratition of snow precipitation
1521                                   !     betas = 1 -> equipartition, betas < 1 -> more on leads
1531   rn_kappa_i  = 1.0               !  radiation attenuation coefficient in sea ice (m-1)
1541   nn_conv_dif = 50                !  maximal number of iterations for heat diffusion computation
1551   rn_terr_dif = 0.0001            !  maximum temperature after heat diffusion (degC)
1561   nn_ice_thcon= 1                 !  sea ice thermal conductivity
1571                                   !     0: k = k0 + beta.S/T (Untersteiner, 1964)
1581                                   !     1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007)
1591   rn_cdsn     = 0.31              !  thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971)
1601                                   !  Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)
1611   nn_monocat  = 0                 !  virtual ITD mono-category parameterizations (1, jpl = 1 only) or not (0)
1621                                   !     2: simple piling instead of ridging --- temporary option
1631                                   !     3: activate G(he) only              --- temporary option
1641                                   !     4: activate lateral melting only    --- temporary option
1651  ln_it_qnsice = .true.            !  iterate the surface non-solar flux with surface temperature (T) or not (F)
1661/
1671!------------------------------------------------------------------------------
1681&namicesal     !   Ice salinity
1691!------------------------------------------------------------------------------
1701   nn_icesal   =  2                !  ice salinity option
1711                                   !     1: constant ice salinity (S=rn_icesal)
1721                                   !     2: varying salinity parameterization S(z,t)
1731                                   !     3: prescribed salinity profile S(z), Schwarzacher, 1959
1741   rn_icesal   =  4.               !  ice salinity (g/kg, nn_icesal = 1 only)
1751   rn_sal_gd   =  5.               !  restoring ice salinity, gravity drainage (g/kg)
1761   rn_time_gd  =  1.73e+6          !  restoring time scale, gravity drainage  (s)
1771   rn_sal_fl   =  2.               !  restoring ice salinity, flushing (g/kg)
1781   rn_time_fl  =  8.64e+5          !  restoring time scale, flushing (s)
1791   rn_simax    = 20.               !  maximum tolerated ice salinity (g/kg)
1801   rn_simin    =  0.1              !  minimum tolerated ice salinity (g/kg)
1811/
1821!------------------------------------------------------------------------------
1831&namiceitdme   !   Ice mechanical redistribution (ridging and rafting)
1841!------------------------------------------------------------------------------
1851   rn_Cs       =   0.5             !  fraction of shearing energy contributing to ridging
1861   rn_fsnowrdg =   0.5             !  snow volume fraction that survives in ridging
1871   rn_fsnowrft =   0.5             !  snow volume fraction that survives in rafting
1881   nn_partfun  =   1               !  type of ridging participation function
1891                                   !     0: linear (Thorndike et al, 1975)
1901                                   !     1: exponential (Lipscomb, 2007
1911   rn_gstar    =   0.15            !  fractional area of thin ice being ridged (nn_partfun = 0)
1921   rn_astar    =   0.05            !  exponential measure of ridging ice fraction (nn_partfun = 1)
1931   rn_hstar    = 100.0             !  determines the maximum thickness of ridged ice (m) (Hibler, 1980)
1941   ln_rafting  =   .true.          !  rafting activated (T) or not (F)
1951   rn_hraft    =   0.75            !  threshold thickness for rafting (m)
1961   rn_craft    =   5.0             !  squeezing coefficient used in the rafting function
1971   rn_por_rdg  =   0.3             !  porosity of newly ridged ice (Lepparanta et al., 1995)
1981/
199
Note: See TracBrowser for help on using the repository browser.