Ignore:
Timestamp:
2015-02-02T11:28:50+01:00 (6 years ago)
Author:
vancop
Message:

new itd boundaries, namelist changes, mono-category and comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r5046 r5048  
    22!! NEMO/LIM3 :  1 - dynamics/advection/thermo          (namicerun) 
    33!! namelists    2 - ice intialisation                  (namiceini) 
    4 !!              3 - ice dynamic                        (namicedyn) 
     4!!              3 - ice dynamics                       (namicedyn) 
    55!!              4 - ice advection                      (namicetrp) 
    6 !!              5 - thermodynamic                      (namicethd) 
     6!!              5 - thermodynamics                     (namicethd) 
    77!!              6 - ice salinity                       (namicesal) 
    88!!              7 - mechanical redistribution of ice   (namiceitdme) 
     
    1717   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output) 
    1818   ln_limdyn     = .true.          !  ice dynamics (T) or thermodynamics only (F) 
    19    amax          = 0.999           !  maximum ice concentration 
    20    ln_nicep      = .false.         !  Ice points output for debug (yes or no) 
    21    ln_limdiahsb  = .false.          !  check the heat and salt budgets (T) or not (F) 
     19   amax          = 0.999           !  maximum tolerated ice concentration 
     20   ln_nicep      = .false.         !  ice points output for debug (yes or no) 
     21   ln_limdiahsb  = .false.         !  check the heat and salt budgets (T) or not (F) 
    2222   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F) 
    2323/ 
    2424!----------------------------------------------------------------------- 
    25 &namiceini     !   ice initialisation 
     25&namiceini     !   ice initialization 
    2626!----------------------------------------------------------------------- 
    2727   ln_limini   = .false.   !  activate ice initialization (T) or not (F) 
    2828   thres_sst   =  0.0      !  threshold water temperature for initial sea ice 
    29    hts_ini_n   =  0.3      !  initial snow thickness in the north 
    30    hts_ini_s   =  0.3      !        "            "          south 
    31    hti_ini_n   =  3.0      !  initial ice thickness in the north 
    32    hti_ini_s   =  1.0      !        "            "         south 
    33    ati_ini_n   =  0.9      !  initial ice concentration in the north 
    34    ati_ini_s   =  0.9      !        "            "             south 
    35    smi_ini_n   =  6.3      !  initial ice salinity in the north 
    36    smi_ini_s   =  6.3      !        "            "    south 
    37    tmi_ini_n   =  270.     !  initial ice/snw temp in the north 
    38    tmi_ini_s   =  270.     !  initial ice/snw temp in the south 
     29   hts_ini_n   =  0.3      !  initial snow thickness (North) 
     30   hts_ini_s   =  0.3      !        "            "   (South) 
     31   hti_ini_n   =  3.0      !  initial ice thickness  (North) 
     32   hti_ini_s   =  1.0      !        "            "   (South) 
     33   ati_ini_n   =  0.9      !  initial ice concentration (North) 
     34   ati_ini_s   =  0.9      !        "            "      (South) 
     35   smi_ini_n   =  6.3      !  initial ice salinity   (North) 
     36   smi_ini_s   =  6.3      !        "            "   (South) 
     37   tmi_ini_n   =  270.     !  initial ice/snw temp   (North) 
     38   tmi_ini_s   =  270.     !  initial ice/snw temp   (South) 
    3939/ 
    4040!----------------------------------------------------------------------- 
    41 &namicedyn     !   ice dynamic 
     41&namicedyn     !   ice dynamics 
    4242!----------------------------------------------------------------------- 
    43    cw          =   5.0e-03 !  drag coefficient for oceanic stress 
    44    pstar       =   2.0e+04 !  1st bulk-rheology parameter 
     43   cw          =   5.0e-03 !  ice-ocean drag coefficient 
     44   pstar       =   2.0e+04 !  1st bulk-rheology parameter (N/m) 
    4545   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    46    creepl      =   1.0e-12 !  creep limit 
     46   creepl      =   1.0e-12 !  creep limit (s-1) 
    4747   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
    48    ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] ! depends on resolution 
    49    nevp        = 120       !  number of iterations for subcycling in EVP 
    50    relast      = 0.333     !  ratio of elastic timescale over ice time step (1/3 if nevp=120 ; 1/9 if nevp=300) 
    51    hminrhg     =   0.001   !  ice volume (a*h in m) below which ice velocity equal ocean velocity 
     48   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice (m2/s), depends on resolution 
     49   nevp        = 120       !  number of EVP subcycles 
     50   relast      =   0.333   !  ratio of elastic timescale to ice time step (1/3 if nevp=120, 1/9 if nevp=300) 
     51   hminrhg     =   0.001   !  ice volume (a*h in m) below which ice velocity equals ocean velocity 
    5252/ 
    5353!----------------------------------------------------------------------- 
    54 &namicethd     !   ice thermodynamic 
     54&namicethd     !   ice thermodynamics 
    5555!----------------------------------------------------------------------- 
    56    hiccrit     = 0.1       !  ice thickness for lateral accretion  
    57                            !         caution 1.0, 1.0 best value to be used!!! (gilles G.)  ???? 
    58    fraz_swi    = 0         !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
     56   hiccrit     = 0.1       !  thickness for new ice formation in open water 
     57   fraz_swi    = 0         !  use frazil ice collection thickness as a function of wind (1.0) or not (0.0) 
    5958   maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom 
    60    vfrazb      = 0.4166667 !  thresold drift speed for frazil ice collecting at the ice bottom 
     59   vfrazb      = 0.417    !  thresold drift speed for frazil ice collecting at the ice bottom 
    6160   Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom 
    6261   hiclim      = 0.10      !  minimum thickness of the first thickness category, must be smaller than hiccrit 
    6362   betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66 
    6463                           !        betas = 1 -> equipartition, betas < 1 -> more on leads 
    65    kappa_i     = 1.0       !  extinction radiation parameter in sea ice (1.0) 
     64   kappa_i     = 1.0       !  extinction radiation parameter in sea ice (m-1) 
    6665   nconv_i_thd = 50        !  maximal number of iterations for heat diffusion computation 
    67    maxer_i_thd = 0.0001    !  maximal error in temperature for heat diffusion computation 
    68    thcon_i_swi = 1         !  switch for computation of thermal conductivity in the ice 
    69                            !        (0) Untersteiner (1964), (1) Pringle et al. (2007) 
     66   maxer_i_thd = 0.0001    !  temperature error tolerance after heat diffusion  
     67   thcon_i_swi = 1         !  sea ice thermal conductivity, Untersteiner_64 (0) or Pringle_et_al_07 (1) 
     68   nn_monocat  = 0         !  virtual ITD mono-category parameterizations (1) or not (0); 1 ok for jpl=1 only 
     69                           !  2=replace ridging by piling, 3=activate G(he) only, 4=activate lateral melting only) --- those are temporary test options 
    7070/ 
    7171!----------------------------------------------------------------------- 
    7272&namicesal     !   ice salinity 
    7373!----------------------------------------------------------------------- 
    74    num_sal     =  2        !  salinity option: 1 -> S = bulk_sal 
    75                            !                   2 -> S = S(z,t) with a simple parameterization 
    76                            !                   3 -> S = S(z) profile of Scwharzacher [1959] 
    77                            !                   4 -> S = S(h) Cox and Weeks [1974] 
     74   num_sal     =  2        !  S=bulk_sal (1), S(z,t) - (2), S(z) Schwarzacher_59 (3)  
    7875   bulk_sal    =  4.0      !  if 1 is used, it represents the ice salinity 
    79    sal_G       =  5.00     !  restoring salinity for GD 
    80    time_G      =  1.728e+6 !  restoring time for GD 
    81    sal_F       =  2.00     !  restoring salinity for flushing 
    82    time_F      =  8.640e+5 !  restoring time for flushing 
    83    s_i_max     = 20.0      !  Maximum salinity  
    84    s_i_min     =  0.1      !  Minimum tolerated ice salinity 
     76   sal_G       =  5.00     !  restoring salinity, gravity drainage (g/kg) 
     77   time_G      =  1.728e+6 !  restoring time scale, gravity drainage  (s) 
     78   sal_F       =  2.00     !  restoring salinity, flushing (g/kg) 
     79   time_F      =  8.640e+5 !  restoring time, flushing (s) 
     80   s_i_max     = 20.0      !  maximum tolerated ice salinity (g/kg) 
     81   s_i_min     =  0.1      !  minimum tolerated ice salinity (g/kg) 
    8582/ 
    8683!----------------------------------------------------------------------- 
    8784&namiceitdme   !   parameters for mechanical redistribution of ice  
    8885!----------------------------------------------------------------------- 
    89    ridge_scheme_swi = 0      !  which ridging scheme using (1=Rothrock,else=Hibler79) 
    90    Cs               =   0.50 !  shearing energy contribution to ridging 
    91    Cf               =  17.0  !  ratio of ridging work to PE change in ridging 
    92    fsnowrdg         =   0.5  !  snow fraction that survives in ridging 
    93    fsnowrft         =   0.5  !  snow fraction that survives in rafting 
    94    Gstar            =   0.15 !  fractional area of thin ice being ridged 
    95    astar            =   0.05 !  equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 
    96    Hstar            = 100.0  !  parameter determining the maximum thickness of ridged ice 
    97    raft_swi         =   1    !  rafting or not 
    98    hparmeter        =   0.75 !  threshold thickness for rafting or not 
     86   ridge_scheme_swi =   0    !  ice strength parameteriztaion, Hibler_79 (0), Rothrock_75 (1) 
     87   Cs               =   0.50 !  ratio of shearing energy contributing to ridging 
     88   Cf               =  17.0  !  ratio of ridging work to potential energy change in ridging 
     89   fsnowrdg         =   0.5  !  snow volume fraction that survives in ridging 
     90   fsnowrft         =   0.5  !  snow volume fraction that survives in rafting 
     91   Gstar            =   0.15 !  fractional area of thin ice being ridged (partfun_swi = 0) 
     92   astar            =   0.05 !  measures fractional area of thin ice being ridged (partfun_swi = 1) 
     93   Hstar            = 100.0  !  determines the maximum thickness of ridged ice (m) 
     94   raft_swi         =   1    !  rafting activated (1) or not (0) 
     95   hparmeter        =   0.75 !  threshold thickness for rafting (m) 
    9996   Craft            =   5.0  !  coefficient used in the rafting function 
    100    ridge_por        =   0.3  !  initial porosity of the ridged ice (typically 0.30) 
    101    partfun_swi      =   1    !  participation function linear, TH75 (0) or exponential Letal07 (1) 
    102    brinstren_swi    =   0    !  (1) use brine volume to diminish ice strength 
     97   ridge_por        =   0.3  !  porosity of newly ridged ice 
     98   partfun_swi      =   1    !  participation function: linear, Thorndike et al 75 (0) or exponential, Lipscomb 07 (1) 
     99   brinstren_swi    =   0    !  ice strength is a function of brine volume fraction (1) or not (0) 
Note: See TracChangeset for help on using the changeset viewer.