source: branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref @ 5060

!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/LIM3 :  1 - dynamics/advection/thermo          (namicerun)
!! namelists    2 - ice intialisation                  (namiceini)
!!              3 - ice dynamics                       (namicedyn)
!!              4 - ice advection                      (namicetrp)
!!              5 - thermodynamics                     (namicethd)
!!              6 - ice salinity                       (namicesal)
!!              7 - mechanical redistribution of ice   (namiceitdme)
!!              8 - ice diagnostics                    (namicedia)
!!              9 - ice outputs                        (namiceout)
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

!-----------------------------------------------------------------------
&namicerun     !   Share parameters for dynamics/advection/thermo
!-----------------------------------------------------------------------
   jpl           = 5               !  number of ice  categories
   nlay_i        = 2               !  number of ice  layers
   nlay_s        = 1               !  number of snow layers
   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input)
   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output)
   ln_limdyn     = .true.          !  ice dynamics (T) or thermodynamics only (F)
   amax          = 0.999           !  maximum tolerated ice concentration
   ln_nicep      = .false.         !  ice points output for debug (yes or no)
   ln_limdiahsb  = .false.         !  check the heat and salt budgets (T) or not (F)
   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F)
/
!-----------------------------------------------------------------------
&namiceitd     !   ice thickness distribution
!-----------------------------------------------------------------------
   nn_itdshp    = 1        !  shape of ice categories distribution following:  
                           !      1 = tanh function
                           !      2 = h^(-alpha)
   rn_itmean    = 2.5      !  mean thickness of the domain (used to compute the distribution, nn_itdshp = 2 only)
/
!-----------------------------------------------------------------------
&namiceini     !   ice initialization
!-----------------------------------------------------------------------
   ln_limini   = .true.   !  activate ice initialization (T) or not (F)
   thres_sst   =  2.0      !  threshold water temperature for initial sea ice
   hts_ini_n   =  0.3      !  initial snow thickness (North)
   hts_ini_s   =  0.3      !        "            "   (South)
   hti_ini_n   =  3.0      !  initial ice thickness  (North)
   hti_ini_s   =  1.0      !        "            "   (South)
   ati_ini_n   =  0.9      !  initial ice concentration (North)
   ati_ini_s   =  0.9      !        "            "      (South)
   smi_ini_n   =  6.3      !  initial ice salinity   (North)
   smi_ini_s   =  6.3      !        "            "   (South)
   tmi_ini_n   =  270.     !  initial ice/snw temp   (North)
   tmi_ini_s   =  270.     !  initial ice/snw temp   (South)
/
!-----------------------------------------------------------------------
&namicedyn     !   ice dynamics
!-----------------------------------------------------------------------
   cw          =   5.0e-03 !  ice-ocean drag coefficient
   pstar       =   2.0e+04 !  1st bulk-rheology parameter (N/m)
   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter
   creepl      =   1.0e-12 !  creep limit (s-1)
   ecc         =   2.0     !  eccentricity of the elliptical yield curve
   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice (m2/s) - DO NOT CHANGE IT
                           !  the value is set up for orca2 and is adapted automatically in the code 
                           !  depending on the resolution of your experiment
   nevp        = 120       !  number of EVP subcycles
   relast      =   0.333   !  ratio of elastic timescale to ice time step (1/3 if nevp=120, 1/9 if nevp=300)
/
!-----------------------------------------------------------------------
&namicethd     !   ice thermodynamics
!-----------------------------------------------------------------------
   hiccrit     = 0.1       !  thickness for new ice formation in open water
   fraz_swi    = 0         !  use frazil ice collection thickness as a function of wind (1.0) or not (0.0)
   maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom
   vfrazb      = 0.417     !  thresold drift speed for frazil ice collecting at the ice bottom
   Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom
   hiclim      = 0.10      !  minimum thickness of the first thickness category, must be smaller than hiccrit
   betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66
                           !        betas = 1 -> equipartition, betas < 1 -> more on leads
   kappa_i     = 1.0       !  extinction radiation parameter in sea ice (m-1)
   nconv_i_thd = 50        !  maximal number of iterations for heat diffusion computation
   maxer_i_thd = 0.0001    !  temperature error tolerance after heat diffusion 
   thcon_i_swi = 1         !  sea ice thermal conductivity, Untersteiner_64 (0) or Pringle_et_al_07 (1)
   nn_monocat  = 0         !  virtual ITD mono-category parameterizations (1) or not (0); 1 ok for jpl=1 only
                           !  2=replace ridging by piling, 3=activate G(he) only, 4=activate lateral melting only) --- those are temporary test options
/
!-----------------------------------------------------------------------
&namicesal     !   ice salinity
!-----------------------------------------------------------------------
   num_sal     =  2        !  S=bulk_sal (1), S(z,t) - (2), S(z) Schwarzacher_59 (3) 
   bulk_sal    =  4.0      !  if 1 is used, it represents the ice salinity
   sal_G       =  5.00     !  restoring salinity, gravity drainage (g/kg)
   time_G      =  1.728e+6 !  restoring time scale, gravity drainage  (s)
   sal_F       =  2.00     !  restoring salinity, flushing (g/kg)
   time_F      =  8.640e+5 !  restoring time, flushing (s)
   s_i_max     = 20.0      !  maximum tolerated ice salinity (g/kg)
   s_i_min     =  0.1      !  minimum tolerated ice salinity (g/kg)
/
!-----------------------------------------------------------------------
&namiceitdme   !   parameters for mechanical redistribution of ice 
!-----------------------------------------------------------------------
   ridge_scheme_swi =   0    !  ice strength parameteriztaion, Hibler_79 (0), Rothrock_75 (1)
   Cs               =   0.50 !  ratio of shearing energy contributing to ridging
   Cf               =  17.0  !  ratio of ridging work to potential energy change in ridging
   fsnowrdg         =   0.5  !  snow volume fraction that survives in ridging
   fsnowrft         =   0.5  !  snow volume fraction that survives in rafting
   Gstar            =   0.15 !  fractional area of thin ice being ridged (partfun_swi = 0)
   astar            =   0.05 !  measures fractional area of thin ice being ridged (partfun_swi = 1)
   Hstar            = 100.0  !  determines the maximum thickness of ridged ice (m)
   raft_swi         =   1    !  rafting activated (1) or not (0)
   hparmeter        =   0.75 !  threshold thickness for rafting (m)
   Craft            =   5.0  !  coefficient used in the rafting function
   ridge_por        =   0.3  !  porosity of newly ridged ice
   partfun_swi      =   1    !  participation function: linear, Thorndike et al 75 (0) or exponential, Lipscomb 07 (1)
   brinstren_swi    =   0    !  ice strength is a function of brine volume fraction (1) or not (0)
/