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 4596 for branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/CONFIG – NEMO

Ignore:
Timestamp:
2014-03-26T12:02:30+01:00 (10 years ago)
Author:
gm
Message:

#1260: LDF simplification + bilap iso-neutral for TRA and GYRE

Location:
branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/CONFIG
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r4370 r4596  
    8282!----------------------------------------------------------------------- 
    8383   nn_fsbc     = 1         !  frequency of surface boundary condition computation 
    84                            !     (also = the frequency of sea-ice model call) 
     84   !                       !     (also = the frequency of sea-ice model call) 
    8585   ln_ana      = .true.    !  analytical formulation                    (T => fill namsbc_ana ) 
    8686   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
     
    212212&namtra_ldf    !   lateral diffusion scheme for tracers 
    213213!---------------------------------------------------------------------------------- 
    214    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    215    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     214   !                       !  Operator type: 
     215   ln_traldf_lap   =  .true.   !  laplacian operator 
     216   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     217   !                       !  Direction of action: 
     218   ln_traldf_lev   =  .false.  !  iso-level 
     219   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     220   ln_traldf_iso   =  .false.  !  iso-neutral 
     221   ln_traldf_triad =  .true.   !  iso-neutral with Griffies triads 
     222   ! 
     223   !                       !  iso-neutral options          
     224   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operator) 
     225   rn_slpmax       =   0.01    !  slope limit                      (both operator) 
     226   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     227   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     228   ! 
     229   !                       !  Coefficients 
     230   nn_aht_ijk_t    =  0        !  space/time variation of eddy coef 
     231   !                                   !   =0 constant ; =10 F(k) ; =20 F(i,j)=F(grid spacing) ; =30 F(i,j,k)  
     232   !                                   !   =21 F(i,jt)=Treguier et al. JPO 1997 formulation 
     233   !                                   !   =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     234   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     235   rn_bht_0        = 5.e+11    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     236/ 
     237!---------------------------------------------------------------------------------- 
     238&namtra_ldfeiv !   eddy induced velocity param. 
     239!---------------------------------------------------------------------------------- 
     240   ln_ldfeiv     =.false.  ! use eddy induced velocity parameterization 
     241   ln_ldfeiv_dia =.false.  ! diagnose eiv stream function and velocities 
     242   rn_aeiv_0     = 1000.   ! eddy induced velocity coefficient   [m2/s] 
     243   nn_aei_ijk_t  =  0      ! space/time variation of the eiv coeficient 
     244   !                              !   =0 constant ; =10 F(k) ; =20 F(i,j) = F(grid spacing) ; =30 F(i,j,k)  
     245   !                              !   =21 F(i,jt)=Treguier et al. JPO 1997 formulation 
     246   !                              !   =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    216247/ 
    217248!----------------------------------------------------------------------- 
     
    242273!namdyn_spg    !   surface pressure gradient   (CPP key only) 
    243274!----------------------------------------------------------------------- 
     275 
    244276!----------------------------------------------------------------------- 
    245277&namdyn_ldf    !   lateral diffusion on momentum 
    246278!----------------------------------------------------------------------- 
    247    rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
     279   !                       !  Type of the operator : 
     280   ln_dynldf_lap =  .true.     !  laplacian operator 
     281   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     282   !                       !  Direction of action  : 
     283   ln_dynldf_lev =  .false.    !  iso-level 
     284   ln_dynldf_hor =  .true.     !  horizontal (geopotential) 
     285   ln_dynldf_iso =  .false.    !  iso-neutral 
     286   !                       !  Coefficient 
     287   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     288   !                                !   =0 constant ; =10 F(k) ; =20 F(i,j)=F(grid spacing) ; =30 F(i,j,k)  
     289   !                                !   =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     290   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     291   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     292   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    248293/ 
    249294!----------------------------------------------------------------------- 
  • branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

    r4230 r4596  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_iomput key_mpp_mpi 
  • branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/CONFIG/SHARED/field_def.xml

    r4341 r4596  
    144144         <field id="precip"       long_name="Total precipitation"                                           unit="kg/m2/s"  /> 
    145145 
    146  
    147146         <field id="qt"           long_name="Net Downward Heat Flux"                                       unit="W/m2"     /> 
    148147         <field id="qns"          long_name="non solar Downward Heat Flux"                                 unit="W/m2"     /> 
     
    274273         <field id="uocet"        long_name="ocean transport along i-axis times temperature" unit="degC.m/s" grid_ref="grid_U_3D" /> 
    275274         <field id="uoces"        long_name="ocean transport along i-axis times salinity"    unit="psu.m/s"  grid_ref="grid_U_3D" /> 
     275         <!-- u-eddy coefficients (ldftra, ldfdyn) --> 
     276         <field id="aeiu_2d"      long_name=" surface u-EIV coefficient"                unit="m2/s" /> 
     277         <field id="aeiu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s"  grid_ref="grid_U_3D"/> 
     278         <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
     279         <field id="ahtu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s or m4/s"  grid_ref="grid_U_3D"/> 
    276280         <!-- variables available with MLE --> 
    277281         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"             unit="m3/s"  grid_ref="grid_U_3D"  /> 
     
    279283         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
    280284         <!-- uoce_eiv: available with key_trabbl --> 
    281          <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
    282     <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"             unit="m3/s" /> 
     285         <field id="uoce_bbl"     long_name="BBL  ocean current along i-axis"             unit="m/s"  grid_ref="grid_U_3D" /> 
     286         <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"             unit="m3/s" /> 
    283287         <!-- variables available with key_diaar5 --> 
    284288         <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"  unit="kg/s" grid_ref="grid_U_3D" /> 
     
    297301         <field id="vocet"        long_name="ocean transport along j-axis times temperature" unit="degC.m/s" grid_ref="grid_V_3D" /> 
    298302         <field id="voces"        long_name="ocean transport along j-axis times salinity"    unit="psu.m/s"  grid_ref="grid_V_3D" /> 
     303         <!-- v-eddy coefficients (ldftra, ldfdyn) --> 
     304         <field id="aeiv_2d"      long_name=" surface v-EIV coefficient"                  unit="m2/s" /> 
     305         <field id="aeiv_3d"      long_name=" 3D v-EIV coefficient"                       unit="m2/s" grid_ref="grid_V_3D" /> 
     306         <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
     307         <field id="ahtv_3d"      long_name=" 3D v-eddy diffusivity coefficient"          unit="m2/s or (m4/s)^1/2" grid_ref="grid_V_3D"/> 
    299308         <!-- variables available with MLE --> 
    300309         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"             unit="m3/s"  grid_ref="grid_V_3D"  /> 
  • branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/CONFIG/SHARED/namelist_ref

    r4384 r4596  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  :  1 - run manager      (namrun, namcfg) 
    3 !! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
    4 !!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    5 !!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    6 !!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
     2!! NEMO/OPA  :  1 - run manager      (namrun) 
     3!! namelists    2 - Domain           (namcfg, namzgr, nam_vvl, namzgr_sco, namdom, namtsd) 
     4!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, 
     5!!                                            namsbc_sas, namsbc_cpl, namtra_qsr , namsbc_rnf , 
     6!!                                            namsbc_apr, namsbc_ssr, namsbc_alb) 
    77!!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    9 !!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     9!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    1010!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    1111!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
     
    3131   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3232   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    33    nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=.true. 
    34    nn_rstctl   =       0   !  restart control => activated only if ln_rstart = T 
     33   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T 
     34   nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T 
    3535                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    3636                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     
    4646   nn_chunksz  =       0   !  chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 
    4747/ 
     48!!====================================================================== 
     49!!                      ***  Domain namelists  *** 
     50!!====================================================================== 
     51!!   namcfg       parameters of the configuration 
     52!!   namzgr       vertical coordinate 
     53!!   nam_vvl      vertical coordinate options (z-star, z-tilde) 
     54!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     55!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     56!!   namtsd       data: temperature & salinity 
     57!!====================================================================== 
    4858! 
    4959!----------------------------------------------------------------------- 
    50 &namcfg     !   default parameters of the configuration       
     60&namcfg     !   parameters of the configuration       
    5161!----------------------------------------------------------------------- 
    5262   cp_cfg      =  "default"            !  name of the configuration 
    53    cp_cfz      =         ''            !  name of the zoom of configuration 
     63   cp_cfz      =  "no zoom"            !  name of the zoom of configuration 
    5464   jp_cfg      =       0               !  resolution of the configuration 
    5565   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
     
    5767   jpkdta      =      31               !  number of levels      ( >= jpk ) 
    5868   jpiglo      =      10               !  1st dimension of global domain --> i =jpidta 
    59    jpjglo      =      12               !  2nd    -                  -    --> j  =jpjdta 
     69   jpjglo      =      12               !  2nd    -                  -    --> j =jpjdta 
    6070   jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
    6171   jpjzoom     =       1               !  in data domain indices 
     
    6777                                       !  = 6 cyclic East-West AND North fold F-point pivot 
    6878/ 
    69 !!====================================================================== 
    70 !!                      ***  Domain namelists  *** 
    71 !!====================================================================== 
    72 !!   namzgr       vertical coordinate 
    73 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    74 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    75 !!   namtsd       data: temperature & salinity 
    76 !!====================================================================== 
    77 ! 
    7879!----------------------------------------------------------------------- 
    7980&namzgr        !   vertical coordinate 
     
    8283   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    8384   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     85/ 
     86!----------------------------------------------------------------------- 
     87&nam_vvl       !   vertical coordinate options 
     88!----------------------------------------------------------------------- 
     89   ln_vvl_zstar  = .true.           !  zstar vertical coordinate                    
     90   ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
     91   ln_vvl_layer  = .false.          !  full layer vertical coordinate 
     92   ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar 
     93   ln_vvl_zstar_at_eqtor  = .false. !  ztilde near the equator 
     94   rn_ahe3       = 0.0e0            !  thickness diffusion coefficient 
     95   rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days] 
     96   rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days] 
     97   rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation 
     98   ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
    8499/ 
    85100!----------------------------------------------------------------------- 
     
    205220!!   namsbc_mfs      MFS  bulk formulae formulation 
    206221!!   namsbc_cpl      CouPLed            formulation                     ("key_coupled") 
    207 !!   namsbc_sas      StAndalone Surface module 
     222!!   namsbc_sas      Stand Alone Surface module 
    208223!!   namtra_qsr      penetrative solar radiation 
    209224!!   namsbc_rnf      river runoffs 
     
    633648!!                        Tracer (T & S ) namelists 
    634649!!====================================================================== 
    635 !!   nameos        equation of state 
    636 !!   namtra_adv    advection scheme 
    637 !!   namtra_ldf    lateral diffusion scheme 
    638 !!   namtra_dmp    T & S newtonian damping 
     650!!   nameos           equation of state 
     651!!   namtra_adv       advection scheme 
     652!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
     653!!   namtra_ldf       lateral diffusion scheme 
     654!!   namtra_ldfeiv    eddy induced velocity param. 
     655!!   namtra_dmp       T & S newtonian damping 
    639656!!====================================================================== 
    640657! 
     
    661678/ 
    662679!----------------------------------------------------------------------- 
    663 &namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
     680&namtra_adv_mle !   mixed layer eddy param. (Fox-Kemper param.) 
    664681!----------------------------------------------------------------------- 
    665682   ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     
    677694!---------------------------------------------------------------------------------- 
    678695   !                       !  Operator type: 
    679    ln_traldf_lap    =  .true.   !  laplacian operator 
    680    ln_traldf_bilap  =  .false.  !  bilaplacian operator 
     696   ln_traldf_lap   =  .true.   !    laplacian operator 
     697   ln_traldf_blp   =  .false.  !  bilaplacian operator 
    681698   !                       !  Direction of action: 
    682    ln_traldf_level  =  .false.  !  iso-level 
    683    ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    684    ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp") 
    685    !                 !  Griffies parameters              (all need "key_ldfslp") 
    686    ln_traldf_grif   =  .false.  !  use griffies triads 
    687    ln_traldf_gdia   =  .false.  !  output griffies eddy velocities 
    688    ln_triad_iso     =  .false.  !  pure lateral mixing in ML 
    689    ln_botmix_grif   =  .false.  !  lateral mixing on bottom 
    690    !                       !  Coefficients 
    691    ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 
    692    ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 
    693    !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 
    694    rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s] 
    695    rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    696    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    697    !                                           (normally=0; not used with Griffies) 
    698    rn_slpmax        =     0.01  !  slope limit 
    699    rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
    700    rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
    701    rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
     699   ln_traldf_lev   =  .false.  !  iso-level 
     700   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     701   ln_traldf_iso   =  .true.   !  iso-neutral 
     702   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     703   ! 
     704   !                       !  iso-neutral options:         
     705   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     706   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     707   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     708   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     709   ! 
     710   !                       !  Coefficients: 
     711   nn_aht_ijk_t    = 21        !  space/time variation of eddy coef 
     712   !                                !   =0 constant ; =10 F(k) ; =20 F(i,j)=F(grid spacing) ; =30 F(i,j,k)  
     713   !                                !   =21 F(i,jt)=Treguier et al. JPO 1997 formulation 
     714   !                                !   =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     715   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     716   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     717 
     718!!gm   rn_chsmag        =    1.    !  multiplicative factor in Smagorinsky diffusivity 
     719!!gm   rn_smsh          =    1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
     720/ 
     721!---------------------------------------------------------------------------------- 
     722&namtra_ldfeiv !   eddy induced velocity param. 
     723!---------------------------------------------------------------------------------- 
     724   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     725   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     726   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     727   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     728   !                           !   =0 constant ; =10 F(k) ; =20 F(i,j) = F(grid spacing) ; =30 F(i,j,k)  
     729   !                           !   =21 F(i,jt)=Treguier et al. JPO 1997 formulation 
     730   !                           !   =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    702731/ 
    703732!----------------------------------------------------------------------- 
     
    735764/ 
    736765!----------------------------------------------------------------------- 
    737 &nam_vvl    !   vertical coordinate options 
    738 !----------------------------------------------------------------------- 
    739    ln_vvl_zstar  = .true.           !  zstar vertical coordinate                    
    740    ln_vvl_ztilde = .false.          !  ztilde vertical coordinate: only high frequency variations 
    741    ln_vvl_layer  = .false.          !  full layer vertical coordinate 
    742    ln_vvl_ztilde_as_zstar = .false. !  ztilde vertical coordinate emulating zstar 
    743    ln_vvl_zstar_at_eqtor = .false.  !  ztilde near the equator 
    744    rn_ahe3       = 0.0e0            !  thickness diffusion coefficient 
    745    rn_rst_e3t    = 30.e0            !  ztilde to zstar restoration timescale [days] 
    746    rn_lf_cutoff  = 5.0e0            !  cutoff frequency for low-pass filter  [days] 
    747    rn_zdef_max   = 0.9e0            !  maximum fractional e3t deformation 
    748    ln_vvl_dbg    = .true.           !  debug prints    (T/F) 
    749 / 
    750 !----------------------------------------------------------------------- 
    751766&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    752767!----------------------------------------------------------------------- 
     
    778793!----------------------------------------------------------------------- 
    779794   !                       !  Type of the operator : 
    780    ln_dynldf_lap    =  .true.   !  laplacian operator 
    781    ln_dynldf_bilap  =  .false.  !  bilaplacian operator 
     795   ln_dynldf_lap =  .true.     !  laplacian operator 
     796   ln_dynldf_blp =  .false.    !  bilaplacian operator 
    782797   !                       !  Direction of action  : 
    783    ln_dynldf_level  =  .false.  !  iso-level 
    784    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    785    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     798   ln_dynldf_lev =  .false.    !  iso-level 
     799   ln_dynldf_hor =  .true.     !  horizontal (geopotential) 
     800   ln_dynldf_iso =  .false.    !  iso-neutral 
    786801   !                       !  Coefficient 
    787    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    788    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    789    rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    790    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    791    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    792    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    793    rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    794    rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
     802   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     803   !                                !   =0 constant ; =10 F(k) ; =20 F(i,j)=F(grid spacing) ; =30 F(i,j,k)  
     804   !                                !   =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     805   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     806   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     807   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     808 
     809 
     810!!gm   rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
     811!!gm   rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
     812!!gm   rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    795813/ 
    796814 
     
    914932!!                  ***  Miscellaneous namelists  *** 
    915933!!====================================================================== 
     934!!   namsol            elliptic solver / island / free surface 
    916935!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    917936!!   namctl            Control prints & Benchmark 
    918 !!   namsol            elliptic solver / island / free surface 
     937!!   namc1d            1D configuration options                         ("key_c1d") 
     938!!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
     939!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
    919940!!====================================================================== 
    920941! 
     
    9811002   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    9821003/ 
     1004 
    9831005!!====================================================================== 
    9841006!!                  ***  Diagnostics namelists  *** 
     
    9891011!!   namptr       Poleward Transport Diagnostics 
    9901012!!   namhsb       Heat and salt budgets 
     1013!!   nam_diaharm  Harmonic analysis of tidal constituents               ('key_diaharm') 
     1014!!   namdct       transports through sections 
    9911015!!====================================================================== 
    9921016! 
     
    10681092!!   namobs       observation and model comparison                      ('key_diaobs') 
    10691093!!   nam_asminc   assimilation increments                               ('key_asminc') 
     1094!!   namsbc_wave  External fields from wave model 
     1095!!   namdyn_nept  Neptune effect  
    10701096!!====================================================================== 
    10711097! 
Note: See TracChangeset for help on using the changeset viewer.