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 1121 for trunk/CONFIG/ORCA2_OFF_PISCES/EXP00 – NEMO

Ignore:
Timestamp:
2008-06-20T17:23:23+02:00 (16 years ago)
Author:
cetlod
Message:

change the style of all top namelist to be consistent with the forthcoming style of OPA namelist, see ticket:196

Location:
trunk/CONFIG/ORCA2_OFF_PISCES/EXP00
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r983 r1121  
    1 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! OPA namelist :  model option and parameter input 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    5 ! 
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/OPA  :  1 - run manager      (namrun) 
     3!! namelists    2 - miscellaneous    (namctl,nammpp) 
     4!!              3 - Domain           (namzgr, namzgr_sco, namdom) 
     5!!              6 - Tracer           (nameos, namcla, namqsr) 
     6!!              7 - Inputs dynamics  (namdyna) 
     7!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     8!  CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. 
     9 
     10!!====================================================================== 
     11!!                   ***  Run management namelists  *** 
     12!!====================================================================== 
     13!!   namrun            parameters of the run 
     14!!====================================================================== 
     15 
    616!----------------------------------------------------------------------- 
    7 !       nam_run   parameters of the run 
     17&namrun        !   parameters of the run 
    818!----------------------------------------------------------------------- 
    9 !  no         job number 
    10 !  cexper     experience name for vairmer format 
    11 !  ln_rstart  boolean term for restart (true or false) 
    12 !  nrstdt     control of the restart timestep: 
    13 !                = 0 restart, do not control nit000 in the restart file. 
    14 !                = 1 restart, control nit000 in the restart file. Do not 
    15 !                    use the date in the restart file (use ndate0 in namelist) 
    16 !                = 2 restart, control nit000 in the restart file, use the date 
    17 !                    in the restart file. ndate0 in the namelist is ignored. 
    18 !  nit000     number of the first time step 
    19 !  nitend     number of the last time step 
    20 !  ndate0     initial calendar date aammjj 
    21 !  nleapy     Leap year calendar (0/1) 
    22 !  ninist     initial state output flag (0/1) 
    23 !  nprint     level of print (0 no print) 
    24 !  nstock     frequency of restart file 
    25 !  nwrite     frequency of OUTPUT file 
    26 !  nrunoff    = 0 no, 1 runoff, 2 runoff+river mouth ups adv 
    27 !  ln_ctl     trend control print (expensive!) 
    28 !  nictls     start i indice to make the control SUM (very usefull to compare mono- 
    29 !  nictle     end   i indice to make the control SUM (-versus multi processor runs) 
    30 !  njctls     start j indice to make the control SUM (very usefull to compare mono- 
    31 !  njctle     end   j indice to make the control SUM (-versus multi processor runs) 
    32 !  isplt      number of processors following i 
    33 !  jsplt      number of processors following j 
    34 !  nbench     Bench parameter (0/1): CAUTION it must be zero except for bench 
    35 !             for which we don't care about physical meaning of the results 
    36 ! 
    37 !  CAUTION: for usual run scripts, logical value of 
    38 !  *******  ln_rstart must be .true. or .false. 
    39 !                     and NOT .TRUE. or .FALSE. 
    40 &nam_run 
    41    no         =       0 
    42    cexper     =  "PISCES" 
    43    ln_rstart  = .false. 
    44    nrstdt     =       0 
    45    nit000     =       1 
    46    nitend     =    6000  
    47    ndate0     =  010101 
    48    nleapy     =       0 
    49    ninist     =       0 
    50    nprint     =       0 
    51    nstock     =    6000  
    52    nwrite     =    1200  
    53    ln_ctl     =  .false. 
    54    nictls     =       2 
    55    nictle     =       720 
    56    njctls     =       2 
    57    njctle     =       510 
    58    isplt      =       1 
    59    jsplt      =       1 
    60    nbench     =       0 
     19   no          =       0   !  job number 
     20   cexper      =  "PISCES"  !  experience name  
     21   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     22   nrstdt      =       0   !  restart control = 0 nit000 is not compared to the restart file value 
     23                           !                  = 1 use ndate0 in namelist (not the value in the restart file) 
     24                           !                  = 2 calendar parameters read in the restart file 
     25   nit000      =       1   !  first time step 
     26   nitend      =     6000  !  last  time step 
     27   ndate0      =  010101   !  initial calendar date yymmdd (used if nrstdt=1) 
     28   nleapy      =       0   !  Leap year calendar (1) or not (0) 
     29   ninist      =       0   !  output the initial state (1) or not (0) 
     30   nstock      =     6000  !  frequency of creation of a restart file 
     31   nwrite      =     6000  !  frequency of write in the output file  
    6132/ 
    6233!----------------------------------------------------------------------- 
    63 !       nam_mpp      Massively Parallel Processing 
     34&namctl       !   Control prints & Benchmark 
    6435!----------------------------------------------------------------------- 
    65 !  c_mpi_send         mpi send/recieve type 
    66 !                      = 'S'  : standard blocking send 
    67 !                      = 'B'  : buffer blocking send 
    68 !                      = 'I'  : immediate non-blocking send 
    69 &nam_mpp 
    70    c_mpi_send =  'S' 
     36   ln_ctl     = .false.    !  trends control print (expensive!) 
     37   nprint     =    0       !  level of print (0 no extra print) 
     38   nictls     =    1       !  start i indice of control sum (use to compare mono versus 
     39   nictle     =    182       !  end   i indice of control sum        multi processor runs 
     40   njctls     =    1       !  start j indice of control               over a subdomain) 
     41   njctle     =    149       !  end   j indice of control  
     42   isplt      =    1       !  number of processors in i-direction 
     43   jsplt      =    1       !  number of processors in j-direction 
     44   nbench     =    0       !  Bench mode (1/0): CAUTION use zero except for bench 
    7145/ 
    7246!----------------------------------------------------------------------- 
    73 !       nam_zgr       vertical coordinate 
     47&nammpp      !   Massively Parallel Processing                         ("key_mpp_mpi) 
    7448!----------------------------------------------------------------------- 
    75 !  ln_zco     z-coordinate - full steps      (T/F) 
    76 !  ln_zps     z-coordinate - partial steps   (T/F) 
    77 !  ln_sco     s- or hybrid z-s-coordinate    (T/F) 
    78 &nam_zgr 
    79    ln_zco   =  .false. 
    80    ln_zps   =  .true. 
    81    ln_sco   =  .false. 
     49   c_mpi_send =  'S'       !  mpi send/recieve type   ='S', 'B', or 'I' for standard send, 
     50                           !  buffer blocking send or immediate non-blocking sends, resp. 
    8251/ 
    8352!----------------------------------------------------------------------- 
    84 !       nam_zgr_sco   s-coordinate or hybrid z-s-coordinate 
     53&namzgr       !   vertical coordinate 
    8554!----------------------------------------------------------------------- 
    86 !  sbot_min   minimum depth of s-bottom surface (>0) (m) 
    87 !  sbot_max   maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    88 !  theta      surface control parameter (0<=theta<=20) 
    89 !  thetb      bottom control parameter  (0<=thetb<= 1) 
    90 !  r_max      maximum cut-off r-value allowed (0<r_max<1) 
    91 &nam_zgr_sco 
    92    sbot_min =  300. 
    93    sbot_max = 5250. 
    94    theta    =    6.0 
    95    thetb    =    0.75 
    96    r_max    =    0.15 
     55   ln_zco      = .false.   !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined) 
     56   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
     57   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
    9758/ 
    9859!----------------------------------------------------------------------- 
    99 !       nam_traadv   advection scheme for tracer (option not control by CPP keys) 
     60&namzgr_sco   !   s-coordinate or hybrid z-s-coordinate 
    10061!----------------------------------------------------------------------- 
    101 !  ln_traadv_cen2     2nd order centered scheme    (default T) 
    102 !  ln_traadv_tvd      TVD scheme                   (default F) 
    103 !  ln_traadv_muscl    MUSCL scheme                 (default F) 
    104 !  ln_traadv_muscl2   MUSCL2 scheme                (default F) 
    105 &nam_traadv 
    106    ln_traadv_cen2   =  .true. 
    107    ln_traadv_tvd    =  .false. 
    108    ln_traadv_muscl  =  .false. 
    109    ln_traadv_muscl2 =  .false. 
     62   sbot_min    =  300.     !  minimum depth of s-bottom surface (>0) (m) 
     63   sbot_max    = 5250.     !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     64   theta       =    6.0    !  surface control parameter (0<=theta<=20) 
     65   thetb       =    0.75   !  bottom control parameter  (0<=thetb<= 1) 
     66   r_max       =    0.15   !  maximum cut-off r-value allowed (0<r_max<1) 
    11067/ 
    11168!----------------------------------------------------------------------- 
    112 !       nam_traldf   lateral diffusion scheme for tracer (option not control by CPP keys) 
     69&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    11370!----------------------------------------------------------------------- 
    114 !  Type of the operator :  
    115 !     ln_traldf_lap    laplacian operator          (default T) 
    116 !     ln_traldf_bilap  bilaplacian operator        (default F) 
    117 !  Direction of action  :  
    118 !     ln_traldf_level  iso-level                   (default F) 
    119 !     ln_traldf_hor    horizontal (geopotential)   (default F)^** 
    120 !     ln_traldf_iso    iso-neutral                 (default T)^* 
    121 !  Coefficient 
    122 !     aht0    horizontal eddy diffusivity for tracers (m2/s) 
    123 !     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s) 
    124 !     aeiv0   eddy induced velocity coefficient (m2/s) 
    125 ! ^* require key_ldfslp to compute the direction of the lateral diffusion 
    126 ! ^** require key_ldfslp in s-coordinate 
    127 &nam_traldf 
    128    ln_traldf_lap    =  .true. 
    129    ln_traldf_bilap  =  .false. 
    130    ln_traldf_level  =  .false. 
    131    ln_traldf_hor    =  .false. 
    132    ln_traldf_iso    =  .true. 
    133    aht0    =  2000. 
    134    ahtb0   =     0. 
    135    aeiv0   =  2000. 
     71   e3zps_min   =    5.     !  the thickness of the partial step is set larger than the minimum 
     72   e3zps_rat   =    0.1    !  of e3zps_min and e3zps_rat * e3t   (N.B. 0<e3zps_rat<1) 
     73   nmsh        =    1      !  create (=1) a mesh file (coordinates, scale factors, masks) or not (=0) 
     74   nacc        =    0      !  =1 acceleration of convergence method used, rdt < rdttra(k) 
     75                           !  =0, no acceleration, rdt = rdttra 
     76   atfp        =    0.1    !  asselin time filter parameter 
     77   rdt         = 26280.    !  time step for the dynamics (and tracer if nacc=0) 
     78   rdtmin      = 26280.    !  minimum time step on tracers (used if nacc=1) 
     79   rdtmax      = 26280.    !  maximum time step on tracers (used if nacc=1) 
     80   rdth        =  800.     !  depth variation of tracer time step  (used if nacc=1) 
    13681/ 
    13782!----------------------------------------------------------------------- 
    138 !       nam_dynldf   lateral diffusion on momentum 
     83&namtraldf    !   lateral diffusion scheme for tracer 
    13984!----------------------------------------------------------------------- 
    140 !  Type of the operator :  
    141 !     ln_dynldf_lap    laplacian operator          (default T) 
    142 !     ln_dynldf_bilap  bilaplacian operator        (default F) 
    143 !  Direction of action  :  
    144 !     ln_dynldf_level  iso-level                   (default F) 
    145 !     ln_dynldf_hor    horizontal (geopotential)   (default F)^** 
    146 !     ln_dynldf_iso    iso-neutral                 (default T)^* 
    147 !  Coefficient 
    148 !  ahm0    horizontal eddy viscosity for the dynamics (m2/s) 
    149 !  ahmb0   background eddy viscosity for isopycnal diffusion (m2/s) 
    150 &nam_dynldf 
    151    ln_dynldf_lap    =  .true. 
    152    ln_dynldf_bilap  =  .false. 
    153    ln_dynldf_level  =  .false. 
    154    ln_dynldf_hor    =  .true. 
    155    ln_dynldf_iso    =  .false. 
    156    ahm0    = 40000. 
    157    ahmb0   =     0. 
     85!                               !  Type of the operator : 
     86   ln_traldf_lap    =  .true.   !     laplacian operator 
     87   ln_traldf_bilap  =  .false.  !     bilaplacian operator 
     88                                !  Direction of action  : 
     89   ln_traldf_level  =  .false.  !     iso-level 
     90   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
     91   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     92!                               !  Coefficient 
     93   aht0        =  2000.         !     horizontal eddy diffusivity for tracers [m2/s] 
     94   ahtb0       =     0.         !     background eddy diffusivity for ldf_iso [m2/s] 
     95   aeiv0       =  2000.         !     eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
    15896/ 
    15997!----------------------------------------------------------------------- 
    160 !       namflg   algorithm flags (algorithm not control by CPP keys) 
     98&namcla        !   cross land advection 
    16199!----------------------------------------------------------------------- 
    162 !  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T) 
    163 !                                                  centered      time scheme  (F) 
    164 &namflg 
    165    ln_dynhpg_imp   =  .false. 
     100   n_cla       =    0      !  advection between 2 ocean pts separates by land 
    166101/ 
    167102!----------------------------------------------------------------------- 
    168 !       nam_dynvor   option of physics/algorithm (not control by CPP keys) 
     103&namqsr        !   penetrative solar radiation 
    169104!----------------------------------------------------------------------- 
    170 !  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T) 
    171 !  ln_dynvor_ene      "         "  : energy conserving scheme    (default F) 
    172 !  ln_dynvor_mix      "         "  : mixed scheme                (default F) 
    173 !  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F) 
    174 &nam_dynvor 
    175    ln_dynvor_ene = .FALSE. 
    176    ln_dynvor_ens = .TRUE. 
    177    ln_dynvor_mix = .FALSE. 
    178    ln_dynvor_een = .FALSE. 
     105   ln_traqsr   = .true.    !  penetrative solar radiation (T) or not (F) 
     106   xsi1        =   0.35    !  first depth of extinction 
    179107/ 
    180108!----------------------------------------------------------------------- 
    181 !       nam_tau   surface wind stress 
     109&nameos        !   ocean physical parameters 
    182110!----------------------------------------------------------------------- 
    183 !  ntau000   gently increase the stress over the first ntau_rst time-steps 
    184 !  tau0x     uniform value used as default surface heat flux 
    185 !  tau0y     uniform value used as default solar radiation flux 
    186 &nam_tau 
    187    ntau000 =      0 
    188    tau0x   =      0.e0 
    189    tau0y   =      0.e0 
     111   neos        =    0      !  type of equation of state and Brunt-Vaisala frequency 
     112                           !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
     113                           !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
     114                           !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
     115   ralpha      =    2.e-4  !  thermal expension coefficient (neos= 1 or 2) 
     116   rbeta       =    0.001  !  saline  expension coefficient (neos= 2) 
    190117/ 
    191118!----------------------------------------------------------------------- 
    192 !       nam_flx   surface fluxes 
     119&namdyn        !   offline parameters  
    193120!----------------------------------------------------------------------- 
    194 !  q0       uniform value used as default surface heat flux 
    195 !  qsr0     uniform value used as default solar radiation flux 
    196 !  emp0     uniform value used as default surface freswater budget (E-P) 
    197 &nam_flx 
    198    q0      =      0.e0 
    199    qsr0    =      0.e0 
    200    emp0    =      0.e0 
    201 / 
    202 !----------------------------------------------------------------------- 
    203 !       nam_alb   albedo parameters 
    204 !----------------------------------------------------------------------- 
    205 !  cgren    correction of the snow or ice albedo to take into account 
    206 !  albice   albedo of melting ice in the arctic and antarctic 
    207 !  alphd    coefficients for linear interpolation used to compute albedo 
    208 !           between two extremes values (Pyane, 1972) 
    209 !  alphc     "                                         " 
    210 !  alphdi    "                                         " 
    211 &nam_alb 
    212    cgren    =      0.06 
    213    albice   =      0.5 
    214    alphd    =      0.80 
    215    alphc    =      0.65 
    216    alphdi   =      0.72 
    217 / 
    218 !----------------------------------------------------------------------- 
    219 !       nam_dom   space and time domain (bathymetry, mesh, timestep) 
    220 !----------------------------------------------------------------------- 
    221 !  e3zps_min  the thickness of the partial step is set larger than the 
    222 !  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t 
    223 !                (N.B. 0<e3zps_rat<1) 
    224 !  nmsh       =1 create a mesh file (coordinates, scale factors, masks) 
    225 !  nacc       the acceleration of convergence method 
    226 !             = 0, no acceleration, rdt = rdttra 
    227 !             = 1, acceleration used, rdt < rdttra(k) 
    228 !  atfp       asselin time filter parameter 
    229 !  rdt        time step for the dynamics (and tracer if nacc=0) 
    230 !  rdtmin     minimum time step on tracers 
    231 !  rdtmax     maximum time step on tracers 
    232 !  rdth       depth variation of tracer time step 
    233 &nam_dom 
    234    e3zps_min =     5. 
    235    e3zps_rat =     0.1 
    236    nmsh      =     1 
    237    nacc      =     0 
    238    atfp      =     0.1 
    239    rdt       =  26280. 
    240    rdtmin    =  26280. 
    241    rdtmax    =  26280. 
    242    rdth      =   800. 
    243 / 
    244 !----------------------------------------------------------------------- 
    245 !       nam_fwb   freshwater budget correction 
    246 !----------------------------------------------------------------------- 
    247 !  ln_fwb     logical flag for freshwater budget correction (0 annual mean) 
    248 &nam_fwb 
    249    ln_fwb    = .true. 
    250 / 
    251 !----------------------------------------------------------------------- 
    252 !       nam_ptr   Poleward Transport Diagnostic 
    253 !----------------------------------------------------------------------- 
    254 !  ln_diaptr  logical flag for Poleward transport computation 
    255 !  nf_ptr     Frequency of computation 
    256 &nam_ptr 
    257    ln_diaptr = .true. 
    258    nf_ptr    = 15 
    259 / 
    260 !----------------------------------------------------------------------- 
    261 !       nam_cro   cross land advection 
    262 !----------------------------------------------------------------------- 
    263 !  n_cla   advection between 2 ocean pts separates by land  
    264 &nam_cla 
    265    n_cla   = 0 
    266 / 
    267 !----------------------------------------------------------------------- 
    268 !       nam_zdf   vertical physics 
    269 !----------------------------------------------------------------------- 
    270 !  avt0       vertical eddy diffusivity for tracers (m2/s) 
    271 !  ln_zdfnpc  Non-Penetrative Convection          (default T) 
    272 &nam_zdf 
    273    avt0     = 1.2e-5 
    274    ln_zdfnpc = .false. 
    275 / 
    276 !----------------------------------------------------------------------- 
    277 !       nam_bbl   bottom boundary layer scheme 
    278 !----------------------------------------------------------------------- 
    279 !  atrbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s) 
    280 &nam_bbl 
    281    atrbbl = 10000. 
    282 / 
    283 !----------------------------------------------------------------------- 
    284 !       nam_ddm   double diffusive mixing parameterization 
    285 !----------------------------------------------------------------------- 
    286 !   avts    maximum avs for dd mixing  
    287 !   hsbfr   heat/salt buoyancy flux ratio 
    288 &nam_ddm 
    289       avts  = 1.e-4 
    290       hsbfr = 1.6 
    291 / 
    292 !----------------------------------------------------------------------- 
    293 !       nam_bbc   bottom temperature boundary condition 
    294 !----------------------------------------------------------------------- 
    295 !  ngeo_flux  = 0 no geothermal heat flux 
    296 !             = 1 constant geothermal heat flux 
    297 !             = 2 variable geothermal heat flux (read in geothermal_heating.nc)  
    298 !                 ( C A U T I O N : flux in mW/m2 in the NetCDF file ) 
    299 !  ngeo_flux_const   Constant value of geothermal heat flux (W/m2)  
    300 &nam_bbc 
    301    ngeo_flux =  2 
    302    ngeo_flux_const = 86.4e-3 
    303 / 
    304 !----------------------------------------------------------------------- 
    305 !       nam_qsr   penetrative solar radiation 
    306 !----------------------------------------------------------------------- 
    307 !  ln_traqsr : penetrative solar radiation (T) or not (F)     (Default=T) 
    308 !  rabs       fraction of qsr associated with xsi1 
    309 !  xsi1       first depth of extinction 
    310 !  xsi2       second depth of extinction 
    311 &nam_qsr 
    312    xsi1     =   0.35 
    313 / 
    314 !----------------------------------------------------------------------- 
    315 !       nam_eos   ocean physical parameters 
    316 !----------------------------------------------------------------------- 
    317 !  neos    type of equation of state and Brunt-Vaisala frequency 
    318 !          = 0, UNESCO (formulation of Jackett and McDougall (1994) 
    319 !                                         and of McDougall (1987) ) 
    320 !          = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    321 !          = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    322 !                               with rau0=1020 set in parcst routine 
    323 !  ralpha  thermal expension coefficient (linear equation of state) 
    324 !  rbeta   saline  expension coefficient (linear equation of state) 
    325 &nam_eos 
    326    neos   =      0 
    327    ralpha =  2.e-4 
    328    rbeta  =  0.001 
    329 / 
    330 !----------------------------------------------------------------------- 
    331 !       nam_offdyn    offline parameters  
    332 !----------------------------------------------------------------------- 
    333 !  ndtadyn   number of period in the file for one year 
    334 !  ndtatot   total number of period in the file 
    335 !  nsptint   indicator for time interpolation 
    336 !  nficdyn    number of file to read 
    337 !  lperdyn  = T periodicity of the unique file 
    338 !             = F  (default)   computed with Blanke' scheme 
    339 &nam_offdyn 
    340     ndtadyn = 73  
    341     ndtatot  =  73 
    342     nsptint  =    1 
    343     nficdyn   = 2 
    344     lperdyn = .true. 
     121    ndtadyn   =  73        ! number of period in the file for one year     
     122    ndtatot   =  73        ! total number of period in the file     
     123    nsptint   =  1         ! indicator for time interpolation     
     124    nficdyn   =  2         ! number of file to read     
     125    lperdyn   = .true.     ! periodicity of the unique file (T) 
     126!                          ! F  (default)   computed with Blanke' scheme     
    345127/ 
    346128 

  • trunk/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_pisces

    r1093 r1121  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! PISCES  :    1  - air-sea exchange                         (nampisext) 
     3!! namelists    2  - biological parameters                    (nampisbio) 
     4!!              3  - parameters for nutrient limitations      (nampislim)     
     5!!              4  - parameters for phytoplankton             (nampisprod,nampismort) 
     6!!              5  - parameters for zooplankton               (nampismes,nampiszoo) 
     7!!              6  - parameters for remineralization          (nampisrem) 
     8!!              7  - parameters for calcite chemistry         (nampiscal) 
     9!!              8  - parameters for inputs deposition         (nampissed) 
     10!!              9  - parameters for Kriest parameterization   (nampiskrp, nampiskrs) 
     11!!              10 - additional 2D/3D  diagnostics            (nampisdia) 
    112!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! General namelist for PISCES model 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     13!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     14&nampisext     !   air-sea exchange 
     15!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     16   atcco2     = 278.    ! atmospheric pCO2 
     17/ 
     18!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     19&nampisbio     !   biological parameters 
     20!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     21   part       =  0.85    ! part of calcite not dissolved in guts 
     22   nrdttrc    =  4       ! time step frequency for biology 
     23   wsbio      =  2.      ! POC sinking speed 
     24   xkmort     =  1.E-7   ! half saturation constant for mortality 
     25   ferat3     =  3.E-6   ! Fe/C in zooplankton  
     26   wsbio2     =  50.     ! Big particles sinking speed 
     27/ 
     28!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     29&nampislim     !   parameters for nutrient limitations 
     30!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     31   conc0      =  2.e-6    ! Phosphate half saturation 
     32   conc1      =  10E-6    ! Phosphate half saturation for diatoms 
     33   conc2      =  0.02E-9  ! Iron half saturation for phyto 
     34   conc2m     =  0.08E-9  ! Max iron half saturation for phyto 
     35   conc3      =  0.1E-9   ! Iron half saturation for diatoms 
     36   conc3m     =  0.4E-9   ! Maxi iron half saturation for diatoms 
     37   concnnh4   =  1.E-7    ! NH4 half saturation for phyto 
     38   concdnh4   =  5.E-7    ! NH4 half saturation for diatoms 
     39   xksi1      =  2.E-6    ! half saturation constant for Si uptake 
     40   xksi2      =  3.33E-6  ! half saturation constant for Si/C 
     41   xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
     42   caco3r     =  0.3      ! mean rain ratio 
     43/ 
     44!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     45&nampisprod     !   parameters for phytoplankton growth 
     46!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     47   pislope    =  4.       ! P-I slope   
     48   pislope2   =  4.       ! P-I slope  for diatoms 
     49   excret     =  0.05     ! excretion ratio of phytoplankton 
     50   excret2    =  0.05     ! excretion ratio of diatoms 
     51   chlcnm     =  0.033    ! Minimum Chl/C in nanophytoplankton 
     52   chlcdm     =  0.05     ! Minimum Chl/C in diatoms 
     53   fecnm      =  10E-6    ! Maximum Fe/C in nanophytoplankton 
     54   fecdm      =  15E-6    ! Minimum Fe/C in diatoms 
     55   grosip     =  0.151    ! mean Si/C ratio 
     56/ 
     57!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     58&nampismort     !   parameters for phytoplankton sinks 
     59!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     60   wchl       =  0.001    ! quadratic mortality of phytoplankton 
     61   wchld      =  0.02     ! maximum quadratic mortality of diatoms 
     62   mprat      =  0.01     ! phytoplankton mortality rate 
     63   mprat2     =  0.01     ! Diatoms mortality rate 
     64   mpratm     =  0.01     ! Phytoplankton minimum mortality rate 
     65/ 
     66!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     67&nampismes     !   parameters for mesozooplankton 
     68!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     69   grazrat2   =  0.7      ! maximal mesozoo grazing rate 
     70   resrat2    =  0.005    ! exsudation rate of mesozooplankton 
     71   mzrat2     =  0.03     ! mesozooplankton mortality rate 
     72   xprefc     =  1.       ! zoo preference for phyto 
     73   xprefp     =  0.2      ! zoo preference for POC 
     74   xprefz     =  1.       ! zoo preference for zoo 
     75   xprefpoc   =  0.2      ! zoo preference for poc 
     76   xkgraz2    =  20.E-6   ! half sturation constant for meso grazing 
     77   epsher2    =  0.33     ! Efficicency of Mesozoo growth  
     78   sigma2     =  0.6      ! Fraction of mesozoo excretion as DOM 
     79   unass2     =  0.3      ! non assimilated fraction of P by mesozoo 
     80   grazflux   =  5.e3     ! flux-feeding rate 
     81/ 
     82!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     83&nampiszoo     !   parameters for microzooplankton 
     84!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     85   grazrat    =  4.0      ! maximal zoo grazing rate    
     86   resrat     =  0.03     ! exsudation rate of zooplankton 
     87   mzrat      =  0.0      ! zooplankton mortality rate 
     88   xpref2c    =  0.0      ! Microzoo preference for POM  
     89   xpref2p    =  0.5      ! Microzoo preference for Nanophyto 
     90   xpref2d    =  0.5      ! Microzoo preference for Diatoms 
     91   xkgraz     =  20.E-6   ! half sturation constant for grazing  
     92   epsher     =  0.33     ! Efficiency of microzoo growth 
     93   sigma1     =  0.6      ! Fraction of microzoo excretion as DOM 
     94   unass      =  0.3      ! non assimilated fraction of phyto by zoo 
     95/ 
     96!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     97&nampisrem     !   parameters for remineralization 
     98!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     99   xremik    =  0.3       ! remineralization rate of DOC 
     100   xremip    =  0.025     ! remineralisation rate of POC 
     101   nitrif    =  0.05      ! NH4 nitrification rate 
     102   xsirem    =  0.015     ! remineralization rate of Si 
     103   xlam1     =  0.005     ! scavenging rate of Iron 
     104   oxymin    =  1.E-6     ! Half-saturation constant for anoxia  
     105/ 
     106!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     107&nampiscal     !   parameters for Calcite chemistry 
     108!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     109   kdca       =  0.327e3  ! calcite dissolution rate constant (1/time) 
     110   nca        =  1.       ! order of dissolution reaction (dimensionless) 
     111/ 
     112!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     113&nampissed     !   parameters for inputs deposition 
     114!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     115   bdustfer   =  .true.   ! boolean for dust input from the atmosphere 
     116   briver     =  .true.  ! boolean for river input of nutrients 
     117   bndepo     =  .true.  ! boolean for atmospheric deposition of N 
     118   bsedinput  =  .true.  ! boolean for Fe input from sediments 
     119   sedfeinput =  1E-9     ! Coastal release of Iron 
     120   dustsolub  =  0.014    ! Solubility of the dust 
     121/ 
     122!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     123&nampiskrp     !   Kriest parameterization : parameters     "key_kriest" 
     124!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     125   xkr_eta      = 1.17    ! Sinking  exponent 
     126   xkr_zeta     = 2.28    ! N content exponent 
     127   xkr_mass_min = 0.0002  ! Minimum mass for Aggregates 
     128   xkr_mass_max = 1.      ! Maximum mass for Aggregates 
     129/ 
     130!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     131&nampiskrs     !   Kriest parameterization : size classes  "key_kriest" 
     132!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     133   xkr_sfact    = 942.    ! Sinking factor 
     134   xkr_stick    = 0.5     ! Stickiness 
     135   xkr_nnano    = 2.337   ! Nbr of cell in nano size class 
     136   xkr_ndiat    = 3.718   ! Nbr of cell in diatoms size class 
     137   xkr_nmeso    = 7.147   ! Nbr of cell in mesozoo size class 
     138   xkr_naggr    = 9.877   ! Nbr of cell in aggregates  size class 
     139/ 
     140!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     141&nampisdia     !   additional 2D/3D tracers diagnostics ("key_trc_diaadd") 
     142!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     143   nwritedia   =  6000     !  time step frequency for tracers diagnostics 
    5144! 
    6 !----------------------------------------------------------------------- 
    7 !   natext   air-sea exchange 
    8 !----------------------------------------------------------------------- 
    9 !   atcco2 : atmospheric pCO2 
    10 ! 
    11 &natext 
    12    atcco2 = 278., 
     145!              !    name   !           title of the field          !     units      ! 
     146!              !           !                                       !                !   
     147   pisdia2d(1)  = 'Cflx     ' , 'DIC flux                          ',  'molC/m2/s    ' 
     148   pisdia2d(2)  = 'Oflx     ' , 'Oxygen flux                       ',  'molC/m2/s    ' 
     149   pisdia2d(3)  = 'Kg       ' , 'Gas transfer                      ',  'mol/m2/s/uatm' 
     150   pisdia2d(4)  = 'Delc     ' , 'Delta CO2                         ',  'uatm         ' 
     151   pisdia2d(5)  = 'PMO      ' , 'POC export                        ',  'molC/m2/s    ' 
     152   pisdia2d(6)  = 'PMO2     ' , 'GOC export                        ',  'molC/m2/s    ' 
     153   pisdia2d(7)  = 'ExpFe1   ' , 'Nano iron export                  ',  'molFe/m2/s   ' 
     154   pisdia2d(8)  = 'ExpFe2   ' , 'Diatoms iron export               ',  'molFe/m2/s   ' 
     155   pisdia2d(9)  = 'ExpSi    ' , 'Silicate export                   ',  'molSi/m2/s   ' 
     156   pisdia2d(10) = 'ExpCaCO3 ' , 'Calcite export                    ',  'molC/m2/s    ' 
     157   pisdia2d(11) = 'heup     ' , 'euphotic layer depth              ',  'm            ' 
     158   pisdia2d(12) = 'Fedep    ' , 'Iron dep                          ',  'molFe/m2/s   ' 
     159   pisdia2d(13) = 'Nfix     ' , 'Nitrogen Fixation                 ',  'molN/m2/s    ' 
     160   pisdia3d(1)  = 'PH       ' , 'PH                                ',  '-            ' 
     161   pisdia3d(2)  = 'CO3      ' , 'Bicarbonates                      ',  'mol/l        ' 
     162   pisdia3d(3)  = 'CO3sat   ' , 'CO3 saturation                    ',  'mol/l        ' 
     163   pisdia3d(4)  = 'PAR      ' , 'light penetration                 ',  'W/m2         ' 
     164   pisdia3d(5)  = 'PPPHY    ' , 'Primary production of nanophyto   ',  'molC/m3/s    ' 
     165   pisdia3d(6)  = 'PPPHY2   ' , 'Primary production of diatoms     ',  'molC/m3/s    ' 
     166   pisdia3d(7)  = 'PPZOO    ' , 'Primary production of microzoo    ',  'molC/m3/s    ' 
     167   pisdia3d(8)  = 'PPZOO2   ' , 'Primary production of mesozoo     ',  'molC/m3/s    ' 
     168   pisdia3d(9)  = 'PBSi     ' , 'Primary production of Si diatoms  ',  'molSi/m3/s   ' 
     169   pisdia3d(10) = 'PFeN     ' , 'Primary production of nano iron   ',  'molFe/m3/s   ' 
     170   pisdia3d(11) = 'PFeD     ' , 'Primary production of diatoms iron',  'molFe/m3/s   ' 
    13171/ 
    14 !----------------------------------------------------------------------- 
    15 !   natbio   biological parameters 
    16 !----------------------------------------------------------------------- 
    17 !  part        part of calcite not dissolved in guts 
    18 !  nrdttrc     time step frequency for biology 
    19 !  wsbio       POC sinking speed 
    20 !  xkmort      half saturation constant for mortality 
    21 !  ferat3      Fe/C in zooplankton 
    22 !  wsbio2      Big particles sinking speed 
    23 ! 
    24 &natbio 
    25    part       =  0.85, 
    26    nrdttrc    =  4, 
    27    wsbio      =  2. , 
    28    xkmort     =  1.E-7 , 
    29    ferat3     =  3.E-6, 
    30    wsbio2     =  50., 
    31 / 
    32 !----------------------------------------------------------------------- 
    33 !   natlim   parameters for nutrient limitations 
    34 !----------------------------------------------------------------------- 
    35 !  conc0       Phosphate half saturation 
    36 !  conc1       Phosphate half saturation for diatoms 
    37 !  conc2       Iron half saturation for phyto 
    38 !  conc3       Iron half saturation for diatoms 
    39 !  concnnh4    NH4 half saturation for phyto 
    40 !  concdnh4    NH4 half saturation for diatoms 
    41 !  xksi1       half saturation constant for Si uptake 
    42 !  xksi2       half saturation constant for Si/C 
    43 !  xkdoc       half-saturation constant of DOC remineralization 
    44 !  caco3r      mean rain ratio 
    45 &natlim 
    46    conc0      =  2.e-6, 
    47    conc1      =  10E-6 , 
    48    conc2      =  0.02E-9 , 
    49    conc2m     =  0.08E-9, 
    50    conc3      =  0.1E-9 , 
    51    conc3m     =  0.4E-9, 
    52    concnnh4   =  1.E-7, 
    53    concdnh4   =  5.E-7, 
    54    xksi1      =  2.E-6 , 
    55    xksi2      =  3.33E-6 , 
    56    xkdoc      =  417.E-6 , 
    57    caco3r     =  0.3, 
    58 / 
    59 !----------------------------------------------------------------------- 
    60 !   natprod   parameters for phytoplankton growth 
    61 !----------------------------------------------------------------------- 
    62 !  pislope     P-I slope 
    63 !  pislope2    P-I slope  for diatoms 
    64 !  excret      excretion ratio of phytoplankton 
    65 !  excret2     excretion ratio of diatoms 
    66 !  chlcnm      Minimum Chl/C in nanophytoplankton 
    67 !  chlcdm      Minimum Chl/C in diatoms 
    68 !  fecnm       Maximum Fe/C in nanophytoplankton 
    69 !  fecdm       Minimum Fe/C in diatoms 
    70 !  grosip      mean Si/C ratio 
    71 &natprod 
    72    pislope    =  4. , 
    73    pislope2   =  4. , 
    74    excret     =  0.05 , 
    75    excret2    =  0.05 , 
    76    chlcnm     =  0.033, 
    77    chlcdm     =  0.05, 
    78    fecnm      =  10E-6, 
    79    fecdm      =  15E-6, 
    80    grosip     =  0.151, 
    81 / 
    82 !----------------------------------------------------------------------- 
    83 !   natmort   parameters for phytoplankton sinks 
    84 !----------------------------------------------------------------------- 
    85 !  wchl        quadratic mortality of phytoplankton 
    86 !  wchld       maximum quadratic mortality of diatoms 
    87 !  mprat       phytoplankton mortality rate 
    88 !  mprat2      Diatoms mortality rate 
    89 !  mpratm      Phytoplankton minimum mortality rate 
    90 &natmort 
    91    wchl       =  0.001, 
    92    wchld      =  0.02, 
    93    mprat      =  0.01, 
    94    mprat2     =  0.01, 
    95    mpratm     =  0.01, 
    96 / 
    97 !----------------------------------------------------------------------- 
    98 !   natmes   parameters for mesozooplankton 
    99 !----------------------------------------------------------------------- 
    100 !  grazrat2    maximal mesozoo grazing rate 
    101 !  resrat2     exsudation rate of mesozooplankton 
    102 !  mzrat2      mesozooplankton mortality rate 
    103 !  xprefc      zoo preference for phyto 
    104 !  xprefp      zoo preference for POC 
    105 !  xprefz      zoo preference for zoo 
    106 !  xprefpoc    zoo preference for poc 
    107 !  xkgraz2     half sturation constant for meso grazing 
    108 !  epsher2     Efficicency of Mesozoo growth 
    109 !  sigma2      Fraction of mesozoo excretion as DOM 
    110 !  unass2      non assimilated fraction of P by mesozoo 
    111 !  grazflux    flux-feeding rate 
    112 &natmes 
    113    grazrat2   =  0.7, 
    114    resrat2    =  0.005, 
    115    mzrat2     =  0.03, 
    116    xprefc     =  1. , 
    117    xprefp     =  0.2 , 
    118    xprefz     =  1 , 
    119    xprefpoc   =  0.2 , 
    120    xkgraz2    =  20.E-6 , 
    121    epsher2    =  0.33 , 
    122    sigma2     =  0.6, 
    123    unass2     =  0.3 , 
    124    grazflux   =  5.e3, 
    125 / 
    126 !----------------------------------------------------------------------- 
    127 !   natbio   biological parameters 
    128 !----------------------------------------------------------------------- 
    129 !  grazrat     maximal zoo grazing rate 
    130 !  resrat      exsudation rate of zooplankton 
    131 !  mzrat       zooplankton mortality rate 
    132 !  xpref2c      Microzoo preference for POM 
    133 !  xpref2p      Microzoo preference for Nanophyto 
    134 !  xpref2d      Microzoo preference for Diatoms 
    135 !  xkgraz      half sturation constant for grazing 
    136 !  epsher      Efficiency of microzoo growth 
    137 !  sigma1      Fraction of microzoo excretion as DOM 
    138 !  unass       non assimilated fraction of phyto by zoo 
    139 &natzoo 
    140    grazrat    =  4.0, 
    141    resrat     =  0.03, 
    142    mzrat      =  0.0, 
    143    xpref2c    =  0.0 , 
    144    xpref2p    =  0.5 , 
    145    xpref2d    =  0.5 , 
    146    xkgraz     =  20.E-6 , 
    147    epsher     =  0.33 , 
    148    sigma1     =  0.6, 
    149    unass      =  0.3 , 
    150 / 
    151 ! 
    152 !----------------------------------------------------------------------- 
    153 !   natrem   parameters for remineralization 
    154 !----------------------------------------------------------------------- 
    155 !  xremik      remineralization rate of DOC 
    156 !  xremip      remineralisation rate of POC 
    157 !  nitrif      NH4 nitrification rate 
    158 !  xsirem      remineralization rate of Si 
    159 !  xlam1       scavenging rate of Iron 
    160 !  oxymin      Half-saturation constant for anoxia 
    161 &natrem 
    162    xremik    =  0.3, 
    163    xremip    =  0.025, 
    164    nitrif    =  0.05, 
    165    xsirem    =  0.015, 
    166    xlam1     =  0.005, 
    167    oxymin    =  1.E-6, 
    168 / 
    169 !----------------------------------------------------------------------- 
    170 !   natcal   Calcite chemistry 
    171 !----------------------------------------------------------------------- 
    172 !  kdca        calcite dissolution rate constant (1/time) 
    173 !  nca         order of dissolution reaction (dimensionless) 
    174 &natcal 
    175    kdca       =  0.327e3, 
    176    nca        =  1., 
    177 / 
    178 !----------------------------------------------------------------------- 
    179 !   natsms   inputs deposition 
    180 !----------------------------------------------------------------------- 
    181 !   bdustfer    boolean for dust input from the atmosphere 
    182 !   briver      boolean for river input of nutrients 
    183 !   bndepo      boolean for atmospheric deposition of N 
    184 !   bsedinput   boolean for Fe input from sediments 
    185 !   sedfeinput  Coastal release of Iron 
    186 ! 
    187 &natsms 
    188    bdustfer   =  .true., 
    189    briver     =  .false., 
    190    bndepo     =  .false., 
    191    bsedinput  =  .false., 
    192    sedfeinput =  1E-9, 
    193    dustsolub  =  0.014, 
    194 / 
    195 !--------------------------------------------------------------------------------- 
    196 !   natkrpar   parameter for Kriest parameterization      (#ifdef key_kriest) 
    197 !--------------------------------------------------------------------------------- 
    198 !  xkr_eta       Sinking  exponent 
    199 !  xkr_zeta      N content exponent 
    200 !  xkr_mass_min  Minimum mass for Aggregates 
    201 !  xkr_mass_max  Maximum mass for Aggregates 
    202 ! 
    203 &natkrpar 
    204    xkr_eta      = 1.17, 
    205    xkr_zeta     = 2.28, 
    206    xkr_mass_min = 0.0002, 
    207    xkr_mass_max = 1., 
    208 / 
    209 !--------------------------------------------------------------------------------- 
    210 !   natkrsize   size calsses for Kriest parameterization      (#ifdef key_kriest) 
    211 !--------------------------------------------------------------------------------- 
    212 !  xkr_sfact     Sinking factor 
    213 !  xkr_stick     Stickiness 
    214 !  xkr_xnnano    Nbr of cell in nano size class 
    215 !  xkr_ndiat     Nbr of cell in diatoms size class 
    216 !  xkr_nmeso     Nbr of cell in mesozoo size class 
    217 !  xkr_naggr     Nbr of cell in aggregates  size class 
    218 ! 
    219 &natkrsize 
    220    xkr_sfact    = 942., 
    221    xkr_stick    = 0.5, 
    222    xkr_nnano    = 2.337, 
    223    xkr_ndiat    = 3.718, 
    224    xkr_nmeso    = 7.147, 
    225    xkr_naggr    = 9.877, 
    226 / 

  • trunk/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top

    r1092 r1121  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! NEMO/TOP1 :  1 - tracer definition                     (namtoptrc) 
     3!! namelists    2 - dynamical tracer trends               (namtoptrd) 
     4!!              6 - tracer advection                      (namtopadv) 
     5!!              7 - tracer bottom boundary                (namtopbbl) 
     6!!              8 - tracer lateral diffusion              (namtopldf) 
     7!!              3 - tracer vertical physics               (namtopzdf) 
     8!!              9 - tracer newtonian damping              (namtopdmp) 
    19!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 ! OPA MODEL general namelist for passive tracers 
    3 ! ------------- 
    4 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     10!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     11&namtoptrc     !   tracers definition 
     12!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     13   ndttrc      =  1        !  time step frequency for passive tracers       
     14   nwritetrc   =  6000     !  time step frequency for tracer outputs 
     15   lrsttr      = .false.   !  start from a restart file (T) or not (F) 
     16   nrsttr      =   0       !  restart control = 0 initial time step is not compared to the restart file value 
     17                           !                  = 1 do not use the value in the restart file 
     18                           !                  = 2 calendar parameters read in the restart file 
    519! 
    6 !----------------------------------------------------------------------- 
    7 !       nattrc   
    8 !----------------------------------------------------------------------- 
    9 ! 
    10 !       NATTRC    
    11 !       nwritetrc time step frequency for tracer outputs 
    12 !       lrsttr    boolean term for tracer model restart (true or false) 
    13 !       nrsttr    control of the time step for tracer model restart (0, 1 or 2) 
    14 !       tracer type defined by : 
    15 !                 *  short name  
    16 !                 *  long_name 
    17 !                 *  units 
    18 !                 *  logical to read initial value from file or not 
    19 !                 *  multiplicative coefficient 
    20 !                 *  logical to save value   
    21 !  
    22 &nattrc 
    23    nwritetrc   = 6000  
    24    lrsttr      = .true. 
    25    nrsttr      = 0 
    26    tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration          ',  'umol/L' ,  .false. ,  .true. 
    27    tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration             ',  'ueq/L ' ,  .false. ,  .true. 
    28    tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration             ',  'umol/L' ,  .false. ,  .true. 
    29    tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                      ',  'umol/L' ,  .false. ,  .true. 
    30    tracer(5)   = 'PO4     ' , 'Phosphate Concentration                    ',  'umol/L' ,  .false. ,  .true. 
    31    tracer(6)   = 'POC     ' , 'Small particle organic carbon Concentration',  'umol/L' ,  .false. ,  .true. 
    32    tracer(7)   = 'Si      ' , 'Silicate Concentration                     ',  'umol/L' ,  .false. ,  .true. 
    33    tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration            ',  'umol/L' ,  .false. ,  .true. 
    34    tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration             ',  'umol/L' ,  .false. ,  .true. 
    35    tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration            ',  'umol/L' ,  .false. ,  .true. 
    36    tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                      ',  'umol/L' ,  .false. ,  .true. 
    37    tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration              ',  'umol/L' ,  .false. ,  .true. 
    38    tracer(13)  = 'BSi     ' , 'Diatoms Silicate Concentration             ',  'umol/L' ,  .false. ,  .true. 
    39    tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration               ',  'umol/L' ,  .false. ,  .true. 
    40    tracer(15)  = 'BFe     ' , 'Big iron particles Concentration           ',  'umol/L' ,  .false. ,  .true. 
    41    tracer(16)  = 'GOC     ' , 'Big particle organic carbon Concentration  ',  'umol/L' ,  .false. ,  .true. 
    42    tracer(17)  = 'SFe     ' , 'Small iron particles Concentration         ',  'umol/L' ,  .false. ,  .true. 
    43    tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration                ',  'umol/L' ,  .false. ,  .true. 
    44    tracer(19)  = 'DSi     ' , 'Sinking biogenic Silicate Concentration    ',  'umol/L' ,  .false. ,  .true. 
    45    tracer(20)  = 'NFe     ' , 'Nano iron Concentration                    ',  'umol/L' ,  .false. ,  .true. 
    46    tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration              ',  'umol/L' ,  .false. ,  .true. 
    47    tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration           ',  'umol/L' ,  .false. ,  .true. 
    48    tracer(23)  = 'NO3     ' , 'Nitrates Concentration                     ',  'umol/L' ,  .false. ,  .true. 
    49    tracer(24)  = 'NH4     ' , 'Ammonium Concentration                     ',  'umol/L' ,  .false. ,  .true. 
     20!              !    name   !           title of the field              !   units   ! initial data ! save   ! 
     21!              !           !                                           !           ! from file    ! or not !  
     22!              !           !                                           !           ! or not       !        ! 
     23   tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'umol/L' ,  .false.     ,  .true. 
     24   tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'ueq/L ' ,  .false.     ,  .true. 
     25   tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'umol/L' ,  .false.     ,  .true. 
     26   tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'umol/L' ,  .false.     ,  .true. 
     27   tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'umol/L' ,  .false.     ,  .true. 
     28   tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'umol/L' ,  .false.     ,  .true. 
     29   tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'umol/L' ,  .false.     ,  .true. 
     30   tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'umol/L' ,  .false.     ,  .true. 
     31   tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'umol/L' ,  .false.     ,  .true. 
     32   tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'umol/L' ,  .false.     ,  .true. 
     33   tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'umol/L' ,  .false.     ,  .true. 
     34   tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'umol/L' ,  .false.     ,  .true. 
     35   tracer(13)  = 'BSi     ' , 'Diatoms Silicate Concentration         ',  'umol/L' ,  .false.     ,  .true. 
     36   tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'umol/L' ,  .false.     ,  .true. 
     37   tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'umol/L' ,  .false.     ,  .true. 
     38   tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'umol/L' ,  .false.     ,  .true. 
     39   tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'umol/L' ,  .false.     ,  .true. 
     40   tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'umol/L' ,  .false.     ,  .true. 
     41   tracer(19)  = 'DSi     ' , 'Sinking biogenic Silicate Concentration',  'umol/L' ,  .false.     ,  .true. 
     42   tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'umol/L' ,  .false.     ,  .true. 
     43   tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'umol/L' ,  .false.     ,  .true. 
     44   tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'umol/L' ,  .false.     ,  .true. 
     45   tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'umol/L' ,  .false.     ,  .true. 
     46   tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'umol/L' ,  .false.     ,  .true. 
    5047/ 
    5148!----------------------------------------------------------------------- 
    52 !       natrtd   dynamical tracers trends (#ifdef key_trc_diatrd)  
     49&namtoptrd     !   dynamical tracers trends ("key_trc_diatrd") 
    5350!----------------------------------------------------------------------- 
    54 !   nwritetrd : time step frequency for dynamical trends outputs 
    55 !   luttrd    : logical to keep large diagnostics with trends or not 
    56 !               one value per tracer 
    57 &natrtd 
    58    nwritetrd  = 6000 
    59    luttrd(1)  = .false. 
    60    luttrd(2)  = .false. 
    61    luttrd(3)  = .false. 
    62    luttrd(4)  = .false. 
    63    luttrd(5)  = .false. 
    64    luttrd(6)  = .false. 
    65    luttrd(7)  = .false. 
    66    luttrd(8)  = .false. 
    67    luttrd(9)  = .false. 
    68    luttrd(10) = .false. 
    69    luttrd(11) = .false. 
    70    luttrd(12) = .false. 
    71    luttrd(13) = .false. 
    72    luttrd(14) = .false. 
    73    luttrd(15) = .false. 
    74    luttrd(16) = .false. 
    75    luttrd(17) = .false. 
    76    luttrd(18) = .false. 
    77    luttrd(19) = .false. 
    78    luttrd(20) = .false. 
    79    luttrd(21) = .false. 
    80    luttrd(22) = .false. 
    81    luttrd(23) = .false. 
    82    luttrd(24) = .false. 
     51   nwritetrd  =    6000   !  time step frequency for dynamical trends outputs 
     52!             !    save trends or not  ! 
     53   luttrd(1)  =   .false. 
     54   luttrd(2)  =   .false. 
     55   luttrd(3)  =   .false. 
     56   luttrd(4)  =   .false. 
     57   luttrd(5)  =   .false. 
     58   luttrd(6)  =   .false. 
     59   luttrd(7)  =   .false. 
     60   luttrd(8)  =   .false. 
     61   luttrd(9)  =   .false. 
     62   luttrd(10) =   .false. 
     63   luttrd(11) =   .false. 
     64   luttrd(12) =   .false. 
     65   luttrd(13) =   .false. 
     66   luttrd(14) =   .false. 
     67   luttrd(15) =   .false. 
     68   luttrd(16) =   .false. 
     69   luttrd(17) =   .false. 
     70   luttrd(18) =   .false. 
     71   luttrd(19) =   .false. 
     72   luttrd(20) =   .false. 
     73   luttrd(21) =   .false. 
     74   luttrd(22) =   .false. 
     75   luttrd(23) =   .false. 
     76   luttrd(24) =   .false. 
    8377/ 
    8478!----------------------------------------------------------------------- 
    85 !       natdia  additional 2D/3D  (#ifdef key_trc_diaadd) 
     79&namtopadv    !   advection scheme for passive tracer  
    8680!----------------------------------------------------------------------- 
    87 !  nwritedia : time step frequency for additional arrays outputs 
    88 !  2D/3D diagnostic type defined by : 
    89 !                 *  short name  
    90 !                 *  long_name 
    91 !                 *  units 
    92 !                 *  logical to save value or not  
    93 !  
    94 &natdia 
    95    nwritedia   = 6000 
    96    diag2d(1)   = 'Cflx     ' , 'DIC flux                          ',  'molC/m2/s    ' 
    97    diag2d(2)   = 'Oflx     ' , 'Oxygen flux                       ',  'molC/m2/s    ' 
    98    diag2d(3)   = 'Kg       ' , 'Gas transfer                      ',  'mol/m2/s/uatm' 
    99    diag2d(4)   = 'Delc     ' , 'Delta CO2                         ',  'uatm         ' 
    100    diag2d(5)   = 'PMO      ' , 'POC export                        ',  'molC/m2/s    ' 
    101    diag2d(6)   = 'PMO2     ' , 'GOC export                        ',  'molC/m2/s    ' 
    102    diag2d(7)   = 'ExpFe1   ' , 'Nano iron export                  ',  'molFe/m2/s   ' 
    103    diag2d(8)   = 'ExpFe2   ' , 'Diatoms iron export               ',  'molFe/m2/s   ' 
    104    diag2d(9)   = 'ExpSi    ' , 'Silicate export                   ',  'molSi/m2/s   ' 
    105    diag2d(10)  = 'ExpCaCO3 ' , 'Calcite export                    ',  'molC/m2/s    ' 
    106    diag2d(11)  = 'heup     ' , 'euphotic layer depth              ',  'm            ' 
    107    diag2d(12)  = 'Fedep    ' , 'Iron dep                          ',  'molFe/m2/s   ' 
    108    diag2d(13)  = 'Nfix     ' , 'Nitrogen Fixation                 ',  'molN/m2/s    ' 
    109    diag3d(1)   = 'PH       ' , 'PH                                ',  '-            ' 
    110    diag3d(2)   = 'CO3      ' , 'Bicarbonates                      ',  'mol/l        ' 
    111    diag3d(3)   = 'CO3sat   ' , 'CO3 saturation                    ',  'mol/l        ' 
    112    diag3d(4)   = 'PAR      ' , 'light penetration                 ',  'W/m2         ' 
    113    diag3d(5)   = 'PPPHY    ' , 'Primary production of nanophyto   ',  'molC/m3/s    ' 
    114    diag3d(6)   = 'PPPHY2   ' , 'Primary production of diatoms     ',  'molC/m3/s    ' 
    115    diag3d(7)   = 'PPZOO    ' , 'Primary production of microzoo    ',  'molC/m3/s    ' 
    116    diag3d(8)   = 'PPZOO2   ' , 'Primary production of mesozoo     ',  'molC/m3/s    ' 
    117    diag3d(9)   = 'PBSi     ' , 'Primary production of Si diatoms  ',  'molSi/m3/s   ' 
    118    diag3d(10)  = 'PFeN     ' , 'Primary production of nano iron   ',  'molFe/m3/s   ' 
    119    diag3d(11)  = 'PFeD     ' , 'Primary production of diatoms iron',  'molFe/m3/s   ' 
     81   ln_trcadv_cen2   =  .false.  !  2nd order centered scheme    
     82   ln_trcadv_tvd    =  .false.  !  TVD scheme    
     83   ln_trcadv_muscl  =  .true.   !  MUSCL scheme 
     84   ln_trcadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
     85   ln_trcadv_smolar =  .false.  !  SMOLAR scheme 
     86   rsc              =  1.       !  tuning coefficient for Smolar. scheme  
     87   ncortrc          =  1        !  number of corrective phases for Smolar. scheme  
     88   crosster         =  .false.  !  computes Smolar crossterms (T) or not (F) 
    12089/ 
    12190!----------------------------------------------------------------------- 
    122 !       natnum  numerical schemes  
     91&namtopbbl    !   bottom boundary layer scheme for passive tracer  
    12392!----------------------------------------------------------------------- 
    124 !  ndttrc time step frequency for passive tracers 
    125 !  lhdf logical if true computes horizontal diffusion 
    126 !  rsc tuning coefficient for Smolar.  
    127 !  rtrn truncation value for Smolar. 
    128 !  ncor number of corrective phases for Smolar.  
    129 !  crosster logical if true computes Smolar crossterms 
    130 &natnum 
    131    ndttrc   = 1 
    132    rsc      = 1. 
    133    rtrn     = 1.e-15 
    134    ncortrc  = 1 
    135    crosster = .false. 
     93   atrcbbl          = 1000.     !  lateral mixing coefficient in the bbl  [m2/s] 
    13694/ 
    13795!----------------------------------------------------------------------- 
    138 !       namtrcadv   advection scheme for tracer (option not control by CPP keys) 
     96&namtopldf    !   lateral diffusion scheme for passive tracer  
    13997!----------------------------------------------------------------------- 
    140 !  ln_trcadv_cen2     2nd order centered scheme    (default F) 
    141 !  ln_trcadv_tvd      TVD scheme                   (default F) 
    142 !  ln_trcadv_muscl    MUSCL scheme                 (default F) 
    143 !  ln_trcadv_muscl2   MUSCL2 scheme                (default F) 
    144 !  ln_trcadv_smolar   SMOLAR scheme                (default T) 
    145 &namtrcadv 
    146    ln_trcadv_cen2   =  .false. 
    147    ln_trcadv_tvd    =  .false. 
    148    ln_trcadv_muscl  =  .true. 
    149    ln_trcadv_muscl2 =  .false. 
    150    ln_trcadv_smolar =  .false. 
    151 / 
    152 ! 
    153 !----------------------------------------------------------------------- 
    154 !       namtrcbbl   bottom boundary layer scheme 
    155 !----------------------------------------------------------------------- 
    156 !  atrcbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s) 
    157 &namtrcbbl 
    158    atrcbbl = 1000. 
     98   ln_trcldf_diff   =  .true.   !  performs lateral diffusion (T) or not (F) 
     99!                               !  Type of the operator :  
     100   ln_trcldf_lap    =  .true.   !     laplacian operator        
     101   ln_trcldf_bilap  =  .false.  !     bilaplacian operator      
     102                                !  Direction of action  : 
     103   ln_trcldf_level  =  .false.  !     iso-level                 
     104   ln_trcldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
     105   ln_trcldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     106!                               !  Coefficient 
     107   ahtrc0      =  2000.         !     horizontal eddy diffusivity for tracers [m2/s] 
     108   ahtrb0      =     0.         !     background eddy diffusivity for ldf_iso [m2/s] 
     109   aeivtr0     =  2000.         !     eddy induced velocity coefficient [m2/s]    (require "key_trcldf_eiv") 
     110   trcrat      =  1.            !  ratio betweeen passive and active tracer diffusion coeff 
    159111/ 
    160112!----------------------------------------------------------------------- 
    161 !       namtrcldf   lateral diffusion scheme for tracer (option not control by CPP keys) 
     113&namtopzdf        !   vertical physics 
    162114!----------------------------------------------------------------------- 
    163 !  Flag to performs lateral diffusion or not : 
    164 !     ln_trcldf_diff 
    165 !  Type of the operator :  
    166 !     ln_trcldf_lap    laplacian operator          (default T) 
    167 !     ln_trcldf_bilap  bilaplacian operator        (default F) 
    168 !  Direction of action  :  
    169 !     ln_trcldf_level  iso-level                   (default F) 
    170 !     ln_trcldf_hor    horizontal (geopotential)   (default F)^** 
    171 !     ln_trcldf_iso    iso-neutral                 (default T)^* 
    172 ! ^* require key_ldfslp to compute the direction of the lateral diffusion 
    173 ! ^** require key_ldfslp in s-coordinate 
    174 !  ahtrb0     background eddy diffusivity for isopycnal diffusion (m2/s) 
    175 !  trcrat     ratio betweeen passive and active tracer diffusion coeff 
    176 !  ahtrc0     horizontal eddy diffus. for passive tracer 
    177 !  aeivtr0    eddy induced veloc. coef. for passive tracer 
    178 &namtrcldf 
    179    ln_trcldf_diff   =  .true. 
    180    ln_trcldf_lap    =  .true. 
    181    ln_trcldf_bilap  =  .false. 
    182    ln_trcldf_level  =  .false. 
    183    ln_trcldf_hor    =  .false. 
    184    ln_trcldf_iso    =  .true. 
    185    ahtrb0           =  0.  
    186    trcrat           =  1. 
    187    ahtrc0           =  2000 
    188    aeivtr0          =  2000. 
     115   ln_trczdf_exp   =  .false.  !  split explicit (T) or implicit (F) time stepping 
     116   n_trczdf_exp    =   3       !  number of sub-timestep for ln_trczdfexp=T 
    189117/ 
    190118!----------------------------------------------------------------------- 
    191 !       namtrczdf   vertical physics 
     119&namtopdmp        !   passive tracer newtonian damping                 ('key_trcdmp') 
    192120!----------------------------------------------------------------------- 
    193 !  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F) 
    194 !                               (=F)  euler backward (F) 
    195 !  n_zdfexp   number of sub-timestep for time splitting scheme 
    196 &namtrczdf 
    197    ln_trczdf_exp =  .false. 
    198    n_trczdf_exp =      3 
     121   ndmptr      =   20      !  type of damping in passive tracers 
     122                           !     ='latitude', damping poleward of 'ndmp' degrees and function  
     123                           !                  of the distance-to-coast. Red and Med Seas as ndmptr=-1 
     124                           !     =-1 damping only in Med and Red Seas 
     125   ndmpftr     =    0      !  create a damping.coeff NetCDF file (=1) or not (=0) 
     126   nmldmptr    =    1      !  type of damping: =0 damping throughout the water column 
     127                           !                   =1 no damping in the mixed layer defined by avt >5cm2/s ) 
     128                           !                   =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
     129   sdmptr      =   50.     !  surface time scale for internal damping (days) 
     130   bdmptr      =  360.     !  bottom  time scale for internal damping (days) 
     131   hdmptr      =  800.     !  depth of transition between sdmptr and bdmptr (meters) 
    199132/ 
    200 !----------------------------------------------------------------------- 
    201 !       namtrcdmp   tracer newtonian damping ('key_trcdmp') 
    202 !----------------------------------------------------------------------- 
    203 !  ndmptr    type of damping in temperature and salinity  
    204 !          (='latitude', damping poleward of 'ndmp' degrees and function  
    205 !             of the distance-to-coast. Red and Med Seas as ndmp=-1) 
    206 !          (=-1 damping only in Med and Red Seas) 
    207 !  ndmpftr   =1 create a damping.coeff NetCDF file (the 3D damping array) 
    208 !  nmldmptr  type of damping in the mixed layer 
    209 !          (=0 damping throughout the water column) 
    210 !     (=1 no damping in the mixed layer defined by avt >5cm2/s ) 
    211 !     (=2 no damping in the mixed layer defined rho<rho(surf)+.01 ) 
    212 !  sdmptr    surface time scale for internal damping (days) 
    213 !  bdmptr    bottom time scale for internal damping (days) 
    214 !  hdmptr    depth of transition between sdmp and bdmp (meters) 
    215 &namtrcdmp 
    216    ndmptr   =  20 
    217    ndmpftr  =   0 
    218    nmldmptr =   1 
    219    sdmptr   =  50. 
    220    bdmptr   = 360. 
    221    hdmptr   = 800. 
    222 / 
Note: See TracChangeset for help on using the changeset viewer.