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_PISCES/EXP00/namelist_ice_lim3 – 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_PISCES/EXP00/namelist_ice_lim3

    r1055 r1113  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     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) 
    111!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     12 
     13!----------------------------------------------------------------------- 
     14&namicerun     !   Share parameters for dynamics/advection/thermo 
     15!----------------------------------------------------------------------- 
     16   ln_limdyn   = .true.    !  ice dynamics (T) or thermodynamics only (F) 
     17   acrit       =  1.0e-02 , 1.0e-02  !  minimum fraction for leads in the Northern (Southern) Hemisphere 
     18   hsndif      =  0.0      !  computation of temperature in snow (=0.0) or not 
     19   hicdif      =  0.0      !  computation of temperature in ice  (=0.0) or not (=9999.0) 
     20   cai         =  1.40e-3  !  atmospheric drag over sea ice 
     21   cao         =  1.00e-3  !  atmospheric drag over ocean 
     22   ln_nicep    = .false.   !  Ice points output for debug (yes or no) 
     23/ 
     24!----------------------------------------------------------------------- 
     25&namiceini     !   ice initialisation 
     26!----------------------------------------------------------------------- 
     27   ttest       =  2.0      !  threshold water temperature for initial sea ice 
     28   hninn       =  0.3      !  initial snow thickness in the north 
     29   hnins       =  0.1      !        "            "          south 
     30   hginn_u     =  3.50     !  initial undeformed ice thickness in the north 
     31   hgins_u     =  1.0      !        "            "              "     south 
     32   aginn_u     =  0.95     !  initial undeformed ice concentration in the north 
     33   agins_u     =  0.9      !        "            "              "         south 
     34   hginn_d     =  0.0      !  initial  deformed  ice thickness     in the north 
     35   hgins_d     =  0.0      !   
     36   aginn_d     =  0.00     !  initial  deformed  ice concentration in the north 
     37   agins_d     =  0.00     !        "            "              "         south 
     38   sinn        =  6.301    !  initial salinity in the north 
     39   sins        =  6.301    !        "            "    south 
     40/ 
     41!----------------------------------------------------------------------- 
     42&namicedyn     !   ice dynamic 
     43!----------------------------------------------------------------------- 
     44   epsd        =   1.0e-20 !  tolerance parameter 
     45   alpha       =   0.5     !  coefficient for semi-implicit coriolis 
     46   dm          =   0.6e+03 !  diffusion constant for dynamics 
     47   nbiter      =   1       !  number of sub-time steps for relaxation 
     48   nbitdr      = 500       !  maximum number of iterations for relaxation 
     49   om          =   0.5     !  relaxation constant  
     50   resl        =   5.0e-05 !  maximum value for the residual of relaxation 
     51   cw          =   5.0e-03 !  drag coefficient for oceanic stress 
     52   angvg       =   0.0     !  turning angle for oceanic stress 
     53   pstar       =   4.0e+04 !  1st bulk-rheology parameter 
     54   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
     55   etamn       =   0.0e+07 !  minimun value for viscosity 
     56   creepl      =   1.0e-12 !  creep limit 
     57   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
     58   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
     59   nevp        = 300       !  number of iterations for subcycling in EVP, SB, 400 
     60   telast      =9600.0     !  timescale for elastic waves, SB, 720.0 
     61   alphaevp    =   1.0     !  coefficient for the solution of internal ice stresses 
     62/ 
     63!----------------------------------------------------------------------- 
     64&namicetrp     !   ice transport 
     65!----------------------------------------------------------------------- 
     66   bound       =   0.      !  boundary conditions (=0.0 no-slip, =1.0 free-slip) 
     67/ 
     68!----------------------------------------------------------------------- 
     69&namicethd     !   ice thermodynamic 
     70!----------------------------------------------------------------------- 
     71   hmelt       = -0.15     !  maximum melting at the bottom 
     72   hiccrit     = 0.1 , 0.1 !  ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
     73                           !         caution 1.0, 1.0 best value to be used!!! (gilles G.)  ???? 
     74   fraz_swi    = 0.0       !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
     75   maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom 
     76   vfrazb      = 0.4166667 !  thresold drift speed for frazil ice collecting at the ice bottom 
     77   Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom 
     78   hicmin      = 0.2       !  ice thickness corr. to max. energy stored in brine pocket 
     79   hiclim      = 0.10      !  minimum ice thickness 
     80   amax        = 0.999     !  maximum lead fraction 
     81   sbeta       = 1.        !  numerical caracteritic of the scheme for diffusion in ice 
     82                           !        Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) 
     83   parlat      = 0.0       !  percentage of energy used for lateral ablation 
     84   hakspl      = 0.25      !  slope of distr. for Hakkinen-Mellor's lateral melting 
     85   hibspl      = 0.5       !  slope of distribution for Hibler's lateral melting 
     86   exld        = 2.0       !  exponent for leads-closure rate 
     87   hakdif      = 1.0       !  coefficient for diffusions of ice and snow 
     88   thth        = 0.2       !  threshold thickness for comp. of eq. thermal conductivity 
     89   hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
     90   parsub      = 1.0       !  switch for snow sublimation or not 
     91   alphs       = 1.0       !  coefficient for snow density when snow ice formation 
     92   betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66 
     93                           !        betas = 1 -> equipartition, betas < 1 -> more on leads 
     94   kappa_i     = 1.0       !  extinction radiation parameter in sea ice (1.0) 
     95   nconv_i_thd = 50        !  maximal number of iterations for heat diffusion computation 
     96   maxer_i_thd = 0.0001    !  maximal error in temperature for heat diffusion computation 
     97   thcon_i_swi = 1         !  switch for computation of thermal conductivity in the ice 
     98                           !        (0) Untersteiner (1964), (1) Pringle et al. (2007) 
     99/ 
     100!----------------------------------------------------------------------- 
     101&namicesal     !   ice salinity 
     102!----------------------------------------------------------------------- 
     103   num_sal     =  2        !  salinity option: 1 -> S = bulk_sal 
     104                           !                   2 -> S = S(z,t) with a simple parameterization 
     105                           !                   3 -> S = S(z) profile of Scwharzacher [1959] 
     106                           !                   4 -> S = S(h) Cox and Weeks [1974] 
     107   bulk_sal    =  4.0      !  if 1 is used, it represents the ice salinity 
     108   sal_G       =  5.00     !  restoring salinity for GD 
     109   time_G      =  1.728e+6 !  restoring time for GD 
     110   sal_F       =  2.00     !  restoring salinity for flushing 
     111   time_F      =  8.640e+5 !  restoring time for flushing 
     112   s_i_max     = 20.0      !  Maximum salinity  
     113   s_i_min     =  0.1      !  Minimum tolerated ice salinity 
     114   s_i_0       =  3.5      !  1st salinity for salinity profile 
     115   s_i_1       =  4.5      !  2nd salinity for salinity profile 
     116/ 
     117!----------------------------------------------------------------------- 
     118&namiceitdme   !   parameters for mechanical redistribution of ice  
     119!----------------------------------------------------------------------- 
     120   ridge_scheme_swi = 0      !  which ridging scheme using (1=Rothrock,else=Hibler79) 
     121   Cs               =   0.50 !  shearing energy contribution to ridging 
     122   Cf               =  17.0  !  ratio of ridging work to PE change in ridging 
     123   fsnowrdg         =   0.5  !  snow fraction that survives in ridging 
     124   fsnowrft         =   0.5  !  snow fraction that survives in rafting 
     125   Gstar            =   0.15 !  fractional area of thin ice being ridged 
     126   astar            =   0.05 !  equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 
     127   Hstar            = 100.0  !  parameter determining the maximum thickness of ridged ice 
     128   raftswi          =   1    !  rafting or not 
     129   hparmeter        =   0.75 !  threshold thickness for rafting or not 
     130   Craft            =   5.0  !  coefficient used in the rafting function 
     131   ridge_por        =   0.3  !  initial porosity of the ridged ice (typically 0.30) 
     132   sal_max_ridge    =  15.0  !  maximum ridged ice salinity 
     133   partfun_swi      =   1    !  participation function linear, TH75 (0) or exponential Letal07 (1) 
     134   transfun_swi     =   0    !  transfer function uniform of H80 (0) or exponential Letal07 (1) 
     135   brinstren_swi    =   0    !  (1) use brine volume to diminish ice strength 
     136/ 
     137!----------------------------------------------------------------------- 
     138&namicedia     !   ice diagnostics 
     139!----------------------------------------------------------------------- 
     140   fmtinf      ='1PE13.5 ' !  format of the output values 
     141   nfrinf      = 4         !  number of variables written in one line 
     142   ntmoy       = 1         !  instantaneous values of ice evolution or averaging 
     143   ninfo       = 1         !  frequency of ouputs on file ice_evolu in case of averaging 
     144/ 
     145!----------------------------------------------------------------------- 
     146&namiceout     !   parameters for outputs 
     147!----------------------------------------------------------------------- 
     148   noumef      =   37      !  number of fields 
     149   add_diag_swi=    1      !  1 -> diagnose distribution in thickness space 
     150                           !  0 -> only simple diagnostics 
    2151! 
    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 !  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere 
    14 !  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0) 
    15 !  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0) 
    16 !  cai         : atmospheric drag over sea ice 
    17 !  cao         : atmospheric drag over ocean 
    18 !  ln_nicep  : Ice points output for debug (yes or no) 
    19 ! 
    20 &namicerun 
    21    ln_limdyn    = .true.  
    22    acrit   = 1.0e-02 , 1.0e-02 
    23    hsndif  = 0.0  
    24    hicdif  = 0.0 
    25    cai     = 1.40e-3 
    26    cao     = 1.00e-3 
    27    ln_nicep  = .false. 
    28 / 
    29 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    30 !       namiceini   parameters for ice initialisation 
    31 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    32 ! 
    33 !  ttest   : threshold water temperature for initial sea ice 
    34 !  hninn   : initial snow thickness in the north 
    35 !  hginn_u : initial undeformed ice thickness in the north 
    36 !  aginn_u : initial undeformed ice concentration in the north 
    37 !  hginn_d : initial  deformed  ice thickness in the north 
    38 !  aginn_d : initial  deformed  ice concentration in the north 
    39 !  hnins   : initial snow thickness in the south 
    40 !  hgins_u : initial undeformed ice thickness in the north 
    41 !  agins_u : initial undeformed ice concentration in the north 
    42 !  hgins_d : initial  deformed  ice thickness in the north 
    43 !  agins_d : initial  deformed  ice concentration in the north 
    44 !  sinn    : initial salinity in the north 
    45 !  sins    : initial salinity in the south 
    46 ! 
    47 &namiceini 
    48    ttest   = 2.0 
    49    hninn   = 0.3 
    50    hginn_u = 3.50 
    51    aginn_u = 0.95 
    52    hginn_d = 0.0 
    53    aginn_d = 0.00 
    54    hnins = 0.1 
    55    hgins_u = 1.0 
    56    agins_u = 0.9 
    57    hgins_d = 0.0 
    58    agins_d = 0.00 
    59    sinn  = 6.301 
    60    sins  = 6.301 
    61 / 
    62 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    63 !       namicedia   parameters for ice diagnostics 
    64 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    65 ! 
    66 !  fmtinf : format of the output values  
    67 !  nfrinf : number of variables written in one line 
    68 !  ntmoy  : instantaneous values of ice evolution or averaging 
    69 !  ninfo  : frequency of ouputs on file ice_evolu in case of averaging 
    70 ! 
    71 &namicedia 
    72    fmtinf = '1PE13.5 ' 
    73    nfrinf = 4 
    74    ntmoy  = 1 
    75    ninfo  = 1 
    76 / 
    77 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    78 !       namicedyn   parameters for ice dynamic 
    79 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    80 ! 
    81 !  epsd   : tolerance parameter 
    82 !  alpha  : coefficient for semi-implicit coriolis 
    83 !  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    84 !  dm     : diffusion constant for dynamics. 
    85 !  nbiter : number of sub-time steps for relaxation 
    86 !  nbitdr : maximum number of iterations for relaxation 
    87 !  om     : relaxation constant  
    88 !  resl   : maximum value for the residual of relaxation  
    89 !  cw     : drag coefficient for oceanic stress 
    90 !  angvg  : turning angle for oceanic stress 
    91 !  pstar  : first bulk-rheology parameter 
    92 !  c_rhg  : second bulk-rhelogy parameter 
    93 !  etamn  : minimun value for viscosity 
    94 !  creepl : creep limit 
    95 !  ecc    : eccentricity of the elliptical yield curve 
    96 !  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s) 
    97 !  nevp   : number of iterations for subcycling in EVP, SB, 400 
    98 !  telest : timescale for elastic waves, SB, 720.0 
    99 !  alphaevp : coefficient for the solution of internal ice stresses 
    100 ! 
    101 &namicedyn  
    102    epsd   = 1.0e-20 
    103    alpha  = 0.5 
    104    dm     = 0.6e+03 
    105    nbiter = 1 
    106    nbitdr = 500 
    107    om     = 0.5 
    108    resl   = 5.0e-05 
    109    cw     = 5.0e-03 
    110    angvg  = 0.0 
    111    pstar  = 4.00e+04  
    112    c_rhg  = 20.0 
    113    etamn  = 0.0e+07 
    114    creepl = 1.0e-12 
    115    ecc    = 2.0 
    116    ahi0   = 350.e0 
    117    nevp   = 300 
    118    telast = 9600.0 
    119    alphaevp = 1.0 
    120 / 
    121 ! 
    122 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    123 !       namicetrp   parameters for ice advection 
    124 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    125 !  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    126 &namicetrp  
    127    bound  = 0. 
    128 / 
    129 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    130 !       namicethd   parameters for thermodynamic computation 
    131 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    132 ! 
    133 !  hmelt       : maximum melting at the bottom 
    134 !  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
    135 !                caution 1.0, 1.0 best value to be used!!! (gilles G.) 
    136 !  fraz_swi    : use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
    137 !  maxfrazb    : maximum portion of frazil ice collecting at the ice bottom 
    138 !  vfrazb      : thresold drift speed for frazil ice collecting at the ice bottom 
    139 !  Cfrazb      : squeezing coefficient for frazil ice collecting at the ice bottom 
    140 !  hicmin      : ice thickness corr. to max. energy stored in brine pocket 
    141 !  hiclim      : minimum ice thickness 
    142 !  amax        : maximum lead fraction 
    143 !  sbeta       : numerical caracteritic of the scheme for diffusion in ice 
    144 !               Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) 
    145 !  parlat      : percentage of energy used for lateral ablation 
    146 !  hakspl      : slope of distr. for Hakkinen-Mellor's lateral melting 
    147 !  hibspl      : slope of distribution for Hibler's lateral melting 
    148 !  exld        : exponent for leads-closure rate 
    149 !  hakdif      : coefficient for diffusions of ice and snow 
    150 !  thth         : threshold thickness for comp. of eq. thermal conductivity 
    151 !  hnzst       : thickness of the surf. layer in temp. computation 
    152 !  parsub      : switch for snow sublimation or not 
    153 !  alphs       : coefficient for snow density when snow ice formation 
    154 !  betas       : exponent in lead-ice fractionation of snow precipitation 0.66 
    155 !                betas = 1 -> equipartition, betas < 1 -> more on leads 
    156 !  kappa_i     : extinction radiation parameter in sea ice (1.0) 
    157 !  nconv_i_thd : maximal number of iterations for heat diffusion computation 
    158 !  maxer_i_thd : maximal error in temperature for heat diffusion computation 
    159 !  thcon_i_swi : switch for computation of thermal conductivity in the ice 
    160 !                (0) Untersteiner (1964), (1) Pringle et al. (2007) 
    161 !  
    162 &namicethd 
    163    hmelt   = -0.15 
    164    hiccrit = 0.1 , 0.1 
    165    fraz_swi= 0.0 
    166    maxfrazb= 0.0 
    167    vfrazb  = 0.4166667 
    168    Cfrazb  = 5.0 
    169    hicmin  = 0.2 
    170    hiclim  = 0.10 
    171    amax    = 0.999  
    172    sbeta   = 1. 
    173    parlat  = 0.0 
    174    hakspl  = 0.25 
    175    hibspl  = 0.5 
    176    exld    = 2.0 
    177    hakdif  = 1.0 
    178    thth    = 0.2 
    179    hnzst   = 0.1 
    180    parsub  = 1.0 
    181    alphs   = 1.0 
    182    betas   = 0.6 
    183    kappa_i = 1.0  
    184    nconv_i_thd = 50 
    185    maxer_i_thd = 0.0001 
    186    thcon_i_swi = 1 
    187 / 
    188 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    189 !       namicesal   parameters for salinity calculation 
    190 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    191 ! 
    192 !   num_sal = number of salinity option 
    193 !      1 -> S = bulk_sal 
    194 !      2 -> S = S(z,t) with a simple parameterization 
    195 !      3 -> S = S(z) profile of Scwharzacher [1959] 
    196 !      4 -> S = S(h) Cox and Weeks [1974] 
    197 !   bulk_sal = if 1 is used, it represents the ice salinity 
    198 !   sal_G  : restoring salinity for GD  
    199 !   time_G : restoring time for GD 
    200 !   sal_F  : restoring salinity for flushing 
    201 !   time_F : restoring time for flushing  
    202 !   s_i_max : Maximum salinity      
    203 !   s_i_min : Minimum tolerated ice salinity    
    204 !   s_i_0   : 1st salinity for salinity profile  
    205 !   s_i_1   : 2nd salinity for salinity profile  
    206  
    207 ! 
    208 &namicesal 
    209    num_sal  = 2 
    210    bulk_sal = 4.0 
    211    sal_G    = 5.00 
    212    time_G   = 1.728e+06 
    213    sal_F    = 2.00 
    214    time_F   = 8.640e+05 
    215    s_i_max  = 20.0 
    216    s_i_min  =  0.1 
    217    s_i_0    =  3.5 
    218    s_i_1    =  4.5 
    219 / 
    220 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    221 !       namiceitdme   parameters for mechanical redistribution of ice  
    222 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    223 ! 
    224 !   ridge_scheme_swi ! which ridging scheme using (1=Rothrock,else=Hibler79) 
    225 !   Cs               ! shearing energy contribution to ridging  
    226 !   Cf               ! ratio of ridging work to PE change in ridging         
    227 !   fsnowrdg         ! snow fraction that survives in ridging 
    228 !   fsnowrdg         ! snow fraction that survives in rafting 
    229 !   Gstar            ! fractional area of thin ice being ridged 
    230 !   astar            ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 
    231 !   Hstar            ! parameter determining the maximum thickness of ridged ice 
    232 !   raftswi          ! rafting or not 
    233 !   hparmeter        ! threshold thickness for rafting or not 
    234 !   Craft            ! coefficient used in the rafting function 
    235 !   ridge_por        ! intial porosity of the ridged ice (typically 0.30) 
    236 !   sal_max_ridge    ! maximum ridged ice salinity 
    237 !   partfun_swi      ! participation function linear, TH75 (0) or exponential Letal07 (1) 
    238 !   transfun_swi     ! transfer function uniform of H80 (0) or exponential Letal07 (1) 
    239 !   brinstren_swi    ! (1) use brine volume to diminish ice strength 
    240 &namiceitdme 
    241    ridge_scheme_swi = 0 
    242    Cs = 0.50 
    243    Cf = 17.0 
    244    fsnowrdg  = 0.5 
    245    fsnowrft  = 0.5 
    246    Gstar     = 0.15 
    247    astar     = 0.05 
    248    Hstar     = 100.0 
    249    raftswi   = 1 
    250    hparmeter = 0.75 
    251    Craft     = 5.0 
    252    ridge_por = 0.3 
    253    sal_max_ridge = 15.0 
    254    partfun_swi   = 1 
    255    transfun_swi  = 0 
    256    brinstren_swi = 0 
    257 / 
    258 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    259 !       namiceout   parameters for outputs 
    260 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    261 ! 
    262 !   noumef      : number of fields 
    263 !   field 
    264 !        title  : title of the field 
    265 !        name   : name of the field 
    266 !        unit   : unit of the field 
    267 !        nc     : switch for saving field ( = 1 ) or not ( = 0 ) 
    268 !        cmulti : multiplicative constant 
    269 !        cadd   : additive constant 
    270 !  add_diag_swi : 1 -> diagnose distribution in thickness space 
    271 !                 0 -> only simple diagnostics 
    272 ! 
    273 &namiceout 
    274    noumef   = 37 
    275    field_1  = 'Ice concentration                  ', 'iiceconc', '%       ', 1, 1.0, 0.0 
    276    field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0 
    277    field_3  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0 
    278    field_4  = 'Daily bottom thermo ice production ', 'iicebopr', 'cm/day  ', 1, 100., 0.0 
    279    field_5  = 'Daily dynamic ice production       ', 'iicedypr', 'cm/day  ', 1, 100., 0.0 
    280    field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0 
    281    field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0 
    282    field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0 
    283    field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15 
    284    field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0 
    285    field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0 
    286    field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0 
    287    field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0 
    288    field_14 = 'Heat flux due to brine release     ', 'iocehebr', 'w/m2    ', 1, 1.0, 0.0 
    289    field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0 
    290    field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0  
    291    field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0 
    292    field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0 
    293    field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0 
    294    field_20 = 'Mean ice salinity                  ', 'iicesali', 'psu     ', 1, 1.0, 0.0 
    295    field_21 = 'Mean ice age                       ', 'iiceages', 'years   ', 1, 0.002739, 0.0 
    296    field_22 = 'Daily lateral thermo ice prod.     ', 'iicelapr', 'cm/day  ', 1, 100., 0.0 
    297    field_23 = 'Daily snowice ice production       ', 'iicesipr', 'cm/day  ', 1, 100., 0.0 
    298    field_24 = 'Mean ice temperature               ', 'iicetemp', 'C       ', 1, 1.0, -273.15 
    299    field_25 = 'Ice total heat content             ', 'iiceheco', '10^9 J  ', 1, 1.0, 0.0 
    300    field_26 = 'Ice surface temperature            ', 'iicesurt', 'C       ', 1, 1.0, -273.15 
    301    field_27 = 'Snow temperature                   ', 'isnotem2', 'C       ', 1, 1.0, -273.15 
    302    field_28 = 'Fsbri - brine salt flux            ', 'iicfsbri', 'kg/m2/s ', 1, 1.0, 0.0 
    303    field_29 = 'Fseqv - equivalent FW salt flux    ', 'iicfseqv', 'kg/m2/s ', 1, 1.0, 0.0 
    304    field_30 = 'Brine volume                       ', 'ibrinvol', '%       ', 1, 100.0, 0.0 
    305    field_31 = 'Frazil ice collection thickness    ', 'iicecolf', 'm       ', 1, 1.0, 0.0 
    306    field_32 = 'Ice strength                       ', 'iicestre', 'N/m     ', 1, 0.001, 0.0 
    307    field_33 = 'Ice velocity                       ', 'iicevelo', 'm/s     ', 1, 1.0, 0.0 
    308    field_34 = 'Surface melt                       ', 'iicesume', 'cm/day  ', 1, 100., 0.0 
    309    field_35 = 'Bottom melt                        ', 'iicebome', 'cm/day  ', 1, 100., 0.0 
    310    field_36 = 'Divergence                         ', 'iicedive', '10-8s-1 ', 1, 1.0e8, 0.0 
    311    field_37 = 'Shear                              ', 'iiceshea', '10-8s-1 ', 1, 1.0e8, 0.0 
    312    add_diag_swi = 1 
     152!           !         title of the field           !  name     !   units   !  save  ! multipl. ! additive ! 
     153!           !                                      !           !           ! or not !  factor  !  factor  ! 
     154   field_1  = 'Ice concentration                  ', 'iiceconc', '%       ',    1   ,  1.0     ,    0.0 
     155   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0 
     156   field_3  = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0 
     157   field_4  = 'Daily bottom thermo ice production ', 'iicebopr', 'cm/day  ',    1   , 100.     ,    0.0 
     158   field_5  = 'Daily dynamic ice production       ', 'iicedypr', 'cm/day  ',    1   , 100.     ,    0.0 
     159   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0 
     160   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0 
     161   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0 
     162   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15 
     163   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0 
     164   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     165   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     166   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0 
     167   field_14 = 'Heat flux due to brine release     ', 'iocehebr', 'w/m2    ',    1   ,  1.0     ,    0.0 
     168   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0 
     169   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0  
     170   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     171   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
     172   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ',    1   ,  1.0     ,    0.0 
     173   field_20 = 'Mean ice salinity                  ', 'iicesali', 'psu     ',    1   ,  1.0     ,    0.0 
     174   field_21 = 'Mean ice age                       ', 'iiceages', 'years   ',    1   ,  0.002739,    0.0 
     175   field_22 = 'Daily lateral thermo ice prod.     ', 'iicelapr', 'cm/day  ',    1   ,100.      ,    0.0 
     176   field_23 = 'Daily snowice ice production       ', 'iicesipr', 'cm/day  ',    1   ,100.      ,    0.0 
     177   field_24 = 'Mean ice temperature               ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15 
     178   field_25 = 'Ice total heat content             ', 'iiceheco', '10^9 J  ',    1   ,  1.0     ,    0.0 
     179   field_26 = 'Ice surface temperature            ', 'iicesurt', 'C       ',    1   ,  1.0     , -273.15 
     180   field_27 = 'Snow temperature                   ', 'isnotem2', 'C       ',    1   ,  1.0     , -273.15 
     181   field_28 = 'Fsbri - brine salt flux            ', 'iicfsbri', 'kg/m2/s ',    1   ,  1.0     ,    0.0 
     182   field_29 = 'Fseqv - equivalent FW salt flux    ', 'iicfseqv', 'kg/m2/s ',    1   ,  1.0     ,    0.0 
     183   field_30 = 'Brine volume                       ', 'ibrinvol', '%       ',    1   ,  100.0   ,    0.0 
     184   field_31 = 'Frazil ice collection thickness    ', 'iicecolf', 'm       ',    1   ,  1.0     ,    0.0 
     185   field_32 = 'Ice strength                       ', 'iicestre', 'N/m     ',    1   ,  0.001   ,    0.0 
     186   field_33 = 'Ice velocity                       ', 'iicevelo', 'm/s     ',    1   ,  1.0     ,    0.0 
     187   field_34 = 'Surface melt                       ', 'iicesume', 'cm/day  ',    1   ,100.      ,    0.0 
     188   field_35 = 'Bottom melt                        ', 'iicebome', 'cm/day  ',    1   ,100.      ,    0.0 
     189   field_36 = 'Divergence                         ', 'iicedive', '10-8s-1 ',    1   ,  1.0e8   ,    0.0 
     190   field_37 = 'Shear                              ', 'iiceshea', '10-8s-1 ',    1   ,  1.0e8   ,    0.0 
    313191/       
Note: See TracChangeset for help on using the changeset viewer.