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 3815 – NEMO

Changeset 3815


Ignore:
Timestamp:
2013-02-18T20:42:54+01:00 (11 years ago)
Author:
vichi
Message:

Put same namelist of GYRE in GYRE_BFM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_MERGE_2012/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist

    r3795 r3815  
    33!! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
    44!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core 
    5 !!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,  
     5!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    77!!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     
    2424&namrun        !   parameters of the run 
    2525!----------------------------------------------------------------------- 
    26    nn_no       =       0   !  job number 
    27    cn_exp      = "GYRE_BFM"!  experience name  
     26   nn_no       =       0   !  job number (no more used...) 
     27   cn_exp      =  "GYRE"   !  experience name 
    2828   nn_it000    =       1   !  first time step 
    2929   nn_itend    =    4320   !  last  time step 
    30    nn_date0    =  010101   !  initial calendar date yymmdd (used if nn_rstctl=1) 
     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) 
    3131   nn_leapy    =      30   !  Leap year calendar (1) or not (0) 
    3232   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    33    nn_rstctl   =       0       !  restart control = 0 nn_it000 is not compared to the restart file value 
    34                                !                  = 1 use nn_date0 in namelist (not the value in the restart file) 
    35                                !                  = 2 calendar parameters read in the restart file 
     33   nn_rstctl   =       0   !  restart control => activated only if ln_rstart = T 
     34                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
     35                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     36                           !    = 2 nn_date0 read in restart  ; nn_it000 : check consistancy between namelist and restart 
    3637   cn_ocerst_in  = "restart"   !  suffix of ocean restart name (input) 
    3738   cn_ocerst_out = "restart"   !  suffix of ocean restart name (output) 
    3839   nn_istate   =       0   !  output the initial state (1) or not (0) 
    3940   nn_stock    =    4320   !  frequency of creation of a restart file (modulo referenced to 1) 
    40    nn_write    =     120   !  frequency of write in the output file   (modulo referenced to nn_it000) 
     41   nn_write    =      60   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    4142   ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4243   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
     
    5152!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    5253!!   namdom       space and time domain (bathymetry, mesh, timestep) 
    53 !!   namtsd       data: temperature & salinity          
     54!!   namtsd       data: temperature & salinity 
    5455!!====================================================================== 
    5556! 
     
    6465&namzgr_sco    !   s-coordinate or hybrid z-s-coordinate 
    6566!----------------------------------------------------------------------- 
    66    rn_sbot_min =  300.     !  minimum depth of s-bottom surface (>0) (m) 
    67    rn_sbot_max = 5250.     !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    68    rn_theta    =    6.0    !  surface control parameter (0<=rn_theta<=20) 
    69    rn_thetb    =    0.75   !  bottom control parameter  (0<=rn_thetb<= 1) 
    70    rn_rmax     =    0.15   !  maximum cut-off r-value allowed (0<rn_max<1) 
    71    ln_s_sigma  = .false.   !  hybrid s-sigma coordinates 
    72    rn_bb       =    0.8    !  stretching with s-sigma 
    73    rn_hc       =  150.0    !  critical depth with s-sigma  
     67   ln_s_sh94   = .true.    !  Song & Haidvogel 1994 hybrid S-sigma   (T)| 
     68   ln_s_sf12   = .false.   !  Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied 
     69   ln_sigcrit  = .false.   !  use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch 
     70                           !  stretching coefficients for all functions 
     71   rn_sbot_min =   10.0    !  minimum depth of s-bottom surface (>0) (m) 
     72   rn_sbot_max = 7000.0    !  maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
     73   rn_hc       =  150.0    !  critical depth for transition to stretched coordinates 
     74                        !!!!!!!  Envelop bathymetry 
     75   rn_rmax     =    0.3    !  maximum cut-off r-value allowed (0<r_max<1) 
     76                        !!!!!!!  SH94 stretching coefficients  (ln_s_sh94 = .true.) 
     77   rn_theta    =    6.0    !  surface control parameter (0<=theta<=20) 
     78   rn_bb       =    0.8    !  stretching with SH94 s-sigma    
     79                        !!!!!!!  SF12 stretching coefficient  (ln_s_sf12 = .true.) 
     80   rn_alpha    =    4.4    !  stretching with SF12 s-sigma 
     81   rn_efold    =    0.0    !  efold length scale for transition to stretched coord 
     82   rn_zs       =    1.0    !  depth of surface grid box 
     83                           !  bottom cell depth (Zb) is a linear function of water depth Zb = H*a + b 
     84   rn_zb_a     =    0.024  !  bathymetry scaling factor for calculating Zb 
     85   rn_zb_b     =   -0.2    !  offset for calculating Zb 
     86                        !!!!!!!! Other stretching (not SH94 or SF12) [also uses rn_theta above] 
     87   rn_thetb    =    1.0    !  bottom control parameter  (0<=thetb<= 1)  
    7488/ 
    7589!----------------------------------------------------------------------- 
     
    93107/ 
    94108!----------------------------------------------------------------------- 
    95 &namtsd    !   data : Temperature  & Salinity                            
     109&namtsd    !   data : Temperature  & Salinity 
    96110!----------------------------------------------------------------------- 
    97111!          ! file name ! frequency (hours)    ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     
    124138&namsbc        !   Surface Boundary Condition (surface module) 
    125139!----------------------------------------------------------------------- 
    126    nn_fsbc     = 1         !  frequency of surface boundary condition computation  
     140   nn_fsbc     = 1         !  frequency of surface boundary condition computation 
    127141                           !     (also = the frequency of sea-ice model call) 
    128    ln_ana      = .true.    !  analytical formulation                    (T => fill namsbc_ana )  
     142   ln_ana      = .true.    !  analytical formulation                    (T => fill namsbc_ana ) 
    129143   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    130    ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio)  
    131    ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core)  
     144   ln_blk_clio = .false.   !  CLIO bulk formulation                     (T => fill namsbc_clio) 
     145   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
    132146   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    133147   ln_cpl      = .false.   !  Coupled formulation                       (T => fill namsbc_cpl ) 
     
    136150                           !  =1 use observed ice-cover      , 
    137151                           !  =2 ice-model used                         ("key_lim3" or "key_lim2) 
     152   nn_ice_embd = 0         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
     153                           !  =1 levitating ice with mass and salt exchange but no presure effect 
     154                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    138155   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    139156   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    140157   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    141    nn_fwb      = 0         !  FreshWater Budget: =0 unchecked  
    142                            !     =1 global mean of e-p-r set to zero at each time step  
     158   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
     159                           !     =1 global mean of e-p-r set to zero at each time step 
    143160                           !     =2 annual global mean of e-p-r set to zero 
    144161                           !     =3 global emp set to zero and spread out over erp area 
     162   ln_wave = .false.       !  Activate coupling with wave (either Stokes Drift or Drag coefficient, or both)  (T => fill namsbc_wave) 
    145163   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave) 
     164   ln_sdw  = .false.       !  Computation of 3D stokes drift                (T => fill namsbc_wave) 
    146165/ 
    147166!----------------------------------------------------------------------- 
     
    167186 
    168187   cn_dir      = './'      !  root directory for the location of the flux files 
    169 /       
     188/ 
    170189!----------------------------------------------------------------------- 
    171190&namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
     
    224243!                    !                       ! categories !  reference  !    orientation       ! grids  ! 
    225244! send 
    226 sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   ''     
    227 sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   ''     
    228 sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   ''     
    229 sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T'        
    230 sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   ''         
     245sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     246sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     247sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     248sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     249sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    231250! receive 
    232 sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   ''     
    233 sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''     
    234 sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V'    
    235 sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''     
    236 sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''     
    237 sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   ''     
    238 sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   ''     
    239 sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''     
    240 sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''     
    241 sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   ''     
     251sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     252sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     253sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     254sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     255sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     256sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     257sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     258sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     259sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     260sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    242261/ 
    243262!----------------------------------------------------------------------- 
     
    247266!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    248267   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , '' 
    249   
     268 
    250269   cn_dir      = './'      !  root directory for the location of the runoff files 
    251270   ln_traqsr   = .true.    !  Light penetration (T) or not (F) 
    252271   ln_qsr_rgb  = .false.   !  RGB (Red-Green-Blue) light penetration 
    253272   ln_qsr_2bd  = .true.    !  2 bands              light penetration 
    254    ln_qsr_bio  = .true.   !  bio-model light penetration 
     273   ln_qsr_bio  = .false.   !  bio-model light penetration 
    255274   nn_chldta   =      0    !  RGB : Chl data (=1) or cst value (=0) 
    256275   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1) 
     
    298317   sn_sst      = 'sst_data'  ,        24         ,  'sst'    ,    .false.   , .false., 'yearly'  , ''       , '' 
    299318   sn_sss      = 'sss_data'  ,        -1         ,  'sss'    ,    .true.    , .true. , 'yearly'  , ''       , '' 
    300      
     319 
    301320   cn_dir      = './'      !  root directory for the location of the runoff files 
    302321   nn_sstr     =     0     !  add a retroaction term in the surface heat       flux (=1) or not (=0) 
     
    307326   ln_sssr_bnd =   .false. !  flag to bound erp term (associated with nn_sssr=2) 
    308327   rn_sssr_bnd =   4.e0    !  ABS(Max/Min) value of the damping erp term [mm/day] 
    309 /       
     328/ 
    310329!----------------------------------------------------------------------- 
    311330&namsbc_alb    !   albedo parameters 
    312331!----------------------------------------------------------------------- 
    313    rn_cloud    =    0.06   !  cloud correction to snow and ice albedo  
     332   rn_cloud    =    0.06   !  cloud correction to snow and ice albedo 
    314333   rn_albice   =    0.53   !  albedo of melting ice in the arctic and antarctic 
    315334   rn_alphd    =    0.80   !  coefficients for linear interpolation used to 
    316    rn_alphc    =    0.65   !  compute albedo between two extremes values  
     335   rn_alphc    =    0.65   !  compute albedo between two extremes values 
    317336   rn_alphdi   =    0.72   !  (Pyane, 1972) 
     337/ 
     338!----------------------------------------------------------------------- 
     339&namberg       !   iceberg parameters 
     340!----------------------------------------------------------------------- 
     341      ln_icebergs              = .false. 
     342      ln_bergdia               = .true.               ! Calculate budgets 
     343      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
     344      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
     345      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     346                                                      ! Initial mass required for an iceberg of each class 
     347      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     348                                                      ! Proportion of calving mass to apportion to each class 
     349      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     350                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     351                                                      ! i.e. number of icebergs represented at a point 
     352      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     353                                                      ! thickness of newly calved bergs (m) 
     354      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     355      rn_rho_bergs             = 850.                 ! Density of icebergs 
     356      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
     357      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
     358      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
     359      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     360      ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
     361      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
     362                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     363      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
     364      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
     365 
     366               ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     367               !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
     368      sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
     369 
     370      cn_dir = './' 
    318371/ 
    319372 
     
    324377!!   namcla        cross land advection 
    325378!!   namobc        open boundaries parameters                           ("key_obc") 
    326 !!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
     379!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
    327380!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
    328381!!   namtide       Tidal forcing at open boundaries                     ("key_bdy_tides") 
     
    346399   ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    347400   ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    348    ln_obc_fla  = .false.   !  Flather open boundary condition  
     401   ln_obc_fla  = .false.   !  Flather open boundary condition 
    349402   nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    350403                           !  = 1 the obc data are read in 'obc.dta' files 
     
    395448&nambdy        !  unstructured open boundaries                          ("key_bdy") 
    396449!----------------------------------------------------------------------- 
    397     nb_bdy = 1                            !  number of open boundary sets        
     450    nb_bdy = 1                            !  number of open boundary sets 
    398451    ln_coords_file = .true.               !  =T : read bdy coordinates from file 
    399452    cn_coords_file = 'coordinates.bdy.nc' !  bdy coordinates files 
     
    431484/ 
    432485!----------------------------------------------------------------------- 
    433 &nambdy_tide     ! tidal forcing at open boundaries               
     486&nambdy_tide     ! tidal forcing at open boundaries 
    434487!----------------------------------------------------------------------- 
    435488   filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     
    470523!!====================================================================== 
    471524!!   nambfr        bottom friction 
    472 !!   nambbc        bottom temperature boundary condition                
     525!!   nambbc        bottom temperature boundary condition 
    473526!!   nambbl        bottom boundary layer scheme                         ("key_trabbl") 
    474527!!====================================================================== 
     
    490543!----------------------------------------------------------------------- 
    491544   ln_trabbc   = .false.   !  Apply a geothermal heating at the ocean bottom 
    492    nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux  
     545   nn_geoflx   =    0      !  geothermal heat flux: = 0 no flux 
    493546                           !     = 1 constant flux 
    494                            !     = 2 variable flux (read in geothermal_heating.nc in mW/m2)  
     547                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    495548   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    496549/ 
     
    510563!!   namtra_adv    advection scheme 
    511564!!   namtra_ldf    lateral diffusion scheme 
    512 !!   namtra_dmp    T & S newtonian damping                         
     565!!   namtra_dmp    T & S newtonian damping 
    513566!!====================================================================== 
    514567! 
     
    524577/ 
    525578!----------------------------------------------------------------------- 
    526 &namtra_adv    !   advection scheme for tracer  
    527 !----------------------------------------------------------------------- 
    528    ln_traadv_cen2   =  .false.  !  2nd order centered scheme    
    529    ln_traadv_tvd    =  .true.   !  TVD scheme                 
    530    ln_traadv_muscl  =  .false.  !  MUSCL scheme              
    531    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries   
    532    ln_traadv_ubs    =  .false.  !  UBS scheme                  
    533    ln_traadv_qck    =  .false.  !  QUICKEST scheme                  
    534 / 
    535 !----------------------------------------------------------------------- 
    536 &namtra_ldf    !   lateral diffusion scheme for tracer  
    537 !----------------------------------------------------------------------- 
    538    !                       !  Type of the operator :  
    539    ln_traldf_lap    =  .true.   !  laplacian operator        
    540    ln_traldf_bilap  =  .false.  !  bilaplacian operator      
    541    !                       !  Direction of action  : 
    542    ln_traldf_level  =  .false.  !  iso-level                 
    543    ln_traldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" when ln_sco=T) 
    544    ln_traldf_iso    =  .true.   !  iso-neutral                          (require "key_ldfslp") 
    545    ln_traldf_grif   =  .false.  !  griffies skew flux formulation       (require "key_ldfslp") 
    546    ln_traldf_gdia   =  .false.  !  griffies operator strfn diagnostics  (require "key_ldfslp") 
    547    ln_triad_iso     =  .false.  !  griffies operator calculates triads twice => pure lateral mixing in ML (require "key_ldfslp") 
    548    ln_botmix_grif   =  .false.  !  griffies operator with lateral mixing on bottom (require "key_ldfslp") 
    549    !                       !  Coefficient 
     579&namtra_adv    !   advection scheme for tracer 
     580!----------------------------------------------------------------------- 
     581   ln_traadv_cen2   =  .false.  !  2nd order centered scheme 
     582   ln_traadv_tvd    =  .true.   !  TVD scheme 
     583   ln_traadv_muscl  =  .false.  !  MUSCL scheme 
     584   ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
     585   ln_traadv_ubs    =  .false.  !  UBS scheme 
     586   ln_traadv_qck    =  .false.  !  QUICKEST scheme 
     587   ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     588/ 
     589!---------------------------------------------------------------------------------- 
     590&namtra_ldf    !   lateral diffusion scheme for tracers 
     591!---------------------------------------------------------------------------------- 
     592   !                       !  Operator type: 
     593   ln_traldf_lap    =  .true.   !  laplacian operator 
     594   ln_traldf_bilap  =  .false.  !  bilaplacian operator 
     595   !                       !  Direction of action: 
     596   ln_traldf_level  =  .false.  !  iso-level 
     597   ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
     598   ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp") 
     599   !                 !  Griffies parameters              (all need "key_ldfslp") 
     600   ln_traldf_grif   =  .false.  !  use griffies triads 
     601   ln_traldf_gdia   =  .false.  !  output griffies eddy velocities 
     602   ln_triad_iso     =  .false.  !  pure lateral mixing in ML 
     603   ln_botmix_grif   =  .false.  !  lateral mixing on bottom 
     604   !                       !  Coefficients 
     605   ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 
     606   ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 
     607   !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 
     608   rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    550609   rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    551610   rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    552    rn_aeiv_0        =     0.    !  eddy induced velocity coef. [m2/s]   (require "key_traldf_eiv") 
    553 / 
    554 !----------------------------------------------------------------------- 
    555 &namtra_dmp    !   tracer: T & S newtonian damping                       
     611   !                                           (normally=0; not used with Griffies) 
     612/ 
     613!----------------------------------------------------------------------- 
     614&namtra_dmp    !   tracer: T & S newtonian damping 
    556615!----------------------------------------------------------------------- 
    557616   ln_tradmp   =  .false.  !  add a damping termn (T) or not (F) 
     
    581640&namdyn_adv    !   formulation of the momentum advection 
    582641!----------------------------------------------------------------------- 
    583    ln_dynadv_vec = .true.  !  vector form (T) or flux form (F)   
     642   ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    584643   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    585    ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme  
    586 /   
     644   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     645/ 
    587646!----------------------------------------------------------------------- 
    588647&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    589648!----------------------------------------------------------------------- 
    590    ln_dynvor_ene = .true.  !  enstrophy conserving scheme   
    591    ln_dynvor_ens = .false. !  energy conserving scheme     
    592    ln_dynvor_mix = .false. !  mixed scheme                
    593    ln_dynvor_een = .false. !  energy & enstrophy scheme   
     649   ln_dynvor_ene = .true.  !  energy    conserving scheme   
     650   ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     651   ln_dynvor_mix = .false. !  mixed scheme 
     652   ln_dynvor_een = .false. !  energy & enstrophy scheme 
    594653/ 
    595654!----------------------------------------------------------------------- 
    596655&namdyn_hpg    !   Hydrostatic pressure gradient option 
    597656!----------------------------------------------------------------------- 
    598    ln_hpg_zco  = .true.    !  z-coordinate - full steps                    
     657   ln_hpg_zco  = .true.    !  z-coordinate - full steps 
    599658   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
    600659   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
     
    614673&namdyn_ldf    !   lateral diffusion on momentum 
    615674!----------------------------------------------------------------------- 
    616    !                       !  Type of the operator :  
    617    ln_dynldf_lap    =  .true.   !  laplacian operator          
    618    ln_dynldf_bilap  =  .false.  !  bilaplacian operator     
    619    !                       !  Direction of action  :  
    620    ln_dynldf_level  =  .false.  !  iso-level                
     675   !                       !  Type of the operator : 
     676   ln_dynldf_lap    =  .true.   !  laplacian operator 
     677   ln_dynldf_bilap  =  .false.  !  bilaplacian operator 
     678   !                       !  Direction of action  : 
     679   ln_dynldf_level  =  .false.  !  iso-level 
    621680   ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    622681   ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     
    697756&namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    698757!------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    699    ln_kpprimix = .true.    !  shear instability mixing  
     758   ln_kpprimix = .true.    !  shear instability mixing 
    700759   rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s] 
    701760   rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
    702761   rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability 
    703762   rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s] 
    704    rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2]  
    705    rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s]  
     763   rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2] 
     764   rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s] 
    706765   nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
    707766   nn_ave      =  1        !  constant (=0) or profile (=1) background on avt 
     
    737796   rn_n2min    = 1.e-8     !  threshold of the Brunt-Vaisala frequency (s-1) 
    738797   rn_tfe      = 0.333     !  tidal dissipation efficiency 
    739    rn_me       = 0.2       !  mixing efficiency  
     798   rn_me       = 0.2       !  mixing efficiency 
    740799   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    741800   rn_tfe_itf  = 1.        !  ITF tidal dissipation efficiency 
     
    747806!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    748807!!   namctl            Control prints & Benchmark 
    749 !!   namsol            elliptic solver / island / free surface  
    750 !!====================================================================== 
    751 ! 
    752 !----------------------------------------------------------------------- 
    753 &namsol        !   elliptic solver / island / free surface  
     808!!   namsol            elliptic solver / island / free surface 
     809!!====================================================================== 
     810! 
     811!----------------------------------------------------------------------- 
     812&namsol        !   elliptic solver / island / free surface 
    754813!----------------------------------------------------------------------- 
    755814   nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
     
    770829   nn_buffer   =   0       !  size in bytes of exported buffer ('B' case), 0 no exportation 
    771830   ln_nnogather=  .false.  !  activate code to avoid mpi_allgather use at the northfold 
    772    jpni        =   2       !  jpni   number of processors following i (set automatically if < 1)      
    773    jpnj        =   2       !  jpnj   number of processors following j (set automatically if < 1)      
    774    jpnij       =   4       !  jpnij  number of local domains (set automatically if < 1)               
     831   jpni        =   0       !  jpni   number of processors following i (set automatically if < 1) 
     832   jpnj        =   0       !  jpnj   number of processors following j (set automatically if < 1) 
     833   jpnij       =   0       !  jpnij  number of local domains (set automatically if < 1) 
    775834/ 
    776835!----------------------------------------------------------------------- 
     
    787846   nn_bench    =    0      !  Bench mode (1/0): CAUTION use zero except for bench 
    788847                           !     (no physical validity of the results) 
    789    nn_timing   =    1      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
     848   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    790849/ 
    791850 
     
    797856!!   namflo       float parameters                                      ("key_float") 
    798857!!   namptr       Poleward Transport Diagnostics 
    799 !!   namhsb       Heat and salt budgets  
     858!!   namhsb       Heat and salt budgets 
    800859!!====================================================================== 
    801860! 
     
    842901   ln_diaptr  = .false.    !  Poleward heat and salt transport (T) or not (F) 
    843902   ln_diaznl  = .true.     !  Add zonal means and meridional stream functions 
    844    ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not  
     903   ln_subbas  = .true.     !  Atlantic/Pacific/Indian basins computation (T) or not 
    845904                           !  (orca configuration only, need input basins mask file named "subbasins.nc" 
    846905   ln_ptrcomp = .true.     !  Add decomposition : overturning 
     
    849908/ 
    850909!----------------------------------------------------------------------- 
    851 &namhsb       !  Heat and salt budgets  
     910&namhsb       !  Heat and salt budgets 
    852911!----------------------------------------------------------------------- 
    853912   ln_diahsb  = .false.    !  check the heat and salt budgets (T) or not (F) 
     
    861920                           !     -1 : debug all section 
    862921                           !  0 < n : debug section number n 
    863 /  
     922/ 
    864923 
    865924!!====================================================================== 
     
    873932&namobs       !  observation usage switch                               ('key_diaobs') 
    874933!----------------------------------------------------------------------- 
    875    ln_t3d     = .false.    ! Logical switch for T profile observations          
    876    ln_s3d     = .false.    ! Logical switch for S profile observations           
    877    ln_ena     = .false.    ! Logical switch for ENACT insitu data set            
    878    !                       !     ln_cor                  Logical switch for Coriolis insitu data set        
    879    ln_profb   = .false.    ! Logical switch for feedback insitu data set      
    880    ln_sla     = .false.    ! Logical switch for SLA observations                
    881  
    882    ln_sladt   = .false.    ! Logical switch for AVISO SLA data               
    883  
    884    ln_slafb   = .false.    ! Logical switch for feedback SLA data             
    885                            !     ln_ssh                  Logical switch for SSH observations               
    886  
    887    ln_sst     = .false.    ! Logical switch for SST observations               
    888                            !     ln_reysst               Logical switch for Reynolds observations        
    889                            !     ln_ghrsst               Logical switch for GHRSST observations           
    890  
    891    ln_sstfb   = .false.    ! Logical switch for feedback SST data           
    892                            !     ln_sss                  Logical switch for SSS observations               
    893                            !     ln_seaice               Logical switch for Sea Ice observations         
    894                            !     ln_vel3d                Logical switch for velocity observations          
    895                            !     ln_velavcur             Logical switch for velocity daily av. cur.     
    896                            !     ln_velhrcur             Logical switch for velocity high freq. cur.    
    897                            !     ln_velavadcp            Logical switch for velocity daily av. ADCP   
     934   ln_t3d     = .false.    ! Logical switch for T profile observations 
     935   ln_s3d     = .false.    ! Logical switch for S profile observations 
     936   ln_ena     = .false.    ! Logical switch for ENACT insitu data set 
     937   !                       !     ln_cor                  Logical switch for Coriolis insitu data set 
     938   ln_profb   = .false.    ! Logical switch for feedback insitu data set 
     939   ln_sla     = .false.    ! Logical switch for SLA observations 
     940 
     941   ln_sladt   = .false.    ! Logical switch for AVISO SLA data 
     942 
     943   ln_slafb   = .false.    ! Logical switch for feedback SLA data 
     944                           !     ln_ssh                  Logical switch for SSH observations 
     945 
     946   ln_sst     = .true.     ! Logical switch for SST observations               
     947   ln_reysst  = .true.     !     ln_reysst               Logical switch for Reynolds observations        
     948   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations           
     949 
     950   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
     951                           !     ln_sss                  Logical switch for SSS observations 
     952                           !     ln_seaice               Logical switch for Sea Ice observations 
     953                           !     ln_vel3d                Logical switch for velocity observations 
     954                           !     ln_velavcur             Logical switch for velocity daily av. cur. 
     955                           !     ln_velhrcur             Logical switch for velocity high freq. cur. 
     956                           !     ln_velavadcp            Logical switch for velocity daily av. ADCP 
    898957                           !     ln_velhradcp            Logical switch for velocity high freq. ADCP 
    899                            !     ln_velfb                Logical switch for feedback velocity data        
    900                            !     ln_grid_global          Global distribtion of observations          
    901                            !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table   
    902                            !     grid_search_file        Grid search lookup file header  
    903                            !     enactfiles              ENACT input observation file names  
    904                            !     coriofiles              Coriolis input observation file name   
    905    !                       ! profbfiles: Profile feedback input observation file name  
     958                           !     ln_velfb                Logical switch for feedback velocity data 
     959                           !     ln_grid_global          Global distribtion of observations 
     960                           !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table 
     961                           !     grid_search_file        Grid search lookup file header 
     962                           !     enactfiles              ENACT input observation file names 
     963                           !     coriofiles              Coriolis input observation file name 
     964   !                       ! profbfiles: Profile feedback input observation file name 
    906965   profbfiles = 'profiles_01.nc' 
    907                            !     ln_profb_enatim         Enact feedback input time setting switch     
     966                           !     ln_profb_enatim         Enact feedback input time setting switch 
    908967                           !     slafilesact             Active SLA input observation file name 
    909                            !     slafilespas             Passive SLA input observation file name  
    910    !                       ! slafbfiles: Feedback SLA input observation file name  
     968                           !     slafilespas             Passive SLA input observation file name 
     969   !                       ! slafbfiles: Feedback SLA input observation file name 
    911970   slafbfiles = 'sla_01.nc' 
    912                            !     sstfiles                GHRSST input observation file name        
    913    !                       ! sstfbfiles: Feedback SST input observation file name  
     971                           !     sstfiles                GHRSST input observation file name 
     972   !                       ! sstfbfiles: Feedback SST input observation file name 
    914973   sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    915                            !     seaicefiles             Sea Ice input observation file name  
    916                            !     velavcurfiles           Vel. cur. daily av. input file name   
    917                            !     velhvcurfiles           Vel. cur. high freq. input file name   
    918                            !     velavadcpfiles          Vel. ADCP daily av. input file name     
    919                            !     velhvadcpfiles          Vel. ADCP high freq. input file name  
    920                            !     velfbfiles              Vel. feedback input observation file name  
    921                            !     dobsini                 Initial date in window YYYYMMDD.HHMMSS        
    922                            !     dobsend                 Final date in window YYYYMMDD.HHMMSS          
    923                            !     n1dint                  Type of vertical interpolation method         
    924                            !     n2dint                  Type of horizontal interpolation method        
    925                            !     ln_nea                  Rejection of observations near land switch     
    926    nmsshc     = 0          ! MSSH correction scheme                          
    927                            !     mdtcorr                 MDT  correction                                
    928                            !     mdtcutoff               MDT cutoff for computed correction           
    929    ln_altbias = .false.    ! Logical switch for alt bias                 
    930    ln_ignmis  = .true.     ! Logical switch for ignoring missing files    
    931                            !     endailyavtypes   ENACT daily average types                     
     974                           !     seaicefiles             Sea Ice input observation file name 
     975                           !     velavcurfiles           Vel. cur. daily av. input file name 
     976                           !     velhvcurfiles           Vel. cur. high freq. input file name 
     977                           !     velavadcpfiles          Vel. ADCP daily av. input file name 
     978                           !     velhvadcpfiles          Vel. ADCP high freq. input file name 
     979                           !     velfbfiles              Vel. feedback input observation file name 
     980                           !     dobsini                 Initial date in window YYYYMMDD.HHMMSS 
     981                           !     dobsend                 Final date in window YYYYMMDD.HHMMSS 
     982                           !     n1dint                  Type of vertical interpolation method 
     983                           !     n2dint                  Type of horizontal interpolation method 
     984                           !     ln_nea                  Rejection of observations near land switch 
     985   nmsshc     = 0          ! MSSH correction scheme 
     986                           !     mdtcorr                 MDT  correction 
     987                           !     mdtcutoff               MDT cutoff for computed correction 
     988   ln_altbias = .false.    ! Logical switch for alt bias 
     989   ln_ignmis  = .true.     ! Logical switch for ignoring missing files 
     990                           !     endailyavtypes   ENACT daily average types 
    932991   ln_grid_global = .true. 
    933992   ln_grid_search_lookup = .false. 
    934 /  
     993/ 
    935994!----------------------------------------------------------------------- 
    936995&nam_asminc   !   assimilation increments                               ('key_asminc') 
    937996!----------------------------------------------------------------------- 
    938     ln_bkgwri = .false.    !  Logical switch for writing out background state  
     997    ln_bkgwri = .false.    !  Logical switch for writing out background state 
    939998    ln_trainc = .false.    !  Logical switch for applying tracer increments 
    940999    ln_dyninc = .false.    !  Logical switch for applying velocity increments 
    941     ln_sshinc = .false.    !  Logical switch for applying SSH increments  
     1000    ln_sshinc = .false.    !  Logical switch for applying SSH increments 
    9421001    ln_asmdin = .false.    !  Logical switch for Direct Initialization (DI) 
    9431002    ln_asmiau = .false.    !  Logical switch for Incremental Analysis Updating (IAU) 
     
    9571016!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    9581017   sn_cdg      =  'cdg_wave' ,        1          , 'drag_coeff' , .true.   , .false. , 'daily'  ,''         , '' 
     1018   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     , .true.   , .false. , 'daily'  ,''         , '' 
     1019   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     , .true.   , .false. , 'daily'  ,''         , '' 
     1020   sn_wn       =  'sdw_wave' ,        1          , 'wave_num'   , .true.   , .false. , 'daily'  ,''         , '' 
    9591021! 
    9601022   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
Note: See TracChangeset for help on using the changeset viewer.