Changeset 5068


Ignore:
Timestamp:
2015-02-06T19:13:17+01:00 (6 years ago)
Author:
clem
Message:

LIM3 change all namelist names to fit with NEMO convention

File:
1 edited

Legend:

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

    r5060 r5068  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    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) 
     2!! LIM3 namelist :   
     3!!              1 - Generic parameters                 (namicerun) 
     4!!              2 - Ice initialization                 (namiceini) 
     5!!              3 - Ice discretization                 (namiceitd) 
     6!!              4 - Ice dynamics and transport         (namicedyn) 
     7!!              5 - Ice thermodynamics                 (namicethd) 
     8!!              6 - Ice salinity                       (namicesal) 
     9!!              7 - Ice mechanical redistribution      (namiceitdme) 
    1110!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    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 
     11! 
     12!------------------------------------------------------------------------------ 
     13&namicerun     !   Generic parameters 
     14!------------------------------------------------------------------------------ 
    1915   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input) 
    2016   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output) 
    2117   ln_limdyn     = .true.          !  ice dynamics (T) or thermodynamics only (F) 
    22    amax          = 0.999           !  maximum tolerated ice concentration 
     18   rn_amax       = 0.999           !  maximum tolerated ice concentration  
    2319   ln_nicep      = .false.         !  ice points output for debug (yes or no) 
    2420   ln_limdiahsb  = .false.         !  check the heat and salt budgets (T) or not (F) 
    2521   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F) 
    2622/ 
    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) 
     23!------------------------------------------------------------------------------ 
     24&namiceini     !   Ice initialization 
     25!------------------------------------------------------------------------------ 
     26   ln_iceini      = .true.         !  activate ice initialization (T) or not (F) 
     27   rn_thres_sst   =  2.0           !  maximum water temperature with initial ice (degC) 
     28   rn_hts_ini_n   =  0.3           !  initial real snow thickness (m), North 
     29   rn_hts_ini_s   =  0.3           !        "            "             South 
     30   rn_hti_ini_n   =  3.0           !  initial real ice thickness  (m), North 
     31   rn_hti_ini_s   =  1.0           !        "            "             South 
     32   rn_ati_ini_n   =  0.9           !  initial ice concentration   (-), North 
     33   rn_ati_ini_s   =  0.9           !        "            "             South 
     34   rn_smi_ini_n   =  6.3           !  initial ice salinity     (g/kg), North 
     35   rn_smi_ini_s   =  6.3           !        "            "             South 
     36   rn_tmi_ini_n   =  270.          !  initial ice/snw temperature (K), North 
     37   rn_tmi_ini_s   =  270.          !        "            "              South 
    3438/ 
    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) 
     39!------------------------------------------------------------------------------ 
     40&namiceitd     !   Ice discretization 
     41!------------------------------------------------------------------------------ 
     42   jpl            =    5           !  number of ice  categories 
     43   nlay_i         =    2           !  number of ice  layers 
     44   nlay_s         =    1           !  number of snow layers (only 1 is working) 
     45   nn_catbnd      =    1           !  computation of ice category boundaries based on 
     46                                   !      1: tanh function 
     47                                   !      2: h^(-alpha), function of rn_himean 
     48   rn_himean      =    2.5         !  expected domain-average ice thickness (m), nn_catbnd = 2 only 
    5049/ 
    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) 
     50!------------------------------------------------------------------------------ 
     51&namicedyn     !   Ice dynamics and transport 
     52!------------------------------------------------------------------------------ 
     53   nn_icestr      =    0           !  ice strength parameteriztaion                          !MV, replaces ridge_scheme_swi 
     54                                   !     0: Hibler_79     P = pstar*<h>*exp(-c_rhg*A) 
     55                                   !     1: Rothrock_75   P = Cf*coeff*integral(wr.h^2)     
     56   nn_icestr_bvf  =    0           !  ice strength function brine volume (1) or not (0)      !MV replaces brinstren_swi 
     57   rn_pstar       =    2.0e+04     !  ice strength thickness parameter (N/m2), nn_icestr = 0  
     58   rn_crhg        =   20.0         !  ice strength conc. parameter (-), nn_icestr = 0        
     59   rn_cio         =    5.0e-03     !  ice-ocean drag coefficient           (-)              
     60   rn_creepl      =    1.0e-12     !  creep limit (s-1)                                    
     61   rn_ecc         =    2.0         !  eccentricity of the elliptical yield curve           
     62   nn_nevp        =  120           !  number of EVP subcycles                              
     63   rn_relast      =    0.333       !  ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
     64                                   !     advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     65!   nn_ahi0        =    2           !  horizontal diffusivity computation ! MV only 2 is coded now 
     66!                                   !     0: use rn_ahi0_ref 
     67!                                   !     1: use rn_ahi0_ref . mean grid cell length / ( 2deg mean grid cell length ) 
     68!                                   !     2: use rn_ahi0_ref . grid cell length      / ( 2deg mean grid cell length ) 
     69   rn_ahi0_ref    = 350.0          !  horizontal sea ice diffusivity (m2/s)  
     70                                   !     if nn_ahi0 > 0, rn_ahi0_ref is the reference value at a nominal 2 deg resolution 
    6471/ 
    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 
     72!------------------------------------------------------------------------------ 
     73&namicethd     !   Ice thermodynamics 
     74!------------------------------------------------------------------------------ 
     75   rn_hnewice  = 0.1               !  thickness for new ice formation in open water (m) 
     76   nn_frazil   = 0                 !  use frazil ice collection thickness as a function of wind (1.0) or not (0.0) 
     77   rn_maxfrazb = 0.0               !  maximum fraction of frazil ice collecting at the ice base 
     78   rn_vfrazb   = 0.417             !  thresold drift speed for frazil ice collecting at the ice bottom (m/s) 
     79   rn_Cfrazb   = 5.0               !  squeezing coefficient for frazil ice collecting at the ice bottom 
     80   rn_himin    = 0.10              !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 
     81   rn_betas    = 0.66              !  exponent in lead-ice repratition of snow precipitation 
     82                                   !     betas = 1 -> equipartition, betas < 1 -> more on leads 
     83   rn_kappa_i  = 1.0               !  radiation attenuation coefficient in sea ice (m-1) 
     84   nn_conv_dif = 50                !  maximal number of iterations for heat diffusion computation 
     85   rn_terr_dif = 0.0001            !  maximum temperature after heat diffusion (degC) 
     86   nn_ice_thcon= 1                 !  sea ice thermal conductivity 
     87                                   !     0: k = k0 + beta.S/T (Untersteiner, 1964) 
     88                                   !     1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 
     89   nn_monocat  = 0                 !  virtual ITD mono-category parameterizations (1, jpl = 1 only) or not (0) 
     90                                   !     2: simple piling instead of ridging --- temporary option 
     91                                   !     3: activate G(he) only              --- temporary option 
     92                                   !     4: activate lateral melting only    --- temporary option 
    8293/ 
    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&namicesal     !   Ice salinity 
     96!------------------------------------------------------------------------------ 
     97   nn_icesal   =  2                !  ice salinity option 
     98                                   !     1: constant ice salinity (S=rn_icesal) 
     99                                   !     2: varying salinity parameterization S(z,t) 
     100                                   !     3: prescribed salinity profile S(z), Schwarzacher, 1959 
     101   rn_icesal   =  4.               !  ice salinity (g/kg, nn_icesal = 1 only) 
     102   rn_sal_gd   =  5.               !  restoring ice salinity, gravity drainage (g/kg) 
     103   rn_time_gd  =  1.73e+6          !  restoring time scale, gravity drainage  (s) 
     104   rn_sal_fl   =  2.               !  restoring ice salinity, flushing (g/kg) 
     105   rn_time_fl  =  8.64e+5          !  restoring time scale, flushing (s) 
     106   rn_simax    = 20.               !  maximum tolerated ice salinity (g/kg) 
     107   rn_simin    =  0.1              !  minimum tolerated ice salinity (g/kg) 
    94108/ 
    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) 
     109!------------------------------------------------------------------------------ 
     110&namiceitdme   !   Ice mechanical redistribution (ridging and rafting) 
     111!------------------------------------------------------------------------------ 
     112   rn_pe_rdg   =   17.0            !  ridging work divided by pot. energy change in ridging, nn_ice str = 1 ! MV replaces Cf 
     113   rn_Cs       =   0.5             !  fraction of shearing energy contributing to ridging 
     114   rn_fsnowrdg =   0.5             !  snow volume fraction that survives in ridging 
     115   rn_fsnowrft =   0.5             !  snow volume fraction that survives in rafting 
     116   nn_partfun  =   1               !  type of ridging participation function 
     117                                   !     0: linear (Thorndike et al, 1975) 
     118                                   !     1: exponential (Lipscomb, 2007 
     119   rn_gstar    =   0.15            !  fractional area of thin ice being ridged (nn_partfun = 0) 
     120   rn_astar    =   0.05            !  exponential measure of ridging ice fraction (nn_partfun = 1) 
     121   rn_hstar    = 100.0             !  determines the maximum thickness of ridged ice (m) (Hibler, 1980) 
     122   nn_rafting  =   1               !  rafting activated (1) or not (0) 
     123   rn_hraft    =   0.75            !  threshold thickness for rafting (m) !MV replaces h_parmeter 
     124   rn_craft    =   5.0             !  squeezing coefficient used in the rafting function 
     125   rn_por_rdg  =   0.3             !  porosity of newly ridged ice (Lepparanta et al., 1995) 
    112126/ 
Note: See TracChangeset for help on using the changeset viewer.