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 1113 for trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist_ice_lim2 – NEMO

Ignore:
Timestamp:
2008-06-14T16:58:05+02:00 (16 years ago)
Author:
rblod
Message:

Namelist style for dynamics and sea-ice, see ticket #197

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/CONFIG/ORCA2_LIM/EXP00/1_namelist_ice_lim2

    r833 r1113  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/LIM2 :  1 - dynamics/advection/thermo          (namicerun) 
     3!! namelists    2 - ice intialisation                  (namiceini) 
     4!!              3 - ice dynamic                        (namicedyn) 
     5!!              5 - ice advection                      (namicetrp) 
     6!!              6 - thermodynamic                      (namicethd) 
     7!!              7 - ice salinity                       (namicesal) 
     8!!              8 - mechanical redistribution of ice   (namiceitdme) 
     9!!              3 - ice diagnostics                    (namicedia) 
     10!!              9 - ice outputs                        (namiceout) 
    111!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! 
    3 ! ICE namelist :  ice model option and parameter input 
    4 ! ------------- 
    5 ! 
    6 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    7 ! 
    8 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    9 !       namicerun   Share parameters for dynamics/advection/thermo 
    10 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    11 ! 
    12 !  ln_limdyn   : switch for ice dynamics (true) or not (false) 
    13 !  ln_limdmp   : restoring Ice thickness and Fraction leads flag 
    14 !  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere 
    15 !  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0) 
    16 !  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0) 
    17 ! 
    18 &namicerun 
    19    ln_limdyn = .true.  
    20    ln_limdmp = .FALSE.  
    21    acrit     = 1.0e-06 , 1.0e-06 
    22    hsndif    = 0.0  
    23    hicdif    = 0.0 
     12 
     13!----------------------------------------------------------------------- 
     14&namicerun     !   Share parameters for dynamics/advection/thermo 
     15!----------------------------------------------------------------------- 
     16   ln_limdyn   = .true.    !  ice dynamics (T) or thermodynamics only (F) 
     17   ln_limdmp   = .false.   !  restoring ice thickness and fraction leads (T) or not (F) 
     18   acrit       = 1.0e-06   , 1.0e-06   ! minimum fraction for leads in the Northern (Southern) Hemisphere 
     19   hsndif      =  0.0      !  computation of temperature in snow (=0.0) or not 
     20   hicdif      =  0.0      !  computation of temperature in ice  (=0.0) or not (=9999.0) 
    2421/ 
    25 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    26 !       namiceini   parameters for ice initialisation 
    27 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    28 ! 
    29 !  ln_limini : Ice initialization state flag 
    30 !  ttest     : threshold water temperature for initial sea ice 
    31 !  hninn     : initial snow thickness in the north 
    32 !  hginn     : initial ice thickness in the north 
    33 !  alinn     : initial leads area in the north 
    34 !  hnins     : initial snow thickness in the south 
    35 !  hgins     : initial ice thickness in the south 
    36 !  alins     : initial leads area in the south 
    37 ! 
    38 &namiceini 
    39    ln_limini = .FALSE.  
    40    ttest = 2.0 
    41    hninn = 0.5 
    42    hginn = 3.0 
    43    alinn = 0.05 
    44    hnins = 0.1 
    45    hgins = 1.0 
    46    alins = 0.1 
     22!----------------------------------------------------------------------- 
     23&namiceini     !   ice initialisation 
     24!----------------------------------------------------------------------- 
     25   ln_limini   = .false.   !  read the ice initial state in the file 'Ice_initialization.nc' (T) or not (F) 
     26   ttest       =  2.0      !  threshold water temperature for initial sea ice 
     27   hninn       =  0.5      !  initial snow thickness in the north 
     28   hginn       =  3.0      !  initial ice  thickness in the north 
     29   alinn       =  0.05     !  initial leads area     in the north 
     30   hnins       =  0.1      !  same  three parameter  in the south 
     31   hgins       =  1.0      !        "                 "     south    
     32   alins       =  0.1      !        "                 "     south 
    4733/ 
    48 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    49 !       namicedia   parameters for ice diagnostics 
    50 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    51 ! 
    52 !  fmtinf : format of the output values  
    53 !  nfrinf : number of variables written in one line 
    54 !  ntmoy  : instantaneous values of ice evolution or averaging 
    55 !  ninfo  : frequency of ouputs on file ice_evolu in case of averaging 
    56 ! 
    57 &namicedia 
    58    fmtinf = '1PE13.5 ' 
    59    nfrinf = 4 
    60    ntmoy  = 1 
    61    ninfo  = 1 
     34!----------------------------------------------------------------------- 
     35&namicedyn     !   ice dynamic 
     36!----------------------------------------------------------------------- 
     37   epsd        =   1.0e-20 !  tolerance parameter 
     38   alpha       =   0.5     !  coefficient for semi-implicit coriolis 
     39   dm          =   0.6e+03 !  diffusion constant for dynamics 
     40   nbiter      =   1       !  number of sub-time steps for relaxation 
     41   nbitdr      = 100       !  maximum number of iterations for relaxation 
     42   om          =   0.5     !  relaxation constant  
     43   resl        =   5.0e-05 !  maximum value for the residual of relaxation 
     44   cw          =   5.0e-03 !  drag coefficient for oceanic stress 
     45   angvg       =   0.0     !  turning angle for oceanic stress 
     46   pstar       =   1.0e+04 !  1st bulk-rheology parameter 
     47   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
     48   etamn       =   0.0e+07 !  minimun value for viscosity 
     49   creepl      =   2.0e-08 !  creep limit 
     50   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
     51   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    6252/ 
    63 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    64 !       namicedyn   parameters for ice dynamic 
    65 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    66 ! 
    67 !  epsd   : tolerance parameter 
    68 !  alpha  : coefficient for semi-implicit coriolis 
    69 !  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    70 !  dm     : diffusion constant for dynamics. 
    71 !  nbiter : number of sub-time steps for relaxation 
    72 !  nbitdr : maximum number of iterations for relaxation 
    73 !  om     : relaxation constant  
    74 !  resl   : maximum value for the residual of relaxation  
    75 !  cw     : drag coefficient for oceanic stress 
    76 !  angvg  : turning angle for oceanic stress 
    77 !  pstar  : first bulk-rheology parameter 
    78 !  c_rhg  : second bulk-rhelogy parameter 
    79 !  etamn  : minimun value for viscosity 
    80 !  creepl : creep limit 
    81 !  ecc    : eccentricity of the elliptical yield curve 
    82 !  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s) 
    83 ! 
    84 &namicedyn  
    85    epsd   = 1.0e-20 
    86    alpha  = 0.5 
    87    dm     = 0.6e+03 
    88    nbiter = 1 
    89    nbitdr = 100 
    90    om     = 0.5 
    91    resl   = 5.0e-05 
    92    cw     = 5.0e-03 
    93    angvg  = 0.0 
    94    pstar  = 1.0e+04  
    95    c_rhg  = 20.0 
    96    etamn  = 0.0e+07 
    97    creepl = 2.0e-08 
    98    ecc    = 2.0 
    99    ahi0   = 350.e0 
     53!----------------------------------------------------------------------- 
     54&namicetrp     !   ice transport 
     55!----------------------------------------------------------------------- 
     56   bound       =   0.      !  boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    10057/ 
    101 ! 
    102 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    103 !       namicetrp   parameters for ice advection 
    104 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    105 !  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    106 &namicetrp  
    107    bound  = 0. 
    108 / 
    109  
    110 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    111 !       namicethd   parameters for thermodynamic computation 
    112 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    113 ! 
     58!----------------------------------------------------------------------- 
     59&namicethd     !   ice thermodynamic 
     60!----------------------------------------------------------------------- 
    11461!  hmelt       : maximum melting at the bottom 
    11562!  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
    … …  
    13279!  
    13380&namicethd 
    134    hmelt   = -0.15 
    135    hiccrit = 0.3 , 0.3 
    136    hicmin  = 0.2 
    137    hiclim  = 0.05 
    138    amax    = 0.999  
    139    swiqst  = 1. 
    140    sbeta   = 1. 
    141    parlat  = 0.0 
    142    hakspl  = 0.5 
    143    hibspl  = 0.5 
    144    exld    = 2.0 
    145    hakdif  = 1.0 
    146    thth    = 0.2 
    147    hnzst   = 0.1 
    148    parsub  = 1.0 
    149    alphs   = 1.0 
     81   hmelt       = -0.15     !  maximum melting at the bottom 
     82   hiccrit     = 0.3 , 0.3 !  ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
     83   hicmin      = 0.2       !  ice thickness corr. to max. energy stored in brine pocket 
     84   hiclim      = 0.05      !  minimum ice thickness 
     85   amax        = 0.999     !  maximum lead fraction 
     86   swiqst      = 1.        !  energy stored in brine pocket (=1) or not (=0)   
     87   sbeta       = 1.        !  numerical caracteritic of the scheme for diffusion in ice 
     88                           !        Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) 
     89   parlat      = 0.0       !  percentage of energy used for lateral ablation 
     90   hakspl      = 0.5       !  slope of distr. for Hakkinen-Mellor's lateral melting 
     91   hibspl      = 0.5       !  slope of distribution for Hibler's lateral melting 
     92   exld        = 2.0       !  exponent for leads-closure rate 
     93   hakdif      = 1.0       !  coefficient for diffusions of ice and snow 
     94   thth        = 0.2       !  threshold thickness for comp. of eq. thermal conductivity 
     95   hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
     96   parsub      = 1.0       !  switch for snow sublimation or not 
     97   alphs       = 1.0       !  coefficient for snow density when snow ice formation 
    15098/ 
     99!----------------------------------------------------------------------- 
     100&namicedia     !   ice diagnostics 
     101!----------------------------------------------------------------------- 
     102   fmtinf      ='1PE13.5 ' !  format of the output values 
     103   nfrinf      = 4         !  number of variables written in one line 
     104   ntmoy       = 1         !  instantaneous values of ice evolution or averaging 
     105   ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging 
     106/ 
     107!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     108&namiceout     !   parameters for outputs 
     109!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     110   noumef      =   19      !  number of fields 
    151111! 
    152 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    153 !       namiceout   parameters for outputs 
    154 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    155 ! 
    156 !   noumef      : number of fields 
    157 !   field 
    158 !        title  : title of the field 
    159 !        name   : name of the field 
    160 !        unit   : unit of the field 
    161 !        nc     : switch for saving field ( = 1 ) or not ( = 0 ) 
    162 !        cmulti : multiplicative constant 
    163 !        cadd   : additive constant 
    164 ! 
    165 &namiceout 
    166    noumef   = 19 
    167    field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0 
    168    field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0 
    169    field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ', 1, 1.0, 0.0 
    170    field_4  = 'Ice concentration                  ', 'ileadfra', '%       ', 1, -1.0, 1.0 
    171    field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ', 1, 1.0, -273.15 
    172    field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0 
    173    field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0 
    174    field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0 
    175    field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15 
    176    field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0 
    177    field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0 
    178    field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0 
    179    field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0 
    180    field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt', 1, 1.0, 0.0 
    181    field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0 
    182    field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0  
    183    field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0 
    184    field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0 
    185    field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0 
     112!           !         title of the field           !  name     !   units   !  save  ! multipl. ! additive ! 
     113!           !                                      !           !           ! or not !  factor  !  factor  ! 
     114   field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0 
     115   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0 
     116   field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ',    1   ,  1.0     ,    0.0 
     117   field_4  = 'Ice concentration                  ', 'ileadfra', '%       ',    1   , -1.0     ,    1.0 
     118   field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15 
     119   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0 
     120   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0 
     121   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0 
     122   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15 
     123   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0 
     124   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     125   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     126   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0 
     127   field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt',    1   ,  1.0     ,    0.0 
     128   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0 
     129   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0  
     130   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     131   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     132   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ',    1   ,  1.0     ,    0.0 
    186133/       
Note: See TracChangeset for help on using the changeset viewer.