Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/1_namelist
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/1_namelist	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/1_namelist	(revision 966)
@@ -0,0 +1,887 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OPA  :       1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl)
+!! namelists         2 - Domain           (nam_zgr, nam_zgr_sco, namdom)
+!!                   3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core
+!!                                         namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb)
+!!                   4 - lateral boundary (namlbc, namcla, namobc)
+!!                   5 - bottom  boundary (nambfr, nambbc, nambbl)
+!!                   6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp)
+!!                   7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)
+!!                   8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm)
+!!                   9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr)
+!!                   9 - miscellaneous    (namsol)
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!  CAUTION: some scripts does not support CAPITALs for logical definition 
+!  *******       use .true. or .false. and NOT .TRUE. or .FALSE.
+
+!!======================================================================
+!!                   ***  Run management namelists  ***
+!!======================================================================
+!!   namrun            parameters of the run
+!!   nam_mpp           Massively Parallel Processing
+!!   nam_mpp_dyndist   
+!!   nam_ctl           Control prints & Benchmark
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namrun   parameters of the run
+!-----------------------------------------------------------------------
+!  no         job number
+!  cexper     experience name for vairmer format
+!  ln_rstart  boolean term for restart (true or false)
+!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file.
+!             !               = 1 restart, control nit000 in the restart file. Do not
+!             !                   use the date in the restart file (use ndate0 in namelist)
+!             !               = 2 restart, control nit000 in the restart file, use the date
+!             !                   in the restart file. ndate0 in the namelist is ignored.
+!  nit000     number of the first time step
+!  nitend     number of the last time step
+!  ndate0     initial calendar date aammjj
+!  nleapy     Leap year calendar (0/1)
+!  ninist     initial state output flag (0/1)
+!  nstock     frequency of restart file
+!  nwrite     frequency of OUTPUT file
+!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc)
+&namrun
+   no         =       0
+   cexper     =  "Agulhas"
+   ln_rstart  = .false.
+   nrstdt     =       0
+   nit000     =       1
+   nitend     =   10950
+   ndate0     =  010101
+   nleapy     =       0
+   ninist     =       0
+   nstock     =   10950 
+   nwrite     =   10950 
+   ln_dimgnnn = .false.
+/
+!-----------------------------------------------------------------------
+!       nam_mpp      Massively Parallel Processing
+!-----------------------------------------------------------------------
+!  c_mpi_send         mpi send/recieve type
+!                      = 'S'  : standard blocking send
+!                      = 'B'  : buffer blocking send
+!                      = 'I'  : immediate non-blocking send
+&nam_mpp
+   c_mpi_send =  'S'
+/
+!-----------------------------------------------------------------------
+!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist")
+!-----------------------------------------------------------------------
+!  jpni   number of processors following i
+!  jpnj   number of processors following j
+!  jpnij  number of local domains
+&nam_mpp_dyndist
+   jpni  = 1
+   jpnj  = 1
+   jpnij = 1
+/
+!-----------------------------------------------------------------------
+!       nam_ctl      Control prints & Benchmark
+!-----------------------------------------------------------------------
+!  ln_ctl     trends control print (expensive!)
+!  nprint     level of print (0 no print)
+!  nictls     start i indice to make the control SUM (very usefull to compare mono-
+!  nictle     end   i indice to make the control SUM (-versus multi processor runs)
+!  njctls     start j indice to make the control SUM (very usefull to compare mono-
+!  njctle     end   j indice to make the control SUM (-versus multi processor runs)
+!  nisplt     number of processors following i
+!  njsplt     number of processors following j
+!  nbench     Bench parameter (0/1): CAUTION it must be zero except for bench
+!             for which we don't care about physical meaning of the results
+!  nbit_cmp   bit comparison mode parameter (0/1): enables bit comparison between 
+!             single and multiple processor runs.
+&namctl
+   ln_ctl =  .false.
+   nprint =       0
+   nictls =       0
+   nictle =       0
+   njctls =       0
+   njctle =       0
+   isplt  =       1
+   jsplt  =       1
+   nbench =       0
+   nbit_cmp =     0
+/
+
+!!======================================================================
+!!                      ***  Domain namelists  ***
+!!======================================================================
+!!   nam_zgr       vertical coordinate
+!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
+!!   namdom        space and time domain (bathymetry, mesh, timestep)
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       nam_zgr       vertical coordinate
+!-----------------------------------------------------------------------
+!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined)
+!  ln_zps     z-coordinate - partial steps   (T/F)
+!  ln_sco     s- or hybrid z-s-coordinate    (T/F)
+&nam_zgr
+   ln_zco   =  .false.
+   ln_zps   =  .true.
+   ln_sco   =  .false.
+/
+!-----------------------------------------------------------------------
+!       nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
+!-----------------------------------------------------------------------
+!  sbot_min   minimum depth of s-bottom surface (>0) (m)
+!  sbot_max   maximum depth of s-bottom surface (= ocean depth) (>0) (m)
+!  theta      surface control parameter (0<=theta<=20)
+!  thetb      bottom control parameter  (0<=thetb<= 1)
+!  r_max      maximum cut-off r-value allowed (0<r_max<1)
+&nam_zgr_sco
+   sbot_min =  300.
+   sbot_max = 5250.
+   theta    =    6.0
+   thetb    =    0.75
+   r_max    =    0.15
+/
+!-----------------------------------------------------------------------
+!       namdom   space and time domain (bathymetry, mesh, timestep)
+!-----------------------------------------------------------------------
+!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels)
+!  e3zps_min  the thickness of the partial step is set larger than the
+!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1)
+!  nmsh       =1 create a mesh file (coordinates, scale factors, masks)
+!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra
+!                                                  = 1,     used, rdt < rdttra(k)
+!  atfp       asselin time filter parameter
+!  rdt        time step for the dynamics (and tracer if nacc=0)
+!  rdtmin     minimum time step on tracers
+!  rdtmax     maximum time step on tracers
+!  rdth       depth variation of tracer time step
+!  rdtbt      barotropic time step (for the time splitting algorithm)
+!  nclosea    =0 remove the closed sea from the global domain (orca configuration)
+!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 
+&namdom
+   ntopo     =     1
+   e3zps_min =     5.
+   e3zps_rat =     0.1
+   nmsh      =     0
+   nacc      =     0
+   atfp      =     0.1
+   rdt       =  2880.
+   rdtmin    =  2880.
+   rdtmax    =  2880.
+   rdth      =   800.
+   rdtbt     =    90.
+   nclosea   =     0
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***
+!!======================================================================
+!!   namsbc        surface boundary condition
+!!   namsbc_ana    analytical         formulation
+!!   namsbc_flx    flux               formulation
+!!   namsbc_clio   CLIO bulk formulea formulation
+!!   namsbc_core   CORE bulk formulea formulation
+!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled")
+!!   namqsr        penetrative solar radiation
+!!   namsbc_rnf    river runoffs
+!!   namsbc_ssr    sea surface restoring term (for T and/or S)
+!!   namalb        albedo parameters
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namsbc       Surface Boundary Condition (surface module)
+!-----------------------------------------------------------------------
+!  nn_fsbc       frequency of surface boundary condition computation 
+!                (= the frequency of sea-ice model call)
+!  ln_ana        analytical formulation (fill namsbc_ana ) 
+!  ln_flx        flux formulation       (fill namsbc_flx ) 
+!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core) 
+!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio) 
+!  ln_cpl        Coupled formulation    (fill namsbc_cpl ) 
+!  nn_ice        =0 no ice boundary condition   ,
+!                =1 observed ice-cover          ,
+!                =2 ice-model used                                      ("key_lim3" or "key_lim2)
+!  nn_ico_cpl    ice-ocean coupling : =0 LIM-3 old case
+!                                     =1 stresses computed using now ocean velocity
+!                                     =2 combination of 0 and 1 cases
+!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr)
+!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr)
+!  nn_fwb       FreshWater Budget: =0 no control                              , 
+!                                  =1 annual global mean of e-p-r set to zero   ,
+!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step
+&namsbc
+   nn_fsbc     = 5
+   ln_ana      = .false.
+   ln_flx      = .false.
+   ln_blk_clio = .true.
+   ln_blk_core = .false.
+   ln_cpl      = .false.
+   nn_ice      = 2
+   nn_ico_cpl  = 0
+   ln_dm2dc    = .false.
+   ln_rnf      = .false.
+   ln_ssr      = .false.
+   nn_fwb      = 0
+/
+!-----------------------------------------------------------------------
+!       namsbc_ana   analytical surface boundary condition
+!-----------------------------------------------------------------------
+!  nn_tau000   gently increase the stress over the first ntau_rst time-steps
+!  rn_utau0    uniform value for the i-stress
+!  rn_vtau0    uniform value for the j-stress
+!  rn_q0       uniform value for the total heat flux
+!  rn_qsr0     uniform value for the solar radiation
+!  rn_emp0     uniform value for the freswater budget (E-P)
+&namsbc_ana
+   nn_tau000  =   0
+   rn_utau0   =   0.5
+   rn_vtau0   =   0.e0
+   rn_q0      =   0.e0
+   rn_qsr0    =   0.e0
+   rn_emp0    =   0.e0
+/
+!-----------------------------------------------------------------------
+!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux")
+!-----------------------------------------------------------------------
+!  cn_dir   directory in which the model is executed
+!
+!  THE ORDER OF THE FILES MATTER:
+!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
+!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
+&namsbc_flx
+   cn_dir       = './'
+   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0
+   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0
+   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0
+   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0
+   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0
+/      
+!-----------------------------------------------------------------------
+!       namsbc_clio  CLIO bulk formulea
+!-----------------------------------------------------------------------
+!  cn_dir   directory in which the model is executed
+!
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting !
+!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  !
+&namsbc_clio
+   cn_dir       = './'
+   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0    
+   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0    
+   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0      
+   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0     
+   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0    
+   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0      
+   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0       
+/
+!-----------------------------------------------------------------------
+!       namsbc_core  CORE bulk formulea
+!-----------------------------------------------------------------------
+!  cn_dir         directory in which the model is executed
+!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m
+!  alpha_precip   multiplication factor for precipitation (total & snow)
+!
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting !
+!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  !
+&namsbc_core
+   cn_dir       = './'
+   ln_2m        = .false.
+   alpha_precip = 1.
+   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0    
+   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0    
+   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0      
+   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0      
+   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0     
+   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0    
+   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0      
+   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0       
+/
+!-----------------------------------------------------------------------
+!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled")
+!-----------------------------------------------------------------------
+&namsbc_cpl
+/
+!-----------------------------------------------------------------------
+!       namqsr   penetrative solar radiation
+!-----------------------------------------------------------------------
+!  ln_traqsr : penetrative solar radiation (T) or not (F)     (Default=T)
+!  rabs       fraction of qsr associated with xsi1
+!  xsi1       first depth of extinction
+!  xsi2       second depth of extinction
+&namqsr
+   ln_traqsr = .true.
+   rabs     =   0.58
+   xsi1     =   0.35
+   xsi2     =   23.0
+/
+!-----------------------------------------------------------------------
+!       namsbc_rnf   runoffs namelist surface boundary condition
+!-----------------------------------------------------------------------
+!cn_dir      Root directory for location of ssr files
+!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv
+!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used
+!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s]
+!sn_rnf      informations about the runoff file to be read
+!sn_cnf      informations about the runoff mouth file to be read
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name        ! frequency (hours) ! variable   ! time interpolation !
+!               !                    !   (if <0  months) !   name     !     (logical)      !
+&namsbc_rnf
+   cn_dir     = './'
+   nn_runoff  = 2
+   rn_hrnf    = 0.e0
+   rn_avt_rnf = 1.e-3
+   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0     
+   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0              
+/
+!-----------------------------------------------------------------------
+!       namsbc_ssr   surface boundary condition : sea surface restoring
+!-----------------------------------------------------------------------
+!  cn_dir    directory in which the model is executed
+!  nn_sstr   =0/1   add a retroaction term in the surface heat flux
+!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux
+!             (=1, salt flux, concentration/dillution emps only)
+!             (=2, volume flux, both emp and emps are updated) 
+!  dqdt       magnitude of the retroaction on temperature [W/m2/K]
+!  deds       magnitude of the damping on salinity [mm/day]
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
+!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
+&namsbc_ssr
+   cn_dir       = './'
+   nn_sstr      = 0
+   nn_sssr      = 0
+   dqdt         = -40.
+   deds         =  27.7
+   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0
+   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0
+/      
+!-----------------------------------------------------------------------
+!       namalb   albedo parameters
+!-----------------------------------------------------------------------
+!  cgren    correction of the snow or ice albedo to take into account
+!  albice   albedo of melting ice in the arctic and antarctic
+!  alphd    coefficients for linear interpolation used to compute albedo
+!           between two extremes values (Pyane, 1972)
+!  alphc     "                                         "
+!  alphdi    "                                         "
+&namalb
+   cgren    =      0.06
+   albice   =      0.5
+   alphd    =      0.80
+   alphc    =      0.65
+   alphdi   =      0.72
+/
+
+!!======================================================================
+!!               ***  Lateral boundary condition  ***
+!!======================================================================
+!!   namlbc        lateral momentum boundary condition
+!!   namcla        cross land advection
+!!   namobc        open boundaries parameters                           ("key_obc")
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namlbc   lateral momentum boundary condition
+!-----------------------------------------------------------------------
+!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip
+!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip
+&namlbc
+   shlat  =      2.
+/
+!-----------------------------------------------------------------------
+!       namcla   cross land advection
+!-----------------------------------------------------------------------
+!  n_cla   advection between 2 ocean pts separates by land 
+&namcla
+   n_cla   = 0
+/
+!-----------------------------------------------------------------------
+!       namobc    open boundaries parameters (#ifdef key_obc)
+!-----------------------------------------------------------------------
+!  nobc_dta   = 0 the obc data are equal to the initial state
+!             = 1 the obc data are read in 'obc   .dta' files
+!  rdpeob  time relaxation (days) for the east open boundary
+!  rdpwob  time relaxation (days) for the west open boundary
+!  rdpnob  time relaxation (days) for the north open boundary
+!  rdpsob  time relaxation (days) for the south open boundary
+!  zbsic1  barotropic stream function on isolated coastline 1
+!  zbsic2  barotropic stream function on isolated coastline 2
+!  zbsic3  barotropic stream function on isolated coastline 3
+!  ln_obc_clim  climatological obc data files (default T)
+!  ln_vol_cst   total volume conserved
+&namobc
+    nobc_dta =    0
+    rdpein   =    1.
+    rdpwin   =    1.
+    rdpnin   =   30.
+    rdpsin   =    1.
+    rdpeob   = 1500.
+    rdpwob   =   15.
+    rdpnob   =  150.
+    rdpsob   =   15.
+    zbsic1   =  140.e+6
+    zbsic2   =    1.e+6
+    zbsic3   =    0.
+    ln_obc_clim = .true.
+    ln_vol_cst  = .false.
+/
+
+!!======================================================================
+!!                 ***  Bottom boundary condition  ***
+!!======================================================================
+!!   nambfr        bottom friction
+!!   nambbc        bottom temperature boundary condition                ("key_trabbc")
+!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv")
+!!======================================================================
+!-----------------------------------------------------------------------
+!       nambfr   bottom friction
+!-----------------------------------------------------------------------
+!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction
+!                                     = 3 : free slip  ;  = 2 : non linear friction
+!  bfri1    bottom drag coefficient (linear case)
+!  bfri2    bottom drag coefficient (non linear case)
+!  bfeb2    bottom turbulent kinetic energy  (m2/s2)
+&nambfr
+   nbotfr =       1
+   bfri1  =   4.e-4
+   bfri2  =   1.e-3
+   bfeb2  =  2.5e-3
+/
+!-----------------------------------------------------------------------
+!       nambbc   bottom temperature boundary condition
+!-----------------------------------------------------------------------
+!  ngeo_flux  = 0 no geothermal heat flux
+!             = 1 constant geothermal heat flux
+!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2)
+!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2) 
+&nambbc
+   ngeo_flux =  2
+   ngeo_flux_const = 86.4e-3
+/
+!-----------------------------------------------------------------------
+!       nambbl   bottom boundary layer scheme
+!-----------------------------------------------------------------------
+!                                !  diffusive bbl                       ("key_trabbl")
+!                                !  advective bbl                       ("key_trabbl_adv")
+!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s)
+&nambbl
+   atrbbl = 10000.
+/
+
+!!======================================================================
+!!                        Tracer (T & S ) namelists
+!!======================================================================
+!!   nameos        equation of state
+!!   nam_traadv    advection scheme
+!!   nam_traldf    lateral diffusion scheme
+!!   namtdp        tracer newtonian damping                             ("key_tradmp")
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       nameos   ocean physical parameters
+!-----------------------------------------------------------------------
+!  neos    type of equation of state and Brunt-Vaisala frequency
+!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
+!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
+!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
+!  ralpha  thermal expension coefficient (neos= 1 or 2)
+!  rbeta   saline  expension coefficient (neos= 2)
+&nameos
+   neos   =      0
+   ralpha =  2.e-4
+   rbeta  =  0.001
+/
+!-----------------------------------------------------------------------
+!       nam_traadv   advection scheme for tracer
+!-----------------------------------------------------------------------
+!  ln_traadv_cen2     2nd order centered scheme          (default T)
+!  ln_traadv_tvd      TVD scheme                         (default F)
+!  ln_traadv_muscl    MUSCL scheme                       (default F)
+!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F)
+!  ln_traadv_ubs      UBS scheme                         (default F)
+!  ln_traadv_qck      QUICKEST scheme                    (default F)
+&nam_traadv
+   ln_traadv_cen2   =  .false.
+   ln_traadv_tvd    =  .true.
+   ln_traadv_muscl  =  .false.
+   ln_traadv_muscl2 =  .false.
+   ln_traadv_ubs    =  .false.
+   ln_traadv_qck    =  .false.
+/
+!-----------------------------------------------------------------------
+!       nam_traldf   lateral diffusion scheme for tracer
+!-----------------------------------------------------------------------
+!  Type of the operator : 
+!     ln_traldf_lap    laplacian operator          (default T)
+!     ln_traldf_bilap  bilaplacian operator        (default F)
+!  Direction of action  : 
+!     ln_traldf_level  iso-level                   (default F)
+!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T)
+!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
+!  Coefficient
+!     aht0    horizontal eddy diffusivity for tracers (m2/s)
+!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s)
+!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv")
+&nam_traldf
+   ln_traldf_lap    =  .true.
+   ln_traldf_bilap  =  .false.
+   ln_traldf_level  =  .false.
+   ln_traldf_hor    =  .false.
+   ln_traldf_iso    =  .true.
+   aht0    =  1000.
+   ahtb0   =     0.
+   aeiv0   =     0.
+/
+!-----------------------------------------------------------------------
+!       namtdp   tracer newtonian damping ('key_tradmp')
+!-----------------------------------------------------------------------
+!  ndmp    type of damping in temperature and salinity 
+!          !   ='latitude', damping poleward of 'ndmp' degrees and function 
+!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1
+!          !   =-1          damping only in Med and Red Seas
+!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array)
+!  nmldmp  type of damping in the mixed layer
+!          !   =0 damping throughout the water column
+!	        !   =1 no damping in the mixed layer defined by avt >5cm2/s 
+!	        !   =2 no damping in the mixed layer defined rho<rho(surf)+.01 
+!  sdmp    surface time scale for internal damping (days)
+!  bdmp    bottom time scale for internal damping (days)
+!  hdmp    depth of transition between sdmp and bdmp (meters)
+&namtdp
+   ndmp   =   90
+   ndmpf  =    1
+   nmldmp =    1
+   sdmp   =  50.
+   bdmp   = 360.
+   hdmp   = 800.
+/
+
+!!======================================================================
+!!                      ***  Dynamics namelists  ***
+!!======================================================================
+!!   nam_dynadv    formulation of the momentum advection
+!!   nam_dynvor    advection scheme
+!!   nam_dynhpg    hydrostatic pressure gradient
+!!   namflg        hydrostatic pressure gradient time stepping
+!!   nam_dynspg    surface pressure gradient                            (CPP key only)
+!!   nam_dynldf    lateral diffusion scheme
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       nam_dynadv   formulation of the momentum advection
+!-----------------------------------------------------------------------
+!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T)
+!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T)
+!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F)
+&nam_dynadv   
+   ln_dynadv_vec   = .true.
+   ln_dynadv_cen2  = .false.
+   ln_dynadv_ubs   = .false.
+/  
+!-----------------------------------------------------------------------
+!       nam_dynvor   option of physics/algorithm (not control by CPP keys)
+!-----------------------------------------------------------------------
+!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T)
+!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F)
+!  ln_dynvor_mix      "         "  : mixed scheme                (default F)
+!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F)
+&nam_dynvor
+   ln_dynvor_ene = .false.
+   ln_dynvor_ens = .false.
+   ln_dynvor_mix = .false.
+   ln_dynvor_een = .true.
+/
+!-----------------------------------------------------------------------
+!       nam_dynhpg   Hydrostatic pressure gradient option
+!-----------------------------------------------------------------------
+!  type of pressure gradient scheme (choose one only!)
+!     ln_hpg_zco    z-coordinate - full steps                   (default T)
+!     ln_hpg_zps    z-coordinate - partial steps (interpolation)
+!     ln_hpg_sco    s-coordinate (standard jacobian formulation)
+!     ln_hpg_hel    s-coordinate (helsinki modification)
+!     ln_hpg_wdj    s-coordinate (weighted density jacobian)
+!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial)
+!     ln_hpg_rot    s-coordinate (ROTated axes scheme)
+!  parameters
+!     gamm          weighting coefficient (wdj scheme)
+&nam_dynhpg
+   ln_hpg_zco = .false.
+   ln_hpg_zps = .true.
+   ln_hpg_sco = .false.
+   ln_hpg_hel = .false.
+   ln_hpg_wdj = .false.
+   ln_hpg_djc = .false.
+   ln_hpg_rot = .false.
+   gamm       = 0.e0
+/
+!-----------------------------------------------------------------------
+!       namflg   algorithm flags (algorithm not control by CPP keys)
+!-----------------------------------------------------------------------
+!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T)
+!                                                  centered     time scheme  (F)
+!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0)
+&namflg
+   ln_dynhpg_imp   =  .false.
+   nn_dynhpg_rst   =  0
+/
+!-----------------------------------------------------------------------
+!       nam_dynspg   surface pressure gradient   (CPP key only)
+!-----------------------------------------------------------------------
+!                                !  explicit free surface               ("key_dynspg_exp")
+!                                !  filtered free surface               ("key_dynspg_flt")
+!                                !  split-explicit free surface         ("key_dynspg_ts")
+!                                !  rigid-lid                           ("key_dynspg_rl")
+
+!-----------------------------------------------------------------------
+!       nam_dynldf   lateral diffusion on momentum
+!-----------------------------------------------------------------------
+!  Type of the operator : 
+!     ln_dynldf_lap    laplacian operator          (default T)
+!     ln_dynldf_bilap  bilaplacian operator        (default F)
+!  Direction of action  : 
+!     ln_dynldf_level  iso-level                   (default F)
+!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T)
+!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
+!  Coefficient
+!     ahm0    horizontal eddy viscosity for the dynamics (m2/s)
+!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s)
+&nam_dynldf
+   ln_dynldf_lap    =  .false.
+   ln_dynldf_bilap  =  .true.
+   ln_dynldf_level  =  .false.
+   ln_dynldf_hor    =  .true.
+   ln_dynldf_iso    =  .false.
+   ahm0    = -8.5e+11
+   ahmb0   =     0.
+/
+
+!!======================================================================
+!!             Tracers & Dynamics vertical physics namelists
+!!======================================================================
+!!       namzdf        vertical physics
+!!       namnpc        non penetrative convection                       
+!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      )
+!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      )
+!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      )
+!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      )
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namzdf   vertical physics
+!-----------------------------------------------------------------------
+!  ln_zdfevd  enhanced vertical diffusion         (default T)
+!  ln_zdfnpc  Non-Penetrative Convection          (default T)
+!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst")
+!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst")
+!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s)
+!  n_evdm     = 0  apply enhanced mixing on tracer only
+!             = 1  apply enhanced mixing on both tracer and momentum
+!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F)
+!                               (=F)  euler backward (F)
+!  n_zdfexp   number of sub-timestep for time splitting scheme
+&namzdf
+   ln_zdfevd = .true.
+   ln_zdfnpc = .false.
+   avm0      = 1.2e-4
+   avt0      = 1.2e-5
+   avevd     =   100.
+   n_evdm    =     1
+   ln_zdfexp =  .false.
+   n_zdfexp  =     3
+/
+!-----------------------------------------------------------------------
+!       namnpc   non penetrative convection
+!-----------------------------------------------------------------------
+!  nnpc1   non penetrative convective scheme frequency   
+!  nnpc2   non penetrative convective scheme print frequency
+&namnpc
+   nnpc1  =      1
+   nnpc2  =    365
+/
+!-----------------------------------------------------------------------
+!       namric   richardson number dependent vertical diffusion         ("key_zdfric" )
+!-----------------------------------------------------------------------
+!  avmri   maximum value of the vertical viscosity
+!  alp     coefficient of the parameterization
+!  nric    coefficient of the parameterization
+&namric
+   avmri = 100.e-4
+   alp   =      5.
+   nric  =       2
+/
+!-----------------------------------------------------------------------
+!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke")
+!-----------------------------------------------------------------------
+!  ln_rstke flag to restart with tke from a run without tke (default F)
+!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) )
+!  ediss    coef. of the Kolmogoroff dissipation  
+!  ebb      coef. of the surface input of tke
+!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm )
+!  emin     minimum value of tke (m^2/s^2)
+!  emin0    surface minimum value of tke (m^2/s^2)
+!  nitke    number of restart iterative loops
+!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded)
+!  nmxl     length used   = 0 bounded by the distance to surface and bottom
+!           !             = 1 bounded by the local vertical scale factor
+!           !             = 2 first vertical derivative of mixing length bounded by 1
+!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm)
+!           !                = 1 prandtl number function of richarson number (avt=pdl*avm)
+!           !                = 2 same as = 1 but a shapiro filter is applied on pdl
+!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1)
+!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb
+&namtke
+   ln_rstke = .false.
+   ediff =       0.1
+   ediss =       0.7
+   ebb   =       60.
+   efave =        1.
+   emin  =     1.e-6
+   emin0 =     1.e-4
+   nitke =        50
+   nmxl  =         2
+   npdl  =         1
+   navb  =         0
+/
+!-----------------------------------------------------------------------
+!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" )
+!-----------------------------------------------------------------------
+!  ln_kpprimix   shear instability mixing  (default T)
+!  difmiw        constant internal wave viscosity (m2/s)
+!  difsiw        constant internal wave diffusivity (m2/s)
+!  Riinfty       local Richardson Number limit for shear instability
+!  difri         maximum shear mixing at Rig = 0    (m2/s)
+!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection
+!  difcon        maximum mixing in interior convection (m2/s)
+!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv
+!  navb          = 0/1 flag for constant or profile background avt
+&namkpp
+   ln_kpprimix  = .true.
+   difmiw       =  1.e-04
+   difsiw       =  0.1e-04
+   Riinfty      =  0.8
+   difri        =  0.0050
+   bvsqcon      = -0.01e-07
+   difcon       =  1.
+   navb         =  0
+   nave         =  1
+/
+!-----------------------------------------------------------------------
+!       namddm   double diffusive mixing parameterization               ("key_zdfddm")
+!-----------------------------------------------------------------------
+!   avts    maximum avs for dd mixing 
+!   hsbfr   heat/salt buoyancy flux ratio
+&namddm
+      avts  = 1.e-4
+      hsbfr = 1.6
+/
+
+!!======================================================================
+!!                  ***  Miscelaneous namelists  ***
+!!======================================================================
+!!   namsol      elliptic solver / island / free surface 
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namsol   elliptic solver / island / free surface 
+!-----------------------------------------------------------------------
+!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg)
+!                            (=2 successive-over-relaxation: sor)
+!                            (=3 FETI currently it does not work!       ("key_feti")
+!  nsol_arp  absolute/relative (0/1) precision convergence test
+!  nmin      minimum of iterations for the SOR solver
+!  nmax      maximum of iterations for the SOR solver
+!  nmod      frequency of test for the SOR solver
+!  eps       absolute precision of the solver
+!  resmax    absolute precision for the SOR solver
+!  sor       optimal coefficient for SOR solver
+!  epsisl    absolute precision on stream function solver
+!  nmisl     maximum pcg iterations for island
+!  rnu       strength of the additional force used in free surface b.c.
+&namsol
+   nsolv     =      1
+   nsol_arp  =      0
+   nmin      =    300
+   nmax      =    800
+   nmod      =     10
+   eps       =  1.E-6
+   resmax    = 1.E-10
+   sor       =   1.92
+   epsisl    = 1.e-10
+   nmisl     =   4000
+   rnu       =     1.
+/
+!!======================================================================
+!!                       ***  Diagnostics namelists  ***
+!!======================================================================
+!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
+!!   namgap       level mean model-data gap                             ("key_diagap")
+!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr")
+!!   namflo       float parameters                                      ("key_float")
+!!   namptr       Poleward Transport Diagnostics
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra")
+!                          or mixed-layer trends                        ('key_trdmld')
+!                          or barotropic vorticity                      ("key_trdvor")
+!-----------------------------------------------------------------------
+!  ntrd                time step frequency dynamics and tracers trends
+!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk)
+!  ln_trdmld_restart   restart for ML diagnostics
+!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day)
+!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S
+&namtrd
+   ntrd              = 365
+   nctls             =   0
+   ln_trdmld_restart = .false.
+   ucf               =  1.
+   ln_trdmld_instant = .false.
+/
+!-----------------------------------------------------------------------
+!       namgap    level mean model-data gap ('key_diagap')
+!-----------------------------------------------------------------------
+!  ngap    time-step frequency of model-data gap computation
+!  nprg    time-step frequency of gap print in model output
+&namgap
+   ngap =  15
+   nprg =  10
+/
+!-----------------------------------------------------------------------
+!       namspr  surface pressure diagnostic
+!-----------------------------------------------------------------------
+!  nmaxp   maximum of iterations for the solver
+!  epsp    absolute precision of the solver
+!  niterp  number of iteration done by the solver
+&namspr
+   nmaxp   =   1000
+   epsp    =  1.e-3
+   niterp  =    400
+/
+!-----------------------------------------------------------------------
+!       namflo    float parameters                                      ("key_float")
+!-----------------------------------------------------------------------
+!  ln_rstflo   boolean term for float restart (true or false)
+!  nwritefl   frequency of float output file 
+!  nstockfl   frequency of float restart file 
+!  ln_argo    Argo type floats (stay at the surface each 10 days)
+!  ln_flork4  = T trajectories computed with a 4th order Runge-Kutta
+!             = F  (default)   computed with Blanke' scheme
+&namflo
+    ln_rstflo = .false.
+    nwritefl  =      75
+    nstockfl  =    5475
+    ln_argo   = .false.
+    ln_flork4 = .false.
+/
+!-----------------------------------------------------------------------
+!       namptr   Poleward Transport Diagnostic
+!-----------------------------------------------------------------------
+!  ln_diaptr  logical flag for Poleward transport computation
+!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation
+!             need input basins mask file named "subbasins.nc"
+!  nf_ptr     Frequency of computation
+&namptr
+   ln_diaptr = .false.
+   ln_subbas = .false.
+   nf_ptr    =  15
+/
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/1_namelist_ice_lim2
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/1_namelist_ice_lim2	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/1_namelist_ice_lim2	(revision 966)
@@ -0,0 +1,186 @@
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+! ICE namelist :  ice model option and parameter input
+! -------------
+!
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicerun   Share parameters for dynamics/advection/thermo
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  ln_limdyn   : switch for ice dynamics (true) or not (false)
+!  ln_limdmp   : restoring Ice thickness and Fraction leads flag
+!  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere
+!  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0)
+!  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0)
+!
+&namicerun
+   ln_limdyn = .true. 
+   ln_limdmp = .FALSE. 
+   acrit     = 1.0e-06 , 1.0e-06
+   hsndif    = 0.0 
+   hicdif    = 0.0
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceini   parameters for ice initialisation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  ln_limini : Ice initialization state flag
+!  ttest     : threshold water temperature for initial sea ice
+!  hninn     : initial snow thickness in the north
+!  hginn     : initial ice thickness in the north
+!  alinn     : initial leads area in the north
+!  hnins     : initial snow thickness in the south
+!  hgins     : initial ice thickness in the south
+!  alins     : initial leads area in the south
+!
+&namiceini
+   ln_limini = .FALSE. 
+   ttest = 2.0
+   hninn = 0.5
+   hginn = 3.0
+   alinn = 0.05
+   hnins = 0.1
+   hgins = 1.0
+   alins = 0.1
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicedia   parameters for ice diagnostics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  fmtinf : format of the output values 
+!  nfrinf : number of variables written in one line
+!  ntmoy  : instantaneous values of ice evolution or averaging
+!  ninfo  : frequency of ouputs on file ice_evolu in case of averaging
+!
+&namicedia
+   fmtinf = '1PE13.5 '
+   nfrinf = 4
+   ntmoy  = 1
+   ninfo  = 1
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicedyn   parameters for ice dynamic
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  epsd   : tolerance parameter
+!  alpha  : coefficient for semi-implicit coriolis
+!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
+!  dm     : diffusion constant for dynamics.
+!  nbiter : number of sub-time steps for relaxation
+!  nbitdr : maximum number of iterations for relaxation
+!  om     : relaxation constant 
+!  resl   : maximum value for the residual of relaxation 
+!  cw     : drag coefficient for oceanic stress
+!  angvg  : turning angle for oceanic stress
+!  pstar  : first bulk-rheology parameter
+!  c_rhg  : second bulk-rhelogy parameter
+!  etamn  : minimun value for viscosity
+!  creepl : creep limit
+!  ecc    : eccentricity of the elliptical yield curve
+!  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s)
+!
+&namicedyn 
+   epsd   = 1.0e-20
+   alpha  = 0.5
+   dm     = 0.6e+03
+   nbiter = 1
+   nbitdr = 100
+   om     = 0.5
+   resl   = 5.0e-05
+   cw     = 5.0e-03
+   angvg  = 0.0
+   pstar  = 1.0e+04 
+   c_rhg  = 20.0
+   etamn  = 0.0e+07
+   creepl = 2.0e-08
+   ecc    = 2.0
+   ahi0   = 350.e0
+/
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicetrp   parameters for ice advection
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
+&namicetrp 
+   bound  = 0.
+/
+
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicethd   parameters for thermodynamic computation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  hmelt       : maximum melting at the bottom
+!  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere
+!                caution 1.0, 1.0 best value to be used!!! (gilles G.)
+!  hicmin      : ice thickness corr. to max. energy stored in brine pocket
+!  hiclim      : minimum ice thickness
+!  amax        : maximum lead fraction
+!  swiqst      : energy stored in brine pocket (=1) or not (=0)	
+!  sbeta       : numerical caracteritic of the scheme for diffusion in ice
+!               Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)
+!  parlat      : percentage of energy used for lateral ablation
+!  hakspl      : slope of distr. for Hakkinen-Mellor's lateral melting
+!  hibspl      : slope of distribution for Hibler's lateral melting
+!  exld        : exponent for leads-closure rate
+!  hakdif      : coefficient for diffusions of ice and snow
+!  thth         : threshold thickness for comp. of eq. thermal conductivity
+!  hnzst       : thickness of the surf. layer in temp. computation
+!  parsub      : switch for snow sublimation or not
+!  alphs       : coefficient for snow density when snow ice formation
+! 
+&namicethd
+   hmelt   = -0.15
+   hiccrit = 0.3 , 0.3
+   hicmin  = 0.2
+   hiclim  = 0.05
+   amax    = 0.999 
+   swiqst  = 1.
+   sbeta   = 1.
+   parlat  = 0.0
+   hakspl  = 0.5
+   hibspl  = 0.5
+   exld    = 2.0
+   hakdif  = 1.0
+   thth    = 0.2
+   hnzst   = 0.1
+   parsub  = 1.0
+   alphs   = 1.0
+/
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceout   parameters for outputs
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!   noumef      : number of fields
+!   field
+!        title  : title of the field
+!        name   : name of the field
+!        unit   : unit of the field
+!        nc     : switch for saving field ( = 1 ) or not ( = 0 )
+!        cmulti : multiplicative constant
+!        cadd   : additive constant
+!
+&namiceout
+   noumef   = 19
+   field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0
+   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0
+   field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ', 1, 1.0, 0.0
+   field_4  = 'Ice concentration                  ', 'ileadfra', '%       ', 1, -1.0, 1.0
+   field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ', 1, 1.0, -273.15
+   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0
+   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0
+   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0
+   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15
+   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0
+   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0
+   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0
+   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0
+   field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt', 1, 1.0, 0.0
+   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0
+   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0 
+   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0
+   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0
+   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0
+/      
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/AA_job
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/AA_job	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/AA_job	(revision 966)
@@ -0,0 +1,177 @@
+#-Q- sxnec ##################
+#-Q- sxnec ## SXNEC IDRIS
+#-Q- sxnec ##################
+#-Q- sxnec #@$-r :XXX:XXX:XXX:           # job name
+#-Q- sxnec #@$-o output_XXXXXXX.1        # outputfile name
+#-Q- sxnec #@$-me                        # sends a mail at the end of execution
+#-Q- sxnec #@$-eo                        # merge standard error and output
+#-Q- sxnec #@$-jo                        # print a short performance analysis
+#-Q- sxnec #@$-lT 3:00:00                # Max total CPU time
+#-Q- sxnec #@$-lt 2:59:00                # Max CPU time per process
+#-Q- sxnec #@$-lM 5Gb                    # Max memory
+#-Q- sxnec ##@$-c 4                     # Max number of simultaneous procs used (here=4)
+#-Q- sxnec ##@$-q multi                 # Use more than 1 proc
+#-Q- aix ##################
+#-Q- aix ## IBM   IDRIS  ##
+#-Q- aix ##################
+#-Q- aix #@job_name= :XXX:XXX:XXX:          
+#-Q- aix #@output= output_XXXXXXX.1          
+#-Q- aix #@error= error_XXXXXXX.1          
+#-Q- aix #@cpu_limit= 10800
+#-Q- aix #@stack_limit= 200Mb        
+#-Q- aix #@data_limit= 5Gb         
+#-Q- aix #@notification= complete
+#-Q- aix #@total_tasks= 1 
+#-Q- aix #@job_type= parallel
+#-Q- aix #@queue
+#-Q- sx8brodie ##################
+#-Q- sx8brodie ## SX8NEC IDRIS ##
+#-Q- sx8brodie ##################
+#-Q- sx8brodie #PBS -N :XXX:XXX:XXX:        # job name
+#-Q- sx8brodie #PBS -o output_XXXXXXX.1     # outputfile name
+#-Q- sx8brodie #PBS -j o                    # merge standard error and output
+#-Q- sx8brodie #PBS -S /usr/bin/ksh         # user shell
+#-Q- sx8brodie #PBS -l memsz_job=5gb      # Max memory
+#-Q- sx8brodie #PBS -l cputim_prc=3:01:00   # Max CPU time/prc
+#-Q- sx8brodie #PBS -l cputim_job=3:01:00   # Max CPU time/job
+#-Q- sx8brodie ##PBS -l cpunum_job=4
+#-Q- sx8brodie ##PBS -q multi
+#-Q- sx6nec ##################
+#-Q- sx6nec ## SX6NEC   CEA ##
+#-Q- sx6nec ##################
+#-Q- sx6nec #PBS -N :XXX:XXX:XXX:        # job name
+#-Q- sx6nec #PBS -o output_XXXXXXX.1     # outputfile name
+#-Q- sx6nec #PBS -j o                    # merge standard error and output
+#-Q- sx6nec #PBS -S /usr/bin/ksh         # user shell
+#-Q- sx6nec #PBS -l memsz_job=5gb      # Max memory
+#-Q- sx6nec #PBS -l cputim_job=3:00:00   # Max CPU time
+#-Q- sx6nec #PBS -q prod
+#-Q- sx6nec ##PBS -l cpunum_prc=4
+#-Q- sx6nec ##PBS -q parallel
+########################################################################
+##
+##      Script to run a short experiment when installing the
+##      ORCA2_LIM configuration 
+##
+##            Rachid Benshila, OPA System Team
+##
+##
+#########################################################################
+##-- print echo of commands before and after shell interpretation
+date
+set -vx
+
+#################################
+# Beginning of user modification
+#################################
+
+##-- define some directories : 
+#- Input file storing 
+#-Q- sxnec R_FORCING="gaya:/u/rech/eee/reee451/FORCINGS/ORCA_R2/NEMO"
+#-Q- sx8brodie R_FORCING="gaya:/u/rech/eee/reee451/FORCINGS/ORCA_R2/NEMO"
+#-Q- aix R_FORCING="gaya:/u/rech/eee/reee451/FORCINGS/ORCA_R2/NEMO"
+#-Q- sx6nec R_FORCING="/dmnfs10/cont003/p48ethe/FORCINGS/ORCA_R2/NEMO"
+#- Launching run repository
+#-Q- sxnec R_EXPER=${QSUB_WORKDIR:-${PWD}}
+#-Q- sx8brodie R_EXPER=${PBS_O_WORKDIR}
+#-Q- aix R_EXPER=${QSUB_WORKDIR:-${PWD}}
+#-Q- sx6nec R_EXPER=${PWD}
+#- modipsl repository
+R_EXEDIR=MMMMMMM     
+#- output files storing 
+R_SORTIE_OCE=AAAAAAA/ORCA2_LIM
+#- execution repository
+R_TMP=$TMPDIR         
+
+##-- define processus number : to change when running in mpp
+NB_PROC=1
+
+#################################
+# End of user modification
+#################################
+
+#--  define a fonction to get the input tarfile
+#-Q- sxnec Rapatrie () { [ -f $1/$2 ] || ( rcp ${R_FORCING}/$2 $1 ; cd $1; tar xvf $2 ;) }
+#-Q- sx8brodie Rapatrie () { [ -f $1/$2 ] || ( rcp ${R_FORCING}/$2 $1 ; cd $1; tar xvf $2 ;) }
+#-Q- aix Rapatrie () { [ -f $1/$2 ] || ( rcp ${R_FORCING}/$2 $1 ; cd $1; tar xvf $2 ;) }
+#-Q- sx6nec Rapatrie () { [ -f $1/$2 ] || ( rcp ${R_FORCING}/$2 $1 ; cd $1; tar xvf $2 ;) }
+
+[ -d ${R_TMP} ] || mkdir ${R_TMP}
+cd ${R_TMP}
+
+#--  get the executable
+#-Q- sxnec rcp ${LOGNAME}@rhodes:${R_EXEDIR}/bin/opa opa.xx
+#-Q- sx8brodie rcp ${LOGNAME}@brodie:${R_EXEDIR}/bin/opa opa.xx
+#-Q- aix cp ${R_EXEDIR}/bin/opa opa.xx
+#-Q- sx6nec cp ${R_EXEDIR}/bin/opa opa.xx
+chmod 777 opa.xx
+
+## --  Copy ancillary files
+#
+#- Namelist for the configuration
+cp ${R_EXPER}/namelist_ice namelist_ice 
+cp ${R_EXPER}/namelist namelist
+cp ${R_EXPER}/namelist.passivetrc namelist.passivetrc
+cp ${R_EXPER}/namelist.trc.sms namelist.trc.sms
+ls -alF
+
+#- Namelist for ocean and ice (agrif fine grid)
+cp ${R_EXPER}/1_namelist_ice 1_namelist_ice 
+cp ${R_EXPER}/1_namelist 1_namelist
+cp ${R_EXPER}/AGRIF_FixedGrids.in AGRIF_FixedGrids.in
+
+#- Files for the configuration and ocean dynamics
+Rapatrie ${R_TMP} ORCA2_LIM_nemo_v2.tar
+Rapatrie ${R_TMP} INPUTS_PISCES.tar
+
+ls -alF
+
+#-Q- sxnec export F_PROGINF=yes
+#-Q- sxnec export F_ERRCNT=5
+#-Q- sx8brodie export F_PROGINF=yes
+#-Q- sx8brodie export F_ERRCNT=5
+#-Q- sx6nec export F_PROGINF=yes
+#-Q- sx6nec export F_ERRCNT=5
+
+#-Q- sxnec #- To be use for a mpp run
+#-Q- sxnec #mpirun -np ${NB_PROC} ./opa.xx
+#-Q- sx8brodie #- To be use for a mpp run
+#-Q- sx8brodie #mpirun -np ${NB_PROC} ./opa.xx
+#-Q- sx6nec #- To be use for a mpp run
+#-Q- sx6nec #mpirun -np ${NB_PROC} ./opa.xx
+
+./opa.xx
+
+## -- Save output files
+ls -alF
+#- to copy on a storage machine
+#-Q- sxnec DPUT=mfput
+#-Q- sx8brodie DPUT=mfput
+#-Q- aix DPUT=mfput
+#-Q- sx6nec DPUT=putfer 
+#-Q- sx6nec [ -d ${R_SORTIE_OCE} ] || mkdir -p ${R_SORTIE_OCE}
+
+#- Save ocean output files
+$DPUT *ocean.output ${R_SORTIE_OCE}/.
+$DPUT *solver.stat ${R_SORTIE_OCE}/.
+[ -f *mpp.output* ] && $DPUT *mpp.output* ${R_SORTIE_OCE}/.
+for file_netcdf in *_grid_*.nc *_diad_T*.nc *_ptrc_T*.nc
+do
+$DPUT ${file_netcdf} ${R_SORTIE_OCE}/${file_netcdf}
+done
+#
+#- Save ice output files
+for file_netcdf in *icemod*
+do
+$DPUT ${file_netcdf} ${R_SORTIE_OCE}/${file_netcdf} 
+done
+$DPUT ice.evolu ${R_SORTIE_OCE}/ice.evolu 
+#
+#- Save restart files
+for restart in *restart*
+do
+$DPUT ${restart} ${R_SORTIE_OCE}/${restart} 
+done
+#
+#- Save ftrace file
+[ -f ftrace.out* ] && $DPUT ftrace.out* ${R_SORTIE_OCE}/.
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/AGRIF_FixedGrids.in
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/AGRIF_FixedGrids.in	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/AGRIF_FixedGrids.in	(revision 966)
@@ -0,0 +1,3 @@
+1
+125 170 36 63 4 4 2
+0
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/EMPave_old.dat
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/EMPave_old.dat	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/EMPave_old.dat	(revision 966)
@@ -0,0 +1,1 @@
+                               0  0.0000000000000000E+00  0.0000000000000000E+00
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist	(revision 966)
@@ -0,0 +1,887 @@
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!! NEMO/OPA  :  1 - run manager      (namrun, nam_ctl, nam_mpp, nam_mpp_dyndist, nam_ctl)
+!! namelists    2 - Domain           (nam_zgr, nam_zgr_sco, namdom)
+!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core
+!!                                    namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb)
+!!              4 - lateral boundary (namlbc, namcla, namobc)
+!!              5 - bottom  boundary (nambfr, nambbc, nambbl)
+!!              6 - Tracer           (nameos, nam_traadv, nam_traldf, namtdp)
+!!              7 - dynamics         (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf)
+!!              8 - Verical physics  (namzdf, namnpc, namric, namtke, namkpp, namddm)
+!!              9 - diagnostics      (namtrd, namgap, namspr, namflo, namptr)
+!!              9 - miscellaneous    (namsol)
+!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!  CAUTION: some scripts does not support CAPITALs for logical definition 
+!  *******       use .true. or .false. and NOT .TRUE. or .FALSE.
+
+!!======================================================================
+!!                   ***  Run management namelists  ***
+!!======================================================================
+!!   namrun            parameters of the run
+!!   nam_mpp           Massively Parallel Processing
+!!   nam_mpp_dyndist   
+!!   nam_ctl           Control prints & Benchmark
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namrun   parameters of the run
+!-----------------------------------------------------------------------
+!  no         job number
+!  cexper     experience name for vairmer format
+!  ln_rstart  boolean term for restart (true or false)
+!  nrstdt     restart control = 0 restart, do not control nit000 in the restart file.
+!             !               = 1 restart, control nit000 in the restart file. Do not
+!             !                   use the date in the restart file (use ndate0 in namelist)
+!             !               = 2 restart, control nit000 in the restart file, use the date
+!             !                   in the restart file. ndate0 in the namelist is ignored.
+!  nit000     number of the first time step
+!  nitend     number of the last time step
+!  ndate0     initial calendar date aammjj
+!  nleapy     Leap year calendar (0/1)
+!  ninist     initial state output flag (0/1)
+!  nstock     frequency of restart file
+!  nwrite     frequency of OUTPUT file
+!  ln_dimgnnn (F/T) 1 DIMG file - (for all proc/per proc)
+&namrun
+   no         =       0
+   cexper     =  "ORCA2"
+   ln_rstart  = .false.
+   nrstdt     =       0
+   nit000     =       1
+   nitend     =    5475
+   ndate0     =  010101
+   nleapy     =       0
+   ninist     =       0
+   nstock     =    5475
+   nwrite     =    5475
+   ln_dimgnnn = .false.
+/
+!-----------------------------------------------------------------------
+!       nam_mpp      Massively Parallel Processing
+!-----------------------------------------------------------------------
+!  c_mpi_send         mpi send/recieve type
+!                      = 'S'  : standard blocking send
+!                      = 'B'  : buffer blocking send
+!                      = 'I'  : immediate non-blocking send
+&nam_mpp
+   c_mpi_send =  'S'
+/
+!-----------------------------------------------------------------------
+!       nam_mpp_dyndist      Massively Parallel Distribution            ("key_agrif" && "key_mpp_dyndist")
+!-----------------------------------------------------------------------
+!  jpni   number of processors following i
+!  jpnj   number of processors following j
+!  jpnij  number of local domains
+&nam_mpp_dyndist
+   jpni  = 1
+   jpnj  = 1
+   jpnij = 1
+/
+!-----------------------------------------------------------------------
+!       nam_ctl      Control prints & Benchmark
+!-----------------------------------------------------------------------
+!  ln_ctl     trends control print (expensive!)
+!  nprint     level of print (0 no print)
+!  nictls     start i indice to make the control SUM (very usefull to compare mono-
+!  nictle     end   i indice to make the control SUM (-versus multi processor runs)
+!  njctls     start j indice to make the control SUM (very usefull to compare mono-
+!  njctle     end   j indice to make the control SUM (-versus multi processor runs)
+!  nisplt     number of processors following i
+!  njsplt     number of processors following j
+!  nbench     Bench parameter (0/1): CAUTION it must be zero except for bench
+!             for which we don't care about physical meaning of the results
+!  nbit_cmp   bit comparison mode parameter (0/1): enables bit comparison between 
+!             single and multiple processor runs.
+&namctl
+   ln_ctl =  .false.
+   nprint =       0
+   nictls =       0
+   nictle =       0
+   njctls =       0
+   njctle =       0
+   isplt  =       1
+   jsplt  =       1
+   nbench =       0
+   nbit_cmp =     0
+/
+
+!!======================================================================
+!!                      ***  Domain namelists  ***
+!!======================================================================
+!!   nam_zgr       vertical coordinate
+!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
+!!   namdom        space and time domain (bathymetry, mesh, timestep)
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       nam_zgr       vertical coordinate
+!-----------------------------------------------------------------------
+!  ln_zco     z-coordinate - full steps      (T/F)                      ("key_zco" may also be defined)
+!  ln_zps     z-coordinate - partial steps   (T/F)
+!  ln_sco     s- or hybrid z-s-coordinate    (T/F)
+&nam_zgr
+   ln_zco   =  .false.
+   ln_zps   =  .true.
+   ln_sco   =  .false.
+/
+!-----------------------------------------------------------------------
+!       nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
+!-----------------------------------------------------------------------
+!  sbot_min   minimum depth of s-bottom surface (>0) (m)
+!  sbot_max   maximum depth of s-bottom surface (= ocean depth) (>0) (m)
+!  theta      surface control parameter (0<=theta<=20)
+!  thetb      bottom control parameter  (0<=thetb<= 1)
+!  r_max      maximum cut-off r-value allowed (0<r_max<1)
+&nam_zgr_sco
+   sbot_min =  300.
+   sbot_max = 5250.
+   theta    =    6.0
+   thetb    =    0.75
+   r_max    =    0.15
+/
+!-----------------------------------------------------------------------
+!       namdom   space and time domain (bathymetry, mesh, timestep)
+!-----------------------------------------------------------------------
+!  ntopo      = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels)
+!  e3zps_min  the thickness of the partial step is set larger than the
+!  e3zps_rat     the minimum of e3zps_min and e3zps_rat * e3t (N.B. 0<e3zps_rat<1)
+!  nmsh       =1 create a mesh file (coordinates, scale factors, masks)
+!  nacc       acceleration of convergence method   = 0, not used, rdt = rdttra
+!                                                  = 1,     used, rdt < rdttra(k)
+!  atfp       asselin time filter parameter
+!  rdt        time step for the dynamics (and tracer if nacc=0)
+!  rdtmin     minimum time step on tracers
+!  rdtmax     maximum time step on tracers
+!  rdth       depth variation of tracer time step
+!  rdtbt      barotropic time step (for the time splitting algorithm)
+!  nclosea    =0 remove the closed sea from the global domain (orca configuration)
+!             =1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 
+&namdom
+   ntopo     =     1
+   e3zps_min =     5.
+   e3zps_rat =     0.1
+   nmsh      =     0
+   nacc      =     0
+   atfp      =     0.1
+   rdt       =  5760.
+   rdtmin    =  5760.
+   rdtmax    =  5760.
+   rdth      =   800.
+   rdtbt     =    90.
+   nclosea   =     0
+/
+
+!!======================================================================
+!!            ***  Surface Boundary Condition namelists  ***
+!!======================================================================
+!!   namsbc        surface boundary condition
+!!   namsbc_ana    analytical         formulation
+!!   namsbc_flx    flux               formulation
+!!   namsbc_clio   CLIO bulk formulea formulation
+!!   namsbc_core   CORE bulk formulea formulation
+!!   namsbc_cpl    CouPLed            formulation                       ("key_coupled")
+!!   namqsr        penetrative solar radiation
+!!   namsbc_rnf    river runoffs
+!!   namsbc_ssr    sea surface restoring term (for T and/or S)
+!!   namalb        albedo parameters
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namsbc       Surface Boundary Condition (surface module)
+!-----------------------------------------------------------------------
+!  nn_fsbc       frequency of surface boundary condition computation 
+!                (= the frequency of sea-ice model call)
+!  ln_ana        analytical formulation (fill namsbc_ana ) 
+!  ln_flx        flux formulation       (fill namsbc_flx ) 
+!  ln_blk_clio   CLIO bulk formulation  (fill namsbc_core) 
+!  ln_blk_core   CORE bulk formulation  (fill namsbc_clio) 
+!  ln_cpl        Coupled formulation    (fill namsbc_cpl ) 
+!  nn_ice        =0 no ice boundary condition   ,
+!                =1 observed ice-cover          ,
+!                =2 ice-model used                                      ("key_lim3" or "key_lim2)
+!  nn_ico_cpl    ice-ocean coupling : =0 LIM-3 old case 
+!                                     =1 stresses computed using now ocean velocity
+!                                     =2 combination of 0 and 1 cases
+!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr)
+!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr)
+!  nn_fwb       FreshWater Budget: =0 no control                              , 
+!                                  =1 annual global mean of e-p-r set to zero   ,
+!                                  =2 global mean of e-p-r set to zero at each nn_fsbc time step
+&namsbc
+   nn_fsbc     = 5
+   ln_ana      = .false.
+   ln_flx      = .false.
+   ln_blk_clio = .true.
+   ln_blk_core = .false.
+   ln_cpl      = .false.
+   nn_ice      = 2
+   nn_ico_cpl  = 0
+   ln_dm2dc    = .false.
+   ln_rnf      = .false.
+   ln_ssr      = .false.
+   nn_fwb      = 0
+/
+!-----------------------------------------------------------------------
+!       namsbc_ana   analytical surface boundary condition
+!-----------------------------------------------------------------------
+!  nn_tau000   gently increase the stress over the first ntau_rst time-steps
+!  rn_utau0    uniform value for the i-stress
+!  rn_vtau0    uniform value for the j-stress
+!  rn_q0       uniform value for the total heat flux
+!  rn_qsr0     uniform value for the solar radiation
+!  rn_emp0     uniform value for the freswater budget (E-P)
+&namsbc_ana
+   nn_tau000  =   0
+   rn_utau0   =   0.5
+   rn_vtau0   =   0.e0
+   rn_q0      =   0.e0
+   rn_qsr0    =   0.e0
+   rn_emp0    =   0.e0
+/
+!-----------------------------------------------------------------------
+!       namsbc_flx   surface boundary condition : flux formulation (#ifdef "key_sbc_flux")
+!-----------------------------------------------------------------------
+!  cn_dir   directory in which the model is executed
+!
+!  THE ORDER OF THE FILES MATTER:
+!               !      file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
+!               !                    !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
+&namsbc_flx
+   cn_dir       = './'
+   sn_utau      =      'utau.nc'     ,        24.        ,  'utau'  ,     .false.       ,    0   ,     0
+   sn_vtau      =      'vtau.nc'     ,        24.        ,  'vtau'  ,     .false.       ,    0   ,     0
+   sn_qtot      =      'qtot.nc'     ,        24.        ,  'qtot'  ,     .false.       ,    0   ,     0
+   sn_qsr       =      'qsr.nc'      ,        24.        ,  'qsr'   ,     .false.       ,    0   ,     0
+   sn_emp       =      'emp.nc'      ,        24.        ,  'emp'   ,     .false.       ,    0   ,     0
+/      
+!-----------------------------------------------------------------------
+!       namsbc_clio  CLIO bulk formulea
+!-----------------------------------------------------------------------
+!  cn_dir   directory in which the model is executed
+!
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting !
+!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  !
+&namsbc_clio
+   cn_dir       = './'
+   sn_utau      =    'taux_1m'    ,      -12.         ,  'sozotaux',    .false.        ,    1   ,    0    
+   sn_vtau      =    'tauy_1m'    ,      -12.         ,  'sometauy',    .false.        ,    1   ,    0    
+   sn_wndm      =    'flx'        ,      -12.         ,  'socliowi',    .false.        ,    1   ,    0      
+   sn_tair      =    'flx'        ,      -12.         ,  'socliot1',    .false.        ,    1   ,    0     
+   sn_humi      =    'flx'        ,      -12.         ,  'socliohu',    .false.        ,    1   ,    0    
+   sn_ccov      =    'flx'        ,      -12.         ,  'socliocl',    .false.        ,    1   ,    0      
+   sn_prec      =    'flx'        ,      -12.         ,  'socliopl',    .false.        ,    1   ,    0       
+/
+!-----------------------------------------------------------------------
+!       namsbc_core  CORE bulk formulea
+!-----------------------------------------------------------------------
+!  cn_dir         directory in which the model is executed
+!  ln_2m          logical flag to use air temp. and hum referenced at 2m instead 10m
+!  alpha_precip   multiplication factor for precipitation (total & snow)
+!
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name        ! frequency (hours) ! variable   ! time intepolation !  clim  ! starting !
+!               !                    !   (if <0  months) !   name     !    (logical)      !  (0/1) !  record  !
+&namsbc_core
+   cn_dir       = './'
+   ln_2m        = .false.
+   alpha_precip = 1.
+   sn_wndi      =    'u10'        ,       24.         ,  'U_10_MOD',    .false.        ,    1   ,    0    
+   sn_wndj      =    'v10'        ,       24.         ,  'V_10_MOD',    .false.        ,    1   ,    0    
+   sn_qsr       =    'rad'        ,       24.         ,  'SWDN_MOD',    .false.        ,    1   ,    0      
+   sn_qlw       =    'rad'        ,       24.         ,  'LWDN_MOD',    .false.        ,    1   ,    0      
+   sn_tair      =    't10'        ,       24.         ,  'T_10_MOD',    .false.        ,    1   ,    0     
+   sn_humi      =    'q10'        ,       24.         ,  'Q_10_MOD',    .false.        ,    1   ,    0    
+   sn_prec      =    'precip'     ,      -12.         ,  'RAIN'    ,    .false.        ,    1   ,    0      
+   sn_snow      =    'precip'     ,      -12.         ,  'SNOW'    ,    .false.        ,    1   ,    0       
+/
+!-----------------------------------------------------------------------
+!       namsbc_cpl    coupled ocean/atmosphere model                    ("key_coupled")
+!-----------------------------------------------------------------------
+&namsbc_cpl
+/
+!-----------------------------------------------------------------------
+!       namqsr   penetrative solar radiation
+!-----------------------------------------------------------------------
+!  ln_traqsr : penetrative solar radiation (T) or not (F)     (Default=T)
+!  rabs       fraction of qsr associated with xsi1
+!  xsi1       first depth of extinction
+!  xsi2       second depth of extinction
+&namqsr
+   ln_traqsr = .true.
+   rabs     =   0.58
+   xsi1     =   0.35
+   xsi2     =   23.0
+/
+!-----------------------------------------------------------------------
+!       namsbc_rnf   runoffs namelist surface boundary condition
+!-----------------------------------------------------------------------
+!cn_dir      Root directory for location of ssr files
+!nn_runoff   =0 no, 1 runoff, 2 runoff+river mouth ups adv
+!rn_hrnf     runoffs, depth over which enhanced vertical mixing is used
+!rn_avt_rnf  runoffs, value of the additional vertical mixing coef. [m2/s]
+!sn_rnf      informations about the runoff file to be read
+!sn_cnf      informations about the runoff mouth file to be read
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name        ! frequency (hours) ! variable   ! time interpolation !
+!               !                    !   (if <0  months) !   name     !     (logical)      !
+&namsbc_rnf
+   cn_dir     = './'
+   nn_runoff  = 2
+   rn_hrnf    = 0.e0
+   rn_avt_rnf = 1.e-3
+   sn_rnf     = 'runoff_1m_nomask.nc' ,      -12.        ,  'sorunoff',    .true.   ,  1  ,  0     
+   sn_cnf     = 'runoff_1m_nomask.nc' ,        0.        ,  'socoefr' ,    .false.  ,  1  ,  0              
+/
+!-----------------------------------------------------------------------
+!       namsbc_ssr   surface boundary condition : sea surface restoring
+!-----------------------------------------------------------------------
+!  cn_dir    directory in which the model is executed
+!  nn_sstr   =0/1   add a retroaction term in the surface heat flux
+!  nn_sssr   =O/1/2 add a damping term in the surface freshwater flux
+!             (=1, salt flux, concentration/dillution emps only)
+!             (=2, volume flux, both emp and emps are updated) 
+!  dqdt       magnitude of the retroaction on temperature [W/m2/K]
+!  deds       magnitude of the damping on salinity [mm/day]
+!  THE ORDER OF THE FILES MATTER:
+!               !   file name     ! frequency (hours) ! variable ! time intepolation !  clim  ! starting !
+!               !                 !   (if <0  months) !   name   !    (logical)      !  (0/1) !  record  !
+&namsbc_ssr
+   cn_dir       = './'
+   nn_sstr      = 0
+   nn_sssr      = 0
+   dqdt         = -40.
+   deds         =  27.7
+   sn_sst       =  'sst_data.nc'  ,         24.       ,  'sst'   ,     .false.       ,    0   ,     0
+   sn_sss       =  'sss_data.nc'  ,        -12.       ,  'sss'   ,     .true.        ,    0   ,     0
+/      
+!-----------------------------------------------------------------------
+!       namalb   albedo parameters
+!-----------------------------------------------------------------------
+!  cgren    correction of the snow or ice albedo to take into account
+!  albice   albedo of melting ice in the arctic and antarctic
+!  alphd    coefficients for linear interpolation used to compute albedo
+!           between two extremes values (Pyane, 1972)
+!  alphc     "                                         "
+!  alphdi    "                                         "
+&namalb
+   cgren    =      0.06
+   albice   =      0.5
+   alphd    =      0.80
+   alphc    =      0.65
+   alphdi   =      0.72
+/
+
+!!======================================================================
+!!               ***  Lateral boundary condition  ***
+!!======================================================================
+!!   namlbc        lateral momentum boundary condition
+!!   namcla        cross land advection
+!!   namobc        open boundaries parameters                           ("key_obc")
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namlbc   lateral momentum boundary condition
+!-----------------------------------------------------------------------
+!  shlat            shlat = 0 , free slip   ;   0 < shlat < 2 , partial slip
+!                   shlat = 2 , no slip     ;   2 < shlat     , strong slip
+&namlbc
+   shlat  =      2.
+/
+!-----------------------------------------------------------------------
+!       namcla   cross land advection
+!-----------------------------------------------------------------------
+!  n_cla   advection between 2 ocean pts separates by land 
+&namcla
+   n_cla   = 1
+/
+!-----------------------------------------------------------------------
+!       namobc    open boundaries parameters (#ifdef key_obc)
+!-----------------------------------------------------------------------
+!  nobc_dta   = 0 the obc data are equal to the initial state
+!             = 1 the obc data are read in 'obc   .dta' files
+!  rdpeob  time relaxation (days) for the east open boundary
+!  rdpwob  time relaxation (days) for the west open boundary
+!  rdpnob  time relaxation (days) for the north open boundary
+!  rdpsob  time relaxation (days) for the south open boundary
+!  zbsic1  barotropic stream function on isolated coastline 1
+!  zbsic2  barotropic stream function on isolated coastline 2
+!  zbsic3  barotropic stream function on isolated coastline 3
+!  ln_obc_clim  climatological obc data files (default T)
+!  ln_vol_cst   total volume conserved
+&namobc
+    nobc_dta =    0
+    rdpein   =    1.
+    rdpwin   =    1.
+    rdpnin   =   30.
+    rdpsin   =    1.
+    rdpeob   = 1500.
+    rdpwob   =   15.
+    rdpnob   =  150.
+    rdpsob   =   15.
+    zbsic1   =  140.e+6
+    zbsic2   =    1.e+6
+    zbsic3   =    0.
+    ln_obc_clim = .true.
+    ln_vol_cst  = .false.
+/
+
+!!======================================================================
+!!                 ***  Bottom boundary condition  ***
+!!======================================================================
+!!   nambfr        bottom friction
+!!   nambbc        bottom temperature boundary condition                ("key_trabbc")
+!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv")
+!!======================================================================
+!-----------------------------------------------------------------------
+!       nambfr   bottom friction
+!-----------------------------------------------------------------------
+!  nbotfr   type of bottom friction : = 0 : no slip    ;  = 1 : linear friction
+!                                     = 3 : free slip  ;  = 2 : non linear friction
+!  bfri1    bottom drag coefficient (linear case)
+!  bfri2    bottom drag coefficient (non linear case)
+!  bfeb2    bottom turbulent kinetic energy  (m2/s2)
+&nambfr
+   nbotfr =       1
+   bfri1  =   4.e-4
+   bfri2  =   1.e-3
+   bfeb2  =  2.5e-3
+/
+!-----------------------------------------------------------------------
+!       nambbc   bottom temperature boundary condition
+!-----------------------------------------------------------------------
+!  ngeo_flux  = 0 no geothermal heat flux
+!             = 1 constant geothermal heat flux
+!             = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2)
+!  ngeo_flux_const   Constant value of geothermal heat flux (W/m2) 
+&nambbc
+   ngeo_flux =  2
+   ngeo_flux_const = 86.4e-3
+/
+!-----------------------------------------------------------------------
+!       nambbl   bottom boundary layer scheme
+!-----------------------------------------------------------------------
+!                                !  diffusive bbl                       ("key_trabbl")
+!                                !  advective bbl                       ("key_trabbl_adv")
+!  atrbbl   mixing coefficient of the bottom boundary layer scheme (m2/s)
+&nambbl
+   atrbbl = 10000.
+/
+
+!!======================================================================
+!!                        Tracer (T & S ) namelists
+!!======================================================================
+!!   nameos        equation of state
+!!   nam_traadv    advection scheme
+!!   nam_traldf    lateral diffusion scheme
+!!   namtdp        tracer newtonian damping                             ("key_tradmp")
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       nameos   ocean physical parameters
+!-----------------------------------------------------------------------
+!  neos    type of equation of state and Brunt-Vaisala frequency
+!          !   = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
+!          !   = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
+!          !   = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
+!  ralpha  thermal expension coefficient (neos= 1 or 2)
+!  rbeta   saline  expension coefficient (neos= 2)
+&nameos
+   neos   =      0
+   ralpha =  2.e-4
+   rbeta  =  0.001
+/
+!-----------------------------------------------------------------------
+!       nam_traadv   advection scheme for tracer
+!-----------------------------------------------------------------------
+!  ln_traadv_cen2     2nd order centered scheme          (default T)
+!  ln_traadv_tvd      TVD scheme                         (default F)
+!  ln_traadv_muscl    MUSCL scheme                       (default F)
+!  ln_traadv_muscl2   MUSCL scheme + cen2 at boundaries  (default F)
+!  ln_traadv_ubs      UBS scheme                         (default F)
+!  ln_traadv_qck      QUICKEST scheme                    (default F)
+&nam_traadv
+   ln_traadv_cen2   =  .true.
+   ln_traadv_tvd    =  .false.
+   ln_traadv_muscl  =  .false.
+   ln_traadv_muscl2 =  .false.
+   ln_traadv_ubs    =  .false.
+   ln_traadv_qck    =  .false.
+/
+!-----------------------------------------------------------------------
+!       nam_traldf   lateral diffusion scheme for tracer
+!-----------------------------------------------------------------------
+!  Type of the operator : 
+!     ln_traldf_lap    laplacian operator          (default T)
+!     ln_traldf_bilap  bilaplacian operator        (default F)
+!  Direction of action  : 
+!     ln_traldf_level  iso-level                   (default F)
+!     ln_traldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T)
+!     ln_traldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
+!  Coefficient
+!     aht0    horizontal eddy diffusivity for tracers (m2/s)
+!     ahtb0   background eddy diffusivity for isopycnal diffusion (m2/s)
+!     aeiv0   eddy induced velocity coefficient (m2/s)           (require "key_traldf_eiv")
+&nam_traldf
+   ln_traldf_lap    =  .true.
+   ln_traldf_bilap  =  .false.
+   ln_traldf_level  =  .false.
+   ln_traldf_hor    =  .false.
+   ln_traldf_iso    =  .true.
+   aht0    =  2000.
+   ahtb0   =     0.
+   aeiv0   =  2000.
+/
+!-----------------------------------------------------------------------
+!       namtdp   tracer newtonian damping ('key_tradmp')
+!-----------------------------------------------------------------------
+!  ndmp    type of damping in temperature and salinity 
+!          !   ='latitude', damping poleward of 'ndmp' degrees and function 
+!          !                of the distance-to-coast. Red and Med Seas as ndmp=-1
+!          !   =-1          damping only in Med and Red Seas
+!  ndmpf   =1 create a damping.coeff NetCDF file (the 3D damping array)
+!  nmldmp  type of damping in the mixed layer
+!          !   =0 damping throughout the water column
+!	        !   =1 no damping in the mixed layer defined by avt >5cm2/s 
+!	        !   =2 no damping in the mixed layer defined rho<rho(surf)+.01 
+!  sdmp    surface time scale for internal damping (days)
+!  bdmp    bottom time scale for internal damping (days)
+!  hdmp    depth of transition between sdmp and bdmp (meters)
+&namtdp
+   ndmp   =   -1
+   ndmpf  =    1
+   nmldmp =    1
+   sdmp   =  50.
+   bdmp   = 360.
+   hdmp   = 800.
+/
+
+!!======================================================================
+!!                      ***  Dynamics namelists  ***
+!!======================================================================
+!!   nam_dynadv    formulation of the momentum advection
+!!   nam_dynvor    advection scheme
+!!   nam_dynhpg    hydrostatic pressure gradient
+!!   namflg        hydrostatic pressure gradient time stepping
+!!   nam_dynspg    surface pressure gradient                            (CPP key only)
+!!   nam_dynldf    lateral diffusion scheme
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       nam_dynadv   formulation of the momentum advection
+!-----------------------------------------------------------------------
+!  ln_dynadv_vec      vector form (T) or flux form (F)         (default T)
+!  ln_dynadv_cen2     flux form - 2nd order centered scheme    (default T)
+!  ln_dynadv_ubs      flux form - 3rd order UBS scheme         (default F)
+&nam_dynadv   
+   ln_dynadv_vec   = .true.
+   ln_dynadv_cen2  = .false.
+   ln_dynadv_ubs   = .false.
+/  
+!-----------------------------------------------------------------------
+!       nam_dynvor   option of physics/algorithm (not control by CPP keys)
+!-----------------------------------------------------------------------
+!  ln_dynvor_ens   vorticity trends: enstrophy conserving scheme (default T)
+!  ln_dynvor_ene      "         "  : energy conserving scheme    (default F)
+!  ln_dynvor_mix      "         "  : mixed scheme                (default F)
+!  ln_dynvor_een      "         "  : energy & enstrophy scheme   (default F)
+&nam_dynvor
+   ln_dynvor_ene = .false.
+   ln_dynvor_ens = .true.
+   ln_dynvor_mix = .false.
+   ln_dynvor_een = .false.
+/
+!-----------------------------------------------------------------------
+!       nam_dynhpg   Hydrostatic pressure gradient option
+!-----------------------------------------------------------------------
+!  type of pressure gradient scheme (choose one only!)
+!     ln_hpg_zco    z-coordinate - full steps                   (default T)
+!     ln_hpg_zps    z-coordinate - partial steps (interpolation)
+!     ln_hpg_sco    s-coordinate (standard jacobian formulation)
+!     ln_hpg_hel    s-coordinate (helsinki modification)
+!     ln_hpg_wdj    s-coordinate (weighted density jacobian)
+!     ln_hpg_djc    s-coordinate (Density Jacobian with Cubic polynomial)
+!     ln_hpg_rot    s-coordinate (ROTated axes scheme)
+!  parameters
+!     gamm          weighting coefficient (wdj scheme)
+&nam_dynhpg
+   ln_hpg_zco = .false.
+   ln_hpg_zps = .true.
+   ln_hpg_sco = .false.
+   ln_hpg_hel = .false.
+   ln_hpg_wdj = .false.
+   ln_hpg_djc = .false.
+   ln_hpg_rot = .false.
+   gamm       = 0.e0
+/
+!-----------------------------------------------------------------------
+!       namflg   algorithm flags (algorithm not control by CPP keys)
+!-----------------------------------------------------------------------
+!  ln_dynhpg_imp   hydrostatic pressure gradient: semi-implicit time scheme  (T)
+!                                                  centered     time scheme  (F)
+!   nn_dynhpg_rst  add dynhpg implicit variables in restart ot not (1/0)
+&namflg
+   ln_dynhpg_imp   =  .false.
+   nn_dynhpg_rst   =  0
+/
+!-----------------------------------------------------------------------
+!       nam_dynspg   surface pressure gradient   (CPP key only)
+!-----------------------------------------------------------------------
+!                                !  explicit free surface               ("key_dynspg_exp")
+!                                !  filtered free surface               ("key_dynspg_flt")
+!                                !  split-explicit free surface         ("key_dynspg_ts")
+!                                !  rigid-lid                           ("key_dynspg_rl")
+
+!-----------------------------------------------------------------------
+!       nam_dynldf   lateral diffusion on momentum
+!-----------------------------------------------------------------------
+!  Type of the operator : 
+!     ln_dynldf_lap    laplacian operator          (default T)
+!     ln_dynldf_bilap  bilaplacian operator        (default F)
+!  Direction of action  : 
+!     ln_dynldf_level  iso-level                   (default F)
+!     ln_dynldf_hor    horizontal (geopotential)   (default F)   (require "key_ldfslp" if ln_sco=T)
+!     ln_dynldf_iso    iso-neutral                 (default T)   (require "key_ldfslp")
+!  Coefficient
+!     ahm0    horizontal eddy viscosity for the dynamics (m2/s)
+!     ahmb0   background eddy viscosity for isopycnal diffusion (m2/s)
+&nam_dynldf
+   ln_dynldf_lap    =  .true.
+   ln_dynldf_bilap  =  .false.
+   ln_dynldf_level  =  .false.
+   ln_dynldf_hor    =  .true.
+   ln_dynldf_iso    =  .false.
+   ahm0    = 40000.
+   ahmb0   =     0.
+/
+
+!!======================================================================
+!!             Tracers & Dynamics vertical physics namelists
+!!======================================================================
+!!       namzdf        vertical physics
+!!       namnpc        non penetrative convection                       
+!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      )
+!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      )
+!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      )
+!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      )
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namzdf   vertical physics
+!-----------------------------------------------------------------------
+!  ln_zdfevd  enhanced vertical diffusion         (default T)
+!  ln_zdfnpc  Non-Penetrative Convection          (default T)
+!  avm0       vertical eddy viscosity for the dynamic (m2/s)   (also background Kz if not "key_zdfcst")
+!  avt0       vertical eddy diffusivity for tracers (m2/s)     (also background Kz if not "key_zdfcst")
+!  avevd      vertical coefficient for enhanced diffusion scheme (m2/s)
+!  n_evdm     = 0  apply enhanced mixing on tracer only
+!             = 1  apply enhanced mixing on both tracer and momentum
+!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F)
+!                               (=F)  euler backward (F)
+!  n_zdfexp   number of sub-timestep for time splitting scheme
+&namzdf
+   ln_zdfevd = .true.
+   ln_zdfnpc = .false.
+   avm0      = 1.2e-4
+   avt0      = 1.2e-5
+   avevd     =   100.
+   n_evdm    =     1
+   ln_zdfexp =  .false.
+   n_zdfexp  =     3
+/
+!-----------------------------------------------------------------------
+!       namnpc   non penetrative convection
+!-----------------------------------------------------------------------
+!  nnpc1   non penetrative convective scheme frequency   
+!  nnpc2   non penetrative convective scheme print frequency
+&namnpc
+   nnpc1  =      1
+   nnpc2  =    365
+/
+!-----------------------------------------------------------------------
+!       namric   richardson number dependent vertical diffusion         ("key_zdfric" )
+!-----------------------------------------------------------------------
+!  avmri   maximum value of the vertical viscosity
+!  alp     coefficient of the parameterization
+!  nric    coefficient of the parameterization
+&namric
+   avmri = 100.e-4
+   alp   =      5.
+   nric  =       2
+/
+!-----------------------------------------------------------------------
+!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke")
+!-----------------------------------------------------------------------
+!  ln_rstke flag to restart with tke from a run without tke (default F)
+!  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) )
+!  ediss    coef. of the Kolmogoroff dissipation  
+!  ebb      coef. of the surface input of tke
+!  efave    coef. to applied to the tke diffusion ( avtke=efave*avm )
+!  emin     minimum value of tke (m^2/s^2)
+!  emin0    surface minimum value of tke (m^2/s^2)
+!  nitke    number of restart iterative loops
+!  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded)
+!  nmxl     length used   = 0 bounded by the distance to surface and bottom
+!           !             = 1 bounded by the local vertical scale factor
+!           !             = 2 first vertical derivative of mixing length bounded by 1
+!  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm)
+!           !                = 1 prandtl number function of richarson number (avt=pdl*avm)
+!           !                = 2 same as = 1 but a shapiro filter is applied on pdl
+!  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1)
+!  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb
+&namtke
+   ln_rstke = .false.
+   ediff =       0.1
+   ediss =       0.7
+   ebb   =      3.75
+   efave =        1.
+   emin  =     1.e-6
+   emin0 =     1.e-4
+   nitke =        50
+   nmxl  =         2
+   npdl  =         1
+   navb  =         0
+/
+!-----------------------------------------------------------------------
+!       namkpp   K-Profile Parameterization dependent vertical mixing   ("key_zdfkpp" )
+!-----------------------------------------------------------------------
+!  ln_kpprimix   shear instability mixing  (default T)
+!  difmiw        constant internal wave viscosity (m2/s)
+!  difsiw        constant internal wave diffusivity (m2/s)
+!  Riinfty       local Richardson Number limit for shear instability
+!  difri         maximum shear mixing at Rig = 0    (m2/s)
+!  bvsqcon       Brunt-Vaisala squared (1/s**2) for maximum convection
+!  difcon        maximum mixing in interior convection (m2/s)
+!  nave          = 0/1 flag for horizontal average on avt, avmu, avmv
+!  navb          = 0/1 flag for constant or profile background avt
+&namkpp
+   ln_kpprimix  = .true.
+   difmiw       =  1.e-04
+   difsiw       =  0.1e-04
+   Riinfty      =  0.8
+   difri        =  0.0050
+   bvsqcon      = -0.01e-07
+   difcon       =  1.
+   navb         =  0
+   nave         =  1
+/
+!-----------------------------------------------------------------------
+!       namddm   double diffusive mixing parameterization               ("key_zdfddm")
+!-----------------------------------------------------------------------
+!   avts    maximum avs for dd mixing 
+!   hsbfr   heat/salt buoyancy flux ratio
+&namddm
+      avts  = 1.e-4
+      hsbfr = 1.6
+/
+
+!!======================================================================
+!!                  ***  Miscelaneous namelists  ***
+!!======================================================================
+!!   namsol      elliptic solver / island / free surface 
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namsol   elliptic solver / island / free surface 
+!-----------------------------------------------------------------------
+!  nsolv     elliptic solver (=1 preconditioned conjugate gradient: pcg)
+!                            (=2 successive-over-relaxation: sor)
+!                            (=3 FETI currently it does not work!       ("key_feti")
+!  nsol_arp  absolute/relative (0/1) precision convergence test
+!  nmin      minimum of iterations for the SOR solver
+!  nmax      maximum of iterations for the SOR solver
+!  nmod      frequency of test for the SOR solver
+!  eps       absolute precision of the solver
+!  resmax    absolute precision for the SOR solver
+!  sor       optimal coefficient for SOR solver
+!  epsisl    absolute precision on stream function solver
+!  nmisl     maximum pcg iterations for island
+!  rnu       strength of the additional force used in free surface b.c.
+&namsol
+   nsolv     =      1
+   nsol_arp  =      0
+   nmin      =    300
+   nmax      =    800
+   nmod      =     10
+   eps       =  1.E-6
+   resmax    = 1.E-10
+   sor       =   1.92
+   epsisl    = 1.e-10
+   nmisl     =   4000
+   rnu       =     1.
+/
+!!======================================================================
+!!                       ***  Diagnostics namelists  ***
+!!======================================================================
+!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
+!!   namgap       level mean model-data gap                             ("key_diagap")
+!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr")
+!!   namflo       float parameters                                      ("key_float")
+!!   namptr       Poleward Transport Diagnostics
+!!======================================================================
+
+!-----------------------------------------------------------------------
+!       namtrd    diagnostics on dynamics and/or tracer trends          ("key_trddyn" and/or "key_trdtra")
+!                          or mixed-layer trends                        ('key_trdmld')
+!                          or barotropic vorticity                      ("key_trdvor")
+!-----------------------------------------------------------------------
+!  ntrd                time step frequency dynamics and tracers trends
+!  nctls               control surface type in mixed-layer trends (0,1 or n<jpk)
+!  ln_trdmld_restart   restart for ML diagnostics
+!  ucf                 unit conversion factor (=1 -> /seconds | =86400. -> /day)
+!  ln_trdmld_instant   flag to diagnose trends of instantantaneous or mean ML T/S
+&namtrd
+   ntrd              = 365
+   nctls             =   0
+   ln_trdmld_restart = .false.
+   ucf               =  1.
+   ln_trdmld_instant = .false.
+/
+!-----------------------------------------------------------------------
+!       namgap    level mean model-data gap ('key_diagap')
+!-----------------------------------------------------------------------
+!  ngap    time-step frequency of model-data gap computation
+!  nprg    time-step frequency of gap print in model output
+&namgap
+   ngap =  15
+   nprg =  10
+/
+!-----------------------------------------------------------------------
+!       namspr  surface pressure diagnostic
+!-----------------------------------------------------------------------
+!  nmaxp   maximum of iterations for the solver
+!  epsp    absolute precision of the solver
+!  niterp  number of iteration done by the solver
+&namspr
+   nmaxp   =   1000
+   epsp    =  1.e-3
+   niterp  =    400
+/
+!-----------------------------------------------------------------------
+!       namflo    float parameters                                      ("key_float")
+!-----------------------------------------------------------------------
+!  ln_rstflo   boolean term for float restart (true or false)
+!  nwritefl   frequency of float output file 
+!  nstockfl   frequency of float restart file 
+!  ln_argo    Argo type floats (stay at the surface each 10 days)
+!  ln_flork4  = T trajectories computed with a 4th order Runge-Kutta
+!             = F  (default)   computed with Blanke' scheme
+&namflo
+    ln_rstflo = .false.
+    nwritefl  =      75
+    nstockfl  =    5475
+    ln_argo   = .false.
+    ln_flork4 = .false.
+/
+!-----------------------------------------------------------------------
+!       namptr   Poleward Transport Diagnostic
+!-----------------------------------------------------------------------
+!  ln_diaptr  logical flag for Poleward transport computation
+!  ln_subbas  logical flag for Atlantic/Pacific/Indian basins computation
+!             need input basins mask file named "subbasins.nc"
+!  nf_ptr     Frequency of computation
+&namptr
+   ln_diaptr = .false.
+   ln_subbas = .false.
+   nf_ptr    =  15
+/
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist.passivetrc
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist.passivetrc	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist.passivetrc	(revision 966)
@@ -0,0 +1,222 @@
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+! OPA MODEL general namelist for passive tracers
+! -------------
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+!-----------------------------------------------------------------------
+!       nattrc  
+!-----------------------------------------------------------------------
+!
+!       NATTRC   
+!       nwritetrc time step frequency for tracer outputs
+!       lrsttr    boolean term for tracer model restart (true or false)
+!       nrsttr    control of the time step for tracer model restart (0, 1 or 2)
+!       tracer type defined by :
+!                 *  short name 
+!                 *  long_name
+!                 *  units
+!                 *  logical to read initial value from file or not
+!                 *  multiplicative coefficient
+!                 *  logical to save value  
+! 
+&nattrc
+   nwritetrc   = 5475 
+   lrsttr      = .true.
+   nrsttr      = 0
+   tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration          ',  'umol/L' ,  .false. ,  .true.
+   tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration             ',  'ueq/L ' ,  .false. ,  .true.
+   tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration             ',  'umol/L' ,  .false. ,  .true.
+   tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                      ',  'umol/L' ,  .false. ,  .true.
+   tracer(5)   = 'PO4     ' , 'Phosphate Concentration                    ',  'umol/L' ,  .false. ,  .true.
+   tracer(6)   = 'POC     ' , 'Small particle organic carbon Concentration',  'umol/L' ,  .false. ,  .true.
+   tracer(7)   = 'Si      ' , 'Silicate Concentration                     ',  'umol/L' ,  .false. ,  .true.
+   tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration            ',  'umol/L' ,  .false. ,  .true.
+   tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration             ',  'umol/L' ,  .false. ,  .true.
+   tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration            ',  'umol/L' ,  .false. ,  .true.
+   tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                      ',  'umol/L' ,  .false. ,  .true.
+   tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration              ',  'umol/L' ,  .false. ,  .true.
+   tracer(13)  = 'BSi     ' , 'Diatoms Silicate Concentration             ',  'umol/L' ,  .false. ,  .true.
+   tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration               ',  'umol/L' ,  .false. ,  .true.
+   tracer(15)  = 'BFe     ' , 'Big iron particles Concentration           ',  'umol/L' ,  .false. ,  .true.
+   tracer(16)  = 'GOC     ' , 'Big particle organic carbon Concentration  ',  'umol/L' ,  .false. ,  .true.
+   tracer(17)  = 'SFe     ' , 'Small iron particles Concentration         ',  'umol/L' ,  .false. ,  .true.
+   tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration                ',  'umol/L' ,  .false. ,  .true.
+   tracer(19)  = 'DSi     ' , 'Sinking biogenic Silicate Concentration    ',  'umol/L' ,  .false. ,  .true.
+   tracer(20)  = 'NFe     ' , 'Nano iron Concentration                    ',  'umol/L' ,  .false. ,  .true.
+   tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration              ',  'umol/L' ,  .false. ,  .true.
+   tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration           ',  'umol/L' ,  .false. ,  .true.
+   tracer(23)  = 'NO3     ' , 'Nitrates Concentration                     ',  'umol/L' ,  .false. ,  .true.
+   tracer(24)  = 'NH4     ' , 'Ammonium Concentration                     ',  'umol/L' ,  .false. ,  .true.
+/
+!-----------------------------------------------------------------------
+!       natrtd   dynamical tracers trends (#ifdef key_trc_diatrd) 
+!-----------------------------------------------------------------------
+!   nwritetrd : time step frequency for dynamical trends outputs
+!   luttrd    : logical to keep large diagnostics with trends or not
+!               one value per tracer
+&natrtd
+   nwritetrd  = 5475
+   luttrd(1)  = .false.
+   luttrd(2)  = .false.
+   luttrd(3)  = .false.
+   luttrd(4)  = .false.
+   luttrd(5)  = .false.
+   luttrd(6)  = .false.
+   luttrd(7)  = .false.
+   luttrd(8)  = .false.
+   luttrd(9)  = .false.
+   luttrd(10) = .false.
+   luttrd(11) = .false.
+   luttrd(12) = .false.
+   luttrd(13) = .false.
+   luttrd(14) = .false.
+   luttrd(15) = .false.
+   luttrd(16) = .false.
+   luttrd(17) = .false.
+   luttrd(18) = .false.
+   luttrd(19) = .false.
+   luttrd(20) = .false.
+   luttrd(21) = .false.
+   luttrd(22) = .false.
+   luttrd(23) = .false.
+   luttrd(24) = .false.
+/
+!-----------------------------------------------------------------------
+!       natdia  additional 2D/3D  (#ifdef key_trc_diaadd)
+!-----------------------------------------------------------------------
+!  nwritedia : time step frequency for additional arrays outputs
+!  2D/3D diagnostic type defined by :
+!                 *  short name 
+!                 *  long_name
+!                 *  units
+!                 *  logical to save value or not 
+! 
+&natdia
+   nwritedia   = 5475
+   diag2d(1)   = 'Cflx     ' , 'DIC flux                          ',  'molC/m2/s    '
+   diag2d(2)   = 'Oflx     ' , 'Oxygen flux                       ',  'molC/m2/s    '
+   diag2d(3)   = 'Kg       ' , 'Gas transfer                      ',  'mol/m2/s/uatm'
+   diag2d(4)   = 'Delc     ' , 'Delta CO2                         ',  'uatm         '
+   diag2d(5)   = 'PMO      ' , 'POC export                        ',  'molC/m2/s    '
+   diag2d(6)   = 'PMO2     ' , 'GOC export                        ',  'molC/m2/s    '
+   diag2d(7)   = 'ExpFe1   ' , 'Nano iron export                  ',  'molFe/m2/s   '
+   diag2d(8)   = 'ExpFe2   ' , 'Diatoms iron export               ',  'molFe/m2/s   '
+   diag2d(9)   = 'ExpSi    ' , 'Silicate export                   ',  'molSi/m2/s   '
+   diag2d(10)  = 'ExpCaCO3 ' , 'Calcite export                    ',  'molC/m2/s    '
+   diag2d(11)  = 'heup     ' , 'euphotic layer depth              ',  'm            '
+   diag2d(12)  = 'Fedep    ' , 'Iron dep                          ',  'molFe/m2/s   '
+   diag2d(13)  = 'Nfix     ' , 'Nitrogen Fixation                 ',  'molN/m2/s    '
+   diag3d(1)   = 'PH       ' , 'PH                                ',  '-            '
+   diag3d(2)   = 'CO3      ' , 'Bicarbonates                      ',  'mol/l        '
+   diag3d(3)   = 'CO3sat   ' , 'CO3 saturation                    ',  'mol/l        '
+   diag3d(4)   = 'PAR      ' , 'light penetration                 ',  'W/m2         '
+   diag3d(5)   = 'PPPHY    ' , 'Primary production of nanophyto   ',  'molC/m3/s    '
+   diag3d(6)   = 'PPPHY2   ' , 'Primary production of diatoms     ',  'molC/m3/s    '
+   diag3d(7)   = 'PPZOO    ' , 'Primary production of microzoo    ',  'molC/m3/s    '
+   diag3d(8)   = 'PPZOO2   ' , 'Primary production of mesozoo     ',  'molC/m3/s    '
+   diag3d(9)   = 'PBSi     ' , 'Primary production of Si diatoms  ',  'molSi/m3/s   '
+   diag3d(10)  = 'PFeN     ' , 'Primary production of nano iron   ',  'molFe/m3/s   '
+   diag3d(11)  = 'PFeD     ' , 'Primary production of diatoms iron',  'molFe/m3/s   '
+/
+!-----------------------------------------------------------------------
+!       natnum  numerical schemes 
+!-----------------------------------------------------------------------
+!  ndttrc time step frequency for passive tracers
+!  lhdf logical if true computes horizontal diffusion
+!  rsc tuning coefficient for Smolar. 
+!  rtrn truncation value for Smolar.
+!  ncor number of corrective phases for Smolar. 
+!  crosster logical if true computes Smolar crossterms
+&natnum
+   ndttrc   = 1
+   rsc      = 1.
+   rtrn     = 1.e-15
+   ncortrc  = 1
+   crosster = .false.
+/
+!-----------------------------------------------------------------------
+!       namtrcadv   advection scheme for tracer (option not control by CPP keys)
+!-----------------------------------------------------------------------
+!  ln_trcadv_cen2     2nd order centered scheme    (default F)
+!  ln_trcadv_tvd      TVD scheme                   (default F)
+!  ln_trcadv_muscl    MUSCL scheme                 (default F)
+!  ln_trcadv_muscl2   MUSCL2 scheme                (default F)
+!  ln_trcadv_smolar   SMOLAR scheme                (default T)
+&namtrcadv
+   ln_trcadv_cen2   =  .false.
+   ln_trcadv_tvd    =  .false.
+   ln_trcadv_muscl  =  .true.
+   ln_trcadv_muscl2 =  .false.
+   ln_trcadv_smolar =  .false.
+/
+!
+!-----------------------------------------------------------------------
+!       namtrcbbl   bottom boundary layer scheme
+!-----------------------------------------------------------------------
+!  atrcbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s)
+&namtrcbbl
+   atrcbbl = 1000.
+/
+!-----------------------------------------------------------------------
+!       namtrcldf   lateral diffusion scheme for tracer (option not control by CPP keys)
+!-----------------------------------------------------------------------
+!  Flag to performs lateral diffusion or not :
+!     ln_trcldf_diff
+!  Type of the operator : 
+!     ln_trcldf_lap    laplacian operator          (default T)
+!     ln_trcldf_bilap  bilaplacian operator        (default F)
+!  Direction of action  : 
+!     ln_trcldf_level  iso-level                   (default F)
+!     ln_trcldf_hor    horizontal (geopotential)   (default F)^**
+!     ln_trcldf_iso    iso-neutral                 (default T)^*
+! ^* require key_ldfslp to compute the direction of the lateral diffusion
+! ^** require key_ldfslp in s-coordinate
+!  ahtrb0     background eddy diffusivity for isopycnal diffusion (m2/s)
+!  trcrat     ratio betweeen passive and active tracer diffusion coeff
+!  ahtrc0     horizontal eddy diffus. for passive tracer
+!  aeivtr0    eddy induced veloc. coef. for passive tracer
+&namtrcldf
+   ln_trcldf_diff   =  .true.
+   ln_trcldf_lap    =  .true.
+   ln_trcldf_bilap  =  .false.
+   ln_trcldf_level  =  .false.
+   ln_trcldf_hor    =  .false.
+   ln_trcldf_iso    =  .true.
+   ahtrb0           =  0. 
+   trcrat           =  1.
+   ahtrc0           =  2000
+   aeivtr0          =  2000.
+/
+!-----------------------------------------------------------------------
+!       namtrczdf   vertical physics
+!-----------------------------------------------------------------------
+!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F)
+!                               (=F)  euler backward (F)
+!  n_zdfexp   number of sub-timestep for time splitting scheme
+&namtrczdf
+   ln_trczdf_exp =  .false.
+   n_trczdf_exp =      3
+/
+!-----------------------------------------------------------------------
+!       namtrcdmp   tracer newtonian damping ('key_trcdmp')
+!-----------------------------------------------------------------------
+!  ndmptr    type of damping in temperature and salinity 
+!          (='latitude', damping poleward of 'ndmp' degrees and function 
+!             of the distance-to-coast. Red and Med Seas as ndmp=-1)
+!          (=-1 damping only in Med and Red Seas)
+!  ndmpftr   =1 create a damping.coeff NetCDF file (the 3D damping array)
+!  nmldmptr  type of damping in the mixed layer
+!          (=0 damping throughout the water column)
+!	   (=1 no damping in the mixed layer defined by avt >5cm2/s )
+!	   (=2 no damping in the mixed layer defined rho<rho(surf)+.01 )
+!  sdmptr    surface time scale for internal damping (days)
+!  bdmptr    bottom time scale for internal damping (days)
+!  hdmptr    depth of transition between sdmp and bdmp (meters)
+&namtrcdmp
+   ndmptr   =  20
+   ndmpftr  =   0
+   nmldmptr =   1
+   sdmptr   =  50.
+   bdmptr   = 360.
+   hdmptr   = 800.
+/
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist.trc.sms
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist.trc.sms	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist.trc.sms	(revision 966)
@@ -0,0 +1,226 @@
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+! General namelist for PISCES model
+! -------------
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+!-----------------------------------------------------------------------
+!   natext   air-sea exchange
+!-----------------------------------------------------------------------
+!   atcco2 : atmospheric pCO2
+!
+&natext
+   atcco2 = 278.,
+/
+!-----------------------------------------------------------------------
+!   natbio   biological parameters
+!-----------------------------------------------------------------------
+!  part        part of calcite not dissolved in guts
+!  nrdttrc     time step frequency for biology
+!  wsbio       POC sinking speed
+!  xkmort      half saturation constant for mortality
+!  ferat3      Fe/C in zooplankton
+!  wsbio2      Big particles sinking speed
+!
+&natbio
+   part       =  0.85,
+   nrdttrc    =  4,
+   wsbio      =  2. ,
+   xkmort     =  1.E-7 ,
+   ferat3     =  3.E-6,
+   wsbio2     =  50.,
+/
+!-----------------------------------------------------------------------
+!   natlim   parameters for nutrient limitations
+!-----------------------------------------------------------------------
+!  conc0       Phosphate half saturation
+!  conc1       Phosphate half saturation for diatoms
+!  conc2       Iron half saturation for phyto
+!  conc3       Iron half saturation for diatoms
+!  concnnh4    NH4 half saturation for phyto
+!  concdnh4    NH4 half saturation for diatoms
+!  xksi1       half saturation constant for Si uptake
+!  xksi2       half saturation constant for Si/C
+!  xkdoc       half-saturation constant of DOC remineralization
+!  caco3r      mean rain ratio
+&natlim
+   conc0      =  2.e-6,
+   conc1      =  10E-6 ,
+   conc2      =  0.02E-9 ,
+   conc2m     =  0.08E-9,
+   conc3      =  0.1E-9 ,
+   conc3m     =  0.4E-9,
+   concnnh4   =  1.E-7,
+   concdnh4   =  5.E-7,
+   xksi1      =  2.E-6 ,
+   xksi2      =  3.33E-6 ,
+   xkdoc      =  417.E-6 ,
+   caco3r     =  0.3,
+/
+!-----------------------------------------------------------------------
+!   natprod   parameters for phytoplankton growth
+!-----------------------------------------------------------------------
+!  pislope     P-I slope
+!  pislope2    P-I slope  for diatoms
+!  excret      excretion ratio of phytoplankton
+!  excret2     excretion ratio of diatoms
+!  chlcnm      Minimum Chl/C in nanophytoplankton
+!  chlcdm      Minimum Chl/C in diatoms
+!  fecnm       Maximum Fe/C in nanophytoplankton
+!  fecdm       Minimum Fe/C in diatoms
+!  grosip      mean Si/C ratio
+&natprod
+   pislope    =  4. ,
+   pislope2   =  4. ,
+   excret     =  0.05 ,
+   excret2    =  0.05 ,
+   chlcnm     =  0.033,
+   chlcdm     =  0.05,
+   fecnm      =  10E-6,
+   fecdm      =  15E-6,
+   grosip     =  0.151,
+/
+!-----------------------------------------------------------------------
+!   natmort   parameters for phytoplankton sinks
+!-----------------------------------------------------------------------
+!  wchl        quadratic mortality of phytoplankton
+!  wchld       maximum quadratic mortality of diatoms
+!  mprat       phytoplankton mortality rate
+!  mprat2      Diatoms mortality rate
+!  mpratm      Phytoplankton minimum mortality rate
+&natmort
+   wchl       =  0.001,
+   wchld      =  0.02,
+   mprat      =  0.01,
+   mprat2     =  0.01,
+   mpratm     =  0.01,
+/
+!-----------------------------------------------------------------------
+!   natmes   parameters for mesozooplankton
+!-----------------------------------------------------------------------
+!  grazrat2    maximal mesozoo grazing rate
+!  resrat2     exsudation rate of mesozooplankton
+!  mzrat2      mesozooplankton mortality rate
+!  xprefc      zoo preference for phyto
+!  xprefp      zoo preference for POC
+!  xprefz      zoo preference for zoo
+!  xprefpoc    zoo preference for poc
+!  xkgraz2     half sturation constant for meso grazing
+!  epsher2     Efficicency of Mesozoo growth
+!  sigma2      Fraction of mesozoo excretion as DOM
+!  unass2      non assimilated fraction of P by mesozoo
+!  grazflux    flux-feeding rate
+&natmes
+   grazrat2   =  0.7,
+   resrat2    =  0.005,
+   mzrat2     =  0.03,
+   xprefc     =  1. ,
+   xprefp     =  0.2 ,
+   xprefz     =  1 ,
+   xprefpoc   =  0.2 ,
+   xkgraz2    =  20.E-6 ,
+   epsher2    =  0.33 ,
+   sigma2     =  0.6,
+   unass2     =  0.3 ,
+   grazflux   =  5.e3,
+/
+!-----------------------------------------------------------------------
+!   natbio   biological parameters
+!-----------------------------------------------------------------------
+!  grazrat     maximal zoo grazing rate
+!  resrat      exsudation rate of zooplankton
+!  mzrat       zooplankton mortality rate
+!  xpref2c      Microzoo preference for POM
+!  xpref2p      Microzoo preference for Nanophyto
+!  xpref2d      Microzoo preference for Diatoms
+!  xkgraz      half sturation constant for grazing
+!  epsher      Efficiency of microzoo growth
+!  sigma1      Fraction of microzoo excretion as DOM
+!  unass       non assimilated fraction of phyto by zoo
+&natzoo
+   grazrat    =  4.0,
+   resrat     =  0.03,
+   mzrat      =  0.0,
+   xpref2c    =  0.0 ,
+   xpref2p    =  0.5 ,
+   xpref2d    =  0.5 ,
+   xkgraz     =  20.E-6 ,
+   epsher     =  0.33 ,
+   sigma1     =  0.6,
+   unass      =  0.3 ,
+/
+!
+!-----------------------------------------------------------------------
+!   natrem   parameters for remineralization
+!-----------------------------------------------------------------------
+!  xremik      remineralization rate of DOC
+!  xremip      remineralisation rate of POC
+!  nitrif      NH4 nitrification rate
+!  xsirem      remineralization rate of Si
+!  xlam1       scavenging rate of Iron
+!  oxymin      Half-saturation constant for anoxia
+&natrem
+   xremik    =  0.3,
+   xremip    =  0.025,
+   nitrif    =  0.05,
+   xsirem    =  0.015,
+   xlam1     =  0.005,
+   oxymin    =  1.E-6,
+/
+!-----------------------------------------------------------------------
+!   natcal   Calcite chemistry
+!-----------------------------------------------------------------------
+!  kdca        calcite dissolution rate constant (1/time)
+!  nca         order of dissolution reaction (dimensionless)
+&natcal
+   kdca       =  0.327e3,
+   nca        =  1.,
+/
+!-----------------------------------------------------------------------
+!   natsms   inputs deposition
+!-----------------------------------------------------------------------
+!   bdustfer    boolean for dust input from the atmosphere
+!   briver      boolean for river input of nutrients
+!   bndepo      boolean for atmospheric deposition of N
+!   bsedinput   boolean for Fe input from sediments
+!   sedfeinput  Coastal release of Iron
+!
+&natsms
+   bdustfer   =  .true.,
+   briver     =  .false.,
+   bndepo     =  .false.,
+   bsedinput  =  .false.,
+   sedfeinput =  1E-9,
+   dustsolub  =  0.014,
+/
+!---------------------------------------------------------------------------------
+!   natkrpar   parameter for Kriest parameterization      (#ifdef key_kriest)
+!---------------------------------------------------------------------------------
+!  xkr_eta       Sinking  exponent
+!  xkr_zeta      N content exponent
+!  xkr_mass_min  Minimum mass for Aggregates
+!  xkr_mass_max  Maximum mass for Aggregates
+!
+&natkrpar
+   xkr_eta      = 1.17,
+   xkr_zeta     = 2.28,
+   xkr_mass_min = 0.0002,
+   xkr_mass_max = 1.,
+/
+!---------------------------------------------------------------------------------
+!   natkrsize   size calsses for Kriest parameterization      (#ifdef key_kriest)
+!---------------------------------------------------------------------------------
+!  xkr_sfact     Sinking factor
+!  xkr_stick     Stickiness
+!  xkr_xnnano    Nbr of cell in nano size class
+!  xkr_ndiat     Nbr of cell in diatoms size class
+!  xkr_nmeso     Nbr of cell in mesozoo size class
+!  xkr_naggr     Nbr of cell in aggregates  size class
+!
+&natkrsize
+   xkr_sfact    = 942.,
+   xkr_stick    = 0.5,
+   xkr_nnano    = 2.337,
+   xkr_ndiat    = 3.718,
+   xkr_nmeso    = 7.147,
+   xkr_naggr    = 9.877,
+/
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_ice_lim2
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_ice_lim2	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_ice_lim2	(revision 966)
@@ -0,0 +1,186 @@
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+! ICE namelist :  ice model option and parameter input
+! -------------
+!
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicerun   Share parameters for dynamics/advection/thermo
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  ln_limdyn   : switch for ice dynamics (true) or not (false)
+!  ln_limdmp   : restoring Ice thickness and Fraction leads flag
+!  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere
+!  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0)
+!  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0)
+!
+&namicerun
+   ln_limdyn = .true. 
+   ln_limdmp = .FALSE. 
+   acrit     = 1.0e-06 , 1.0e-06
+   hsndif    = 0.0 
+   hicdif    = 0.0
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceini   parameters for ice initialisation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  ln_limini : Ice initialization state flag
+!  ttest     : threshold water temperature for initial sea ice
+!  hninn     : initial snow thickness in the north
+!  hginn     : initial ice thickness in the north
+!  alinn     : initial leads area in the north
+!  hnins     : initial snow thickness in the south
+!  hgins     : initial ice thickness in the south
+!  alins     : initial leads area in the south
+!
+&namiceini
+   ln_limini = .FALSE. 
+   ttest = 2.0
+   hninn = 0.5
+   hginn = 3.0
+   alinn = 0.05
+   hnins = 0.1
+   hgins = 1.0
+   alins = 0.1
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicedia   parameters for ice diagnostics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  fmtinf : format of the output values 
+!  nfrinf : number of variables written in one line
+!  ntmoy  : instantaneous values of ice evolution or averaging
+!  ninfo  : frequency of ouputs on file ice_evolu in case of averaging
+!
+&namicedia
+   fmtinf = '1PE13.5 '
+   nfrinf = 4
+   ntmoy  = 1
+   ninfo  = 1
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicedyn   parameters for ice dynamic
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  epsd   : tolerance parameter
+!  alpha  : coefficient for semi-implicit coriolis
+!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
+!  dm     : diffusion constant for dynamics.
+!  nbiter : number of sub-time steps for relaxation
+!  nbitdr : maximum number of iterations for relaxation
+!  om     : relaxation constant 
+!  resl   : maximum value for the residual of relaxation 
+!  cw     : drag coefficient for oceanic stress
+!  angvg  : turning angle for oceanic stress
+!  pstar  : first bulk-rheology parameter
+!  c_rhg  : second bulk-rhelogy parameter
+!  etamn  : minimun value for viscosity
+!  creepl : creep limit
+!  ecc    : eccentricity of the elliptical yield curve
+!  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s)
+!
+&namicedyn 
+   epsd   = 1.0e-20
+   alpha  = 0.5
+   dm     = 0.6e+03
+   nbiter = 1
+   nbitdr = 100
+   om     = 0.5
+   resl   = 5.0e-05
+   cw     = 5.0e-03
+   angvg  = 0.0
+   pstar  = 1.0e+04 
+   c_rhg  = 20.0
+   etamn  = 0.0e+07
+   creepl = 2.0e-08
+   ecc    = 2.0
+   ahi0   = 350.e0
+/
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicetrp   parameters for ice advection
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
+&namicetrp 
+   bound  = 0.
+/
+
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicethd   parameters for thermodynamic computation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  hmelt       : maximum melting at the bottom
+!  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere
+!                caution 1.0, 1.0 best value to be used!!! (gilles G.)
+!  hicmin      : ice thickness corr. to max. energy stored in brine pocket
+!  hiclim      : minimum ice thickness
+!  amax        : maximum lead fraction
+!  swiqst      : energy stored in brine pocket (=1) or not (=0)	
+!  sbeta       : numerical caracteritic of the scheme for diffusion in ice
+!               Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)
+!  parlat      : percentage of energy used for lateral ablation
+!  hakspl      : slope of distr. for Hakkinen-Mellor's lateral melting
+!  hibspl      : slope of distribution for Hibler's lateral melting
+!  exld        : exponent for leads-closure rate
+!  hakdif      : coefficient for diffusions of ice and snow
+!  thth         : threshold thickness for comp. of eq. thermal conductivity
+!  hnzst       : thickness of the surf. layer in temp. computation
+!  parsub      : switch for snow sublimation or not
+!  alphs       : coefficient for snow density when snow ice formation
+! 
+&namicethd
+   hmelt   = -0.15
+   hiccrit = 0.3 , 0.3
+   hicmin  = 0.2
+   hiclim  = 0.05
+   amax    = 0.999 
+   swiqst  = 1.
+   sbeta   = 1.
+   parlat  = 0.0
+   hakspl  = 0.5
+   hibspl  = 0.5
+   exld    = 2.0
+   hakdif  = 1.0
+   thth    = 0.2
+   hnzst   = 0.1
+   parsub  = 1.0
+   alphs   = 1.0
+/
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceout   parameters for outputs
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!   noumef      : number of fields
+!   field
+!        title  : title of the field
+!        name   : name of the field
+!        unit   : unit of the field
+!        nc     : switch for saving field ( = 1 ) or not ( = 0 )
+!        cmulti : multiplicative constant
+!        cadd   : additive constant
+!
+&namiceout
+   noumef   = 19
+   field_1  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0
+   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0
+   field_3  = 'Ice produced                       ', 'iiceprod', 'm/kt    ', 1, 1.0, 0.0
+   field_4  = 'Ice concentration                  ', 'ileadfra', '%       ', 1, -1.0, 1.0
+   field_5  = 'Ice temperature                    ', 'iicetemp', 'C       ', 1, 1.0, -273.15
+   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0
+   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0
+   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0
+   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15
+   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0
+   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0
+   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0
+   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0
+   field_14 = 'Salt flux at ocean surface         ', 'iocesafl', 'kg/m2/kt', 1, 1.0, 0.0
+   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0
+   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0 
+   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0
+   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0
+   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0
+/      
Index: /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_ice_lim3
===================================================================
--- /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_ice_lim3	(revision 966)
+++ /trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_ice_lim3	(revision 966)
@@ -0,0 +1,313 @@
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+! ICE namelist :  ice model option and parameter input
+! -------------
+!
+!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicerun   Share parameters for dynamics/advection/thermo
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  ln_limdyn        : switch for ice dynamics (true) or not (false)
+!  acrit(1/2)  : minimum fraction for leads in the Northern (Southern) Hemisphere
+!  hnsdif      : computation of temperature in snow (=0.0) or not (=9999.0)
+!  hicdif      : computation of temperature in ice  (=0.0) or not (=9999.0)
+!  cai         : atmospheric drag over sea ice
+!  cao         : atmospheric drag over ocean
+!  ln_nicep  : Ice points output for debug (yes or no)
+!
+&namicerun
+   ln_limdyn    = .true. 
+   acrit   = 1.0e-02 , 1.0e-02
+   hsndif  = 0.0 
+   hicdif  = 0.0
+   cai     = 1.40e-3
+   cao     = 1.00e-3
+   ln_nicep  = .true.
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceini   parameters for ice initialisation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  ttest   : threshold water temperature for initial sea ice
+!  hninn   : initial snow thickness in the north
+!  hginn_u : initial undeformed ice thickness in the north
+!  aginn_u : initial undeformed ice concentration in the north
+!  hginn_d : initial  deformed  ice thickness in the north
+!  aginn_d : initial  deformed  ice concentration in the north
+!  hnins   : initial snow thickness in the south
+!  hgins_u : initial undeformed ice thickness in the north
+!  agins_u : initial undeformed ice concentration in the north
+!  hgins_d : initial  deformed  ice thickness in the north
+!  agins_d : initial  deformed  ice concentration in the north
+!  sinn    : initial salinity in the north
+!  sins    : initial salinity in the south
+!
+&namiceini
+   ttest   = 2.0
+   hninn   = 0.3
+   hginn_u = 3.50
+   aginn_u = 0.95
+   hginn_d = 0.0
+   aginn_d = 0.00
+   hnins = 0.1
+   hgins_u = 1.0
+   agins_u = 0.9
+   hgins_d = 0.0
+   agins_d = 0.00
+   sinn  = 6.301
+   sins  = 6.301
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicedia   parameters for ice diagnostics
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  fmtinf : format of the output values 
+!  nfrinf : number of variables written in one line
+!  ntmoy  : instantaneous values of ice evolution or averaging
+!  ninfo  : frequency of ouputs on file ice_evolu in case of averaging
+!
+&namicedia
+   fmtinf = '1PE13.5 '
+   nfrinf = 4
+   ntmoy  = 1
+   ninfo  = 1
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicedyn   parameters for ice dynamic
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  epsd   : tolerance parameter
+!  alpha  : coefficient for semi-implicit coriolis
+!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
+!  dm     : diffusion constant for dynamics.
+!  nbiter : number of sub-time steps for relaxation
+!  nbitdr : maximum number of iterations for relaxation
+!  om     : relaxation constant 
+!  resl   : maximum value for the residual of relaxation 
+!  cw     : drag coefficient for oceanic stress
+!  angvg  : turning angle for oceanic stress
+!  pstar  : first bulk-rheology parameter
+!  c_rhg  : second bulk-rhelogy parameter
+!  etamn  : minimun value for viscosity
+!  creepl : creep limit
+!  ecc    : eccentricity of the elliptical yield curve
+!  ahi0   : horizontal eddy diffusivity coefficient for sea-ice (m2/s)
+!  nevp   : number of iterations for subcycling in EVP, SB, 400
+!  telest : timescale for elastic waves, SB, 720.0
+!  alphaevp : coefficient for the solution of internal ice stresses
+!
+&namicedyn 
+   epsd   = 1.0e-20
+   alpha  = 0.5
+   dm     = 0.6e+03
+   nbiter = 1
+   nbitdr = 500
+   om     = 0.5
+   resl   = 5.0e-05
+   cw     = 5.0e-03
+   angvg  = 0.0
+   pstar  = 4.00e+04 
+   c_rhg  = 20.0
+   etamn  = 0.0e+07
+   creepl = 1.0e-12
+   ecc    = 2.0
+   ahi0   = 350.e0
+   nevp   = 300
+   telast = 9600.0
+   alphaevp = 1.0
+/
+!
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicetrp   parameters for ice advection
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!  bound  : boundary conditions (=0.0 no-slip, =1.0 free-slip)
+&namicetrp 
+   bound  = 0.
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicethd   parameters for thermodynamic computation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!  hmelt       : maximum melting at the bottom
+!  hiccrit(1/2): ice thickness for lateral accretion in the Northern (Southern) Hemisphere
+!                caution 1.0, 1.0 best value to be used!!! (gilles G.)
+!  fraz_swi    : use of frazil ice collection thickness in function of wind (1.0) or not (0.0)
+!  maxfrazb    : maximum portion of frazil ice collecting at the ice bottom
+!  vfrazb      : thresold drift speed for frazil ice collecting at the ice bottom
+!  Cfrazb      : squeezing coefficient for frazil ice collecting at the ice bottom
+!  hicmin      : ice thickness corr. to max. energy stored in brine pocket
+!  hiclim      : minimum ice thickness
+!  amax        : maximum lead fraction
+!  sbeta       : numerical caracteritic of the scheme for diffusion in ice
+!               Cranck-Nicholson (=0.5), implicit (=1), explicit (=0)
+!  parlat      : percentage of energy used for lateral ablation
+!  hakspl      : slope of distr. for Hakkinen-Mellor's lateral melting
+!  hibspl      : slope of distribution for Hibler's lateral melting
+!  exld        : exponent for leads-closure rate
+!  hakdif      : coefficient for diffusions of ice and snow
+!  thth         : threshold thickness for comp. of eq. thermal conductivity
+!  hnzst       : thickness of the surf. layer in temp. computation
+!  parsub      : switch for snow sublimation or not
+!  alphs       : coefficient for snow density when snow ice formation
+!  betas       : exponent in lead-ice fractionation of snow precipitation 0.66
+!                betas = 1 -> equipartition, betas < 1 -> more on leads
+!  kappa_i     : extinction radiation parameter in sea ice (1.0)
+!  nconv_i_thd : maximal number of iterations for heat diffusion computation
+!  maxer_i_thd : maximal error in temperature for heat diffusion computation
+!  thcon_i_swi : switch for computation of thermal conductivity in the ice
+!                (0) Untersteiner (1964), (1) Pringle et al. (2007)
+! 
+&namicethd
+   hmelt   = -0.15
+   hiccrit = 0.1 , 0.1
+   fraz_swi= 0.0
+   maxfrazb= 0.0
+   vfrazb  = 0.4166667
+   Cfrazb  = 5.0
+   hicmin  = 0.2
+   hiclim  = 0.10
+   amax    = 0.999 
+   sbeta   = 1.
+   parlat  = 0.0
+   hakspl  = 0.25
+   hibspl  = 0.5
+   exld    = 2.0
+   hakdif  = 1.0
+   thth    = 0.2
+   hnzst   = 0.1
+   parsub  = 1.0
+   alphs   = 1.0
+   betas   = 0.6
+   kappa_i = 1.0 
+   nconv_i_thd = 50
+   maxer_i_thd = 0.0001
+   thcon_i_swi = 1
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namicesal   parameters for salinity calculation
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!   num_sal = number of salinity option
+!      1 -> S = bulk_sal
+!      2 -> S = S(z,t) with a simple parameterization
+!      3 -> S = S(z) profile of Scwharzacher [1959]
+!      4 -> S = S(h) Cox and Weeks [1974]
+!   bulk_sal = if 1 is used, it represents the ice salinity
+!   sal_G  : restoring salinity for GD 
+!   time_G : restoring time for GD
+!   sal_F  : restoring salinity for flushing
+!   time_F : restoring time for flushing 
+!   s_i_max : Maximum salinity     
+!   s_i_min : Minimum tolerated ice salinity   
+!   s_i_0   : 1st salinity for salinity profile 
+!   s_i_1   : 2nd salinity for salinity profile 
+
+!
+&namicesal
+   num_sal  = 2
+   bulk_sal = 4.0
+   sal_G    = 5.00
+   time_G   = 1.728e+06
+   sal_F    = 2.00
+   time_F   = 8.640e+05
+   s_i_max  = 20.0
+   s_i_min  =  0.1
+   s_i_0    =  3.5
+   s_i_1    =  4.5
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceitdme   parameters for mechanical redistribution of ice 
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!   ridge_scheme_swi ! which ridging scheme using (1=Rothrock,else=Hibler79)
+!   Cs               ! shearing energy contribution to ridging 
+!   Cf               ! ratio of ridging work to PE change in ridging        
+!   fsnowrdg         ! snow fraction that survives in ridging
+!   fsnowrdg         ! snow fraction that survives in rafting
+!   Gstar            ! fractional area of thin ice being ridged
+!   astar            ! equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice)
+!   Hstar            ! parameter determining the maximum thickness of ridged ice
+!   raftswi          ! rafting or not
+!   hparmeter        ! threshold thickness for rafting or not
+!   Craft            ! coefficient used in the rafting function
+!   ridge_por        ! intial porosity of the ridged ice (typically 0.30)
+!   sal_max_ridge    ! maximum ridged ice salinity
+!   partfun_swi      ! participation function linear, TH75 (0) or exponential Letal07 (1)
+!   transfun_swi     ! transfer function uniform of H80 (0) or exponential Letal07 (1)
+!   brinstren_swi    ! (1) use brine volume to diminish ice strength
+&namiceitdme
+   ridge_scheme_swi = 0
+   Cs = 0.50
+   Cf = 17.0
+   fsnowrdg  = 0.5
+   fsnowrft  = 0.5
+   Gstar     = 0.15
+   astar     = 0.05
+   Hstar     = 100.0
+   raftswi   = 1
+   hparmeter = 0.75
+   Craft     = 5.0
+   ridge_por = 0.3
+   sal_max_ridge = 15.0
+   partfun_swi   = 1
+   transfun_swi  = 0
+   brinstren_swi = 0
+/
+!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
+!       namiceout   parameters for outputs
+!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
+!
+!   noumef      : number of fields
+!   field
+!        title  : title of the field
+!        name   : name of the field
+!        unit   : unit of the field
+!        nc     : switch for saving field ( = 1 ) or not ( = 0 )
+!        cmulti : multiplicative constant
+!        cadd   : additive constant
+!  add_diag_swi : 1 -> diagnose distribution in thickness space
+!                 0 -> only simple diagnostics
+!
+&namiceout
+   noumef   = 37
+   field_1  = 'Ice concentration                  ', 'iiceconc', '%       ', 1, 1.0, 0.0
+   field_2  = 'Ice thickness                      ', 'iicethic', 'm       ', 1, 1.0, 0.0
+   field_3  = 'Snow thickness                     ', 'isnowthi', 'm       ', 1, 1.0, 0.0
+   field_4  = 'Daily bottom thermo ice production ', 'iicebopr', 'cm/day  ', 1, 100., 0.0
+   field_5  = 'Daily dynamic ice production       ', 'iicedypr', 'cm/day  ', 1, 100., 0.0
+   field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ', 1, 1.0, 0.0
+   field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ', 1, 1.0, 0.0
+   field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ', 1, 1.0, 0.0
+   field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ', 1, 1.0, -273.15
+   field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ', 1, 1.0, 0.0
+   field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ', 1, 1.0, 0.0
+   field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ', 1, 1.0, 0.0
+   field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ', 1, 1.0, 0.0
+   field_14 = 'Heat flux due to brine release     ', 'iocehebr', 'w/m2    ', 1, 1.0, 0.0
+   field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ', 1, 1.0, 0.0
+   field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ', 1, 1.0, 0.0 
+   field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ', 1, 1.0, 0.0
+   field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ', 1, 1.0, 0.0
+   field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/day  ', 1, 1.0, 0.0
+   field_20 = 'Mean ice salinity                  ', 'iicesali', 'psu     ', 1, 1.0, 0.0
+   field_21 = 'Mean ice age                       ', 'iiceages', 'years   ', 1, 0.002739, 0.0
+   field_22 = 'Daily lateral thermo ice prod.     ', 'iicelapr', 'cm/day  ', 1, 100., 0.0
+   field_23 = 'Daily snowice ice production       ', 'iicesipr', 'cm/day  ', 1, 100., 0.0
+   field_24 = 'Mean ice temperature               ', 'iicetemp', 'C       ', 1, 1.0, -273.15
+   field_25 = 'Ice total heat content             ', 'iiceheco', '10^9 J  ', 1, 1.0, 0.0
+   field_26 = 'Ice surface temperature            ', 'iicesurt', 'C       ', 1, 1.0, -273.15
+   field_27 = 'Snow temperature                   ', 'isnotem2', 'C       ', 1, 1.0, -273.15
+   field_28 = 'Fsbri - brine salt flux            ', 'iicfsbri', 'kg/m2/s ', 1, 1.0, 0.0
+   field_29 = 'Fseqv - equivalent FW salt flux    ', 'iicfseqv', 'kg/m2/s ', 1, 1.0, 0.0
+   field_30 = 'Brine volume                       ', 'ibrinvol', '%       ', 1, 100.0, 0.0
+   field_31 = 'Frazil ice collection thickness    ', 'iicecolf', 'm       ', 1, 1.0, 0.0
+   field_32 = 'Ice strength                       ', 'iicestre', 'N/m     ', 1, 0.001, 0.0
+   field_33 = 'Ice velocity                       ', 'iicevelo', 'm/s     ', 1, 1.0, 0.0
+   field_34 = 'Surface melt                       ', 'iicesume', 'cm/day  ', 1, 100., 0.0
+   field_35 = 'Bottom melt                        ', 'iicebome', 'cm/day  ', 1, 100., 0.0
+   field_36 = 'Divergence                         ', 'iicedive', '10-8s-1 ', 1, 1.0e8, 0.0
+   field_37 = 'Shear                              ', 'iiceshea', '10-8s-1 ', 1, 1.0e8, 0.0
+   add_diag_swi = 1
+/      
