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 6043 for branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG – NEMO

Ignore:
Timestamp:
2015-12-14T10:27:28+01:00 (8 years ago)
Author:
timgraham
Message:

Merged head of trunk into branch

Location:
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG
Files:
37 edited
1 copied

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r5600 r6043  
    7272/ 
    7373!----------------------------------------------------------------------- 
    74 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    75 !----------------------------------------------------------------------- 
    76    ln_bt_nn_auto =    .FALSE.          !  Set nn_baro automatically to be just below 
    77                                        !  a user defined maximum courant number (rn_bt_cmax) 
    78    nn_baro       =    30               !  Number of iterations of barotropic mode 
    79 / 
    80 !----------------------------------------------------------------------- 
    8174&namcrs        !   Grid coarsening for dynamics output and/or 
    8275               !   passive tracer coarsened online simulations 
     
    194187   rn_shlat    =     0     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    195188                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    196 / 
    197 !----------------------------------------------------------------------- 
    198 &namcla        !   cross land advection 
    199 !----------------------------------------------------------------------- 
    200 / 
    201 !----------------------------------------------------------------------- 
    202 &namobc        !   open boundaries parameters                           ("key_obc") 
    203 !----------------------------------------------------------------------- 
    204189/ 
    205190!----------------------------------------------------------------------- 
     
    281266&namtra_adv    !   advection scheme for tracer 
    282267!----------------------------------------------------------------------- 
     268   ln_traadv_fct =  .true.   !  FCT scheme 
     269      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     270      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     271      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     272      !                             !        (number of sub-timestep = nn_fct_zts) 
    283273/ 
    284274!----------------------------------------------------------------------- 
     
    289279&namtra_ldf    !   lateral diffusion scheme for tracers 
    290280!---------------------------------------------------------------------------------- 
    291    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    292    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    293    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    294    rn_aht_0         =    50.    !  horizontal eddy diffusivity for tracers [m2/s] 
     281   !                       !  Operator type: 
     282   ln_traldf_lap   =  .true.   !    laplacian operator 
     283   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     284   !                       !  Direction of action: 
     285   ln_traldf_lev   =  .false.  !  iso-level 
     286   ln_traldf_hor   =  .true.   !  horizontal (geopotential) 
     287   ln_traldf_iso   =  .false.  !  iso-neutral 
     288   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     289   ! 
     290   !                       !  iso-neutral options:         
     291   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     292   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     293   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     294   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     295   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     296   ! 
     297   !                       !  Coefficients: 
     298   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     299   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     300   !                                !   =  0           constant  
     301   !                                !   = 10 F(k)      =ldf_c1d  
     302   !                                !   = 20 F(i,j)    =ldf_c2d  
     303   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     304   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     305   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     306   rn_aht_0        =   50.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     307   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     308/ 
     309!---------------------------------------------------------------------------------- 
     310&namtra_ldfeiv !   eddy induced velocity param. 
     311!---------------------------------------------------------------------------------- 
     312   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    295313/ 
    296314!----------------------------------------------------------------------- 
     
    306324&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    307325!----------------------------------------------------------------------- 
     326   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     327   ln_dynvor_ens = .false. !  energy conserving scheme 
     328   ln_dynvor_mix = .false. !  mixed scheme 
     329   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     330      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    308331/ 
    309332!----------------------------------------------------------------------- 
     
    314337/ 
    315338!----------------------------------------------------------------------- 
    316 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    317 !----------------------------------------------------------------------- 
    318 !                          !  explicit free surface                     ("key_dynspg_exp") 
    319 !                          !  filtered free surface                     ("key_dynspg_flt") 
    320 !                          !  split-explicit free surface               ("key_dynspg_ts") 
    321  
     339&namdyn_spg    !   surface pressure gradient  
     340!----------------------------------------------------------------------- 
     341   ln_dynspg_ts   = .true.    ! split-explicit free surface 
     342      ln_bt_auto    = .false.   ! Number of sub-step defined from: 
     343         nn_baro      = 30        ! =F : the number of sub-step in rn_rdt seconds 
     344/ 
    322345!----------------------------------------------------------------------- 
    323346&namdyn_ldf    !   lateral diffusion on momentum 
    324347!----------------------------------------------------------------------- 
    325348   !                       !  Type of the operator : 
    326    ln_dynldf_bilap  =  .true.   !  bilaplacian operator 
    327    ln_dynldf_lap    =  .false.  !  bilaplacian operator 
     349   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     350   ln_dynldf_lap =  .false.    !    laplacian operator 
     351   ln_dynldf_blp =  .true.     !  bilaplacian operator 
    328352   !                       !  Direction of action  : 
    329    ln_dynldf_level  =  .true.   !  iso-level 
    330    ln_dynldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    331                            !  Coefficient 
    332    rn_ahm_0_lap     = 60.0      !  horizontal laplacian eddy viscosity   [m2/s] 
    333    rn_ahm_0_blp     = -1.0e+10  !  horizontal bilaplacian eddy viscosity [m4/s] 
     353   ln_dynldf_lev =  .true.     !  iso-level 
     354   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     355   ln_dynldf_iso =  .false.    !  iso-neutral 
     356   !                       !  Coefficient 
     357   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     358   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     359   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     360   !                                !  =  0  constant  
     361   !                                !  = 10  F(k)=c1d 
     362   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     363   !                                !  = 30  F(i,j,k)=c2d*c1d 
     364   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     365   rn_ahm_0      =     60.     !  horizontal laplacian eddy viscosity   [m2/s] 
     366   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     367   rn_bhm_0      = 1.0e+10     !  horizontal bilaplacian eddy viscosity [m4/s] 
     368   ! 
     369   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    334370/ 
    335371!----------------------------------------------------------------------- 
     
    349385!----------------------------------------------------------------------- 
    350386/ 
    351 !------------------------------------------------------------------------ 
    352 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    353 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    354 / 
    355387!----------------------------------------------------------------------- 
    356388&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     
    369401/ 
    370402!----------------------------------------------------------------------- 
    371 &namsol        !   elliptic solver / island / free surface 
    372 !----------------------------------------------------------------------- 
    373 / 
    374 !----------------------------------------------------------------------- 
    375403&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    376404!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm

    r4245 r6043  
    1  bld::tool::fppkeys  key_bdy key_tide key_dynspg_ts key_ldfslp key_zdfgls  key_vvl key_diainstant key_mpp_mpi key_iomput 
     1 bld::tool::fppkeys  key_bdy key_tide key_zdfgls  key_vvl key_diainstant key_mpp_mpi key_iomput 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r5600 r6043  
    6464/ 
    6565!----------------------------------------------------------------------- 
    66 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    67 !----------------------------------------------------------------------- 
    68 / 
    69 !----------------------------------------------------------------------- 
    7066&namcrs        !   Grid coarsening for dynamics output and/or 
    7167               !   passive tracer coarsened online simulations 
     
    141137&namtra_qsr    !   penetrative solar radiation 
    142138!----------------------------------------------------------------------- 
     139!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     140!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     141   sn_chl      ='chlorophyll_PAPASTATION', -1    , 'CHLA'    ,   .true.     , .true. , 'yearly'  , ''       , ''       , '' 
    143142/ 
    144143!----------------------------------------------------------------------- 
     
    170169!----------------------------------------------------------------------- 
    171170   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    172 / 
    173 !----------------------------------------------------------------------- 
    174 &namcla        !   cross land advection 
    175 !----------------------------------------------------------------------- 
    176 / 
    177 !----------------------------------------------------------------------- 
    178 &namobc        !   open boundaries parameters                           ("key_obc") 
    179 !----------------------------------------------------------------------- 
    180171/ 
    181172!----------------------------------------------------------------------- 
     
    226217&namtra_adv    !   advection scheme for tracer 
    227218!----------------------------------------------------------------------- 
     219! C1D : no advection scheme  
    228220/ 
    229221!----------------------------------------------------------------------- 
     
    231223!----------------------------------------------------------------------- 
    232224/ 
    233 !---------------------------------------------------------------------------------- 
     225!----------------------------------------------------------------------- 
    234226&namtra_ldf    !   lateral diffusion scheme for tracers 
    235 !---------------------------------------------------------------------------------- 
    236 !---------------------------------------------------------------------------------- 
    237    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    238    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    239    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    240    rn_aht_0         =     0.    !  horizontal eddy diffusivity for tracers [m2/s] 
     227!----------------------------------------------------------------------- 
     228! C1D : no lateral diffusion   
     229/ 
     230!----------------------------------------------------------------------- 
     231&namtra_ldfeiv !   eddy induced velocity param. 
     232!----------------------------------------------------------------------- 
     233! C1D : no eiv   
    241234/ 
    242235!----------------------------------------------------------------------- 
     
    248241&namdyn_adv    !   formulation of the momentum advection 
    249242!----------------------------------------------------------------------- 
     243! C1D : no advection scheme  
    250244/ 
    251245!----------------------------------------------------------------------- 
     
    266260/ 
    267261!----------------------------------------------------------------------- 
    268 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    269 !----------------------------------------------------------------------- 
     262&namdyn_spg    !   surface pressure gradient 
     263!----------------------------------------------------------------------- 
     264/ 
    270265!----------------------------------------------------------------------- 
    271266&namdyn_ldf    !   lateral diffusion on momentum 
     
    286281!----------------------------------------------------------------------- 
    287282/ 
    288 !------------------------------------------------------------------------ 
    289 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    291 / 
    292283!----------------------------------------------------------------------- 
    293284&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     
    302293!----------------------------------------------------------------------- 
    303294   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    304 / 
    305 !----------------------------------------------------------------------- 
    306 &namsol        !   elliptic solver / island / free surface 
    307 !----------------------------------------------------------------------- 
    308    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    309    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    310    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    311295/ 
    312296!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r5600 r6043  
    6666/ 
    6767!----------------------------------------------------------------------- 
    68 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    69 !----------------------------------------------------------------------- 
    70 / 
    71 !----------------------------------------------------------------------- 
    7268&namcrs        !   Grid coarsening for dynamics output and/or 
    7369               !   passive tracer coarsened online simulations 
     
    154150!----------------------------------------------------------------------- 
    155151   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    156 / 
    157 !----------------------------------------------------------------------- 
    158 &namcla        !   cross land advection 
    159 !----------------------------------------------------------------------- 
    160 / 
    161 !----------------------------------------------------------------------- 
    162 &namobc        !   open boundaries parameters                           ("key_obc") 
    163 !----------------------------------------------------------------------- 
    164152/ 
    165153!----------------------------------------------------------------------- 
     
    223211&namtra_adv    !   advection scheme for tracer 
    224212!----------------------------------------------------------------------- 
    225    ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
    226    ln_traadv_tvd    =  .true.    !  TVD scheme 
    227    ln_traadv_muscl  =  .false.   !  MUSCL scheme 
    228    ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
    229    ln_traadv_ubs    =  .false.   !  UBS scheme 
    230    ln_traadv_qck    =  .false.   !  QUICKEST scheme 
    231    ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
     213   ln_traadv_fct =  .true.   !  FCT scheme 
     214      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     215      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     216      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     217      !                             !        (number of sub-timestep = nn_fct_zts) 
    232218/ 
    233219!----------------------------------------------------------------------- 
     
    238224&namtra_ldf    !   lateral diffusion scheme for tracers 
    239225!---------------------------------------------------------------------------------- 
    240    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    241    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     226   !                       !  Operator type: 
     227   ln_traldf_lap   =  .true.   !    laplacian operator 
     228   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     229   !                       !  Direction of action: 
     230   ln_traldf_lev   =  .false.  !  iso-level 
     231   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     232   ln_traldf_iso   =  .true.   !  iso-neutral 
     233   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     234   ! 
     235   !                       !  iso-neutral options:         
     236   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
     237   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     238   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     239   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     240   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     241   ! 
     242   !                       !  Coefficients: 
     243   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     244   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     245   !                                !   =  0           constant  
     246   !                                !   = 10 F(k)      =ldf_c1d  
     247   !                                !   = 20 F(i,j)    =ldf_c2d  
     248   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     249   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     250   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     251   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     252   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     253/ 
     254!---------------------------------------------------------------------------------- 
     255&namtra_ldfeiv !   eddy induced velocity param. 
     256!---------------------------------------------------------------------------------- 
     257   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    242258/ 
    243259!----------------------------------------------------------------------- 
     
    253269&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    254270!----------------------------------------------------------------------- 
    255    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    256    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     271   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     272   ln_dynvor_ens = .false. !  energy conserving scheme 
     273   ln_dynvor_mix = .false. !  mixed scheme 
    257274   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     275      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    258276/ 
    259277!----------------------------------------------------------------------- 
     
    264282/ 
    265283!----------------------------------------------------------------------- 
    266 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    267 !----------------------------------------------------------------------- 
    268  
     284&namdyn_spg    !   surface pressure gradient  
     285!----------------------------------------------------------------------- 
     286   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     287/ 
    269288!----------------------------------------------------------------------- 
    270289&namdyn_ldf    !   lateral diffusion on momentum 
    271290!----------------------------------------------------------------------- 
    272    rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
     291   !                       !  Type of the operator : 
     292   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     293   ln_dynldf_lap =  .true.     !    laplacian operator 
     294   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     295   !                       !  Direction of action  : 
     296   ln_dynldf_lev =  .true.     !  iso-level 
     297   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     298   ln_dynldf_iso =  .false.    !  iso-neutral 
     299   !                       !  Coefficient 
     300   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     301   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     302   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     303   !                                !  =  0  constant  
     304   !                                !  = 10  F(k)=c1d 
     305   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     306   !                                !  = 30  F(i,j,k)=c2d*c1d 
     307   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     308   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     309   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     310   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
    273311/ 
    274312!----------------------------------------------------------------------- 
     
    286324   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    287325/ 
    288 !------------------------------------------------------------------------ 
    289 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    291 / 
    292326!----------------------------------------------------------------------- 
    293327&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     
    302336!----------------------------------------------------------------------- 
    303337   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    304 / 
    305 !----------------------------------------------------------------------- 
    306 &namsol        !   elliptic solver / island / free surface 
    307 !----------------------------------------------------------------------- 
    308    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    309    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    310    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    311338/ 
    312339!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

    r5034 r6043  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi key_nosignedzero 
     1 bld::tool::fppkeys key_zdftke key_iomput key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r5600 r6043  
    6969/ 
    7070!----------------------------------------------------------------------- 
    71 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    72 !----------------------------------------------------------------------- 
    73 / 
    74 !----------------------------------------------------------------------- 
    7571&namcrs        !   Grid coarsening for dynamics output and/or 
    7672               !   passive tracer coarsened online simulations 
     
    159155!----------------------------------------------------------------------- 
    160156   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    161 / 
    162 !----------------------------------------------------------------------- 
    163 &namcla        !   cross land advection 
    164 !----------------------------------------------------------------------- 
    165 / 
    166 !----------------------------------------------------------------------- 
    167 &namobc        !   open boundaries parameters                           ("key_obc") 
    168 !----------------------------------------------------------------------- 
    169157/ 
    170158!----------------------------------------------------------------------- 
     
    228216&namtra_adv    !   advection scheme for tracer 
    229217!----------------------------------------------------------------------- 
    230    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     218   ln_traadv_fct =  .true.   !  FCT scheme 
     219      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     220      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     221      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     222      !                             !        (number of sub-timestep = nn_fct_zts) 
    231223/ 
    232224!----------------------------------------------------------------------- 
     
    237229&namtra_ldf    !   lateral diffusion scheme for tracers 
    238230!---------------------------------------------------------------------------------- 
    239    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    240    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     231   !                       !  Operator type: 
     232   ln_traldf_lap   =  .true.   !    laplacian operator 
     233   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     234   !                       !  Direction of action: 
     235   ln_traldf_lev   =  .false.  !  iso-level 
     236   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     237   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     238   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     239   ! 
     240   !                       !  iso-neutral options:         
     241   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     242   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     243   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     244   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     245   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     246   ! 
     247   !                       !  Coefficients: 
     248   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     249   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     250   !                                !   =  0           constant  
     251   !                                !   = 10 F(k)      =ldf_c1d  
     252   !                                !   = 20 F(i,j)    =ldf_c2d  
     253   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     254   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     255   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     256   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     257   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     258/ 
     259!---------------------------------------------------------------------------------- 
     260&namtra_ldfeiv !   eddy induced velocity param. 
     261!---------------------------------------------------------------------------------- 
     262   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    241263/ 
    242264!----------------------------------------------------------------------- 
     
    252274&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    253275!----------------------------------------------------------------------- 
    254    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    255    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
    256    ln_dynvor_een = .false. !  energy & enstrophy scheme 
     276   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
    257277/ 
    258278!----------------------------------------------------------------------- 
     
    263283/ 
    264284!----------------------------------------------------------------------- 
    265 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    266 !----------------------------------------------------------------------- 
     285&namdyn_spg    !   surface pressure gradient 
     286!----------------------------------------------------------------------- 
     287   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     288/ 
    267289!----------------------------------------------------------------------- 
    268290&namdyn_ldf    !   lateral diffusion on momentum 
    269291!----------------------------------------------------------------------- 
     292   !                       !  Type of the operator : 
     293   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     294   ln_dynldf_lap =  .true.     !    laplacian operator 
     295   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     296   !                       !  Direction of action  : 
     297   ln_dynldf_lev =  .true.     !  iso-level 
     298   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     299   ln_dynldf_iso =  .false.    !  iso-neutral 
     300   !                       !  Coefficient 
     301   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     302   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     303   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     304   !                                !  =  0  constant  
     305   !                                !  = 10  F(k)=c1d 
     306   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     307   !                                !  = 30  F(i,j,k)=c2d*c1d 
     308   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     309   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     310   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     311   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
     312   ! 
     313   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     314/ 
    270315   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    271316/ 
     
    284329   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    285330/ 
    286 !------------------------------------------------------------------------ 
    287 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    288 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    289 / 
    290331!----------------------------------------------------------------------- 
    291332&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     
    300341!----------------------------------------------------------------------- 
    301342   ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    302 / 
    303 !----------------------------------------------------------------------- 
    304 &namsol        !   elliptic solver / island / free surface 
    305 !----------------------------------------------------------------------- 
    306    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    307    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    308    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    309343/ 
    310344!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_top_cfg

    r4152 r6043  
    2323!----------------------------------------------------------------------- 
    2424&namtrc_adv    !   advection scheme for passive tracer  
    25 !-----------------------------------------------------------------------   
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     25!----------------------------------------------------------------------- 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

    r4230 r6043  
    1 bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_vectopt_loop key_top key_my_trc key_mpp_mpi key_iomput 
     1bld::tool::fppkeys key_zdftke key_top key_my_trc key_mpp_mpi key_iomput 
    22inc $BFMDIR/src/nemo/bfm.fcm 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r5600 r6043  
    5959/ 
    6060!----------------------------------------------------------------------- 
    61 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    62 !----------------------------------------------------------------------- 
    63 / 
    64 !----------------------------------------------------------------------- 
    6561&namcrs        !   Grid coarsening for dynamics output and/or 
    6662               !   passive tracer coarsened online simulations 
     
    10298!----------------------------------------------------------------------- 
    10399   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    104 / 
    105 !----------------------------------------------------------------------- 
    106 &namcla        !   cross land advection 
    107 !----------------------------------------------------------------------- 
    108100/ 
    109101!----------------------------------------------------------------------- 
     
    143135&namtra_adv    !   advection scheme for tracer 
    144136!----------------------------------------------------------------------- 
    145    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     137   ln_traadv_fct =  .true.   !  FCT scheme 
     138      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     139      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     140      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     141      !                             !        (number of sub-timestep = nn_fct_zts) 
    146142/ 
    147143!---------------------------------------------------------------------------------- 
    148144&namtra_ldf    !   lateral diffusion scheme for tracers 
    149145!---------------------------------------------------------------------------------- 
    150    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    151    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     146   !                       !  Operator type: 
     147   ln_traldf_lap   =  .true.   !    laplacian operator 
     148   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     149   !                       !  Direction of action: 
     150   ln_traldf_lev   =  .false.  !  iso-level 
     151   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     152   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     153   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     154   ! 
     155   !                       !  iso-neutral options:         
     156   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     157   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     158   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     159   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     160   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     161   ! 
     162   !                       !  Coefficients: 
     163   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     164   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     165   !                                !   =  0           constant  
     166   !                                !   = 10 F(k)      =ldf_c1d  
     167   !                                !   = 20 F(i,j)    =ldf_c2d  
     168   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     169   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     170   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     171   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     172   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     173/ 
     174!---------------------------------------------------------------------------------- 
     175&namtra_ldfeiv !   eddy induced velocity param. 
     176!---------------------------------------------------------------------------------- 
     177   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    152178/ 
    153179!----------------------------------------------------------------------- 
     
    163189&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    164190!----------------------------------------------------------------------- 
    165    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    166    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     191   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     192   ln_dynvor_ens = .false. !  energy conserving scheme 
     193   ln_dynvor_mix = .false. !  mixed scheme 
    167194   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     195      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    168196/ 
    169197!----------------------------------------------------------------------- 
     
    174202/ 
    175203!----------------------------------------------------------------------- 
     204&namdyn_spg    !   surface pressure gradient 
     205!----------------------------------------------------------------------- 
     206   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     207/ 
     208!----------------------------------------------------------------------- 
    176209&namdyn_ldf    !   lateral diffusion on momentum 
    177210!----------------------------------------------------------------------- 
     211   !                       !  Type of the operator : 
     212   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     213   ln_dynldf_lap =  .true.    !    laplacian operator 
     214   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     215   !                       !  Direction of action  : 
     216   ln_dynldf_lev =  .true.    !  iso-level 
     217   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     218   ln_dynldf_iso =  .false.    !  iso-neutral 
     219   !                       !  Coefficient 
     220   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     221   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     222   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     223   !                                !  =  0  constant  
     224   !                                !  = 10  F(k)=c1d 
     225   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     226   !                                !  = 30  F(i,j,k)=c2d*c1d 
     227   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     228   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     229   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     230   rn_bhm_0      =      0.      !  horizontal bilaplacian eddy viscosity [m4/s] 
     231   ! 
     232   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     233/ 
    178234   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    179235/ 
     
    189245/ 
    190246!----------------------------------------------------------------------- 
    191 &namsol        !   elliptic solver / island / free surface 
    192 !----------------------------------------------------------------------- 
    193    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    194    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    195    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    196 / 
    197 !----------------------------------------------------------------------- 
    198247&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    199248!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_top_cfg

    r4340 r6043  
    2424&namtrc_adv    !   advection scheme for passive tracer  
    2525!-----------------------------------------------------------------------    
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r5034 r6043  
    1 bld::tool::fppkeys  key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_zdftke key_top key_pisces_reduced key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_XIOS/EXP00/namelist_cfg

    r5600 r6043  
    4141   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    4242   rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    43 !   nn_baro     =   60      !  number of barotropic time step            ("key_dynspg_ts") 
    4443   rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    4544   rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     
    150149/ 
    151150!----------------------------------------------------------------------- 
    152 &namcla        !   cross land advection 
    153 !----------------------------------------------------------------------- 
    154 / 
    155 !----------------------------------------------------------------------- 
    156 &namobc        !   open boundaries parameters                           ("key_obc") 
    157 !----------------------------------------------------------------------- 
    158 / 
    159 !----------------------------------------------------------------------- 
    160151&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    161152!----------------------------------------------------------------------- 
     
    200191&namtra_adv    !   advection scheme for tracer 
    201192!----------------------------------------------------------------------- 
    202    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     193   ln_traadv_fct =  .true.    !  FCT scheme 
     194      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     195      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     196      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     197      !                             !        (number of sub-timestep = nn_fct_zts) 
    203198/ 
    204199!----------------------------------------------------------------------- 
     
    209204&namtra_ldf    !   lateral diffusion scheme for tracers 
    210205!---------------------------------------------------------------------------------- 
    211    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    212    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     206   !                       !  Operator type: 
     207   ln_traldf_lap   =  .true.   !    laplacian operator 
     208   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     209   !                       !  Direction of action: 
     210   ln_traldf_lev   =  .false.  !  iso-level 
     211   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     212   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     213   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     214   ! 
     215   !                       !  iso-neutral options:         
     216   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     217   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     218   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     219   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     220   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     221   ! 
     222   !                       !  Coefficients: 
     223   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     224   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     225   !                                !   =  0           constant  
     226   !                                !   = 10 F(k)      =ldf_c1d  
     227   !                                !   = 20 F(i,j)    =ldf_c2d  
     228   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     229   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     230   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     231   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     232   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     233/ 
     234!---------------------------------------------------------------------------------- 
     235&namtra_ldfeiv !   eddy induced velocity param. 
     236!---------------------------------------------------------------------------------- 
     237   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    213238/ 
    214239!----------------------------------------------------------------------- 
     
    235260/ 
    236261!----------------------------------------------------------------------- 
    237 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    238 !----------------------------------------------------------------------- 
     262&namdyn_spg    !   surface pressure gradient 
     263!----------------------------------------------------------------------- 
     264   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     265/ 
    239266!----------------------------------------------------------------------- 
    240267&namdyn_ldf    !   lateral diffusion on momentum 
     
    274301/ 
    275302!----------------------------------------------------------------------- 
    276 &namsol        !   elliptic solver / island / free surface 
    277 !----------------------------------------------------------------------- 
    278    nn_solv     =      2    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    279    nn_nmin     =    210    !  minimum of iterations for the SOR solver 
    280    rn_sor      =  1.96     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    281 / 
    282 !----------------------------------------------------------------------- 
    283303&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    284304!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_XIOS/cpp_GYRE_XIOS.fcm

    r4373 r6043  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_zdftke key_iomput key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r5600 r6043  
    44!----------------------------------------------------------------------- 
    55&namrun        !   parameters of the run 
     6 nn_it000=1 
    67!----------------------------------------------------------------------- 
    78   cn_exp      = "Agulhas" !  experience name  
    8    nn_itend    =     480   !  last  time step 
     9   nn_itend    =   10950 
    910   nn_stock    =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
    1011   nn_write    =   10950   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    11    ln_clobber  = .true.    !  clobber (overwrite) an existing file 
     12   ln_clobber  =  .true. 
    1213/ 
    1314!----------------------------------------------------------------------- 
     
    5556/ 
    5657!----------------------------------------------------------------------- 
    57 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    58 !----------------------------------------------------------------------- 
    59 / 
    60 !----------------------------------------------------------------------- 
    6158&namcrs        !   Grid coarsening for dynamics output and/or 
    6259               !   passive tracer coarsened online simulations 
     
    116113/ 
    117114!----------------------------------------------------------------------- 
    118 &namcla        !   cross land advection 
    119 !----------------------------------------------------------------------- 
    120 / 
    121 !----------------------------------------------------------------------- 
    122115&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    123116!----------------------------------------------------------------------- 
     
    145138/ 
    146139!----------------------------------------------------------------------- 
    147 &namtra_adv    !   advection scheme for tracer  
    148 !----------------------------------------------------------------------- 
    149 / 
    150 !----------------------------------------------------------------------- 
    151 &namtra_ldf    !   lateral diffusion scheme for tracers  
    152 !----------------------------------------------------------------------- 
    153  
    154    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    155    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     140&namtra_adv    !   advection scheme for tracer 
     141!----------------------------------------------------------------------- 
     142   ln_traadv_fct =  .true.    !  FCT scheme 
     143      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order 
     144      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order 
     145      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     146      !                             !        (number of sub-timestep = nn_fct_zts) 
     147/ 
     148!----------------------------------------------------------------------- 
     149&namtra_ldf    !   lateral diffusion scheme for tracers 
     150!----------------------------------------------------------------------- 
     151   !                       !  Operator type: 
     152   ln_traldf_lap   =  .true.   !    laplacian operator 
     153   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     154   !                       !  Direction of action: 
     155   ln_traldf_lev   =  .false.  !  iso-level 
     156   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     157   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     158   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     159   ! 
     160   !                       !  iso-neutral options:         
     161   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     162   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     163   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     164   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     165   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     166   ! 
     167   !                       !  Coefficients: 
     168   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     169   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     170   !                                !   =  0           constant 
     171   !                                !   = 10 F(k)      =ldf_c1d 
     172   !                                !   = 20 F(i,j)    =ldf_c2d 
     173   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     174   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     175   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     176   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     177   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     178/ 
     179!---------------------------------------------------------------------------------- 
     180&namtra_ldfeiv !   eddy induced velocity param. 
     181!---------------------------------------------------------------------------------- 
     182   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    156183/ 
    157184!----------------------------------------------------------------------- 
     
    162189&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    163190!----------------------------------------------------------------------- 
     191   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     192   ln_dynvor_ens = .false. !  energy conserving scheme 
     193   ln_dynvor_mix = .false. !  mixed scheme 
     194   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     195      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    164196/ 
    165197!----------------------------------------------------------------------- 
     
    168200/ 
    169201!----------------------------------------------------------------------- 
     202&namdyn_spg    !   surface pressure gradient 
     203!----------------------------------------------------------------------- 
     204   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     205/ 
     206!----------------------------------------------------------------------- 
    170207&namdyn_ldf    !   lateral diffusion on momentum 
    171208!----------------------------------------------------------------------- 
    172    !                       !  Type of the operator :  
    173    ln_dynldf_lap    =  .false.  !  laplacian operator          
    174    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    175    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     209   !                       !  Type of the operator : 
     210   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F 
     211   ln_dynldf_lap =  .false.    !    laplacian operator 
     212   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     213   !                       !  Direction of action  : 
     214   ln_dynldf_lev =  .true.     !  iso-level 
     215   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     216   ln_dynldf_iso =  .false.    !  iso-neutral 
     217   !                       !  Coefficient 
     218   nn_ahm_ijk_t  =   0         !  space/time variation of eddy coef 
     219   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     220   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     221   !                                !  =  0  constant 
     222   !                                !  = 10  F(k)=c1d 
     223   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     224   !                                !  = 30  F(i,j,k)=c2d*c1d 
     225   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     226   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     227   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     228   rn_bhm_0      = 8.5e+11      !  horizontal bilaplacian eddy viscosity [m4/s] 
    176229/ 
    177230!----------------------------------------------------------------------- 
     
    193246/ 
    194247!----------------------------------------------------------------------- 
    195 &namsol        !   elliptic solver / island / free surface  
    196 !----------------------------------------------------------------------- 
    197 / 
    198 !----------------------------------------------------------------------- 
    199248&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    200249!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r5600 r6043  
    153153  </context> 
    154154   
     155 <context id="1_nemo" time_origin="1950-01-01 00:00:00" > 
     156     
     157    <!-- $id$ --> 
     158     
     159    <!--  
     160============================================================================================================ 
     161=                                  definition of all existing variables                                    = 
     162=                                            DO NOT CHANGE                                                 = 
     163============================================================================================================ 
     164    --> 
     165    <field_definition src="./field_def.xml"/> 
     166    <!--  
     167============================================================================================================ 
     168=                                           output files definition                                        = 
     169=                                            Define your own files                                         = 
     170=                                         put the variables you want...                                    = 
     171============================================================================================================ 
     172    --> 
     173     
     174    <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="10d" min_digits="4"> 
     175     
     176      <file_group id="1ts" output_freq="1ts"  output_level="10" enabled=".TRUE."/> <!-- 1 time step files --> 
     177      <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."/> <!-- 1h files --> 
     178      <file_group id="2h" output_freq="2h"  output_level="10" enabled=".TRUE."/> <!-- 2h files --> 
     179      <file_group id="3h" output_freq="3h"  output_level="10" enabled=".TRUE."/> <!-- 3h files -->      
     180      <file_group id="4h" output_freq="4h"  output_level="10" enabled=".TRUE."/> <!-- 4h files --> 
     181      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files -->      
     182      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."/> <!-- 1d files --> 
     183      <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files -->     
     184 
     185      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE." >  <!-- 5d files -->   
     186  
     187   <file id="file1" name_suffix="_grid_T" description="ocean T grid variables" > 
     188     <field field_ref="sst"          name="tos"      long_name="sea_surface_temperature"                       /> 
     189     <field field_ref="sss"          name="sos"      long_name="sea_surface_salinity"                          /> 
     190     <field field_ref="ssh"          name="zos"      long_name="sea_surface_height_above_geoid"                /> 
     191     <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"               /> 
     192     <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                            /> 
     193     <field field_ref="sst2"         name="tossq"    long_name="square_of_sea_surface_temperature"             /> 
     194     <field field_ref="ssh2"         name="zossq"    long_name="square_of_sea_surface_height_above_geoid"      /> 
     195     <field field_ref="mldkz5"       /> 
     196     <field field_ref="mldr10_1"     /> 
     197     <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     /> 
     198     <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           /> 
     199     <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          /> 
     200     <field field_ref="saltflx"      name="sosflxdo"  /> 
     201     <field field_ref="taum"         name="taum" /> 
     202     <field field_ref="wspd"         name="sowindsp"  /> 
     203          <field field_ref="precip"       name="soprecip" /> 
     204   </file> 
     205 
     206   <file id="file3" name_suffix="_grid_U" description="ocean U grid variables" > 
     207     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
     208     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
     209     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
     210          <!-- variables available with MLE 
     211          <field field_ref="psiu_mle"     name="psiu_mle"  long_name="MLE_streamfunction_along_i-axis" /> 
     212     --> 
     213   </file> 
     214    
     215   <file id="file4" name_suffix="_grid_V" description="ocean V grid variables" > 
     216     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
     217     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
     218     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
     219          <!-- variables available with MLE 
     220          <field field_ref="psiv_mle"     name="psiv_mle"  long_name="MLE_streamfunction_along_j-axis" /> 
     221     --> 
     222   </file> 
     223    
     224   <file id="file5" name_suffix="_grid_W" description="ocean W grid variables" > 
     225     <field field_ref="woce"         name="wo"      long_name="ocean vertical velocity"         /> 
     226     <field field_ref="avt"          name="difvho"  long_name="ocean_vertical_heat_diffusivity" /> 
     227   </file> 
     228<!-- 
     229   <file id="file6" name_suffix="_icemod" description="ice variables" > 
     230     <field field_ref="ice_pres"                     /> 
     231     <field field_ref="snowthic_cea" name="snd"     long_name="surface_snow_thickness"   /> 
     232     <field field_ref="icethic_cea"  name="sit"     long_name="sea_ice_thickness"        /> 
     233     <field field_ref="iceprod_cea"  name="sip"     long_name="sea_ice_thickness"        /> 
     234     <field field_ref="ist_ipa"      /> 
     235     <field field_ref="uice_ipa"     /> 
     236     <field field_ref="vice_ipa"     /> 
     237     <field field_ref="utau_ice"     /> 
     238     <field field_ref="vtau_ice"     /> 
     239     <field field_ref="qsr_io_cea"   /> 
     240     <field field_ref="qns_io_cea"   /> 
     241     <field field_ref="snowpre"      /> 
     242   </file> 
     243 
     244   <file id="file8" name_suffix="_Tides" description="tidal harmonics" > 
     245     <field field_ref="M2x"          name="M2x"      long_name="M2 Elevation harmonic real part"                       /> 
     246     <field field_ref="M2y"          name="M2y"      long_name="M2 Elevation harmonic imaginary part"                  /> 
     247     <field field_ref="M2x_u"        name="M2x_u"    long_name="M2 current barotrope along i-axis harmonic real part "       /> 
     248     <field field_ref="M2y_u"        name="M2y_u"    long_name="M2 current barotrope along i-axis harmonic imaginary part "  /> 
     249     <field field_ref="M2x_v"        name="M2x_v"    long_name="M2 current barotrope along j-axis harmonic real part "       /> 
     250     <field field_ref="M2y_v"        name="M2y_v"    long_name="M2 current barotrope along j-axis harmonic imaginary part "  /> 
     251   </file> 
     252   --> 
     253      </file_group> 
     254 
     255 
     256      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."/> <!-- real monthly files --> 
     257 
     258 
     259      <file_group id="2m" output_freq="2mo" output_level="10" enabled=".TRUE."/> <!-- real 2m files --> 
     260      <file_group id="3m" output_freq="3mo" output_level="10" enabled=".TRUE."/> <!-- real 3m files --> 
     261      <file_group id="4m" output_freq="4mo" output_level="10" enabled=".TRUE."/> <!-- real 4m files --> 
     262      <file_group id="6m" output_freq="6mo" output_level="10" enabled=".TRUE."/> <!-- real 6m files --> 
     263 
     264      <file_group id="1y"  output_freq="1y" output_level="10" enabled=".TRUE."/> <!-- real yearly files --> 
     265      <file_group id="2y"  output_freq="2y" output_level="10" enabled=".TRUE."/> <!-- real 2y files --> 
     266      <file_group id="5y"  output_freq="5y" output_level="10" enabled=".TRUE."/> <!-- real 5y files --> 
     267      <file_group id="10y" output_freq="10y" output_level="10" enabled=".TRUE."/> <!-- real 10y files --> 
     268 
     269   </file_definition> 
     270     
     271    <!--  
     272============================================================================================================ 
     273= grid definition = = DO NOT CHANGE = 
     274============================================================================================================ 
     275    --> 
     276     
     277   <axis_definition>   
     278      <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" /> 
     279      <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" /> 
     280      <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" /> 
     281      <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" /> 
     282      <axis id="nfloat" long_name="Float number"      unit="-"  /> 
     283      <axis id="icbcla" long_name="Iceberg class"     unit="-"  /> 
     284   </axis_definition>  
     285     
     286   <domain_definition src="./domain_def.xml"/> 
     287    
     288   <grid_definition>     
     289     <grid id="grid_T_2D" domain_ref="grid_T"/> 
     290     <grid id="grid_T_3D" domain_ref="grid_T" axis_ref="deptht"/> 
     291     <grid id="grid_U_2D" domain_ref="grid_U"/> 
     292     <grid id="grid_U_3D" domain_ref="grid_U" axis_ref="depthu"/> 
     293     <grid id="grid_V_2D" domain_ref="grid_V"/> 
     294     <grid id="grid_V_3D" domain_ref="grid_V" axis_ref="depthv"/> 
     295     <grid id="grid_W_2D" domain_ref="grid_W"/> 
     296     <grid id="grid_W_3D" domain_ref="grid_W" axis_ref="depthw"/> 
     297    </grid_definition>    
     298  </context> 
    155299 
    156300  <context id="xios"> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r5034 r6043  
    4848/ 
    4949!----------------------------------------------------------------------- 
    50 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    51 !----------------------------------------------------------------------- 
    52 / 
    53 !----------------------------------------------------------------------- 
    5450&namcrs        !   Grid coarsening for dynamics output and/or 
    5551               !   passive tracer coarsened online simulations 
     
    9389/ 
    9490!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9991&nambfr        !   bottom friction 
    10092!----------------------------------------------------------------------- 
     
    115107&namtra_adv    !   advection scheme for tracer 
    116108!----------------------------------------------------------------------- 
     109   ln_traadv_fct =  .true.    !  FCT scheme 
     110      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     111      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     112      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     113      !                             !        (number of sub-timestep = nn_fct_zts) 
    117114/ 
    118115!----------------------------------------------------------------------- 
     
    120117!----------------------------------------------------------------------- 
    121118/ 
    122 !----------------------------------------------------------------------- 
    123 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    124 !----------------------------------------------------------------------- 
    125 / 
    126119!---------------------------------------------------------------------------------- 
    127120&namtra_ldf    !   lateral diffusion scheme for tracers 
    128121!---------------------------------------------------------------------------------- 
     122   !                       !  Operator type: 
     123   ln_traldf_lap   =  .true.   !    laplacian operator 
     124   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     125   !                       !  Direction of action: 
     126   ln_traldf_lev   =  .false.  !  iso-level 
     127   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     128   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     129   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     130   ! 
     131   !                       !  iso-neutral options:         
     132   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     133   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     134   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     135   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     136   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     137   ! 
     138   !                       !  Coefficients: 
     139   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     140   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     141   !                                !   =  0           constant 
     142   !                                !   = 10 F(k)      =ldf_c1d 
     143   !                                !   = 20 F(i,j)    =ldf_c2d 
     144   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     145   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     146   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     147   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     148   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     149/ 
     150!---------------------------------------------------------------------------------- 
     151&namtra_ldfeiv !   eddy induced velocity param. 
     152!---------------------------------------------------------------------------------- 
     153   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     154   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     155   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     156   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     157   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     158   !                                !   =  0           constant 
     159   !                                !   = 10 F(k)      =ldf_c1d 
     160   !                                !   = 20 F(i,j)    =ldf_c2d 
     161   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     162   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    129163/ 
    130164!----------------------------------------------------------------------- 
     
    139173&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    140174!----------------------------------------------------------------------- 
     175   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     176   ln_dynvor_ens = .false. !  energy conserving scheme 
     177   ln_dynvor_mix = .false. !  mixed scheme 
     178   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     179      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    141180/ 
    142181!----------------------------------------------------------------------- 
     
    145184/ 
    146185!----------------------------------------------------------------------- 
     186&namdyn_spg    !   surface pressure gradient 
     187!----------------------------------------------------------------------- 
     188   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     189/ 
     190!----------------------------------------------------------------------- 
    147191&namdyn_ldf    !   lateral diffusion on momentum 
    148192!----------------------------------------------------------------------- 
     193   !                       !  Type of the operator : 
     194   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     195   ln_dynldf_lap =  .true.     !    laplacian operator 
     196   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     197   !                       !  Direction of action  : 
     198   ln_dynldf_lev =  .true.     !  iso-level 
     199   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     200   ln_dynldf_iso =  .false.    !  iso-neutral 
     201   !                       !  Coefficient 
     202   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     203   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     204   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     205   !                                !  =  0  constant  
     206   !                                !  = 10  F(k)=c1d 
     207   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     208   !                                !  = 30  F(i,j,k)=c2d*c1d 
     209   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     210   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     211   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     212   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    149213/ 
    150214!----------------------------------------------------------------------- 
     
    165229/ 
    166230!----------------------------------------------------------------------- 
    167 &namsol        !   elliptic solver / island / free surface 
    168 !----------------------------------------------------------------------- 
    169 / 
    170 !----------------------------------------------------------------------- 
    171231&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    172232!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r5600 r6043  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc  
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_iomput 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r5600 r6043  
    5555/ 
    5656!----------------------------------------------------------------------- 
    57 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    58 !----------------------------------------------------------------------- 
    59 / 
    60 !----------------------------------------------------------------------- 
    6157&namcrs        !   Grid coarsening for dynamics output and/or 
    6258               !   passive tracer coarsened online simulations 
     
    116112/ 
    117113!----------------------------------------------------------------------- 
    118 &namcla        !   cross land advection 
    119 !----------------------------------------------------------------------- 
    120 / 
    121 !----------------------------------------------------------------------- 
    122114&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    123115!----------------------------------------------------------------------- 
     
    145137/ 
    146138!----------------------------------------------------------------------- 
    147 &namtra_adv    !   advection scheme for tracer  
    148 !----------------------------------------------------------------------- 
    149 / 
    150 !----------------------------------------------------------------------- 
    151 &namtra_ldf    !   lateral diffusion scheme for tracers  
    152 !----------------------------------------------------------------------- 
    153  
    154    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    155    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     139&namtra_adv    !   advection scheme for tracer 
     140!----------------------------------------------------------------------- 
     141   ln_traadv_fct =  .true.    !  FCT scheme 
     142      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     143      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     144      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     145      !                             !        (number of sub-timestep = nn_fct_zts) 
     146/ 
     147!----------------------------------------------------------------------- 
     148&namtra_ldf    !   lateral diffusion scheme for tracers 
     149!---------------------------------------------------------------------------------- 
     150   !                       !  Operator type: 
     151   ln_traldf_lap   =  .true.   !    laplacian operator 
     152   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     153   !                       !  Direction of action: 
     154   ln_traldf_lev   =  .false.  !  iso-level 
     155   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     156   ln_traldf_iso   =  .true.   !  iso-neutral (Standard operator) 
     157   ln_traldf_triad =  .false.  !  iso-neutral (Triads   operator) 
     158   ! 
     159   !                       !  iso-neutral options:         
     160   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     161   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     162   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     163   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     164   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     165   ! 
     166   !                       !  Coefficients: 
     167   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     168   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     169   !                                !   =  0           constant  
     170   !                                !   = 10 F(k)      =ldf_c1d  
     171   !                                !   = 20 F(i,j)    =ldf_c2d  
     172   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     173   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     174   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     175   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     176   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     177/ 
     178!---------------------------------------------------------------------------------- 
     179&namtra_ldfeiv !   eddy induced velocity param. 
     180!---------------------------------------------------------------------------------- 
     181   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    156182/ 
    157183!----------------------------------------------------------------------- 
     
    162188&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    163189!----------------------------------------------------------------------- 
     190   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     191   ln_dynvor_ens = .false. !  energy conserving scheme 
     192   ln_dynvor_mix = .false. !  mixed scheme 
     193   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     194      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     195   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    164196/ 
    165197!----------------------------------------------------------------------- 
     
    168200/ 
    169201!----------------------------------------------------------------------- 
     202&namdyn_spg    !   surface pressure gradient 
     203!----------------------------------------------------------------------- 
     204   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     205/ 
     206!----------------------------------------------------------------------- 
    170207&namdyn_ldf    !   lateral diffusion on momentum 
    171208!----------------------------------------------------------------------- 
    172    !                       !  Type of the operator :  
    173    ln_dynldf_lap    =  .false.  !  laplacian operator          
    174    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    175    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     209   !                       !  Type of the operator : 
     210   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     211   ln_dynldf_lap =  .false.    !    laplacian operator 
     212   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     213   !                       !  Direction of action  : 
     214   ln_dynldf_lev =  .true.     !  iso-level 
     215   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     216   ln_dynldf_iso =  .false.    !  iso-neutral 
     217   !                       !  Coefficient 
     218   nn_ahm_ijk_t  = 20          !  space/time variation of eddy coef 
     219   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     220   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     221   !                                !  =  0  constant  
     222   !                                !  = 10  F(k)=c1d 
     223   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     224   !                                !  = 30  F(i,j,k)=c2d*c1d 
     225   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     226   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     227   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     228   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    176229/ 
    177230!----------------------------------------------------------------------- 
     
    193246/ 
    194247!----------------------------------------------------------------------- 
    195 &namsol        !   elliptic solver / island / free surface  
    196 !----------------------------------------------------------------------- 
    197 / 
    198 !----------------------------------------------------------------------- 
    199248&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    200249!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r5034 r6043  
    4848/ 
    4949!----------------------------------------------------------------------- 
    50 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    51 !----------------------------------------------------------------------- 
    52 / 
    53 !----------------------------------------------------------------------- 
    5450&namcrs        !   Grid coarsening for dynamics output and/or 
    5551               !   passive tracer coarsened online simulations 
     
    9389/ 
    9490!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9991&nambfr        !   bottom friction 
    10092!----------------------------------------------------------------------- 
     
    115107&namtra_adv    !   advection scheme for tracer 
    116108!----------------------------------------------------------------------- 
     109   ln_traadv_fct =  .true.    !  FCT scheme 
     110      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     111      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     112      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     113      !                             !        (number of sub-timestep = nn_fct_zts) 
    117114/ 
    118115!----------------------------------------------------------------------- 
     
    123120&namtra_ldf    !   lateral diffusion scheme for tracers 
    124121!---------------------------------------------------------------------------------- 
    125 / 
    126 !----------------------------------------------------------------------- 
     122   !                       !  Operator type: 
     123   ln_traldf_lap   =  .true.   !    laplacian operator 
     124   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     125   !                       !  Direction of action: 
     126   ln_traldf_lev   =  .false.  !  iso-level 
     127   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     128   ln_traldf_iso   =  .true.   !  iso-neutral (Standard operator) 
     129   ln_traldf_triad =  .false.  !  iso-neutral (Triads   operator) 
     130   ! 
     131   !                       !  iso-neutral options:         
     132   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     133   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     134   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     135   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     136   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     137   ! 
     138   !                       !  Coefficients: 
     139   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     140   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     141   !                                !   =  0           constant  
     142   !                                !   = 10 F(k)      =ldf_c1d  
     143   !                                !   = 20 F(i,j)    =ldf_c2d  
     144   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     145   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     146   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     147   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     148   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     149/ 
     150!---------------------------------------------------------------------------------- 
     151&namtra_ldfeiv !   eddy induced velocity param. 
     152!---------------------------------------------------------------------------------- 
     153   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     154   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     155   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     156   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     157   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     158   !                                !   =  0           constant  
     159   !                                !   = 10 F(k)      =ldf_c1d  
     160   !                                !   = 20 F(i,j)    =ldf_c2d  
     161   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     162   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     163/!----------------------------------------------------------------------- 
    127164&namtra_dmp    !   tracer: T & S newtonian damping 
    128165!----------------------------------------------------------------------- 
     
    135172&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    136173!----------------------------------------------------------------------- 
     174   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     175   ln_dynvor_ens = .false. !  energy conserving scheme 
     176   ln_dynvor_mix = .false. !  mixed scheme 
     177   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     178      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    137179/ 
    138180!----------------------------------------------------------------------- 
     
    144186   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    145187   ln_hpg_prj  = .true.    !  s-coordinate (Pressure Jacobian scheme) 
    146    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    147                                  !           centered      time scheme  (F) 
     188/ 
     189!----------------------------------------------------------------------- 
     190&namdyn_spg    !   surface pressure gradient 
     191!----------------------------------------------------------------------- 
     192   ln_dynspg_ts  = .true.  !  split-explicit free surface 
    148193/ 
    149194!----------------------------------------------------------------------- 
    150195&namdyn_ldf    !   lateral diffusion on momentum 
    151196!----------------------------------------------------------------------- 
     197   !                       !  Type of the operator : 
     198   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     199   ln_dynldf_lap =  .true.     !    laplacian operator 
     200   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     201   !                       !  Direction of action  : 
     202   ln_dynldf_lev =  .true.     !  iso-level 
     203   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     204   ln_dynldf_iso =  .false.    !  iso-neutral 
     205   !                       !  Coefficient 
     206   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     207   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     208   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     209   !                                !  =  0  constant  
     210   !                                !  = 10  F(k)=c1d 
     211   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     212   !                                !  = 30  F(i,j,k)=c2d*c1d 
     213   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     214   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     215   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     216   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     217   ! 
     218   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    152219/ 
    153220!----------------------------------------------------------------------- 
     
    168235/ 
    169236!----------------------------------------------------------------------- 
    170 &namsol        !   elliptic solver / island / free surface 
    171 !----------------------------------------------------------------------- 
    172 / 
    173 !----------------------------------------------------------------------- 
    174237&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    175238!----------------------------------------------------------------------- 
     
    191254!----------------------------------------------------------------------- 
    192255/ 
     256!----------------------------------------------------------------------- 
     257&namobs       !  observation usage                                      ('key_diaobs') 
     258!----------------------------------------------------------------------- 
     259/ 
     260!----------------------------------------------------------------------- 
     261&nam_asminc   !   assimilation increments                               ('key_asminc') 
     262!----------------------------------------------------------------------- 
     263/ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm

    r5600 r6043  
    1  bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc  
     1 bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_diaobs key_asminc key_iomput key_nosignedzero 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg

    r5600 r6043  
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,  
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     7!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    99!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     
    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    (nammpp, nammpp_dyndist, namctl) 
    1414!!             11 - Obs & Assim      (namobs, nam_asminc) 
    1515!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    8686                           ! 
    8787   rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    88    nn_baro     =   64      !  number of barotropic time step            ("key_dynspg_ts") 
    8988   rn_atfp     =    0.1    !  asselin time filter parameter 
    9089   nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
     
    302301!!====================================================================== 
    303302!!   namlbc        lateral momentum boundary condition 
    304 !!   namcla        cross land advection 
    305 !!   namobc        open boundaries parameters                           ("key_obc") 
    306303!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
    307304!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
     
    314311   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    315312                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 &namcla        !   cross land advection 
    319 !----------------------------------------------------------------------- 
    320    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    321 / 
    322 !----------------------------------------------------------------------- 
    323 &namobc        !   open boundaries parameters                           ("key_obc") 
    324 !----------------------------------------------------------------------- 
    325    ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    326    ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    327    ln_obc_fla  = .false.   !  Flather open boundary condition  
    328    nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    329                            !  = 1 the obc data are read in 'obc.dta' files 
    330    cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    331                            !  set to monthly if obc datafile hold 1 month of data 
    332    rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
    333    rn_dpwin    =    1.     !     -           -         -       west    -      - 
    334    rn_dpnin    =    1.     !     -           -         -       north   -      - 
    335    rn_dpsin    =    1.     !     -           -         -       south   -      - 
    336    rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
    337    rn_dpwob    =   15.     !     -           -         -     west    -      - 
    338    rn_dpnob    = 3000.     !     -           -         -     north   -      - 
    339    rn_dpsob    =   15.     !     -           -         -     south   -      - 
    340    rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    341                            !  = 1 the total volume remains constant 
    342313/ 
    343314!----------------------------------------------------------------------- 
     
    394365&nambdy_tide     ! tidal forcing at open boundaries               
    395366!----------------------------------------------------------------------- 
    396    filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     367    filtide       = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
    397368    tide_cpt(1)   ='Q1'  !  names of tidal components used 
    398369    tide_cpt(2)   ='O1'  !  names of tidal components used 
     
    485456/ 
    486457!----------------------------------------------------------------------- 
    487 &namtra_adv    !   advection scheme for tracer  
    488 !----------------------------------------------------------------------- 
    489    ln_traadv_cen2   =  .false.  !  2nd order centered scheme    
    490    ln_traadv_tvd    =  .true.   !  TVD scheme                 
    491    ln_traadv_muscl  =  .false.  !  MUSCL scheme              
    492    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries   
    493    ln_traadv_ubs    =  .false.  !  UBS scheme                  
    494    ln_traadv_qck    =  .false.  !  QUCIKEST scheme                  
     458&namtra_adv    !   advection scheme for tracer 
     459!----------------------------------------------------------------------- 
     460   ln_traadv_fct =  .true.    !  FCT scheme 
     461      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     462      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     463      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     464      !                             !        (number of sub-timestep = nn_fct_zts) 
    495465/ 
    496466!----------------------------------------------------------------------- 
    497467&namtra_ldf    !   lateral diffusion scheme for tracer  
    498468!----------------------------------------------------------------------- 
    499    !                       !  Type of the operator :  
    500    ln_traldf_lap    =  .true.   !  laplacian operator        
    501    ln_traldf_bilap  =  .false.  !  bilaplacian operator      
    502    !                       !  Direction of action  : 
    503    ln_traldf_level  =  .false.  !  iso-level                 
    504    ln_traldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" when ln_sco=T) 
    505    ln_traldf_iso    =  .true.   !  iso-neutral                          (require "key_ldfslp") 
    506    ln_traldf_grif   =  .false.  !  griffies skew flux formulation       (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    507    ln_traldf_gdia   =  .false.  !  griffies operator strfn diagnostics  (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    508    !                       !  Coefficient 
    509    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    510    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    511    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     469   !                       !  Operator type: 
     470   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     471   ln_traldf_lap   =  .true.   !    laplacian operator 
     472   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     473   !                       !  Direction of action: 
     474   ln_traldf_lev   =  .false.  !  iso-level 
     475   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     476   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     477   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     478   ! 
     479   !                       !  iso-neutral options:         
     480   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     481   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     482   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     483   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     484   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     485   ! 
     486   !                       !  Coefficients: 
     487   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     488   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     489   !                                !   =  0           constant  
     490   !                                !   = 10 F(k)      =ldf_c1d  
     491   !                                !   = 20 F(i,j)    =ldf_c2d  
     492   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     493   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     494   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     495   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     496   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     497/ 
     498!---------------------------------------------------------------------------------- 
     499&namtra_ldfeiv !   eddy induced velocity param. 
     500!---------------------------------------------------------------------------------- 
     501   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    512502/ 
    513503!----------------------------------------------------------------------- 
     
    546536&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    547537!----------------------------------------------------------------------- 
    548    ln_dynvor_ene = .false. !  enstrophy conserving scheme   
    549    ln_dynvor_ens = .false. !  energy conserving scheme     
    550    ln_dynvor_mix = .false. !  mixed scheme                
    551    ln_dynvor_een = .true.  !  energy & enstrophy scheme   
     538   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     539   ln_dynvor_ens = .false. !  energy conserving scheme 
     540   ln_dynvor_mix = .false. !  mixed scheme 
     541   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     542      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     543   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    552544/ 
    553545!----------------------------------------------------------------------- 
     
    559551   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    560552   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    561    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    562                                  !           centered      time scheme  (F) 
    563 / 
    564 !----------------------------------------------------------------------- 
    565 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    566 !----------------------------------------------------------------------- 
    567 !                          !  explicit free surface                     ("key_dynspg_exp") 
    568 !                          !  filtered free surface                     ("key_dynspg_flt") 
    569 !                          !  split-explicit free surface               ("key_dynspg_ts") 
    570  
     553/ 
     554!----------------------------------------------------------------------- 
     555&namdyn_spg    !   surface pressure gradient 
     556!----------------------------------------------------------------------- 
     557   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     558/ 
    571559!----------------------------------------------------------------------- 
    572560&namdyn_ldf    !   lateral diffusion on momentum 
    573561!----------------------------------------------------------------------- 
    574    !                       !  Type of the operator :  
    575    ln_dynldf_lap    =  .false.  !  laplacian operator          
    576    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    577    !                       !  Direction of action  :  
    578    ln_dynldf_level  =  .false.  !  iso-level                
    579    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    580    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     562   !                       !  Type of the operator : 
     563   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     564   ln_dynldf_lap =  .false.    !    laplacian operator 
     565   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     566   !                       !  Direction of action  : 
     567   ln_dynldf_lev =  .true.     !  iso-level 
     568   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     569   ln_dynldf_iso =  .false.    !  iso-neutral 
    581570   !                       !  Coefficient 
    582    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    583    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    584    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     571   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     572   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     573   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     574   !                                !  =  0  constant  
     575   !                                !  = 10  F(k)=c1d 
     576   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     577   !                                !  = 30  F(i,j,k)=c2d*c1d 
     578   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     579   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     580   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     581   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    585582/ 
    586583 
     
    700697!!   nammpp_dyndist    Massively Parallel domain decomposition          ("key_agrif" && "key_mpp_dyndist") 
    701698!!   namctl            Control prints & Benchmark 
    702 !!   namsol            elliptic solver / island / free surface  
    703 !!====================================================================== 
    704 ! 
    705 !----------------------------------------------------------------------- 
    706 &namsol        !   elliptic solver / island / free surface  
    707 !----------------------------------------------------------------------- 
    708    nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    709                            !                   =2 successive-over-relaxation (sor) 
    710    nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test 
    711    rn_eps      =  1.e-6    !  absolute precision of the solver 
    712    nn_nmin     =    300    !  minimum of iterations for the SOR solver 
    713    nn_nmax     =    800    !  maximum of iterations for the SOR solver 
    714    nn_nmod     =     10    !  frequency of test for the SOR solver 
    715    rn_resmax   =  1.e-10   !  absolute precision for the SOR solver 
    716    rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    717 / 
     699!!====================================================================== 
     700! 
    718701!----------------------------------------------------------------------- 
    719702&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_top_cfg

    r4147 r6043  
    3939&namtrc_adv    !   advection scheme for passive tracer  
    4040!----------------------------------------------------------------------- 
    41    ln_trcadv_cen2   =  .false.  !  2nd order centered scheme    
    42    ln_trcadv_tvd    =  .true.  !  TVD scheme 
    43    ln_trcadv_muscl  =  .false.   !  MUSCL scheme 
    44    ln_trcadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    45    ln_trcadv_ubs    =  .false.  !  UBS scheme 
    46    ln_trcadv_qck    =  .false.  !  QUICKEST scheme 
     41   ln_trcadv_fct =  .true.   !  FCT scheme 
     42      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     43      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     44      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     45      !                             !        (number of sub-timestep = nn_fct_zts) 
    4746/ 
    4847!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg

    r5600 r6043  
    6767/ 
    6868!----------------------------------------------------------------------- 
    69 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    70 !----------------------------------------------------------------------- 
    71 / 
    72 !----------------------------------------------------------------------- 
    7369&namcrs        !   Grid coarsening for dynamics output and/or 
    7470               !   passive tracer coarsened online simulations 
     
    132128/ 
    133129!----------------------------------------------------------------------- 
    134 &namcla        !   cross land advection 
    135 !----------------------------------------------------------------------- 
    136 / 
    137 !----------------------------------------------------------------------- 
    138 &namobc        !   open boundaries parameters                           ("key_obc") 
    139 !----------------------------------------------------------------------- 
    140 / 
    141 !----------------------------------------------------------------------- 
    142130&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    143131!----------------------------------------------------------------------- 
     
    176164/ 
    177165!----------------------------------------------------------------------- 
    178 &namtra_adv    !   advection scheme for tracer  
    179 !-----------------------------------------------------------------------      
     166&namtra_adv    !   advection scheme for tracer 
     167!----------------------------------------------------------------------- 
     168   ln_traadv_fct =  .true.    !  FCT scheme 
     169      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     170      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     171      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     172      !                             !        (number of sub-timestep = nn_fct_zts) 
    180173/ 
    181174!---------------------------------------------------------------------------------- 
    182175&namtra_ldf    !   lateral diffusion scheme for tracers 
    183176!---------------------------------------------------------------------------------- 
     177   !                       !  Operator type: 
     178   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     179   ln_traldf_lap   =  .true.   !    laplacian operator 
     180   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     181   !                       !  Direction of action: 
     182   ln_traldf_lev   =  .false.  !  iso-level 
     183   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     184   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     185   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     186   ! 
     187   !                       !  iso-neutral options:         
     188   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     189   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     190   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     191   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     192   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     193   ! 
     194   !                       !  Coefficients: 
     195   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     196   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     197   !                                !   =  0           constant  
     198   !                                !   = 10 F(k)      =ldf_c1d  
     199   !                                !   = 20 F(i,j)    =ldf_c2d  
     200   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     201   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     202   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     203   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     204   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     205/ 
     206!---------------------------------------------------------------------------------- 
     207&namtra_ldfeiv !   eddy induced velocity param. 
     208!---------------------------------------------------------------------------------- 
     209   ln_ldfeiv     =.true.    ! use eddy induced velocity parameterization 
     210   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     211   rn_aeiv_0     = 2000.    ! eddy induced velocity coefficient   [m2/s] 
     212   nn_aei_ijk_t  = 21       ! space/time variation of the eiv coeficient 
     213   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     214   !                                !   =  0           constant  
     215   !                                !   = 10 F(k)      =ldf_c1d  
     216   !                                !   = 20 F(i,j)    =ldf_c2d  
     217   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     218   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    184219/ 
    185220!----------------------------------------------------------------------- 
     
    188223/ 
    189224!----------------------------------------------------------------------- 
     225&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
     226!----------------------------------------------------------------------- 
     227/ 
     228!----------------------------------------------------------------------- 
    190229&namdyn_adv    !   formulation of the momentum advection 
    191230!----------------------------------------------------------------------- 
    192231 
    193232!----------------------------------------------------------------------- 
    194 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    195 !----------------------------------------------------------------------- 
    196 / 
    197 !----------------------------------------------------------------------- 
    198233&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    199234!----------------------------------------------------------------------- 
     235   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     236   ln_dynvor_ens = .false. !  energy conserving scheme 
     237   ln_dynvor_mix = .false. !  mixed scheme 
     238   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     239      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    200240/ 
    201241!----------------------------------------------------------------------- 
     
    204244/ 
    205245!----------------------------------------------------------------------- 
    206 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    207 !----------------------------------------------------------------------- 
     246&namdyn_spg    !   surface pressure gradient 
     247!----------------------------------------------------------------------- 
     248   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     249/ 
    208250!----------------------------------------------------------------------- 
    209251&namdyn_ldf    !   lateral diffusion on momentum 
    210252!-----------------------------------------------------------------------       
     253   !                       !  Type of the operator : 
     254   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     255   ln_dynldf_lap =  .true.     !    laplacian operator 
     256   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     257   !                       !  Direction of action  : 
     258   ln_dynldf_lev =  .true.     !  iso-level 
     259   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     260   ln_dynldf_iso =  .false.    !  iso-neutral 
     261   !                       !  Coefficient 
     262   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     263   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     264   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     265   !                                !  =  0  constant  
     266   !                                !  = 10  F(k)=c1d 
     267   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     268   !                                !  = 30  F(i,j,k)=c2d*c1d 
     269   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     270   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     271   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     272   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    211273/ 
    212274!----------------------------------------------------------------------- 
     
    222284!----------------------------------------------------------------------- 
    223285/ 
    224 !------------------------------------------------------------------------ 
    225 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    226 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    227 / 
    228286!----------------------------------------------------------------------- 
    229287&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     
    236294!----------------------------------------------------------------------- 
    237295&namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    238 !----------------------------------------------------------------------- 
    239 / 
    240 !----------------------------------------------------------------------- 
    241 &namsol        !   elliptic solver / island / free surface  
    242296!----------------------------------------------------------------------- 
    243297/ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_top_cfg

    r4340 r6043  
    2424&namtrc_adv    !   advection scheme for passive tracer  
    2525!-----------------------------------------------------------------------    
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm

    r4523 r6043  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg

    r4370 r6043  
    4848/ 
    4949!----------------------------------------------------------------------- 
    50 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    51 !----------------------------------------------------------------------- 
    52 / 
    53 !----------------------------------------------------------------------- 
    5450&namcrs        !   Grid coarsening for dynamics output and/or 
    5551               !   passive tracer coarsened online simulations 
     
    9389/ 
    9490!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9991&nambfr        !   bottom friction 
    10092!----------------------------------------------------------------------- 
     
    115107&namtra_adv    !   advection scheme for tracer 
    116108!----------------------------------------------------------------------- 
     109   ln_traadv_fct =  .true.    !  FCT scheme 
     110      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     111      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     112      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     113      !                             !        (number of sub-timestep = nn_fct_zts) 
    117114/ 
    118115!----------------------------------------------------------------------- 
     
    120117!----------------------------------------------------------------------- 
    121118/ 
    122 !----------------------------------------------------------------------- 
    123 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    124 !----------------------------------------------------------------------- 
    125 / 
    126119!---------------------------------------------------------------------------------- 
    127120&namtra_ldf    !   lateral diffusion scheme for tracers 
    128121!---------------------------------------------------------------------------------- 
     122   !                       !  Operator type: 
     123   ln_traldf_lap   =  .true.   !    laplacian operator 
     124   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     125   !                       !  Direction of action: 
     126   ln_traldf_lev   =  .false.  !  iso-level 
     127   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     128   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     129   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     130   ! 
     131   !                       !  iso-neutral options:         
     132   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     133   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     134   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     135   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     136   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     137   ! 
     138   !                       !  Coefficients: 
     139   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     140   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     141   !                                !   =  0           constant 
     142   !                                !   = 10 F(k)      =ldf_c1d 
     143   !                                !   = 20 F(i,j)    =ldf_c2d 
     144   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     145   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     146   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     147   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     148   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     149/ 
     150!---------------------------------------------------------------------------------- 
     151&namtra_ldfeiv !   eddy induced velocity param. 
     152!---------------------------------------------------------------------------------- 
     153   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     154   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     155   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     156   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     157   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     158   !                                !   =  0           constant 
     159   !                                !   = 10 F(k)      =ldf_c1d 
     160   !                                !   = 20 F(i,j)    =ldf_c2d 
     161   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     162   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    129163/ 
    130164!----------------------------------------------------------------------- 
     
    139173&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    140174!----------------------------------------------------------------------- 
     175   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     176   ln_dynvor_ens = .false. !  energy conserving scheme 
     177   ln_dynvor_mix = .false. !  mixed scheme 
     178   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     179      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    141180/ 
    142181!----------------------------------------------------------------------- 
     
    145184/ 
    146185!----------------------------------------------------------------------- 
     186&namdyn_spg    !   surface pressure gradient 
     187!----------------------------------------------------------------------- 
     188   ln_dynspg_ts  = .true.  !  split-explicit free surface 
     189/ 
     190!----------------------------------------------------------------------- 
    147191&namdyn_ldf    !   lateral diffusion on momentum 
    148192!----------------------------------------------------------------------- 
     193   !                       !  Type of the operator : 
     194   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     195   ln_dynldf_lap =  .true.     !    laplacian operator 
     196   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     197   !                       !  Direction of action  : 
     198   ln_dynldf_lev =  .true.     !  iso-level 
     199   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     200   ln_dynldf_iso =  .false.    !  iso-neutral 
     201   !                       !  Coefficient 
     202   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     203   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     204   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     205   !                                !  =  0  constant  
     206   !                                !  = 10  F(k)=c1d 
     207   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     208   !                                !  = 30  F(i,j,k)=c2d*c1d 
     209   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     210   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     211   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     212   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    149213/ 
    150214!----------------------------------------------------------------------- 
     
    165229/ 
    166230!----------------------------------------------------------------------- 
    167 &namsol        !   elliptic solver / island / free surface 
    168 !----------------------------------------------------------------------- 
    169 / 
    170 !----------------------------------------------------------------------- 
    171231&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    172232!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_top_cfg

    r5600 r6043  
    6262&namtrc_adv    !   advection scheme for passive tracer  
    6363!----------------------------------------------------------------------- 
    64    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    65    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     64   ln_trcadv_mus  =  .true.  !  MUSCL scheme 
     65      ln_mus_ups  =  .false.         !  use upstream scheme near river mouths 
    6666/ 
    6767!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

    r4523 r6043  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_mpp_mpi key_iomput 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r5600 r6043  
    5555/ 
    5656!----------------------------------------------------------------------- 
    57 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    58 !----------------------------------------------------------------------- 
    59 / 
    60 !----------------------------------------------------------------------- 
    6157&namcrs        !   Grid coarsening for dynamics output and/or 
    6258               !   passive tracer coarsened online simulations 
     
    7874!----------------------------------------------------------------------- 
    7975/       
    80  
    81 !----------------------------------------------------------------------- 
    82 &namcla        !   cross land advection 
    83 !----------------------------------------------------------------------- 
    84 / 
    8576!----------------------------------------------------------------------- 
    8677&nambbl        !   bottom boundary layer scheme 
     
    9182!----------------------------------------------------------------------- 
    9283/ 
    93 !----------------------------------------------------------------------- 
    94 &namtra_ldf    !   lateral diffusion scheme for tracer  
    95 !----------------------------------------------------------------------- 
    96    ln_triad_iso     =  .true.  !  griffies operator calculates triads twice => pure lateral mixing in ML (require "key_ldfslp") 
     84!---------------------------------------------------------------------------------- 
     85&namtra_ldf    !   lateral diffusion scheme for tracers 
     86!---------------------------------------------------------------------------------- 
     87   !                       !  Operator type: 
     88   ln_traldf_lap   =  .true.   !    laplacian operator 
     89   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     90   !                       !  Direction of action: 
     91   ln_traldf_lev   =  .false.  !  iso-level 
     92   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     93   ln_traldf_iso   =  .true.   !  iso-neutral 
     94   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     95   ! 
     96   !                       !  iso-neutral options:         
     97   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     98   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     99   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     100   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     101   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     102   ! 
     103   !                       !  Coefficients: 
     104   nn_aht_ijk_t    = 0        !  space/time variation of eddy coef 
     105   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     106   !                                !   =  0           constant  
     107   !                                !   = 10 F(k)      =ldf_c1d  
     108   !                                !   = 20 F(i,j)    =ldf_c2d  
     109   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     110   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     111   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     112   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     113   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     114/ 
     115!---------------------------------------------------------------------------------- 
     116&namtra_ldfeiv !   eddy induced velocity param. 
     117!---------------------------------------------------------------------------------- 
     118   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     119   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     120   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     121   nn_aei_ijk_t  = 0      ! space/time variation of the eiv coeficient 
     122   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     123   !                                !   =  0           constant 
     124   !                                !   = 10 F(k)      =ldf_c1d 
     125   !                                !   = 20 F(i,j)    =ldf_c2d 
     126   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     127   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    97128/ 
    98129!----------------------------------------------------------------------- 
     
    115146   sn_ubl  = 'dyna_grid_U' ,    120            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    116147   sn_vbl  = 'dyna_grid_V' ,    120            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    117    sn_ahu  = 'dyna_grid_U' ,    120            , 'vozoahtu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    118    sn_ahv  = 'dyna_grid_V' ,    120            , 'vomeahtv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    119    sn_ahw  = 'dyna_grid_W' ,    120            , 'voveahtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    120    sn_eiu  = 'dyna_grid_U' ,    120            , 'vozoaeiu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    121    sn_eiv  = 'dyna_grid_V' ,    120            , 'vomeaeiv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    122    sn_eiw  = 'dyna_grid_W' ,    120            , 'soleaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    123148! 
    124149   cn_dir      = './'      !  root directory for the location of the dynamical files 
    125    ln_degrad   =  .false.  !  flag for degradation -                requires ("key_degrad") 
    126150   ln_dynwzv   =  .true.   !  computation of vertical velocity instead of using the one read in file 
    127151   ln_dynbbl   =  .true.   !  bbl coef are in files, so read them - requires ("key_trabbl") 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top_cfg

    r5600 r6043  
    6565&namtrc_adv    !   advection scheme for passive tracer  
    6666!----------------------------------------------------------------------- 
    67    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    68    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     67   ln_trcadv_mus =  .true.   !  MUSCL scheme 
     68      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
    6969/ 
    7070!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r4523 r6043  
    1 bld::tool::fppkeys key_trabbl key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_top key_offline key_pisces key_iomput key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist_cfg

    r4370 r6043  
    5050/ 
    5151!----------------------------------------------------------------------- 
    52 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    53 !----------------------------------------------------------------------- 
    54 / 
    55 !----------------------------------------------------------------------- 
    5652&namcrs        !   Grid coarsening for dynamics output and/or 
    5753               !   passive tracer coarsened online simulations 
     
    8783/ 
    8884!----------------------------------------------------------------------- 
    89 &namcla        !   cross land advection 
    90 !----------------------------------------------------------------------- 
    91 / 
    92 !----------------------------------------------------------------------- 
    9385&nameos        !   ocean physical parameters 
    9486!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_SAS_LIM/cpp_ORCA2_SAS_LIM.fcm

    r4523 r6043  
    1  bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/domain_def.xml

    r5600 r6043  
    163163 
    164164      <domain id="grid_V" long_name="grid V"/> 
     165       
    165166      <domain_group id="grid_W"> 
    166167         <domain id="grid_W" long_name="grid W"/> 
     
    168169         <domain id="EqW" zoom_ibegin="1" zoom_jbegin="0000" zoom_ni="0000" zoom_nj="1" /> 
    169170      </domain_group> 
     171 
     172      <domain id="grid_F" long_name="grid F"/> 
    170173 
    171174      <domain_group id="scalarpoint">  
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/field_def.xml

    r5600 r6043  
    2525 
    2626         <field id="toce"         long_name="temperature"         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
    27          <field id="toce_e3t"     long_name="temperature * e3t"                                                     unit="degC*m"   grid_ref="grid_T_3D" > toce * e3t </field > 
     27         <field id="toce_e3t"     long_name="temperature (thickness weighted)"                                      unit="degC"     grid_ref="grid_T_3D" > toce * e3t </field > 
    2828         <field id="soce"         long_name="salinity"            standard_name="sea_water_practical_salinity"      unit="1e-3"     grid_ref="grid_T_3D"/> 
    29          <field id="soce_e3t"     long_name="salinity * e3t"                                                        unit="1e-3*m"   grid_ref="grid_T_3D" > soce * e3t </field > 
     29         <field id="soce_e3t"     long_name="salinity    (thickness weighted)"                                      unit="1e-3"     grid_ref="grid_T_3D" > soce * e3t </field > 
     30 
     31         <!-- t-eddy viscosity coefficients (ldfdyn) --> 
     32         <field id="ahmt_2d"      long_name=" surface t-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     33         <field id="ahmt_3d"      long_name=" 3D      t-eddy viscosity coefficient"   unit="m2/s or m4/s"                           grid_ref="grid_T_3D"/> 
    3034 
    3135         <field id="sst"          long_name="sea surface temperature"             standard_name="sea_surface_temperature"             unit="degC"     /> 
     
    367371         <field id="utau"         long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
    368372         <field id="uoce"         long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
    369          <field id="uoce_e3u"     long_name="ocean current along i-axis * e3u"                                                                   unit="m2/s"       grid_ref="grid_U_3D"  > uoce * e3u </field> 
     373         <field id="uoce_e3u"     long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
    370374         <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             /> 
    371375         <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                             /> 
     
    375379         <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
    376380 
     381         <!-- u-eddy coefficients (ldftra) --> 
     382         <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
     383         <field id="ahtu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s or m4/s"      grid_ref="grid_U_3D"/> 
     384         <field id="aeiu_2d"      long_name=" surface u-EIV coefficient"                unit="m2/s" /> 
     385         <field id="aeiu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s"              grid_ref="grid_U_3D"/> 
     386 
    377387         <!-- variables available with MLE --> 
    378388         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" /> 
    379389 
    380          <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     390         <!-- uoce_eiv: available EIV --> 
    381391         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"   standard_name="bolus_sea_water_x_velocity"   unit="m/s"   grid_ref="grid_U_3D" /> 
    382392 
     
    402412         <field id="vtau"         long_name="Wind Stress along j-axis"                               standard_name="surface_downward_y_stress"   unit="N/m2"                            /> 
    403413         <field id="voce"         long_name="ocean current along j-axis"                             standard_name="sea_water_y_velocity"        unit="m/s"        grid_ref="grid_V_3D" /> 
    404          <field id="voce_e3v"     long_name="ocean current along j-axis * e3v"                                                                   unit="m2/s"       grid_ref="grid_V_3D"  > voce * e3v </field> 
     414         <field id="voce_e3v"     long_name="ocean current along j-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_V_3D"  > voce * e3v </field> 
    405415         <field id="ssv"          long_name="ocean surface current along j-axis"                                                                 unit="m/s"                             /> 
    406416         <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"                             /> 
     
    410420         <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
    411421 
     422         <!-- v-eddy coefficients (ldftra, ldfdyn) --> 
     423         <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
     424         <field id="ahtv_3d"      long_name=" 3D v-eddy diffusivity coefficient"          unit="m2/s or (m4/s)^1/2"           grid_ref="grid_V_3D"/> 
     425         <field id="aeiv_2d"      long_name=" surface v-EIV coefficient"                  unit="m2/s" /> 
     426         <field id="aeiv_3d"      long_name=" 3D v-EIV coefficient"                       unit="m2/s"                         grid_ref="grid_V_3D" /> 
     427 
    412428         <!-- variables available with MLE --> 
    413429         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"   unit="m3/s"   grid_ref="grid_V_3D" /> 
    414430 
    415          <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     431         <!-- voce_eiv: available with EIV --> 
    416432         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"   standard_name="bolus_sea_water_y_velocity"   unit="m/s"   grid_ref="grid_V_3D" /> 
    417433 
     
    438454        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                               unit="m3/s" /> 
    439455 
    440         <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     456        <!-- woce_eiv: available with EIV --> 
    441457        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"   standard_name="bolus_upward_sea_water_velocity"   unit="m/s" /> 
    442458 
     
    449465 
    450466        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    451         <field id="avt_evd"      long_name="convective enhancement to vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    452         <field id="avm_evd"      long_name="convective enhancement to vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
     467        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
     468        <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    453469 
    454470        <!-- avt_tide: available with key_zdftmx --> 
     
    459475        <field id="w_masstr2"    long_name="square of vertical mass transport"   standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" /> 
    460476 
    461         <!-- aht2d and  aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d --> 
     477        <!-- aht2d and  aht2d_eiv --> 
    462478        <field id="aht2d"        long_name="lateral eddy diffusivity"       standard_name="ocean_tracer_xy_laplacian_diffusivity"      unit="m2/s"   grid_ref="grid_W_2D" /> 
    463479        <field id="aht2d_eiv"    long_name="EIV lateral eddy diffusivity"   standard_name="ocean_tracer_bolus_laplacian_diffusivity"   unit="m2/s"   grid_ref="grid_W_2D" /> 
    464480      </field_group> 
    465            
     481         
     482      <!-- F grid --> 
     483         <!-- f-eddy viscosity coefficients (ldfdyn) --> 
     484         <field id="ahmf_2d"      long_name=" surface f-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     485         <field id="ahmf_3d"      long_name=" 3D      f-eddy viscosity coefficient"   unit="m2/s or m4/s"                           grid_ref="grid_T_3D"/> 
     486 
    466487      <!-- scalar variables available with key_diaar5 --> 
    467488 
     
    651672      <field id="ketrd_hpg"     long_name="ke-trend: hydrostatic pressure gradient"          unit="W/s^3"                        /> 
    652673      <field id="ketrd_spg"     long_name="ke-trend: surface     pressure gradient"          unit="W/s^3"                        /> 
    653       <field id="ketrd_spgexp"  long_name="ke-trend: surface pressure gradient (explicit)"   unit="W/s^3"                        /> 
    654       <field id="ketrd_spgflt"  long_name="ke-trend: surface pressure gradient (filter)"     unit="W/s^3"                        /> 
    655       <field id="ssh_flt"       long_name="filtered contribution to ssh (dynspg_flt)"        unit="m"       grid_ref="grid_T_2D" /> 
    656674      <field id="w0"            long_name="surface vertical velocity"                        unit="m/s"     grid_ref="grid_T_2D" /> 
    657675      <field id="pw0_exp"       long_name="surface pressure flux due to ssh"                 unit="W/s^2"   grid_ref="grid_T_2D" /> 
     
    695713     <field id="utrd_hpg"       long_name="i-trend: hydrostatic pressure gradient"          unit="m/s^2"                        /> 
    696714     <field id="utrd_spg"       long_name="i-trend: surface     pressure gradient"          unit="m/s^2"                        /> 
    697      <field id="utrd_spgexp"    long_name="i-trend: surface pressure gradient (explicit)"   unit="m/s^2"                        /> 
    698      <field id="utrd_spgflt"    long_name="i-trend: surface pressure gradient (filtered)"   unit="m/s^2"                        /> 
    699715     <field id="utrd_keg"       long_name="i-trend: KE gradient         or hor. adv."       unit="m/s^2"                        /> 
    700716     <field id="utrd_rvo"       long_name="i-trend: relative  vorticity or metric term"     unit="m/s^2"                        /> 
     
    715731     <field id="vtrd_hpg"       long_name="j-trend: hydrostatic pressure gradient"          unit="m/s^2"                        /> 
    716732     <field id="vtrd_spg"       long_name="j-trend: surface     pressure gradient"          unit="m/s^2"                        /> 
    717      <field id="vtrd_spgexp"    long_name="j-trend: surface pressure gradient (explicit)"   unit="m/s^2"                        /> 
    718      <field id="vtrd_spgflt"    long_name="j-trend: surface pressure gradient (filtered)"   unit="m/s^2"                        /> 
    719733     <field id="vtrd_keg"       long_name="j-trend: KE gradient         or hor. adv."       unit="m/s^2"                        /> 
    720734     <field id="vtrd_rvo"       long_name="j-trend: relative  vorticity or metric term"     unit="m/s^2"                        /> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5600 r6043  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!!                            namelist_ref 
    13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    24!! NEMO/OPA  :  1 - run manager      (namrun) 
     
    57!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    68!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     9!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    810!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    9 !!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     11!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    1012!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    11 !!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
     13!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 
    1214!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    13 !!             10 - miscellaneous    (namsol, nammpp, namctl) 
     15!!             10 - miscellaneous    (nammpp, namctl) 
    1416!!             11 - Obs & Assim      (namobs, nam_asminc) 
    1517!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    4850   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4951   ln_cfmeta   = .false.   !  output additional data to netCDF files required for compliance with the CF metadata standard 
    50    ln_clobber  = .false.   !  clobber (overwrite) an existing file 
     52   ln_clobber  = .true.   !  clobber (overwrite) an existing file 
    5153   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    5254/ 
     
    161163/ 
    162164!----------------------------------------------------------------------- 
    163 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    164 !----------------------------------------------------------------------- 
    165    ln_bt_fw      =    .TRUE.           !  Forward integration of barotropic equations 
    166    ln_bt_av      =    .TRUE.           !  Time filtering of barotropic variables 
    167    ln_bt_nn_auto =    .TRUE.           !  Set nn_baro automatically to be just below 
    168                                        !  a user defined maximum courant number (rn_bt_cmax) 
    169    nn_baro       =    30               !  Number of iterations of barotropic mode 
    170                                        !  during rn_rdt seconds. Only used if ln_bt_nn_auto=F 
    171    rn_bt_cmax    =    0.8              !  Maximum courant number allowed if ln_bt_nn_auto=T 
    172    nn_bt_flt     =    1                !  Time filter choice 
    173                                        !  = 0 None 
    174                                        !  = 1 Boxcar over   nn_baro barotropic steps 
    175                                        !  = 2 Boxcar over 2*nn_baro     "        " 
    176 / 
    177 !----------------------------------------------------------------------- 
    178165&namcrs        !   Grid coarsening for dynamics output and/or 
    179166               !   passive tracer coarsened online simulations 
     
    200187!----------------------------------------------------------------------- 
    201188&namtsd    !   data : Temperature  & Salinity 
    202 !----------------------------------------------------------------------- 
    203189!----------------------------------------------------------------------- 
    204190!          !  file name                            ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    545531!!====================================================================== 
    546532!!   namlbc        lateral momentum boundary condition 
    547 !!   namcla        cross land advection 
    548533!!   namobc        open boundaries parameters                           ("key_obc") 
    549534!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     
    560545/ 
    561546!----------------------------------------------------------------------- 
    562 &namcla        !   cross land advection 
    563 !----------------------------------------------------------------------- 
    564    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    565 / 
    566 !----------------------------------------------------------------------- 
    567 &namobc        !   open boundaries parameters                           ("key_obc") 
    568 !----------------------------------------------------------------------- 
    569    ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    570    ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    571    ln_obc_fla  = .false.   !  Flather open boundary condition 
    572    nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    573                            !  = 1 the obc data are read in 'obc.dta' files 
    574    cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    575                            !  set to monthly if obc datafile hold 1 month of data 
    576    rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
    577    rn_dpwin    =    1.     !     -           -         -       west    -      - 
    578    rn_dpnin    =    1.     !     -           -         -       north   -      - 
    579    rn_dpsin    =    1.     !     -           -         -       south   -      - 
    580    rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
    581    rn_dpwob    =   15.     !     -           -         -     west    -      - 
    582    rn_dpnob    = 3000.     !     -           -         -     north   -      - 
    583    rn_dpsob    =   15.     !     -           -         -     south   -      - 
    584    rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    585                            !  = 1 the total volume remains constant 
    586 / 
    587 !----------------------------------------------------------------------- 
    588547&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    589548!----------------------------------------------------------------------- 
     
    592551   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    593552   rn_sponge_dyn = 2880.   !  coefficient for dynamics sponge layer [m2/s] 
     553   ln_chk_bathy  = .FALSE. ! 
    594554/ 
    595555!----------------------------------------------------------------------- 
     
    710670                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    711671   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    712  
    713672/ 
    714673!----------------------------------------------------------------------- 
     
    724683!!                        Tracer (T & S ) namelists 
    725684!!====================================================================== 
    726 !!   nameos        equation of state 
    727 !!   namtra_adv    advection scheme 
     685!!   nameos           equation of state 
     686!!   namtra_adv       advection scheme 
    728687!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    729 !!   namtra_ldf    lateral diffusion scheme 
    730 !!   namtra_dmp    T & S newtonian damping 
     688!!   namtra_ldf       lateral diffusion scheme 
     689!!   namtra_ldfeiv    eddy induced velocity param. 
     690!!   namtra_dmp       T & S newtonian damping 
    731691!!====================================================================== 
    732692! 
     
    739699                                 !  = 1, S-EOS   (simplified eos) 
    740700   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
    741    !                             ! 
     701                                 ! 
    742702   !                     ! S-EOS coefficients : 
    743    !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     703                                 !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    744704   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    745705   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     
    753713&namtra_adv    !   advection scheme for tracer 
    754714!----------------------------------------------------------------------- 
    755    ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
    756    ln_traadv_tvd    =  .true.    !  TVD scheme 
    757    ln_traadv_muscl  =  .false.   !  MUSCL scheme 
    758    ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
    759    ln_traadv_ubs    =  .false.   !  UBS scheme 
    760    ln_traadv_qck    =  .false.   !  QUICKEST scheme 
    761    ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
    762    ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection 
     715   ln_traadv_cen =  .false.  !  2nd order centered scheme 
     716      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     717      nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     718   ln_traadv_fct =  .false.  !  FCT scheme 
     719      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     720      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     721      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     722      !                             !        (number of sub-timestep = nn_fct_zts) 
     723   ln_traadv_mus =  .false.  !  MUSCL scheme 
     724      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     725   ln_traadv_ubs =  .false.  !  UBS scheme 
     726      nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT 
     727   ln_traadv_qck =  .false.  !  QUICKEST scheme 
    763728/ 
    764729!----------------------------------------------------------------------- 
    765730&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
    766731!----------------------------------------------------------------------- 
    767    ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     732   ln_mle    = .false.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    768733   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
    769734   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     
    779744!---------------------------------------------------------------------------------- 
    780745   !                       !  Operator type: 
    781    ln_traldf_lap    =  .true.   !  laplacian operator 
    782    ln_traldf_bilap  =  .false.  !  bilaplacian operator 
     746   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     747   ln_traldf_lap   =  .false.  !    laplacian operator 
     748   ln_traldf_blp   =  .false.  !  bilaplacian operator 
    783749   !                       !  Direction of action: 
    784    ln_traldf_level  =  .false.  !  iso-level 
    785    ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    786    ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp") 
    787    !                 !  Griffies parameters              (all need "key_ldfslp") 
    788    ln_traldf_grif   =  .false.  !  use griffies triads 
    789    ln_traldf_gdia   =  .false.  !  output griffies eddy velocities 
    790    ln_triad_iso     =  .false.  !  pure lateral mixing in ML 
    791    ln_botmix_grif   =  .false.  !  lateral mixing on bottom 
    792    !                       !  Coefficients 
    793    ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 
    794    ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 
    795    !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 
    796    rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s] 
    797    rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    798    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    799    !                                           (normally=0; not used with Griffies) 
    800    rn_slpmax        =     0.01  !  slope limit 
    801    rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
    802    rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
    803    rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
     750   ln_traldf_lev   =  .false.  !  iso-level 
     751   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     752   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator) 
     753   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     754   ! 
     755   !                       !  iso-neutral options:         
     756   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
     757   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     758   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     759   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     760   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     761   ! 
     762   !                       !  Coefficients: 
     763   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     764   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     765   !                                !   =  0           constant  
     766   !                                !   = 10 F(k)      =ldf_c1d  
     767   !                                !   = 20 F(i,j)    =ldf_c2d  
     768   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     769   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     770   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     771   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     772   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     773/ 
     774!---------------------------------------------------------------------------------- 
     775&namtra_ldfeiv !   eddy induced velocity param. 
     776!---------------------------------------------------------------------------------- 
     777   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
     778   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     779   rn_aeiv_0     = 2000.    ! eddy induced velocity coefficient   [m2/s] 
     780   nn_aei_ijk_t  = 21       ! space/time variation of the eiv coeficient 
     781   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     782   !                                !   =  0           constant  
     783   !                                !   = 10 F(k)      =ldf_c1d  
     784   !                                !   = 20 F(i,j)    =ldf_c2d  
     785   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     786   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    804787/ 
    805788!----------------------------------------------------------------------- 
     
    819802!!   namdyn_vor    advection scheme 
    820803!!   namdyn_hpg    hydrostatic pressure gradient 
    821 !!   namdyn_spg    surface pressure gradient                            (CPP key only) 
     804!!   namdyn_spg    surface pressure gradient 
    822805!!   namdyn_ldf    lateral diffusion scheme 
    823806!!====================================================================== 
     
    852835   ln_dynvor_ens = .false. !  energy conserving scheme 
    853836   ln_dynvor_mix = .false. !  mixed scheme 
    854    ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    855    ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation 
     837   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     838      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     839   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT USE 
    856840/ 
    857841!----------------------------------------------------------------------- 
     
    864848   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    865849   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    866    ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    867                                  !           centered      time scheme  (F) 
    868 / 
    869 !----------------------------------------------------------------------- 
    870 !namdyn_spg    !   surface pressure gradient   (CPP key only) 
    871 !----------------------------------------------------------------------- 
    872 !                          !  explicit free surface                     ("key_dynspg_exp") 
    873 !                          !  filtered free surface                     ("key_dynspg_flt") 
    874 !                          !  split-explicit free surface               ("key_dynspg_ts") 
    875  
     850/ 
     851!----------------------------------------------------------------------- 
     852&namdyn_spg    !   surface pressure gradient 
     853!----------------------------------------------------------------------- 
     854   ln_dynspg_exp  = .false.   ! explicit free surface 
     855   ln_dynspg_ts   = .false.   ! split-explicit free surface 
     856      ln_bt_fw      = .true.    ! Forward integration of barotropic Eqs. 
     857      ln_bt_av      = .true.    ! Time filtering of barotropic variables 
     858         nn_bt_flt    =  1        ! Time filter choice  = 0 None 
     859         !                        !                     = 1 Boxcar over   nn_baro sub-steps 
     860         !                        !                     = 2 Boxcar over 2*nn_baro  "    " 
     861      ln_bt_auto    = .true.    ! Number of sub-step defined from: 
     862         rn_bt_cmax   =  0.8      ! =T : the Maximum Courant Number allowed 
     863         nn_baro      = 30        ! =F : the number of sub-step in rn_rdt seconds 
     864/ 
    876865!----------------------------------------------------------------------- 
    877866&namdyn_ldf    !   lateral diffusion on momentum 
    878867!----------------------------------------------------------------------- 
    879868   !                       !  Type of the operator : 
    880    ln_dynldf_lap    =  .true.   !  laplacian operator 
    881    ln_dynldf_bilap  =  .false.  !  bilaplacian operator 
     869   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     870   ln_dynldf_lap =  .false.    !    laplacian operator 
     871   ln_dynldf_blp =  .false.    !  bilaplacian operator 
    882872   !                       !  Direction of action  : 
    883    ln_dynldf_level  =  .false.  !  iso-level 
    884    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    885    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     873   ln_dynldf_lev =  .false.    !  iso-level 
     874   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     875   ln_dynldf_iso =  .false.    !  iso-neutral 
    886876   !                       !  Coefficient 
    887    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    888    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    889    rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    890    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    891    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    892    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    893    rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    894    rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
     877   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     878   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     879   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     880   !                                !  =  0  constant  
     881   !                                !  = 10  F(k)=c1d 
     882   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     883   !                                !  = 30  F(i,j,k)=c2d*c1d 
     884   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     885   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     886   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     887   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     888   ! 
     889   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    895890/ 
    896891 
     
    901896!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric") 
    902897!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke") 
    903 !!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp") 
    904898!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm") 
    905899!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx") 
     
    963957                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    964958/ 
    965 !------------------------------------------------------------------------ 
    966 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    967 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    968    ln_kpprimix = .true.    !  shear instability mixing 
    969    rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s] 
    970    rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
    971    rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability 
    972    rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s] 
    973    rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2] 
    974    rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s] 
    975    nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
    976    nn_ave      =  1        !  constant (=0) or profile (=1) background on avt 
    977 / 
    978959!----------------------------------------------------------------------- 
    979960&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     
    1014995!!                  ***  Miscellaneous namelists  *** 
    1015996!!====================================================================== 
    1016 !!   namsol            elliptic solver / island / free surface 
    1017997!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    1018998!!   namctl            Control prints & Benchmark 
     
    10231003!!====================================================================== 
    10241004! 
    1025 !----------------------------------------------------------------------- 
    1026 &namsol        !   elliptic solver / island / free surface 
    1027 !----------------------------------------------------------------------- 
    1028    nn_solv     =      1    !  elliptic solver: =1 preconditioned conjugate gradient (pcg) 
    1029                            !                   =2 successive-over-relaxation (sor) 
    1030    nn_sol_arp  =      0    !  absolute/relative (0/1) precision convergence test 
    1031    rn_eps      =  1.e-6    !  absolute precision of the solver 
    1032    nn_nmin     =    300    !  minimum of iterations for the SOR solver 
    1033    nn_nmax     =    800    !  maximum of iterations for the SOR solver 
    1034    nn_nmod     =     10    !  frequency of test for the SOR solver 
    1035    rn_resmax   =  1.e-10   !  absolute precision for the SOR solver 
    1036    rn_sor      =  1.92     !  optimal coefficient for SOR solver (to be adjusted with the domain) 
    1037 / 
    10381005!----------------------------------------------------------------------- 
    10391006&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_top_ref

    r5600 r6043  
    4141&namtrc_adv    !   advection scheme for passive tracer  
    4242!----------------------------------------------------------------------- 
    43    ln_trcadv_cen2    =  .false.  !  2nd order centered scheme    
    44    ln_trcadv_tvd     =  .true.  !  TVD scheme 
    45    ln_trcadv_muscl   =  .false.   !  MUSCL scheme 
    46    ln_trcadv_muscl2  =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    47    ln_trcadv_ubs     =  .false.  !  UBS scheme 
    48    ln_trcadv_qck     =  .false.  !  QUICKEST scheme 
    49    ln_trcadv_msc_ups =  .false.  !  use upstream scheme within muscl 
     43   ln_trcadv_cen =  .false.  !  2nd order centered scheme 
     44      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     45      nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     46   ln_trcadv_fct =  .false.  !  FCT scheme 
     47      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     48      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     49      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     50      !                             !        (number of sub-timestep = nn_fct_zts) 
     51   ln_trcadv_mus =  .false.  !  MUSCL scheme 
     52      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     53   ln_trcadv_ubs =  .false.  !  UBS scheme 
     54      nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT 
     55   ln_trcadv_qck =  .false.  !  QUICKEST scheme 
    5056/ 
    5157!----------------------------------------------------------------------- 
    5258&namtrc_ldf    !   lateral diffusion scheme for passive tracer  
    5359!----------------------------------------------------------------------- 
    54 !                               !  Type of the operator :  
    55    ln_trcldf_lap    =  .true.   !     laplacian operator        
    56    ln_trcldf_bilap  =  .false.  !     bilaplacian operator      
    57                                 !  Direction of action  : 
    58    ln_trcldf_level  =  .false.  !     iso-level                 
    59    ln_trcldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    60    ln_trcldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
    61 !                               !  Coefficient 
    62    rn_ahtrc_0       =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    63    rn_ahtrb_0       =     0.    !     background eddy diffusivity for ldf_iso [m2/s] 
     60!                          !  Type of the operator:   
     61   ln_trcldf_lap   =  .true.   !    laplacian operator 
     62   ln_trcldf_blp   =  .false.  !  bilaplacian operator 
     63   !                       !  Direction of action: 
     64   ln_trcldf_lev   =  .false.  !  iso-level 
     65   ln_trcldf_hor   =  .false.  !  horizontal (geopotential) 
     66   ln_trcldf_iso   =  .true.   !  iso-neutral (standard operator) 
     67   ln_trcldf_triad =  .false.  !  iso-neutral (triad    operator) 
     68   !                       !  Coefficient  
     69   rn_ahtrc_0      = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     70   rn_bhtrc_0      = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    6471/ 
    6572!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.