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.
Changeset 5709 for branches/UKMO/DEV_r5107_dynvor_updates/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref – NEMO

Ignore:
Timestamp:
2015-08-26T16:14:53+02:00 (9 years ago)
Author:
davestorkey
Message:

Updating UKMO/DEV_r5107_dynvor_updates branch to be relative to revision 5518 of the trunk
(= branching point for NEMO 3.6_STABLE).

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/DEV_r5107_dynvor_updates/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r4990 r5709  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    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) 
     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 !----------------------------------------------------------------------- 
     11! 
     12!------------------------------------------------------------------------------ 
     13&namicerun     !   Generic parameters 
     14!------------------------------------------------------------------------------ 
     15   jpl            =    5           !  number of ice  categories 
     16   nlay_i         =    2           !  number of ice  layers 
     17   nlay_s         =    1           !  number of snow layers (only 1 is working) 
    1618   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input) 
     19   cn_icerst_indir = "."           !  directory from which to read input ice restarts 
    1720   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output) 
     21   cn_icerst_outdir = "."          !  directory in which to write output ice restarts 
    1822   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 (clio) 
    21    cao           = 1.00e-3         !  atmospheric drag over ocean   (clio) 
    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) 
     23   rn_amax       = 0.999           !  maximum tolerated ice concentration  
     24   ln_limdiahsb  = .false.         !  check the heat and salt budgets (T) or not (F) 
    2425   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F) 
     26   ln_icectl     = .false.         !  ice points output for debug (T or F) 
     27   iiceprt       = 10              !  i-index for debug 
     28   jiceprt       = 10              !  j-index for debug 
    2529/ 
    26 !----------------------------------------------------------------------- 
    27 &namiceini     !   ice initialisation 
    28 !----------------------------------------------------------------------- 
    29    ln_limini   = .false.   !  activate ice initialization (T) or not (F) 
    30    thres_sst   =  0.0      !  threshold water temperature for initial sea ice 
    31    hts_ini_n   =  0.3      !  initial snow thickness in the north 
    32    hts_ini_s   =  0.3      !        "            "          south 
    33    hti_ini_n   =  1.0      !  initial ice thickness in the north 
    34    hti_ini_s   =  1.0      !        "            "         south 
    35    ati_ini_n   =  0.9      !  initial ice concentration in the north 
    36    ati_ini_s   =  0.9      !        "            "             south 
    37    smi_ini_n   =  6.301    !  initial ice salinity in the north 
    38    smi_ini_s   =  6.301    !        "            "    south 
    39    tmi_ini_n   =  270.     !  initial ice/snw temp in the north 
    40    tmi_ini_s   =  270.     !  initial ice/snw temp in the south 
     30!------------------------------------------------------------------------------ 
     31&namiceini     !   Ice initialization 
     32!------------------------------------------------------------------------------ 
     33   ln_iceini      = .true.         !  activate ice initialization (T) or not (F) 
     34   rn_thres_sst   =  2.0           !  maximum water temperature with initial ice (degC) 
     35   rn_hts_ini_n   =  0.3           !  initial real snow thickness (m), North 
     36   rn_hts_ini_s   =  0.3           !        "            "             South 
     37   rn_hti_ini_n   =  3.0           !  initial real ice thickness  (m), North 
     38   rn_hti_ini_s   =  1.0           !        "            "             South 
     39   rn_ati_ini_n   =  0.9           !  initial ice concentration   (-), North 
     40   rn_ati_ini_s   =  0.9           !        "            "             South 
     41   rn_smi_ini_n   =  6.3           !  initial ice salinity     (g/kg), North 
     42   rn_smi_ini_s   =  6.3           !        "            "             South 
     43   rn_tmi_ini_n   =  270.          !  initial ice/snw temperature (K), North 
     44   rn_tmi_ini_s   =  270.          !        "            "             South 
    4145/ 
    42 !----------------------------------------------------------------------- 
    43 &namicedyn     !   ice dynamic 
    44 !----------------------------------------------------------------------- 
    45    epsd        =   1.0e-20 !  tolerance parameter 
    46    om          =   0.5     !  relaxation constant  
    47    cw          =   5.0e-03 !  drag coefficient for oceanic stress 
    48    pstar       =   2.0e+04 !  1st bulk-rheology parameter 
    49    c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    50    creepl      =   1.0e-12 !  creep limit 
    51    ecc         =   2.0     !  eccentricity of the elliptical yield curve 
    52    ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    53    nevp        = 120       !  number of iterations for subcycling in EVP 
    54    relast      = 0.333     !  ratio of elastic timescale over ice time step (1/3 if nevp=120 ; 1/9 if nevp=300) 
    55    alphaevp    =   1.0     !  coefficient for the solution of internal ice stresses 
    56    hminrhg     =   0.001   !  ice volume (a*h in m) below which ice velocity equal ocean velocity 
     46!------------------------------------------------------------------------------ 
     47&namiceitd     !   Ice discretization 
     48!------------------------------------------------------------------------------ 
     49   nn_catbnd      =    2           !  computation of ice category boundaries based on 
     50                                   !      1: tanh function 
     51                                   !      2: h^(-alpha), function of rn_himean 
     52   rn_himean      =    2.0         !  expected domain-average ice thickness (m), nn_catbnd = 2 only 
    5753/ 
    58 !----------------------------------------------------------------------- 
    59 &namicethd     !   ice thermodynamic 
    60 !----------------------------------------------------------------------- 
    61    hmelt       = -0.15     !  maximum melting at the bottom 
    62    hiccrit     = 0.1       !  ice thickness for lateral accretion  
    63                            !         caution 1.0, 1.0 best value to be used!!! (gilles G.)  ???? 
    64    fraz_swi    = 0         !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
    65    maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom 
    66    vfrazb      = 0.4166667 !  thresold drift speed for frazil ice collecting at the ice bottom 
    67    Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom 
    68    hiclim      = 0.10      !  minimum ice thickness 
    69    hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
    70    parsub      = 1.0       !  switch for snow sublimation or not 
    71    betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66 
    72                            !        betas = 1 -> equipartition, betas < 1 -> more on leads 
    73    kappa_i     = 1.0       !  extinction radiation parameter in sea ice (1.0) 
    74    nconv_i_thd = 50        !  maximal number of iterations for heat diffusion computation 
    75    maxer_i_thd = 0.0001    !  maximal error in temperature for heat diffusion computation 
    76    thcon_i_swi = 1         !  switch for computation of thermal conductivity in the ice 
    77                            !        (0) Untersteiner (1964), (1) Pringle et al. (2007) 
     54!------------------------------------------------------------------------------ 
     55&namicedyn     !   Ice dynamics and transport 
     56!------------------------------------------------------------------------------ 
     57   nn_icestr      =    0           !  ice strength parameteriztaion                       
     58                                   !     0: Hibler_79     P = pstar*<h>*exp(-c_rhg*A) 
     59                                   !     1: Rothrock_75   P = Cf*coeff*integral(wr.h^2)     
     60   ln_icestr_bvf  =    .false.     !  ice strength function brine volume (T) or not (F)      
     61   rn_pe_rdg      =   17.0         !  ridging work divided by pot. energy change in ridging, if nn_icestr = 1 
     62   rn_pstar       =    2.0e+04     !  ice strength thickness parameter (N/m2), nn_icestr = 0  
     63   rn_crhg        =   20.0         !  ice strength conc. parameter (-), nn_icestr = 0        
     64   rn_cio         =    5.0e-03     !  ice-ocean drag coefficient           (-)              
     65   rn_creepl      =    1.0e-12     !  creep limit (s-1)                                    
     66   rn_ecc         =    2.0         !  eccentricity of the elliptical yield curve           
     67   nn_nevp        =  120           !  number of EVP subcycles                              
     68   rn_relast      =    0.333       !  ratio of elastic timescale to ice time step: Telast = dt_ice * rn_relast  
     69                                   !     advised value: 1/3 (rn_nevp=120) or 1/9 (rn_nevp=300) 
     70   nn_ahi0        =    2           !  horizontal diffusivity computation 
     71                                   !     0: use rn_ahi0_ref 
     72                                   !     1: use rn_ahi0_ref x mean grid cell length / ( 2deg mean grid cell length ) 
     73                                   !     2: use rn_ahi0_ref x grid cell length      / ( 2deg mean grid cell length ) 
     74   rn_ahi0_ref    = 350.0          !  horizontal sea ice diffusivity (m2/s)  
     75                                   !     if nn_ahi0 > 0, rn_ahi0_ref is the reference value at a nominal 2 deg resolution 
    7876/ 
    79 !----------------------------------------------------------------------- 
    80 &namicesal     !   ice salinity 
    81 !----------------------------------------------------------------------- 
    82    num_sal     =  2        !  salinity option: 1 -> S = bulk_sal 
    83                            !                   2 -> S = S(z,t) with a simple parameterization 
    84                            !                   3 -> S = S(z) profile of Scwharzacher [1959] 
    85                            !                   4 -> S = S(h) Cox and Weeks [1974] 
    86    bulk_sal    =  4.0      !  if 1 is used, it represents the ice salinity 
    87    sal_G       =  5.00     !  restoring salinity for GD 
    88    time_G      =  1.728e+6 !  restoring time for GD 
    89    sal_F       =  2.00     !  restoring salinity for flushing 
    90    time_F      =  8.640e+5 !  restoring time for flushing 
    91    s_i_max     = 20.0      !  Maximum salinity  
    92    s_i_min     =  0.1      !  Minimum tolerated ice salinity 
    93    s_i_0       =  3.5      !  1st salinity for salinity profile 
    94    s_i_1       =  4.5      !  2nd salinity for salinity profile 
     77!------------------------------------------------------------------------------ 
     78&namicehdf     !   Ice horizontal diffusion 
     79!------------------------------------------------------------------------------ 
     80   nn_convfrq     = 5              !  convergence check frequency of the Crant-Nicholson scheme (perf. optimization) 
    9581/ 
    96 !----------------------------------------------------------------------- 
    97 &namiceitdme   !   parameters for mechanical redistribution of ice  
    98 !----------------------------------------------------------------------- 
    99    ridge_scheme_swi = 0      !  which ridging scheme using (1=Rothrock,else=Hibler79) 
    100    Cs               =   0.50 !  shearing energy contribution to ridging 
    101    Cf               =  17.0  !  ratio of ridging work to PE change in ridging 
    102    fsnowrdg         =   0.5  !  snow fraction that survives in ridging 
    103    fsnowrft         =   0.5  !  snow fraction that survives in rafting 
    104    Gstar            =   0.15 !  fractional area of thin ice being ridged 
    105    astar            =   0.05 !  equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 
    106    Hstar            = 100.0  !  parameter determining the maximum thickness of ridged ice 
    107    raft_swi         =   1    !  rafting or not 
    108    hparmeter        =   0.75 !  threshold thickness for rafting or not 
    109    Craft            =   5.0  !  coefficient used in the rafting function 
    110    ridge_por        =   0.3  !  initial porosity of the ridged ice (typically 0.30) 
    111    partfun_swi      =   1    !  participation function linear, TH75 (0) or exponential Letal07 (1) 
    112    brinstren_swi    =   0    !  (1) use brine volume to diminish ice strength 
     82!------------------------------------------------------------------------------ 
     83&namicethd     !   Ice thermodynamics 
     84!------------------------------------------------------------------------------ 
     85   rn_hnewice  = 0.1               !  thickness for new ice formation in open water (m) 
     86   ln_frazil   = .false.           !  use frazil ice collection thickness as a function of wind (T) or not (F) 
     87   rn_maxfrazb = 0.0               !  maximum fraction of frazil ice collecting at the ice base 
     88   rn_vfrazb   = 0.417             !  thresold drift speed for frazil ice collecting at the ice bottom (m/s) 
     89   rn_Cfrazb   = 5.0               !  squeezing coefficient for frazil ice collecting at the ice bottom 
     90   rn_himin    = 0.10              !  minimum ice thickness (m) used in remapping, must be smaller than rn_hnewice 
     91   rn_betas    = 0.66              !  exponent in lead-ice repratition of snow precipitation 
     92                                   !     betas = 1 -> equipartition, betas < 1 -> more on leads 
     93   rn_kappa_i  = 1.0               !  radiation attenuation coefficient in sea ice (m-1) 
     94   nn_conv_dif = 50                !  maximal number of iterations for heat diffusion computation 
     95   rn_terr_dif = 0.0001            !  maximum temperature after heat diffusion (degC) 
     96   nn_ice_thcon= 1                 !  sea ice thermal conductivity 
     97                                   !     0: k = k0 + beta.S/T (Untersteiner, 1964) 
     98                                   !     1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 
     99   nn_monocat  = 0                 !  virtual ITD mono-category parameterizations (1, jpl = 1 only) or not (0) 
     100                                   !     2: simple piling instead of ridging --- temporary option 
     101                                   !     3: activate G(he) only              --- temporary option 
     102                                   !     4: activate lateral melting only    --- temporary option 
     103  ln_it_qnsice = .true.            !  iterate the surface non-solar flux with surface temperature (T) or not (F) 
    113104/ 
    114 !----------------------------------------------------------------------- 
    115 &namicedia     !   ice diagnostics 
    116 !----------------------------------------------------------------------- 
    117    fmtinf      ='1PE13.5 ' !  format of the output values 
    118    nfrinf      = 4         !  number of variables written in one line 
    119    ntmoy       = 1         !  instantaneous values of ice evolution or averaging 
    120    ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging 
     105!------------------------------------------------------------------------------ 
     106&namicesal     !   Ice salinity 
     107!------------------------------------------------------------------------------ 
     108   nn_icesal   =  2                !  ice salinity option 
     109                                   !     1: constant ice salinity (S=rn_icesal) 
     110                                   !     2: varying salinity parameterization S(z,t) 
     111                                   !     3: prescribed salinity profile S(z), Schwarzacher, 1959 
     112   rn_icesal   =  4.               !  ice salinity (g/kg, nn_icesal = 1 only) 
     113   rn_sal_gd   =  5.               !  restoring ice salinity, gravity drainage (g/kg) 
     114   rn_time_gd  =  1.73e+6          !  restoring time scale, gravity drainage  (s) 
     115   rn_sal_fl   =  2.               !  restoring ice salinity, flushing (g/kg) 
     116   rn_time_fl  =  8.64e+5          !  restoring time scale, flushing (s) 
     117   rn_simax    = 20.               !  maximum tolerated ice salinity (g/kg) 
     118   rn_simin    =  0.1              !  minimum tolerated ice salinity (g/kg) 
    121119/ 
    122 !!----------------------------------------------------------------------- 
    123 !&namicehsb       !  Heat and salt budgets  
    124 !!----------------------------------------------------------------------- 
    125 !/ 
    126  
     120!------------------------------------------------------------------------------ 
     121&namiceitdme   !   Ice mechanical redistribution (ridging and rafting) 
     122!------------------------------------------------------------------------------ 
     123   rn_Cs       =   0.5             !  fraction of shearing energy contributing to ridging 
     124   rn_fsnowrdg =   0.5             !  snow volume fraction that survives in ridging 
     125   rn_fsnowrft =   0.5             !  snow volume fraction that survives in rafting 
     126   nn_partfun  =   1               !  type of ridging participation function 
     127                                   !     0: linear (Thorndike et al, 1975) 
     128                                   !     1: exponential (Lipscomb, 2007 
     129   rn_gstar    =   0.15            !  fractional area of thin ice being ridged (nn_partfun = 0) 
     130   rn_astar    =   0.05            !  exponential measure of ridging ice fraction (nn_partfun = 1) 
     131   rn_hstar    = 100.0             !  determines the maximum thickness of ridged ice (m) (Hibler, 1980) 
     132   ln_rafting  =   .true.          !  rafting activated (T) or not (F) 
     133   rn_hraft    =   0.75            !  threshold thickness for rafting (m) 
     134   rn_craft    =   5.0             !  squeezing coefficient used in the rafting function 
     135   rn_por_rdg  =   0.3             !  porosity of newly ridged ice (Lepparanta et al., 1995) 
     136/ 
Note: See TracChangeset for help on using the changeset viewer.