Ignore:
Timestamp:
03/28/14 17:02:32 (10 years ago)
Author:
aclsce
Message:

IPSLCM6_rc0 : modified to fit with new NEMO version (trunk rev 4565 instead of branch dev_MERGE_2013)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/PARAM/namelist_ORCA2_ref

    r2188 r2221  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  :  1 - run manager      (namrun) 
     2!! NEMO/OPA  :  1 - run manager      (namrun, namcfg) 
    33!! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
    4 !!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
     4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
     
    1111!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
    1212!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namptr, namhsb) 
    13 !!             10 - miscellaneous    (namsol, nammpp, nammpp_dyndist, namctl) 
     13!!             10 - miscellaneous    (namsol, nammpp, namctl) 
    1414!!             11 - Obs & Assim      (namobs, nam_asminc) 
    1515!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    2525!----------------------------------------------------------------------- 
    2626   nn_no       =       0   !  job number (no more used...) 
    27    cn_exp      =  "ORCA2"  !  AUTO - experience name 
    28    nn_it000    =       1   !  AUTO - first time step 
    29    nn_itend    =    5475   !  AUTO - last  time step (std 5475) 
    30    nn_date0    =  010101   !  AUTO - date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    31    nn_leapy    =       0   !  AUTO - Leap year calendar (1) or not (0) 
    32    ln_rstart   = .false.   !  AUTO - start from rest (F) or from a restart file (T) 
    33    nn_rstctl   =       0   !  AUTO - restart control => activated only if ln_rstart = T 
     27   cn_exp      =  "ORCA2"  !  experience name 
     28   nn_it000    =       1   !  first time step 
     29   nn_itend    =    5475   !  last  time step (std 5475) 
     30   nn_date0    =  010101   !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
     31   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
     32   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     33   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=.true. 
     34   nn_rstctl   =       0   !  restart control => activated only if ln_rstart = T 
    3435                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    3536                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
    3637                           !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    37    cn_ocerst_in  = "restartopa"   !  suffix of ocean restart name (input) 
     38   cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input) 
    3839   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
    3940   nn_istate   =       0   !  output the initial state (1) or not (0) 
    40    nn_stock    =    5475   !  AUTO - frequency of creation of a restart file (modulo referenced to 1) 
    41    nn_write    =    5475   !  AUTO - frequency of write in the output file   (modulo referenced to nn_it000) 
     41   nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
     42   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    4243   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    43    ln_mskland  = .true.   !  mask land points in NetCDF outputs (costly: + ~15%) 
     44   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4445   ln_clobber  = .false.   !  clobber (overwrite) an existing file 
    4546   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    4647/ 
    47  
     48! 
     49!----------------------------------------------------------------------- 
     50&namcfg     !   default parameters of the configuration       
     51!----------------------------------------------------------------------- 
     52   cp_cfg      =  "default"            !  name of the configuration 
     53   cp_cfz      =         ''            !  name of the zoom of configuration 
     54   jp_cfg      =       0               !  resolution of the configuration 
     55   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
     56   jpjdta      =      12               !  2nd    "         "    ( >= jpj ) 
     57   jpkdta      =      31               !  number of levels      ( >= jpk ) 
     58   jpiglo      =      10               !  1st dimension of global domain --> i =jpidta 
     59   jpjglo      =      12               !  2nd    -                  -    --> j  =jpjdta 
     60   jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
     61   jpjzoom     =       1               !  in data domain indices 
     62   jperio      =       0               !  lateral cond. type (between 0 and 6) 
     63                                       !  = 0 closed                 ;   = 1 cyclic East-West 
     64                                       !  = 2 equatorial symmetric   ;   = 3 North fold T-point pivot  
     65                                       !  = 4 cyclic East-West AND North fold T-point pivot 
     66                                       !  = 5 North fold F-point pivot 
     67                                       !  = 6 cyclic East-West AND North fold F-point pivot 
     68/ 
    4869!!====================================================================== 
    4970!!                      ***  Domain namelists  *** 
     
    91112!----------------------------------------------------------------------- 
    92113   nn_bathy    =    1      !  compute (=0) or read (=1) the bathymetry file 
    93    nn_closea    =   1      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
    94    nn_msh      =    0      !  AUTO - create (=1) a mesh file or not (=0) 
     114   rn_bathy    =    0.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
     115   nn_closea   =    0      !  remove (=0) or keep (=1) closed seas and lakes (ORCA) 
     116   nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
    95117   rn_hmin     =   -3.     !  min depth of the ocean (>0) or min number of ocean level (<0) 
    96118   rn_e3zps_min=   20.     !  partial step thickness is set larger than the minimum of 
     
    98120                           ! 
    99121   rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
    100    nn_baro     =   64      !  number of barotropic time step            ("key_dynspg_ts") 
    101122   rn_atfp     =    0.1    !  asselin time filter parameter 
    102123   nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
     
    105126   rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1) 
    106127   rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
     128   ln_crs      = .false.      !  Logical switch for coarsening module 
     129   jphgr_msh   =       0               !  type of horizontal mesh 
     130                                       !  = 0 curvilinear coordinate on the sphere read in coordinate.nc 
     131                                       !  = 1 geographical mesh on the sphere with regular grid-spacing 
     132                                       !  = 2 f-plane with regular grid-spacing 
     133                                       !  = 3 beta-plane with regular grid-spacing 
     134                                       !  = 4 Mercator grid with T/U point at the equator 
     135   ppglam0     =       0.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
     136   ppgphi0     =     -35.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
     137   ppe1_deg    =       1.0             !  zonal      grid-spacing (degrees) 
     138   ppe2_deg    =       0.5             !  meridional grid-spacing (degrees) 
     139   ppe1_m      =    5000.0             !  zonal      grid-spacing (degrees) 
     140   ppe2_m      =    5000.0             !  meridional grid-spacing (degrees) 
     141   ppsur       =    -4762.96143546300  !  ORCA r4, r2 and r05 coefficients 
     142   ppa0        =      255.58049070440  ! (default coefficients) 
     143   ppa1        =      245.58132232490  ! 
     144   ppkth       =       21.43336197938  ! 
     145   ppacr       =        3.0            ! 
     146   ppdzmin     =       10.             !  Minimum vertical spacing 
     147   pphmax      =     5000.             !  Maximum depth 
     148   ldbletanh   =    .TRUE.             !  Use/do not use double tanf function for vertical coordinates 
     149   ppa2        =      100.760928500000 !  Double tanh function parameters 
     150   ppkth2      =       48.029893720000 ! 
     151   ppacr2      =       13.000000000000 ! 
     152/ 
     153!----------------------------------------------------------------------- 
     154&namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
     155!----------------------------------------------------------------------- 
     156   ln_bt_fw      =    .TRUE.           !  Forward integration of barotropic equations 
     157   ln_bt_av      =    .TRUE.           !  Time filtering of barotropic variables 
     158   ln_bt_nn_auto =    .TRUE.           !  Set nn_baro automatically to be just below 
     159                                       !  a user defined maximum courant number (rn_bt_cmax) 
     160   nn_baro       =    30               !  Number of iterations of barotropic mode 
     161                                       !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F 
     162   rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T  
     163   nn_bt_flt     =    1                !  Time filter choice 
     164                                       !  = 0 None 
     165                                       !  = 1 Boxcar over   nn_baro barotropic steps 
     166                                       !  = 2 Boxcar over 2*nn_baro     "        "   
     167/ 
     168!----------------------------------------------------------------------- 
     169&namcrs        !   Grid coarsening for dynamics output and/or 
     170               !   passive tracer coarsened online simulations 
     171!----------------------------------------------------------------------- 
     172   nn_factx    = 3         !  Reduction factor of x-direction 
     173   nn_facty    = 3         !  Reduction factor of y-direction 
     174   nn_binref   = 0         !  Bin centering preference: NORTH or EQUAT 
     175                           !  0, coarse grid is binned with preferential treatment of the north fold 
     176                           !  1, coarse grid is binned with centering at the equator 
     177                           !    Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty. 
     178   nn_msh_crs  = 1         !  create (=1) a mesh file or not (=0) 
     179   nn_crs_kz   = 0         ! 0, MEAN of volume boxes 
     180                           ! 1, MAX of boxes 
     181                           ! 2, MIN of boxes 
     182   ln_crs_wn   = .true.    ! wn coarsened (T) or computed using horizontal divergence ( F ) 
    107183/ 
    108184!----------------------------------------------------------------------- 
    109185&namtsd    !   data : Temperature  & Salinity 
    110186!----------------------------------------------------------------------- 
    111 !          ! file name ! frequency (hours)    ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
    112 !          !           !  (if <0  months)     !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
    113    sn_tem  = 'data_1m_potential_temperature_nomask', -1,'votemper',  .true.  , .true., 'yearly'   , ' '      , ' ' 
    114    sn_sal  = 'data_1m_salinity_nomask'             , -1,'vosaline',  .true.  , .true., 'yearly'   , ''       , ' ' 
     187!----------------------------------------------------------------------- 
     188!          !  file name                            ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     189!          !                                       !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     190   sn_tem  = 'data_1m_potential_temperature_nomask',         -1        ,'votemper' ,    .true.    , .true. , 'yearly'   , ''       ,   ''    ,    '' 
     191   sn_sal  = 'data_1m_salinity_nomask'             ,         -1        ,'vosaline' ,    .true.    , .true. , 'yearly'   , ''       ,   ''    ,    '' 
    115192   ! 
    116193   cn_dir        = './'     !  root directory for the location of the runoff files 
    117194   ln_tsd_init   = .true.   !  Initialisation of ocean T & S with T &S input data (T) or not (F) 
    118    ln_tsd_tradmp = .false.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
     195   ln_tsd_tradmp = .true.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
    119196/ 
    120197!!====================================================================== 
     
    128205!!   namsbc_mfs      MFS  bulk formulae formulation 
    129206!!   namsbc_cpl      CouPLed            formulation                     ("key_coupled") 
     207!!   namsbc_sas      StAndalone Surface module 
    130208!!   namtra_qsr      penetrative solar radiation 
    131209!!   namsbc_rnf      river runoffs 
     
    143221   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    144222   ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio) 
    145    ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
     223   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    146224   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    147    ln_cpl      = .true.    !  Coupled formulation                       (T => fill namsbc_cpl ) 
     225   ln_cpl      = .false.   !  Coupled formulation                       (T => fill namsbc_cpl ) 
    148226   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    149227   nn_ice      = 2         !  =0 no ice boundary condition   , 
    150228                           !  =1 use observed ice-cover      , 
    151                            !  =2 ice-model used                         ("key_lim3" or "key_lim2") 
    152    cn_iceflx   = 'linear'  !  AUTO -  
    153                            !  ="ave"    - Average flux over all ice categories ("key_lim3") 
    154                            !  ="linear" - Linear repartition of fluxes over ice categories ("key_lim3") 
    155    nn_ice_embd = 0         !  AUTO -  
    156                            !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
     229                           !  =2 ice-model used                         ("key_lim3" or "key_lim2) 
     230   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    157231                           !  =1 levitating ice with mass and salt exchange but no presure effect 
    158232                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    159233   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    160    ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    161    ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    162    nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
     234   ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     235   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
     236   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
    163237                           !     =1 global mean of e-p-r set to zero at each time step 
    164238                           !     =2 annual global mean of e-p-r set to zero 
    165                            !     =3 global emp set to zero and spread out over erp area 
    166239   ln_wave = .false.       !  Activate coupling with wave (either Stokes Drift or Drag coefficient, or both)  (T => fill namsbc_wave) 
    167240   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave) 
    168241   ln_sdw  = .false.       !  Computation of 3D stokes drift                (T => fill namsbc_wave) 
     242   cn_iceflx = 'linear'    !  redistribution of solar input into ice categories during coupling ice/atm. 
    169243/ 
    170244!----------------------------------------------------------------------- 
     
    181255&namsbc_flx    !   surface boundary condition : flux formulation 
    182256!----------------------------------------------------------------------- 
    183 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    184 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    185    sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , '' 
    186    sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , '' 
    187    sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , '' 
    188    sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , '' 
    189    sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , '' 
     257!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !  
     258!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     259   sn_utau     = 'utau'      ,        24         , 'utau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
     260   sn_vtau     = 'vtau'      ,        24         , 'vtau'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
     261   sn_qtot     = 'qtot'      ,        24         , 'qtot'    , .false.      , .false., 'yearly'  , ''       , ''       , '' 
     262   sn_qsr      = 'qsr'       ,        24         , 'qsr'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
     263   sn_emp      = 'emp'       ,        24         , 'emp'     , .false.      , .false., 'yearly'  , ''       , ''       , '' 
    190264 
    191265   cn_dir      = './'      !  root directory for the location of the flux files 
     
    194268&namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
    195269!----------------------------------------------------------------------- 
    196 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    197 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    198    sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , '' 
    199    sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , '' 
    200    sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , '' 
    201    sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , '' 
    202    sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , '' 
    203    sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , '' 
    204    sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , '' 
     270!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     271!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     272   sn_utau     = 'taux_1m'   ,       -1          , 'sozotaux',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     273   sn_vtau     = 'tauy_1m'   ,       -1          , 'sometauy',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     274   sn_wndm     = 'flx'       ,       -1          , 'socliowi',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     275   sn_tair     = 'flx'       ,       -1          , 'socliot2',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     276   sn_humi     = 'flx'       ,       -1          , 'socliohu',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     277   sn_ccov     = 'flx'       ,       -1          , 'socliocl',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
     278   sn_prec     = 'flx'       ,       -1          , 'socliopl',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    205279 
    206280   cn_dir      = './'      !  root directory for the location of the bulk files are 
     
    209283&namsbc_core   !   namsbc_core  CORE bulk formulae 
    210284!----------------------------------------------------------------------- 
    211 !              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    212 !              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    213    sn_wndi     = 'u_10.15JUNE2009_orca2'       ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , ''       , 'Uwnd' 
    214    sn_wndj     = 'v_10.15JUNE2009_orca2'       ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , ''       , 'Vwnd' 
    215    sn_qsr      = 'ncar_rad.15JUNE2009_orca2'   ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
    216    sn_qlw      = 'ncar_rad.15JUNE2009_orca2'   ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
    217    sn_tair     = 't_10.15JUNE2009_orca2'       ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
    218    sn_humi     = 'q_10.15JUNE2009_orca2'       ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , ''       , '' 
    219    sn_prec     = 'ncar_precip.15JUNE2009_orca2',        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , ''       , '' 
    220    sn_snow     = 'ncar_precip.15JUNE2009_orca2',        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , ''       , '' 
    221    sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , ''       , '' 
     285!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
     286!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
     287   sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Uwnd'   , '' 
     288   sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bicubic_noc.nc'   , 'Vwnd'   , '' 
     289   sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     290   sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     291   sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     292   sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     293   sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     294   sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
     295   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core_orca2_bilinear_noc.nc'  , ''       , '' 
    222296 
    223297   cn_dir      = './'      !  root directory for the location of the bulk files 
    224298   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
    225299   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
     300   ln_bulk2z   = .false.   !  Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 
     301   rn_zqt      = 3.        !  Air temperature and humidity reference height (m) (ln_bulk2z) 
     302   rn_zu       = 4.        !  Wind vector reference height (m)                  (ln_bulk2z) 
    226303   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
     304   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     305   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity  
     306                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    227307/ 
    228308!----------------------------------------------------------------------- 
    229309&namsbc_mfs   !   namsbc_mfs  MFS bulk formulae 
    230310!----------------------------------------------------------------------- 
    231 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    232 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    233    sn_wndi     =   'ecmwf'   ,        6          , 'u10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
    234    sn_wndj     =   'ecmwf'   ,        6          , 'v10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
    235    sn_clc      =   'ecmwf'   ,        6          , 'clc'     ,    .true.    , .false. , 'daily'  ,'bilinear.nc', '' 
    236    sn_msl      =   'ecmwf'   ,        6          , 'msl'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
    237    sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , '' 
    238    sn_rhm      =   'ecmwf'   ,        6          , 'rh'      ,    .true.    , .false. , 'daily'  ,'bilinear.nc', '' 
    239    sn_prec     =   'ecmwf'   ,        6          , 'precip'  ,    .true.    , .true.  , 'daily'  ,'bicubic.nc' , '' 
     311!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights     ! rotation ! land/sea mask ! 
     312!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename    ! pairing  ! filename      ! 
     313   sn_wndi     =   'ecmwf'   ,        6          , 'u10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
     314   sn_wndj     =   'ecmwf'   ,        6          , 'v10'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
     315   sn_clc      =   'ecmwf'   ,        6          , 'clc'     ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , '' 
     316   sn_msl      =   'ecmwf'   ,        6          , 'msl'     ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
     317   sn_tair     =   'ecmwf'   ,        6          , 't2'      ,    .true.    , .false. , 'daily'  ,'bicubic.nc' , ''       , '' 
     318   sn_rhm      =   'ecmwf'   ,        6          , 'rh'      ,    .true.    , .false. , 'daily'  ,'bilinear.nc', ''       , '' 
     319   sn_prec     =   'ecmwf'   ,        6          , 'precip'  ,    .true.    , .true.  , 'daily'  ,'bicubic.nc' , ''       , '' 
    240320 
    241321   cn_dir      = './ECMWF/'      !  root directory for the location of the bulk files 
     
    249329sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    250330sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    251 sn_snd_thick  =       'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
    252 sn_snd_crt    =       'mixed oce-ice'        ,    'no'    , 'cartesian' , 'eastward-northward' ,  'T' 
     331sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     332sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    253333sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    254334! receive 
    255 sn_rcv_w10m   =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    256 sn_rcv_taumod =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    257 sn_rcv_tau    =       'mixed oce-ice'        ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     335sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     336sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     337sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    258338sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    259 sn_rcv_qsr    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    260 sn_rcv_qns    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     339sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     340sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    261341sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    262342sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     
    265345/ 
    266346!----------------------------------------------------------------------- 
     347&namsbc_sas    !   analytical surface boundary condition 
     348!----------------------------------------------------------------------- 
     349!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     350!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     351   sn_usp      = 'sas_grid_U' ,    120           , 'vozocrtx' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , ''       
     352   sn_vsp      = 'sas_grid_V' ,    120           , 'vomecrty' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     353   sn_tem      = 'sas_grid_T' ,    120           , 'sosstsst' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     354   sn_sal      = 'sas_grid_T' ,    120           , 'sosaline' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     355   sn_ssh      = 'sas_grid_T' ,    120           , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''             , '' 
     356 
     357   ln_3d_uv    = .true.    !  specify whether we are supplying a 3D u,v field 
     358   cn_dir      = './'      !  root directory for the location of the bulk files are 
     359/ 
     360!----------------------------------------------------------------------- 
    267361&namtra_qsr    !   penetrative solar radiation 
    268362!----------------------------------------------------------------------- 
    269 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    270 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    271    sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , '' 
     363!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     364!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     365   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    272366 
    273367   cn_dir      = './'      !  root directory for the location of the runoff files 
    274368   ln_traqsr   = .true.    !  Light penetration (T) or not (F) 
    275    ln_qsr_rgb  = .false.   !  RGB (Red-Green-Blue) light penetration 
     369   ln_qsr_rgb  = .true.    !  RGB (Red-Green-Blue) light penetration 
    276370   ln_qsr_2bd  = .false.   !  2 bands              light penetration 
    277    ln_qsr_bio  = .true.    !  bio-model light penetration 
    278    nn_chldta   =      0    !  RGB : Chl data (=1) or cst value (=0) 
     371   ln_qsr_bio  = .false.   !  bio-model light penetration 
     372   nn_chldta   =      1    !  RGB : Chl data (=1) or cst value (=0) 
    279373   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1) 
    280374   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction 
    281375   rn_si1      =   23.0    !  2 bands: longest depth of extinction 
     376   ln_qsr_ice  = .true.    !  light penetration for ice-model LIM3 
    282377/ 
    283378!----------------------------------------------------------------------- 
    284379&namsbc_rnf    !   runoffs namelist surface boundary condition 
    285380!----------------------------------------------------------------------- 
    286 !              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    287 !              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    288    sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , '' 
    289    sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , '' 
    290    sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , '' 
    291    sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , '' 
    292    sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , '' 
     381!              !  file name           ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     382!              !                      !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     383   sn_rnf      = 'runoff_core_monthly',        -1         , 'sorunoff',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     384   sn_cnf      = 'runoff_core_monthly',         0         , 'socoefr0',   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
     385   sn_s_rnf    = 'runoffs'            ,        24         , 'rosaline',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     386   sn_t_rnf    = 'runoffs'            ,        24         , 'rotemper',   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
     387   sn_dep_rnf  = 'runoffs'            ,         0         , 'rodepth' ,   .false.    , .true. , 'yearly'  , ''       , ''       , '' 
    293388 
    294389   cn_dir       = './'      !  root directory for the location of the runoff files 
     
    305400&namsbc_apr    !   Atmospheric pressure used as ocean forcing or in bulk 
    306401!----------------------------------------------------------------------- 
    307 !              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    308 !              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
    309    sn_apr      = 'patm'    ,         -1        ,'somslpre',    .true.      , .true.  , 'yearly'  ,  ''      ,  '' 
     402!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     403!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     404   sn_apr      = 'patm'      ,         -1        ,'somslpre',    .true.     , .true. , 'yearly'  ,  ''      ,   ''     , '' 
    310405 
    311406   cn_dir      = './'       !  root directory for the location of the bulk files 
     
    317412&namsbc_ssr    !   surface boundary condition : sea surface restoring 
    318413!----------------------------------------------------------------------- 
    319 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    320 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    321    sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , '' 
    322    sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , '' 
     414!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     415!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     416   sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , ''       , '' 
     417   sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , ''       , '' 
    323418 
    324419   cn_dir      = './'      !  root directory for the location of the runoff files 
     
    334429&namsbc_alb    !   albedo parameters 
    335430!----------------------------------------------------------------------- 
    336    rn_cloud    =    0.00   !  cloud correction to snow and ice albedo 
    337    rn_albice   =    0.5    !  albedo of melting ice in the arctic and antarctic 
     431   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo 
     432   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic 
    338433   rn_alphd    =    0.80   !  coefficients for linear interpolation used to 
    339434   rn_alphc    =    0.65   !  compute albedo between two extremes values 
     
    368463      rn_speed_limit           = 0.                   ! CFL speed limit for a berg    
    369464 
    370                ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
    371                !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
    372       sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
     465!              ! file name ! frequency (hours) !   variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     466!              !           !  (if <0  months)  !     name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     467      sn_icb =  'calving' ,       -1           , 'calvingmask',  .true.        , .true.  , 'yearly'  , ''       , ''       , '' 
    373468    
    374469      cn_dir = './'  
     
    431526!----------------------------------------------------------------------- 
    432527   ln_tide_pot   = .true.   !  use tidal potential forcing 
    433    nb_harmo      =    11    !  number of constituents used 
     528   ln_tide_ramp  = .false.  ! 
     529   rdttideramp   =    0.    ! 
    434530   clname(1)     =   'M2'   !  name of constituent 
    435531   clname(2)     =   'S2' 
     
    447543&nambdy        !  unstructured open boundaries                          ("key_bdy") 
    448544!----------------------------------------------------------------------- 
    449     nb_bdy = 1                            !  number of open boundary sets 
     545    nb_bdy         = 0                    !  number of open boundary sets 
    450546    ln_coords_file = .true.               !  =T : read bdy coordinates from file 
    451547    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files 
    452     ln_mask_file = .false.                !  =T : read mask from file 
    453     cn_mask_file = ''                     !  name of mask file (if ln_mask_file=.TRUE.) 
    454     nn_dyn2d      =  2                    !  boundary conditions for barotropic fields 
    455     nn_dyn2d_dta  =  3                    !  = 0, bdy data are equal to the initial state 
     548    ln_mask_file   = .false.              !  =T : read mask from file 
     549    cn_mask_file   = ''                   !  name of mask file (if ln_mask_file=.TRUE.) 
     550    cn_dyn2d       = 'none'               ! 
     551    nn_dyn2d_dta   =  0                   !  = 0, bdy data are equal to the initial state 
    456552                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
    457553                                          !  = 2, use tidal harmonic forcing data from files 
    458554                                          !  = 3, use external data AND tidal harmonic forcing 
    459     nn_dyn3d      =  0                    !  boundary conditions for baroclinic velocities 
     555    cn_dyn3d      =  'none'               !   
    460556    nn_dyn3d_dta  =  0                    !  = 0, bdy data are equal to the initial state 
    461                            !  = 1, bdy data are read in 'bdydata   .nc' files 
    462     nn_tra        =  1                    !  boundary conditions for T and S 
    463     nn_tra_dta    =  1                    !  = 0, bdy data are equal to the initial state 
    464                            !  = 1, bdy data are read in 'bdydata   .nc' files 
    465     nn_rimwidth  = 10                      !  width of the relaxation zone 
    466     ln_vol     = .false.                  !  total volume correction (see nn_volctl parameter) 
    467     nn_volctl  = 1                        !  = 0, the total water flux across open boundaries is zero 
     557                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     558    cn_tra        =  'none'               !  
     559    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
     560                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     561    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
     562    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
     563    rn_time_dmp   =  1.                   ! Damping time scale in days  
     564    rn_time_dmp_out =  1.                 ! Outflow damping time scale 
     565    nn_rimwidth   = 10                    !  width of the relaxation zone 
     566    ln_vol        = .false.               !  total volume correction (see nn_volctl parameter) 
     567    nn_volctl     = 1                     !  = 0, the total water flux across open boundaries is zero 
    468568/ 
    469569!----------------------------------------------------------------------- 
    470570&nambdy_dta      !  open boundaries - external data           ("key_bdy") 
    471571!----------------------------------------------------------------------- 
    472 !              !   file name    ! frequency (hours) !  variable  ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
    473 !              !                !  (if <0  months)  !    name    !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
    474    bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
    475    bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
    476    bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
    477    bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
    478    bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
    479    bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
    480    bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   '' 
     572!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     573!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     574   bn_ssh =     'amm12_bdyT_u2d' ,         24        , 'sossheig' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     575   bn_u2d =     'amm12_bdyU_u2d' ,         24        , 'vobtcrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     576   bn_v2d =     'amm12_bdyV_u2d' ,         24        , 'vobtcrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     577   bn_u3d  =    'amm12_bdyU_u3d' ,         24        , 'vozocrtx' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     578   bn_v3d  =    'amm12_bdyV_u3d' ,         24        , 'vomecrty' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     579   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     580   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    481581   cn_dir  =    'bdydta/' 
    482582   ln_full_vel = .false. 
     
    485585&nambdy_tide     ! tidal forcing at open boundaries 
    486586!----------------------------------------------------------------------- 
    487    filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
    488     tide_cpt(1)   ='Q1'  !  names of tidal components used 
    489     tide_cpt(2)   ='O1'  !  names of tidal components used 
    490     tide_cpt(3)   ='P1'  !  names of tidal components used 
    491     tide_cpt(4)   ='S1'  !  names of tidal components used 
    492     tide_cpt(5)   ='K1'  !  names of tidal components used 
    493     tide_cpt(6)   ='2N2' !  names of tidal components used 
    494     tide_cpt(7)   ='MU2' !  names of tidal components used 
    495     tide_cpt(8)   ='N2'  !  names of tidal components used 
    496     tide_cpt(9)   ='NU2' !  names of tidal components used 
    497     tide_cpt(10)   ='M2'  !  names of tidal components used 
    498     tide_cpt(11)   ='L2'  !  names of tidal components used 
    499     tide_cpt(12)   ='T2'  !  names of tidal components used 
    500     tide_cpt(13)   ='S2'  !  names of tidal components used 
    501     tide_cpt(14)   ='K2'  !  names of tidal components used 
    502     tide_cpt(15)   ='M4'  !  names of tidal components used 
    503     tide_speed(1)   = 13.398661 !  phase speeds of tidal components (deg/hour) 
    504     tide_speed(2)   = 13.943036 !  phase speeds of tidal components (deg/hour) 
    505     tide_speed(3)   = 14.958932 !  phase speeds of tidal components (deg/hour) 
    506     tide_speed(4)   = 15.000001 !  phase speeds of tidal components (deg/hour) 
    507     tide_speed(5)   = 15.041069 !  phase speeds of tidal components (deg/hour) 
    508     tide_speed(6)   = 27.895355 !  phase speeds of tidal components (deg/hour) 
    509     tide_speed(7)   = 27.968210 !  phase speeds of tidal components (deg/hour) 
    510     tide_speed(8)   = 28.439730 !  phase speeds of tidal components (deg/hour) 
    511     tide_speed(9)   = 28.512585 !  phase speeds of tidal components (deg/hour) 
    512     tide_speed(10)   = 28.984106 !  phase speeds of tidal components (deg/hour) 
    513     tide_speed(11)   = 29.528479 !  phase speeds of tidal components (deg/hour) 
    514     tide_speed(12)   = 29.958935 !  phase speeds of tidal components (deg/hour) 
    515     tide_speed(13)   = 30.000002 !  phase speeds of tidal components (deg/hour) 
    516     tide_speed(14)   = 30.082138 !  phase speeds of tidal components (deg/hour) 
    517     tide_speed(15)   = 57.968212 !  phase speeds of tidal components (deg/hour) 
    518     ln_tide_date = .true.               !  adjust tidal harmonics for start date of run 
     587   filtide          = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     588   ln_bdytide_2ddta = .false. 
     589   ln_bdytide_conj  = .false. 
    519590/ 
    520591!!====================================================================== 
     
    532603                           !                              = 2 : nonlinear friction 
    533604   rn_bfri1    =    4.e-4  !  bottom drag coefficient (linear case) 
    534    rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case) 
     605   rn_bfri2    =    1.e-3  !  bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
     606   rn_bfri2_max =   1.e-1  !  max. bottom drag coefficient (non linear case and ln_loglayer=T) 
    535607   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
     608   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T  
     609   ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    536610   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    537611   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     
    552626   nn_bbl_ldf  =  1      !  diffusive bbl (=1)   or not (=0) 
    553627   nn_bbl_adv  =  0      !  advective bbl (=1/2) or not (=0) 
    554    rn_ahtbbl   =  10000.  !  lateral mixing coefficient in the bbl  [m2/s] 
     628   rn_ahtbbl   =  1000.  !  lateral mixing coefficient in the bbl  [m2/s] 
    555629   rn_gambbl   =  10.    !  advective bbl coefficient                 [s] 
    556630/ 
     
    585659   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    586660   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     661/ 
     662!----------------------------------------------------------------------- 
     663&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
     664!----------------------------------------------------------------------- 
     665   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     666   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
     667   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     668   rn_lf     = 5.e+3       ! typical scale of mixed layer front (meters)                      (case rn_mle=0) 
     669   rn_time   = 172800.     ! time scale for mixing momentum across the mixed layer (seconds)  (case rn_mle=0) 
     670   rn_lat    = 20.         ! reference latitude (degrees) of MLE coef.                        (case rn_mle=1) 
     671   nn_mld_uv = 0           ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) 
     672   nn_conv   = 0           ! =1 no MLE in case of convection ; =0 always MLE 
     673   rn_rho_c_mle  = 0.01    ! delta rho criterion used to calculate MLD for FK 
    587674/ 
    588675!---------------------------------------------------------------------------------- 
     
    609696   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    610697   !                                           (normally=0; not used with Griffies) 
     698   rn_slpmax        =     0.01  !  slope limit 
     699   rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
     700   rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
     701   rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
    611702/ 
    612703!----------------------------------------------------------------------- 
    613704&namtra_dmp    !   tracer: T & S newtonian damping 
    614705!----------------------------------------------------------------------- 
    615    ln_tradmp   =  .false.   !  add a damping termn (T) or not (F) 
     706   ln_tradmp   =  .true.   !  add a damping termn (T) or not (F) 
    616707   nn_hdmp     =   -1      !  horizontal shape =-1, damping in Med and Red Seas only 
    617708                           !                   =XX, damping poleward of XX degrees (XX>0) 
     
    642733   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    643734   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     735/ 
     736!----------------------------------------------------------------------- 
     737&nam_vvl    !   vertical coordinate options 
     738!----------------------------------------------------------------------- 
     739   ln_vvl_zstar  = .true.           !  zstar vertical coordinate                    
     740   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
     741   ln_vvl_layer  = .false.          !  full layer vertical coordinate 
     742   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar 
     743   ln_vvl_zstar_at_eqtor = .false.  !  ztilde near the equator 
     744   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient 
     745   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days] 
     746   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days] 
     747   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation 
     748   ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
    644749/ 
    645750!----------------------------------------------------------------------- 
     
    683788   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    684789   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
     790   rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
     791   rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
     792   rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
     793   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
     794   rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
    685795/ 
    686796 
     
    730840   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
    731841   rn_ediss    =   0.7     !  coef. of the Kolmogoroff dissipation 
    732    rn_ebb      =  60.0    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 
     842   rn_ebb      =  67.83    !  coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) 
    733843   rn_emin     =   1.e-6   !  minimum value of tke [m2/s2] 
    734844   rn_emin0    =   1.e-4   !  surface minimum value of tke [m2/s2] 
     845   rn_bshear   =   1.e-20  ! background shear (>0) currently a numerical threshold (do not change it) 
    735846   nn_mxl      =   2       !  mixing length: = 0 bounded by the distance to surface and bottom 
    736847                           !                 = 1 bounded by the local vertical scale factor 
     
    739850   nn_pdl      =   1       !  Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) 
    740851   ln_mxl0     = .true.    !  surface mixing length scale = F(wind stress) (T) or not (F) 
    741    rn_mxl0     =   0.01    !  surface  buoyancy lenght scale minimum value 
     852   rn_mxl0     =   0.04    !  surface  buoyancy lenght scale minimum value 
    742853   ln_lc       = .true.    !  Langmuir cell parameterisation (Axell 2002) 
    743854   rn_lc       =   0.15    !  coef. associated to Langmuir cells 
    744    nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
     855   nn_etau     =   1       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    745856                           !        = 0 no penetration 
    746857                           !        = 1 add a tke source below the ML 
     
    828939   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
    829940   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold 
    830    jpni        =   1       !  jpni   number of processors following i (set automatically if < 1) 
    831    jpnj        =   5       !  jpnj   number of processors following j (set automatically if < 1) 
    832    jpnij       =   5       !  jpnij  number of local domains (set automatically if < 1) 
     941   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1) 
     942   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1) 
     943   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1) 
    833944/ 
    834945!----------------------------------------------------------------------- 
     
    847958   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    848959/ 
    849  
     960!----------------------------------------------------------------------- 
     961&namc1d        !   1D configuration options                             ("key_c1d") 
     962!----------------------------------------------------------------------- 
     963   rn_lat      =    50     !  Column latitude 
     964   rn_lon      =    -145   !  Column longitude 
     965/ 
     966!----------------------------------------------------------------------- 
     967&namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
     968!----------------------------------------------------------------------- 
     969!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     970!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     971   sn_ucur     = 'ucurrent'  ,         -1        ,'u_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Ume'   , '' 
     972   sn_vcur     = 'vcurrent'  ,         -1        ,'v_current',   .false.    , .true. , 'monthly' ,  ''      ,  'Vme'   , '' 
     973! 
     974   cn_dir        = './'    !  root directory for the location of the files 
     975   ln_uvd_init   = .false. !  Initialisation of ocean U & V with U & V input data (T) or not (F) 
     976   ln_uvd_dyndmp = .false. !  damping of ocean U & V toward U & V input data (T) or not (F) 
     977/ 
     978!----------------------------------------------------------------------- 
     979&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d") 
     980!----------------------------------------------------------------------- 
     981   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
     982/ 
    850983!!====================================================================== 
    851984!!                  ***  Diagnostics namelists  *** 
     
    8981031&namptr       !   Poleward Transport Diagnostic 
    8991032!----------------------------------------------------------------------- 
    900    ln_diaptr  = .true.     !  Poleward heat and salt transport (T) or not (F) 
     1033   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    9011034   ln_diaznl  = .true.     !  Add zonal means and meridional stream functions 
    9021035   ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not 
     
    10211154&namsbc_wave   ! External fields from wave model 
    10221155!----------------------------------------------------------------------- 
    1023 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    1024 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    1025    sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' , .true.   , .false. , 'daily'  ,''         , '' 
    1026    sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     , .true.   , .false. , 'daily'  ,''         , '' 
    1027    sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     , .true.   , .false. , 'daily'  ,''         , '' 
    1028    sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   , .true.   , .false. , 'daily'  ,''         , '' 
     1156!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     1157!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     1158   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     1159   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     1160   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     1161   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    10291162! 
    10301163   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
Note: See TracChangeset for help on using the changeset viewer.