New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5034 for branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_ref – NEMO

Ignore:
Timestamp:
2015-01-15T14:48:42+01:00 (9 years ago)
Author:
andrewryan
Message:

merge with trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_ref

    r4384 r5034  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  :  1 - run manager      (namrun, namcfg) 
    3 !! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
     2!! NEMO/OPA  :  1 - run manager      (namrun) 
     3!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd) 
    44!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
     
    1818!!                   ***  Run management namelists  *** 
    1919!!====================================================================== 
    20 !!   namrun        parameters of the run 
     20!!   namrun       parameters of the run 
    2121!!====================================================================== 
    2222! 
     
    3131   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3232   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 
     33   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T 
     34   nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T 
    3535                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    3636                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     
    4747/ 
    4848! 
    49 !----------------------------------------------------------------------- 
    50 &namcfg     !   default parameters of the configuration       
     49!!====================================================================== 
     50!!                      ***  Domain namelists  *** 
     51!!====================================================================== 
     52!!   namcfg       parameters of the configuration       
     53!!   namzgr       vertical coordinate 
     54!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     55!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     56!!   namtsd       data: temperature & salinity 
     57!!====================================================================== 
     58! 
     59!----------------------------------------------------------------------- 
     60&namcfg     !   parameters of the configuration       
    5161!----------------------------------------------------------------------- 
    5262   cp_cfg      =  "default"            !  name of the configuration 
    53    cp_cfz      =         ''            !  name of the zoom of configuration 
     63   cp_cfz      =  "no zoom"            !  name of the zoom of configuration 
    5464   jp_cfg      =       0               !  resolution of the configuration 
    5565   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
     
    5767   jpkdta      =      31               !  number of levels      ( >= jpk ) 
    5868   jpiglo      =      10               !  1st dimension of global domain --> i =jpidta 
    59    jpjglo      =      12               !  2nd    -                  -    --> j  =jpjdta 
     69   jpjglo      =      12               !  2nd    -                  -    --> j =jpjdta 
    6070   jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
    6171   jpjzoom     =       1               !  in data domain indices 
     
    6777                                       !  = 6 cyclic East-West AND North fold F-point pivot 
    6878/ 
    69 !!====================================================================== 
    70 !!                      ***  Domain namelists  *** 
    71 !!====================================================================== 
    72 !!   namzgr       vertical coordinate 
    73 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    74 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    75 !!   namtsd       data: temperature & salinity 
    76 !!====================================================================== 
    77 ! 
    7879!----------------------------------------------------------------------- 
    7980&namzgr        !   vertical coordinate 
     
    8283   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    8384   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     85   ln_isfcav   = .false.   !  ice shelf cavity 
    8486/ 
    8587!----------------------------------------------------------------------- 
     
    183185/ 
    184186!----------------------------------------------------------------------- 
     187&namc1d        !   1D configuration options                             ("key_c1d") 
     188!----------------------------------------------------------------------- 
     189   rn_lat1d    =      50   !  Column latitude (default at PAPA station) 
     190   rn_lon1d    =    -145   !  Column longitude (default at PAPA station) 
     191   ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F) 
     192/ 
     193!----------------------------------------------------------------------- 
    185194&namtsd    !   data : Temperature  & Salinity 
    186195!----------------------------------------------------------------------- 
     
    208217!!   namtra_qsr      penetrative solar radiation 
    209218!!   namsbc_rnf      river runoffs 
     219!!   namsbc_isf      ice shelf melting/freezing 
    210220!!   namsbc_apr      Atmospheric Pressure 
    211221!!   namsbc_ssr      sea surface restoring term (for T and/or S) 
     
    223233   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    224234   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    225    ln_cpl      = .false.   !  Coupled formulation                       (T => fill namsbc_cpl ) 
    226235   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    227236   nn_ice      = 2         !  =0 no ice boundary condition   , 
    228237                           !  =1 use observed ice-cover      , 
    229                            !  =2 ice-model used                         ("key_lim3" or "key_lim2) 
     238                           !  =2 ice-model used                         ("key_lim3" or "key_lim2") 
    230239   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    231240                           !  =1 levitating ice with mass and salt exchange but no presure effect 
    232241                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    233242   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    234    ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     243   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf) 
     244   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
     245                           !  0 =no isf                  1 = presence of ISF  
     246                           !  2 = bg03 parametrisation   3 = rnf file for isf    
     247                           !  4 = ISF fwf specified 
     248                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    235249   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    236250   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     
    240254   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave) 
    241255   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. 
     256   nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
     257                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
     258   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
     259                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
     260                           !  = 0  Average per-category fluxes (forced and coupled mode) 
     261                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
     262                           !  = 2  Redistribute a single flux over categories (coupled mode only) 
    243263/ 
    244264!----------------------------------------------------------------------- 
     
    296316 
    297317   cn_dir      = './'      !  root directory for the location of the bulk files 
    298    ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
    299318   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) 
     319   rn_zqt      = 10.        !  Air temperature and humidity reference height (m) 
     320   rn_zu       = 10.        !  Wind vector reference height (m)                  
    303321   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    304322   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     
    327345!                    !                       ! categories !  reference  !    orientation       ! grids  ! 
    328346! send 
    329 sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    330 sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    331 sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
    332 sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    333 sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     347   sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     348   sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     349   sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     350   sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     351   sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    334352! receive 
    335 sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    336 sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    337 sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    338 sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    339 sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    340 sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    341 sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    342 sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    343 sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    344 sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     353   sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     354   sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     355   sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     356   sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     357   sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     358   sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     359   sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     360   sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     361   sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     362   sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     363! 
     364   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     365   ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models 
     366                             !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    345367/ 
    346368!----------------------------------------------------------------------- 
     
    396418   ln_rnf_tem   = .false.   !  read in temperature information for runoff 
    397419   ln_rnf_sal   = .false.   !  read in salinity information for runoff 
     420/ 
     421!----------------------------------------------------------------------- 
     422&namsbc_isf    !  Top boundary layer (ISF)  
     423!----------------------------------------------------------------------- 
     424!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     425!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
     426! nn_isf == 4 
     427   sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     428   sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     429! nn_isf == 3 
     430   sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     431! nn_isf == 2 and 3 
     432   sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     433   sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     434! nn_isf == 2 
     435   sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     436! for all case 
     437   ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 
     438! only for nn_isf = 1 or 2 
     439   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula 
     440   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula 
     441! only for nn_isf = 1 
     442   nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 
     443   rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008) 
     444                          ! 0 => thickness of the tbl = thickness of the first wet cell 
     445   ln_conserve = .true.   ! conservative case (take into account meltwater advection) 
     446   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s) 
     447                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     448                          !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 
     449                          ! 2 = velocity and stability dependent Gamma    Holland et al. 1999 
    398450/ 
    399451!----------------------------------------------------------------------- 
     
    559611    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
    560612                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     613    cn_ice_lim      =  'none'             !   
     614    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state 
     615                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     616    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice 
     617    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           -- 
     618    rn_ice_age      = 30.                 !  lim3 only:      --   age                -- 
     619 
    561620    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
    562621    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
     
    579638   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    580639   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     640! for lim2 
     641!   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     642!   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     643!   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     644! for lim3 
     645!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     646!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     647!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    581648   cn_dir  =    'bdydta/' 
    582649   ln_full_vel = .false. 
     
    607674   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    608675   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T  
    609    ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    610676   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    611677   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     678   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case) 
     679   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
     680   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
     681   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2) 
     682   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T 
     683   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file ) 
     684   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T) 
     685 
    612686   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
     687   ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    613688/ 
    614689!----------------------------------------------------------------------- 
     
    635710!!   nameos        equation of state 
    636711!!   namtra_adv    advection scheme 
     712!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    637713!!   namtra_ldf    lateral diffusion scheme 
    638714!!   namtra_dmp    T & S newtonian damping 
     
    642718&nameos        !   ocean physical parameters 
    643719!----------------------------------------------------------------------- 
    644    nn_eos      =   0       !  type of equation of state and Brunt-Vaisala frequency 
    645                            !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
    646                            !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    647                            !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    648    rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
    649    rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     720   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency 
     721                                 !  =-1, TEOS-10  
     722                                 !  = 0, EOS-80  
     723                                 !  = 1, S-EOS   (simplified eos) 
     724   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     725   !                             ! 
     726   !                     ! S-EOS coefficients : 
     727   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     728   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
     729   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     730   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
     731   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     732   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
     733   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
     734   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
    650735/ 
    651736!----------------------------------------------------------------------- 
    652737&namtra_adv    !   advection scheme for tracer 
    653738!----------------------------------------------------------------------- 
    654    ln_traadv_cen2   =  .false.  !  2nd order centered scheme 
    655    ln_traadv_tvd    =  .true.   !  TVD scheme 
    656    ln_traadv_muscl  =  .false.  !  MUSCL scheme 
    657    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    658    ln_traadv_ubs    =  .false.  !  UBS scheme 
    659    ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    660    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     739   ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
     740   ln_traadv_tvd    =  .true.    !  TVD scheme 
     741   ln_traadv_muscl  =  .false.   !  MUSCL scheme 
     742   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
     743   ln_traadv_ubs    =  .false.   !  UBS scheme 
     744   ln_traadv_qck    =  .false.   !  QUICKEST scheme 
     745   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
     746   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection 
    661747/ 
    662748!----------------------------------------------------------------------- 
     
    733819   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    734820   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     821   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection 
    735822/ 
    736823!----------------------------------------------------------------------- 
     
    755842   ln_dynvor_mix = .false. !  mixed scheme 
    756843   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     844   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation 
    757845/ 
    758846!----------------------------------------------------------------------- 
     
    9141002!!                  ***  Miscellaneous namelists  *** 
    9151003!!====================================================================== 
     1004!!   namsol            elliptic solver / island / free surface 
    9161005!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    9171006!!   namctl            Control prints & Benchmark 
    918 !!   namsol            elliptic solver / island / free surface 
     1007!!   namc1d            1D configuration options                         ("key_c1d") 
     1008!!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
     1009!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
    9191010!!====================================================================== 
    9201011! 
     
    9591050/ 
    9601051!----------------------------------------------------------------------- 
    961 &namc1d        !   1D configuration options                             ("key_c1d") 
    962 !----------------------------------------------------------------------- 
    963    rn_lat      =    50     !  Column latitude 
    964    rn_lon      =    -145   !  Column longitude 
    965 / 
    966 !----------------------------------------------------------------------- 
    9671052&namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
    9681053!----------------------------------------------------------------------- 
     
    9811066   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    9821067/ 
     1068 
    9831069!!====================================================================== 
    9841070!!                  ***  Diagnostics namelists  *** 
    9851071!!====================================================================== 
    9861072!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    987 !!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
     1073!!   namtrd       dynamics and/or tracer trends 
    9881074!!   namflo       float parameters                                      ("key_float") 
    9891075!!   namptr       Poleward Transport Diagnostics 
     
    10031089/ 
    10041090!----------------------------------------------------------------------- 
    1005 &namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra") 
    1006 !              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor") 
    1007 !----------------------------------------------------------------------- 
    1008    nn_trd      = 365       !  time step frequency dynamics and tracers trends 
    1009    nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    1010    rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    1011    cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
    1012    cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
    1013    ln_trdmld_restart = .false.         !  restart for ML diagnostics 
    1014    ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    1015 / 
     1091&namtrd        !   diagnostics on dynamics and/or tracer trends 
     1092!              !       and/or mixed-layer trends and/or barotropic vorticity 
     1093!----------------------------------------------------------------------- 
     1094   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
     1095   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
     1096   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
     1097   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet) 
     1098   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
     1099   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
     1100   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output 
     1101   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
     1102   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
     1103/ 
     1104!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     1105!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     1106!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
     1107!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
     1108!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
     1109!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
     1110!!gm 
    10161111!----------------------------------------------------------------------- 
    10171112&namflo       !   float parameters                                      ("key_float") 
     
    10851180                           !     ln_ssh                  Logical switch for SSH observations 
    10861181 
    1087    ln_sst     = .true.     ! Logical switch for SST observations 
    1088    ln_reysst  = .true.     !     ln_reysst               Logical switch for Reynolds observations 
     1182   ln_sst     = .false.     ! Logical switch for SST observations 
     1183   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations 
    10891184   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations       
    10901185 
    10911186   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
    10921187                           !     ln_sss                  Logical switch for SSS observations 
    1093                            !     ln_seaice              Logical switch for Sea Ice observations 
     1188   ln_seaice  = .false.    ! Logical switch for Sea Ice observations 
    10941189                           !     ln_vel3d                Logical switch for velocity observations 
    10951190                           !     ln_velavcur             Logical switch for velocity daily av. cur. 
     
    11121207                           !     sstfiles                GHRSST input observation file name 
    11131208   !                       ! sstfbfiles: Feedback SST input observation file name 
    1114    sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    1115                            !     seaicefiles             Sea Ice input observation file name 
     1209   sstfbfiles = 'sst_01.nc' 
     1210                           !     seaicefiles             Sea Ice input observation file names 
     1211   seaicefiles = 'seaice_01.nc'   
    11161212                           !     velavcurfiles           Vel. cur. daily av. input file name 
    11171213                           !     velhvcurfiles           Vel. cur. high freq. input file name 
Note: See TracChangeset for help on using the changeset viewer.