Changeset 7954


Ignore:
Timestamp:
2017-04-24T09:19:00+02:00 (3 years ago)
Author:
gm
Message:

#1880 (HPC-09): OPA & TOP replace key_trabbl by ln_trabbl + in CONFIG, remove reference to key_zdfXXX and key_trabbl

Location:
branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM
Files:
58 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r7953 r7954  
    189189/ 
    190190!----------------------------------------------------------------------- 
    191 &nambbl        !   bottom boundary layer scheme 
    192 !----------------------------------------------------------------------- 
    193    nn_bbl_ldf  =  0      !  diffusive bbl (=1)   or not (=0) 
     191&nambbl        !   bottom boundary layer scheme                         (default: NO) 
     192!----------------------------------------------------------------------- 
    194193/ 
    195194!----------------------------------------------------------------------- 
     
    305304/ 
    306305!----------------------------------------------------------------------- 
    307 &namzdf        !   vertical physics 
    308 !----------------------------------------------------------------------- 
    309    rn_avm0     =   0.1e-6  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    310    rn_avt0     =   0.1e-6  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    311    ln_zdfevd   = .false.   !  enhanced vertical diffusion (evd) (T) or not (F) 
    312    nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    313 / 
    314 !----------------------------------------------------------------------- 
    315 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    316 !----------------------------------------------------------------------- 
    317 / 
    318 !----------------------------------------------------------------------- 
    319 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    320 !----------------------------------------------------------------------- 
    321 / 
    322 !----------------------------------------------------------------------- 
    323 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
     306&namzdf        !   vertical physics                                     (default: NO selection) 
     307!----------------------------------------------------------------------- 
     308   !                       ! type of vertical closure 
     309   ln_zdfcst   = .false.      !  constant mixing 
     310   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     311   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     312   ln_zdfgls   = .true.       !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     313   ! 
     314   !                       ! convection 
     315   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     316      nn_evdm     =    0         ! apply on tracer (=0) or on tracer and momentum (=1) 
     317      rn_evd      =  100.        ! mixing coefficient [m2/s] 
     318   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     319      nn_npc      =    1         ! frequency of application of npc 
     320      nn_npcp     =  365         ! npc control print frequency 
     321   ! 
     322   ln_zdfddm   = .false.   ! double diffusive mixing 
     323      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     324      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     325   ! 
     326   ln_zdftmx   = .false.   ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
     327   ! 
     328   ln_zdfqiao  = .false.   ! surface wave-induced mixing                (T => ln_wave=ln_sdw=T ) 
     329   ! 
     330   !                       ! time-stepping 
     331   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     332      nn_zdfexp=    3            !  number of sub-timestep for ln_zdfexp=T 
     333   ! 
     334   !                       ! coefficients 
     335   rn_avm0     =   0.1e-6     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     336   rn_avt0     =   0.1e-6     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     337   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     338   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
     339/ 
     340!----------------------------------------------------------------------- 
     341&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T) 
     342!----------------------------------------------------------------------- 
     343/ 
     344!----------------------------------------------------------------------- 
     345&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     346!----------------------------------------------------------------------- 
     347/ 
     348!----------------------------------------------------------------------- 
     349&namzdf_gls                !   GLS vertical diffusion                   (ln_zdfgls =T) 
    324350!----------------------------------------------------------------------- 
    325351   rn_charn =  100000.     !  Charnock constant for wb induced roughness length 
     
    327353/ 
    328354!----------------------------------------------------------------------- 
    329 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    330 !----------------------------------------------------------------------- 
    331 / 
    332 !----------------------------------------------------------------------- 
    333 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    334 !----------------------------------------------------------------------- 
    335    ln_tmx_itf  = .FALSE.   !  ITF specific parameterisation 
     355&namzdf_tmx    !   internal wave-driven mixing parameterization         (ln_zdftmx =T) 
     356!----------------------------------------------------------------------- 
    336357/ 
    337358!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm

    r7646 r7954  
    1  bld::tool::fppkeys key_zdfgls key_diainstant key_mpp_mpi key_iomput 
     1 bld::tool::fppkeys   key_diainstant key_mpp_mpi key_iomput 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r7646 r7954  
    255255/ 
    256256!----------------------------------------------------------------------- 
    257 &namzdf        !   vertical physics 
    258 !----------------------------------------------------------------------- 
    259    ln_zdfevd   = .false.        !  enhanced vertical diffusion (evd) (T) or not (F) 
    260 / 
    261 !----------------------------------------------------------------------- 
    262 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    263 !----------------------------------------------------------------------- 
    264 / 
    265 !----------------------------------------------------------------------- 
    266 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    267 !----------------------------------------------------------------------- 
    268 / 
    269 !----------------------------------------------------------------------- 
    270 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    271 !----------------------------------------------------------------------- 
    272 / 
    273 !----------------------------------------------------------------------- 
    274 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    275 !----------------------------------------------------------------------- 
    276 / 
    277 !----------------------------------------------------------------------- 
    278 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    279 !----------------------------------------------------------------------- 
    280    ln_tmx_itf  = .false.   !  ITF specific parameterisation 
     257&namzdf        !   vertical physics                                     (default: NO selection) 
     258!----------------------------------------------------------------------- 
     259   !                       ! type of vertical closure 
     260   ln_zdfcst   = .false.      !  constant mixing 
     261   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     262   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     263   ln_zdfgls   = .true.       !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     264   ! 
     265   !                       ! convection 
     266   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     267      nn_evdm     =    0         ! apply on tracer (=0) or on tracer and momentum (=1) 
     268      rn_evd      =  100.        ! mixing coefficient [m2/s] 
     269   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     270      nn_npc      =    1         ! frequency of application of npc 
     271      nn_npcp     =  365         ! npc control print frequency 
     272   ! 
     273   ln_zdfddm   = .false.   ! double diffusive mixing 
     274      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     275      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     276   ! 
     277   ln_zdftmx   = .false.   ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
     278   ! 
     279   ln_zdfqiao  = .false.   ! surface wave-induced mixing                (T => ln_wave=ln_sdw=T ) 
     280   ! 
     281   !                       ! time-stepping 
     282   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     283   ! 
     284   !                       ! coefficients 
     285   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     286   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     287   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     288   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
     289/ 
     290!----------------------------------------------------------------------- 
     291&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T) 
     292!----------------------------------------------------------------------- 
     293/ 
     294!----------------------------------------------------------------------- 
     295&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     296!----------------------------------------------------------------------- 
     297/ 
     298!----------------------------------------------------------------------- 
     299&namzdf_gls                !   GLS vertical diffusion                   (ln_zdfgls =T) 
     300!----------------------------------------------------------------------- 
     301/ 
     302!----------------------------------------------------------------------- 
     303&namzdf_tmx    !    internal wave-driven mixing parameterization        (ln_zdftmx =T) 
     304!----------------------------------------------------------------------- 
    281305/ 
    282306!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/C1D_PAPA/cpp_C1D_PAPA.fcm

    r4667 r7954  
    1  bld::tool::fppkeys key_c1d key_zdfgls 
     1 bld::tool::fppkeys   key_c1d 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r7715 r7954  
    4343/ 
    4444!----------------------------------------------------------------------- 
    45 &namcrs        !   Grid coarsening for dynamics output and/or 
    46                !   passive tracer coarsened online simulations 
     45&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    4746!----------------------------------------------------------------------- 
    4847/ 
     
    245244/ 
    246245!----------------------------------------------------------------------- 
    247 &namzdf        !   vertical physics 
    248 !----------------------------------------------------------------------- 
    249    nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    250 / 
    251 !----------------------------------------------------------------------- 
    252 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    253 !----------------------------------------------------------------------- 
    254 / 
    255 !----------------------------------------------------------------------- 
    256 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
     246&namzdf        !   vertical physics                                     (default: NO selection) 
     247!----------------------------------------------------------------------- 
     248   !                       ! type of vertical closure 
     249   ln_zdfcst   = .false.      !  constant mixing 
     250   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     251   ln_zdftke   = .true.       !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     252   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     253   ! 
     254   !                       ! convection 
     255   ln_zdfevd   = .true.       !  enhanced vertical diffusion 
     256      nn_evdm     =    1         ! apply on tracer (=0) or on tracer and momentum (=1) 
     257      rn_evd      =  100.        ! mixing coefficient [m2/s] 
     258   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     259      nn_npc      =    1         ! frequency of application of npc 
     260      nn_npcp     =  365         ! npc control print frequency 
     261   ! 
     262   ln_zdfddm   = .false.   ! double diffusive mixing 
     263      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     264      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     265   ! 
     266   ln_zdftmx   = .false.   ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
     267   ! 
     268   ln_zdfqiao  = .false.   ! surface wave-induced mixing                (T => ln_wave=ln_sdw=T ) 
     269   ! 
     270   !                       ! time-stepping 
     271   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     272      nn_zdfexp=    3            !  number of sub-timestep for ln_zdfexp=T 
     273   ! 
     274   !                       ! coefficients 
     275   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     276   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     277   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     278   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
     279/ 
     280!----------------------------------------------------------------------- 
     281&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T) 
     282!----------------------------------------------------------------------- 
     283/ 
     284!----------------------------------------------------------------------- 
     285&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke") 
    257286!----------------------------------------------------------------------- 
    258287   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    259288/ 
    260289!----------------------------------------------------------------------- 
    261 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    262 !----------------------------------------------------------------------- 
    263 / 
    264 !----------------------------------------------------------------------- 
    265 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    266 !----------------------------------------------------------------------- 
    267 / 
    268 !----------------------------------------------------------------------- 
    269 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    270 !----------------------------------------------------------------------- 
    271    ln_tmx_itf  = .false.   !  ITF specific parameterisation 
     290&namzdf_gls                !   GLS vertical diffusion                   (ln_zdfgls =T) 
     291!----------------------------------------------------------------------- 
     292/ 
     293!----------------------------------------------------------------------- 
     294&namzdf_tmx    !    internal wave-driven mixing parameterization        (ln_zdftmx =T) 
     295!----------------------------------------------------------------------- 
    272296/ 
    273297!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

    r5930 r7954  
    1 bld::tool::fppkeys key_zdftke key_top key_my_trc key_mpp_mpi key_iomput 
     1bld::tool::fppkeys   key_top key_my_trc key_mpp_mpi key_iomput 
    22inc $BFMDIR/src/nemo/bfm.fcm 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r7953 r7954  
    8383/ 
    8484!----------------------------------------------------------------------- 
    85 &nambbl        !   bottom boundary layer scheme 
     85&nambbl        !   bottom boundary layer scheme                         (default: NO) 
    8686!----------------------------------------------------------------------- 
    8787/ 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r7953 r7954  
    1 bld::tool::fppkeys key_zdftke key_top key_mpp_mpi key_nosignedzero 
     1bld::tool::fppkeys   key_top key_mpp_mpi key_nosignedzero 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/1_namelist_cfg

    r7715 r7954  
    2020/ 
    2121!----------------------------------------------------------------------- 
    22 &namzgr        !   vertical coordinate 
    23 !----------------------------------------------------------------------- 
    24    ln_zps      = .true.    !  z-coordinate - partial steps 
    25 / 
    26 !----------------------------------------------------------------------- 
    2722&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2823!-----------------------------------------------------------------------   
     
    3429/ 
    3530!----------------------------------------------------------------------- 
    36 &namcrs        !   Grid coarsening for dynamics output and/or 
    37 !              !   passive tracer coarsened online simulations 
     31&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    3832!----------------------------------------------------------------------- 
    3933/ 
     
    10397&nambbl        !   bottom boundary layer scheme 
    10498!----------------------------------------------------------------------- 
     99   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
     100   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
     101   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     102   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
     103   rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     104/ 
    105105/ 
    106106!----------------------------------------------------------------------- 
     
    205205   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    206206/ 
    207 !----------------------------------------------------------------------- 
    208 &namzdf        !   vertical physics 
    209 !----------------------------------------------------------------------- 
    210 / 
    211 !----------------------------------------------------------------------- 
    212 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    213 !----------------------------------------------------------------------- 
    214 / 
    215 !----------------------------------------------------------------------- 
    216 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    217 !----------------------------------------------------------------------- 
    218 / 
    219 !----------------------------------------------------------------------- 
    220 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    221 !----------------------------------------------------------------------- 
    222    ln_tmx_itf  = .false.   !  ITF specific parameterisation 
     207!!====================================================================== 
     208!!                     vertical physics namelists                     !! 
     209!!====================================================================== 
     210!----------------------------------------------------------------------- 
     211&namzdf        !   vertical physics                                     (default: NO selection) 
     212!----------------------------------------------------------------------- 
     213   !                       ! type of vertical closure 
     214   ln_zdfcst   = .false.      !  constant mixing 
     215   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     216   ln_zdftke   = .true.       !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     217   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     218   ! 
     219   !                       ! convection 
     220   ln_zdfevd   = .true.       !  Enhanced Vertical Diffusion scheme 
     221      nn_evdm  =    0            !  evd apply on tracer (=0) or on tracer and momentum (=1) 
     222      rn_evd   =  100.           !  evd mixing coefficient [m2/s] 
     223   ! 
     224   ln_zdfddm   = .true.    ! double diffusive mixing 
     225      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     226      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     227   ! 
     228   ln_zdftmx   = .true.    ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
     229   ! 
     230   ln_zdfqiao  = .false.   ! enhanced wave vertical mixing Qiao (2010) (T => ln_wave=T & ln_sdw=T & fill namsbc_wave) 
     231   ! 
     232   !                       ! time-stepping 
     233   ln_zdfexp   = .false.   ! split-explicit (T) or implicit (F) time stepping scheme 
     234      nn_zdfexp=    3         !  number of sub-timestep for ln_zdfexp=T 
     235   ! 
     236   !                       !  Coefficients 
     237   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     238   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     239   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     240   nn_havtb    =    1         !  horizontal shape for avtb (=1) or not (=0) 
     241/ 
     242!----------------------------------------------------------------------- 
     243&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     244!----------------------------------------------------------------------- 
     245/ 
     246!----------------------------------------------------------------------- 
     247&namzdf_tmx    !    internal wave-driven mixing parameterization        (ln_zdftmx =T) 
     248!----------------------------------------------------------------------- 
     249   nn_zpyc     = 1         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 
     250   ln_mevar    = .true.    !  variable (T) or constant (F) mixing efficiency 
     251   ln_tsdiff   = .true.    !  account for differential T/S mixing (T) or not (F) 
    223252/ 
    224253!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/file_def_nemo.xml

    r7828 r7954  
    9595     <field field_ref="avt"          name="difvho"   /> 
    9696     <field field_ref="w_masstr"     name="vovematr" /> 
    97           <!-- variables available with key_zdftmx_new --> 
     97          <!-- variables available when ln_zdftmx =T --> 
    9898          <field field_ref="av_wave"      name="av_wave"    /> 
    9999          <field field_ref="bn2"          name="bn2"        /> 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/EXP00/namelist_cfg

    r7953 r7954  
    1717      !                    !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules 
    1818      cn_domcfg = "ORCA_R2_zps_domcfg"    ! domain configuration filename 
    19 / 
    20 !----------------------------------------------------------------------- 
    21 &namzgr        !   vertical coordinate 
    22 !----------------------------------------------------------------------- 
    23    ln_zps      = .true.    !  z-coordinate - partial steps 
    2419/ 
    2520!----------------------------------------------------------------------- 
     
    114109/ 
    115110!----------------------------------------------------------------------- 
    116 &nambbl        !   bottom boundary layer scheme 
    117 !----------------------------------------------------------------------- 
     111&nambbl        !   bottom boundary layer scheme                         (default: NO) 
     112!----------------------------------------------------------------------- 
     113   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
     114   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
     115   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     116   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
     117   rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
    118118/ 
    119119!----------------------------------------------------------------------- 
     
    234234   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    235235/ 
     236!!====================================================================== 
     237!!                     vertical physics namelists                     !! 
     238!!====================================================================== 
    236239!----------------------------------------------------------------------- 
    237240&namzdf        !   vertical physics                                     (default: NO selection) 
     
    271274/ 
    272275!----------------------------------------------------------------------- 
    273 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    274 !----------------------------------------------------------------------- 
    275 / 
    276 !----------------------------------------------------------------------- 
    277 &namzdf_tmx_new !   internal wave-driven mixing parameterization        ("key_zdftmx_new" & "key_zdfddm") 
     276&namzdf_tmx    !   tidal mixing parameterization                        (ln_zdftmx =T) 
    278277!----------------------------------------------------------------------- 
    279278   nn_zpyc     = 2         !  pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_LIM3_PISCES/cpp_ORCA2_LIM3_PISCES.fcm

    r7953 r7954  
    1 bld::tool::fppkeys key_trabbl key_lim3 key_zdftke key_zdftmx_new key_iomput key_mpp_mpi key_top key_nosignedzero 
     1bld::tool::fppkeys   key_lim3 key_top key_iomput key_mpp_mpi key_nosignedzero 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r7646 r7954  
    1818/ 
    1919!----------------------------------------------------------------------- 
    20 &namzgr        !   vertical coordinate 
    21 !----------------------------------------------------------------------- 
    22    ln_zps      = .true.    !  z-coordinate - partial steps 
    23 / 
    24 !----------------------------------------------------------------------- 
    2520&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2621!----------------------------------------------------------------------- 
     
    3530   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    3631   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs. 
    37 / 
    38 !----------------------------------------------------------------------- 
    39 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    40 !----------------------------------------------------------------------- 
    41 / 
    42 !----------------------------------------------------------------------- 
    43 &namcrs        !   Grid coarsening for dynamics output and/or 
    44                !   passive tracer coarsened online simulations 
    45 !----------------------------------------------------------------------- 
    4632/ 
    4733!----------------------------------------------------------------------- 
     
    6955&nambbl        !   bottom boundary layer scheme 
    7056!----------------------------------------------------------------------- 
     57   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
     58   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
     59   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     60   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
     61   rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
     62/ 
    7163/ 
    7264!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r7646 r7954  
    1 bld::tool::fppkeys key_trabbl key_top key_iomput key_mpp_mpi 
     1bld::tool::fppkeys   key_top key_iomput key_mpp_mpi 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_OFF_TRC/EXP00/namelist_cfg

    r7445 r7954  
    2020/ 
    2121!----------------------------------------------------------------------- 
    22 &namzgr        !   vertical coordinate 
    23 !----------------------------------------------------------------------- 
    24    ln_zps      = .true.    !  z-coordinate - partial steps 
    25 / 
    26 !----------------------------------------------------------------------- 
    2722&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2823!----------------------------------------------------------------------- 
     
    3732   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    3833   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical Eqs. 
    39 / 
    40 !----------------------------------------------------------------------- 
    41 &namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
    42 !----------------------------------------------------------------------- 
    43 / 
    44 !----------------------------------------------------------------------- 
    45 &namcrs        !   Grid coarsening for dynamics output and/or 
    46                !   passive tracer coarsened online simulations 
    47 !----------------------------------------------------------------------- 
    4834/ 
    4935!----------------------------------------------------------------------- 
     
    7157&nambbl        !   bottom boundary layer scheme 
    7258!----------------------------------------------------------------------- 
     59   ln_trabbl   = .true.    !  Bottom Boundary Layer parameterisation flag 
     60   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
     61   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     62   rn_ahtbbl   =  1000.    !  lateral mixing coefficient in the bbl  [m2/s] 
     63   rn_gambbl   =  10.      !  advective bbl coefficient                 [s] 
    7364/ 
    7465!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_OFF_TRC/cpp_ORCA2_OFF_TRC.fcm

    r7485 r7954  
    1 bld::tool::fppkeys key_trabbl key_top key_iomput key_mpp_mpi 
     1bld::tool::fppkeys   key_top key_iomput key_mpp_mpi 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_SAS_LIM3/EXP00/namelist_cfg

    r7404 r7954  
    1919/ 
    2020!----------------------------------------------------------------------- 
    21 &namzgr        !   vertical coordinate 
    22 !----------------------------------------------------------------------- 
    23    ln_zps      = .true.    !  z-coordinate - partial steps 
    24 / 
    25 !----------------------------------------------------------------------- 
    2621&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2722!----------------------------------------------------------------------- 
     
    3227/ 
    3328!----------------------------------------------------------------------- 
    34 &namcrs        !   Grid coarsening for dynamics output and/or 
    35                !   passive tracer coarsened online simulations 
     29&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    3630!----------------------------------------------------------------------- 
    3731/ 
     
    7569/ 
    7670!----------------------------------------------------------------------- 
    77 &nambfr        !   bottom friction 
    78 !----------------------------------------------------------------------- 
    79 / 
    80 !----------------------------------------------------------------------- 
    81 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    82 !----------------------------------------------------------------------- 
    83    ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    84 / 
    85 !----------------------------------------------------------------------- 
    86 &nambbl        !   bottom boundary layer scheme 
    87 !----------------------------------------------------------------------- 
    88 / 
    89 !----------------------------------------------------------------------- 
    9071&nameos        !   ocean physical parameters 
    9172!----------------------------------------------------------------------- 
     
    10283/ 
    10384!----------------------------------------------------------------------- 
    104 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    105 !----------------------------------------------------------------------- 
    106 / 
    107 !---------------------------------------------------------------------------------- 
    108 &namtra_ldf    !   lateral diffusion scheme for tracers 
    109 !---------------------------------------------------------------------------------- 
    110    !                       !  Operator type: 
    111    ln_traldf_lap   =  .true.   !    laplacian operator 
    112    ln_traldf_blp   =  .false.  !  bilaplacian operator 
    113    !                       !  Direction of action: 
    114    ln_traldf_lev   =  .false.  !  iso-level 
    115    ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    116    ln_traldf_iso   =  .true.   !  iso-neutral (Standard operator) 
    117    ln_traldf_triad =  .false.  !  iso-neutral (Triads   operator) 
    118    ! 
    119    !                       !  iso-neutral options:         
    120    ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
    121    rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    122    ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
    123    rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
    124    ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
    125    ! 
    126    !                       !  Coefficients: 
    127    nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
    128    !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    129    !                                !   =  0           constant  
    130    !                                !   = 10 F(k)      =ldf_c1d  
    131    !                                !   = 20 F(i,j)    =ldf_c2d  
    132    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    133    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    134    !                                !   = 31 F(i,j,k,t)=F(local velocity) 
    135    rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
    136    rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    137 / 
    138 !---------------------------------------------------------------------------------- 
    139 &namtra_ldfeiv !   eddy induced velocity param. 
    140 !---------------------------------------------------------------------------------- 
    141    ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
    142    ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
    143    rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    144    nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
    145    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    146    !                                !   =  0           constant  
    147    !                                !   = 10 F(k)      =ldf_c1d  
    148    !                                !   = 20 F(i,j)    =ldf_c2d  
    149    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    150    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    151 / 
    152 !----------------------------------------------------------------------- 
    153 &namtra_dmp    !   tracer: T & S newtonian damping                      (default: NO) 
    154 !----------------------------------------------------------------------- 
    155 !----------------------------------------------------------------------- 
    156 &namdyn_adv    !   formulation of the momentum advection 
    157 !----------------------------------------------------------------------- 
    158 / 
    159 !----------------------------------------------------------------------- 
    160 &namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    161 !----------------------------------------------------------------------- 
    162    ln_dynvor_ene = .false. !  enstrophy conserving scheme 
    163    ln_dynvor_ens = .false. !  energy conserving scheme 
    164    ln_dynvor_mix = .false. !  mixed scheme 
    165    ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    166       nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    167 / 
    168 !----------------------------------------------------------------------- 
    169 &namdyn_hpg    !   Hydrostatic pressure gradient option 
    170 !----------------------------------------------------------------------- 
    171 / 
    172 !----------------------------------------------------------------------- 
    173 &namdyn_spg    !   surface pressure gradient 
    174 !----------------------------------------------------------------------- 
    175    ln_dynspg_ts  = .true.  !  split-explicit free surface 
    176 / 
    177 !----------------------------------------------------------------------- 
    178 &namdyn_ldf    !   lateral diffusion on momentum 
    179 !----------------------------------------------------------------------- 
    180    !                       !  Type of the operator : 
    181    !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
    182    ln_dynldf_lap =  .true.     !    laplacian operator 
    183    ln_dynldf_blp =  .false.    !  bilaplacian operator 
    184    !                       !  Direction of action  : 
    185    ln_dynldf_lev =  .true.     !  iso-level 
    186    ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
    187    ln_dynldf_iso =  .false.    !  iso-neutral 
    188    !                       !  Coefficient 
    189    nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
    190    !                                !  =-30  read in eddy_viscosity_3D.nc file 
    191    !                                !  =-20  read in eddy_viscosity_2D.nc file 
    192    !                                !  =  0  constant  
    193    !                                !  = 10  F(k)=c1d 
    194    !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
    195    !                                !  = 30  F(i,j,k)=c2d*c1d 
    196    !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
    197    rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
    198    rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
    199    rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    200    ! 
    201    ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    202 / 
    203 !----------------------------------------------------------------------- 
    204 &namzdf        !   vertical physics 
    205 !----------------------------------------------------------------------- 
    206 / 
    207 !----------------------------------------------------------------------- 
    208 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    209 !----------------------------------------------------------------------- 
    210 / 
    211 !----------------------------------------------------------------------- 
    212 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    213 !----------------------------------------------------------------------- 
    214 / 
    215 !----------------------------------------------------------------------- 
    216 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    217 !----------------------------------------------------------------------- 
    218 / 
    219 !----------------------------------------------------------------------- 
    22085&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    22186!----------------------------------------------------------------------- 
     
    22590!----------------------------------------------------------------------- 
    22691/ 
    227 !----------------------------------------------------------------------- 
    228 &namptr       !   Poleward Transport Diagnostic 
    229 !----------------------------------------------------------------------- 
    230 / 
    231 !----------------------------------------------------------------------- 
    232 &namhsb       !  Heat and salt budgets                                  (default F) 
    233 !----------------------------------------------------------------------- 
    234 / 
    235 !----------------------------------------------------------------------- 
    236 &namobs       !  observation usage                                      ('key_diaobs') 
    237 !----------------------------------------------------------------------- 
    238 / 
    239 !----------------------------------------------------------------------- 
    240 &nam_asminc   !   assimilation increments                               ('key_asminc') 
    241 !----------------------------------------------------------------------- 
    242 / 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/ORCA2_SAS_LIM3/cpp_ORCA2_SAS_LIM3.fcm

    r7423 r7954  
    1  bld::tool::fppkeys key_trabbl key_lim3  key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys   key_lim3 key_iomput key_mpp_mpi 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/SHARED/field_def_nemo-opa.xml

    r7953 r7954  
    307307         <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" /> 
    308308 
    309          <!-- uoce_eiv: available with key_trabbl --> 
     309         <!-- variables available when ln_trabbl = T --> 
    310310         <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"    unit="m/s"  /> 
    311311         <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"   unit="m3/s" /> 
     
    355355         <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" /> 
    356356 
    357          <!-- voce_eiv: available with key_trabbl --> 
     357         <!-- variables available when ln_trabbl = T --> 
    358358         <field id="voce_bbl"     long_name="BBL ocean current along j-axis"    unit="m/s"  /> 
    359359         <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"   unit="m3/s" /> 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/SHARED/namelist_ref

    r7953 r7954  
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!! NEMO/OPA  :  1 - run manager      (namrun) 
    5 !! namelists    2 - Domain           (namcfg, namzgr, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
     5!! namelists    2 - Domain           (namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd) 
    66!!              3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_sas) 
    77!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
     
    5959!!====================================================================== 
    6060!!   namcfg       parameters of the configuration 
    61 !!   namzgr       vertical coordinate                                   (default: NO selection) 
    6261!!   namdom       space and time domain (bathymetry, mesh, timestep) 
    6362!!   namwad       Wetting and drying                                    (default F) 
    6463!!   namtsd       data: temperature & salinity 
    65 !!   namcrs       coarsened grid (for outputs and/or TOP)               ("key_crs") 
     64!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T) 
    6665!!   namc1d       1D configuration options                              ("key_c1d") 
    6766!!   namc1d_dyndmp 1D newtonian damping applied on currents             ("key_c1d") 
     
    9493   rn_atfp     =    0.1    !  asselin time filter parameter 
    9594   ! 
    96    ln_crs      = .false.   !  Logical switch for coarsening module 
     95   ln_crs      = .false.   !  Logical switch for coarsening module      (T => fill namcrs) 
    9796/ 
    9897!----------------------------------------------------------------------- 
     
    118117/ 
    119118!----------------------------------------------------------------------- 
    120 &namcrs        !   coarsened grid (for outputs and/or TOP)              ("key_crs") 
     119&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    121120!----------------------------------------------------------------------- 
    122121   nn_factx    = 3         !  Reduction factor of x-direction 
     
    606605!!   nambfr        bottom friction 
    607606!!   nambbc        bottom temperature boundary condition 
    608 !!   nambbl        bottom boundary layer scheme                         ("key_trabbl") 
     607!!   nambbl        bottom boundary layer scheme                         (default: NO) 
    609608!!====================================================================== 
    610609! 
     
    647646/ 
    648647!----------------------------------------------------------------------- 
    649 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    650 !----------------------------------------------------------------------- 
     648&nambbl        !   bottom boundary layer scheme                         (default: NO) 
     649!----------------------------------------------------------------------- 
     650   ln_trabbl   = .false.   !  Bottom Boundary Layer parameterisation flag 
    651651   nn_bbl_ldf  =  1        !  diffusive bbl (=1)   or not (=0) 
    652652   nn_bbl_adv  =  0        !  advective bbl (=1/2) or not (=0) 
     
    901901   ln_zdftmx   = .false.   ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
    902902   ! 
    903    ln_zdfqiao  = .false.   ! surface wave-induced mixing (Qiao et al. 2010) (T =>   ln_wave=ln_sdw=T. & fill namsbc_wave) 
     903   ln_zdfqiao  = .false.   ! surface wave-induced mixing                (T => ln_wave=ln_sdw=T ) 
    904904   ! 
    905905   !                       ! time-stepping 
    906    ln_zdfexp   = .false.   ! split-explicit (T) or implicit (F) scheme 
    907       nn_zdfexp=    3         !  number of sub-timestep for ln_zdfexp=T 
     906   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     907      nn_zdfexp=    3            !  number of sub-timestep for ln_zdfexp=T 
    908908   ! 
    909909   !                       ! coefficients 
     
    912912   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
    913913   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    914    ! 
    915 / 
    916 !----------------------------------------------------------------------- 
    917 &namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric=T) 
     914/ 
     915!----------------------------------------------------------------------- 
     916&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T) 
    918917!----------------------------------------------------------------------- 
    919918   rn_avmri    =  100.e-4  !  maximum value of the vertical viscosity 
     
    928927/ 
    929928!----------------------------------------------------------------------- 
    930 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke=T) 
     929&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
    931930!----------------------------------------------------------------------- 
    932931   rn_ediff    =   0.1     !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     
    956955/ 
    957956!----------------------------------------------------------------------- 
    958 &namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls) 
     957&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T) 
    959958!----------------------------------------------------------------------- 
    960959   rn_emin       = 1.e-7   !  minimum value of e   [m2/s2] 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/EXP00/namelist_cfg

    r7715 r7954  
    326326/ 
    327327!----------------------------------------------------------------------- 
    328 &namzdf        !   vertical physics 
    329 !----------------------------------------------------------------------- 
    330    rn_avm0     =   1.0e-3  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    331    rn_avt0     =   5.0e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    332    nn_avb      =    0      !  profile for background avt & avm (=1) or not (=0) 
    333    nn_havtb    =    0      !  horizontal shape for avtb (=1) or not (=0) 
    334    ln_zdfevd   = .true.    !  enhanced vertical diffusion (evd) (T) or not (F) 
    335    nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    336    rn_avevd    =   0.1     !  evd mixing coefficient [m2/s] 
    337    ln_zdfnpc   = .false.   !  Non-Penetrative Convective algorithm (T) or not (F) 
    338    nn_npc      =    1            !  frequency of application of npc 
    339    nn_npcp     =  365            !  npc control print frequency 
    340    ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping 
    341    nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T 
    342 / 
    343 !----------------------------------------------------------------------- 
    344 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    345 !----------------------------------------------------------------------- 
    346 / 
    347 !----------------------------------------------------------------------- 
    348 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    349 !----------------------------------------------------------------------- 
    350 / 
    351 !----------------------------------------------------------------------- 
    352 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    353 !----------------------------------------------------------------------- 
    354 / 
    355 !----------------------------------------------------------------------- 
    356 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    357 !----------------------------------------------------------------------- 
    358 / 
    359 !----------------------------------------------------------------------- 
    360 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    361 !----------------------------------------------------------------------- 
    362    ln_tmx_itf  = .false.   !  ITF specific parameterisation 
     328&namzdf        !   vertical physics                                     (default: NO selection) 
     329!----------------------------------------------------------------------- 
     330   !                       ! type of vertical closure 
     331   ln_zdfcst   = .true.       !  constant mixing 
     332   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     333   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     334   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     335   ! 
     336   !                       ! convection 
     337   ln_zdfevd   = .true.       !  enhanced vertical diffusion 
     338      nn_evdm     =    1         ! apply on tracer (=0) or on tracer and momentum (=1) 
     339      rn_evd      =   0.1        ! mixing coefficient [m2/s] 
     340   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     341      nn_npc      =    1         ! frequency of application of npc 
     342      nn_npcp     =  365         ! npc control print frequency 
     343   ! 
     344   ln_zdfddm   = .false.   ! double diffusive mixing 
     345   ! 
     346   ln_zdftmx   = .false.   ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
     347   ! 
     348   ln_zdfqiao  = .false.   ! surface wave-induced mixing                (T => ln_wave=ln_sdw=T ) 
     349   ! 
     350   !                       ! time-stepping 
     351   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     352   ! 
     353   !                       ! coefficients 
     354   rn_avm0     =   1.e-3     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     355   rn_avt0     =   5.e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     356   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     357   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    363358/ 
    364359!----------------------------------------------------------------------- 
     
    380375/ 
    381376!----------------------------------------------------------------------- 
    382 &namflo       !   float parameters                                      ("key_float") 
    383 !----------------------------------------------------------------------- 
    384 / 
    385 !----------------------------------------------------------------------- 
    386 &namptr       !   Poleward Transport Diagnostic 
    387 !----------------------------------------------------------------------- 
    388 / 
    389 !----------------------------------------------------------------------- 
    390377&namhsb       !  Heat and salt budgets 
    391378!----------------------------------------------------------------------- 
    392379/ 
    393380!----------------------------------------------------------------------- 
    394 &namdct        ! transports through sections 
    395 !----------------------------------------------------------------------- 
    396 / 
    397 !----------------------------------------------------------------------- 
    398 &namobs       !  observation usage switch                               ('key_diaobs') 
     381&namobs       !  observation usage switch                               (ln_diaobs =T) 
    399382!----------------------------------------------------------------------- 
    400383/ 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/cpp_ISOMIP.fcm

    r7715 r7954  
    1  bld::tool::fppkeys key_zdfcst key_iomput key_mpp_mpi key_nosignedzero  
     1 bld::tool::fppkeys   key_iomput key_mpp_mpi key_nosignedzero  
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_flux_cen2_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_flux_ubs_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_eenH_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_een_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_ene_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT2_vect_ens_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_flux_cen2_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_flux_ubs_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_eenH_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_een_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_ene_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_FCT4_vect_ens_cfg

    r7640 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/EXP00/namelist_cfg

    r7623 r7954  
    7272   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    7373                           !                              = 2 : nonlinear friction 
    74 / 
    75 !----------------------------------------------------------------------- 
    76 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    77 !----------------------------------------------------------------------- 
    78 / 
    79 !----------------------------------------------------------------------- 
    80 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    81 !----------------------------------------------------------------------- 
    8274/ 
    8375!----------------------------------------------------------------------- 
     
    197189   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    198190/ 
    199 !----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    203    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    204    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    205    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     191!!====================================================================== 
     192!!                     vertical physics namelists                     !! 
     193!!====================================================================== 
     194!----------------------------------------------------------------------- 
     195&namzdf        !   vertical physics                                     (default: NO selection) 
     196!----------------------------------------------------------------------- 
     197   !                       ! type of vertical closure 
     198   ln_zdfcst   = .true.       !  constant mixing 
     199   ! 
     200   !                       ! convection 
     201   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     202   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     203   ! 
     204   !                       ! time-stepping 
     205   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     206   ! 
     207   !                       ! coefficients 
     208   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     209   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     210   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     211   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    206212/ 
    207213!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/LOCK_EXCHANGE/cpp_LOCK_EXCHANGE.fcm

    r7423 r7954  
    1  bld::tool::fppkeys key_zdfcst key_iomput key_mpp_mpi key_nosignedzero 
     1 bld::tool::fppkeys   key_iomput key_mpp_mpi key_nosignedzero 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/OVERFLOW/EXP00/namelist_cfg

    r7624 r7954  
    6666   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    6767                           !                              = 2 : nonlinear friction 
    68 / 
    69 !----------------------------------------------------------------------- 
    70 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    71 !----------------------------------------------------------------------- 
    72 / 
    73 !----------------------------------------------------------------------- 
    74 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    75 !----------------------------------------------------------------------- 
    7668/ 
    7769!----------------------------------------------------------------------- 
     
    191183   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    192184/ 
    193 !----------------------------------------------------------------------- 
    194 &namzdf        !   vertical physics 
    195 !----------------------------------------------------------------------- 
    196    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    197    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    198    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    199    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     185!!====================================================================== 
     186!!                     vertical physics namelists                     !! 
     187!!====================================================================== 
     188!----------------------------------------------------------------------- 
     189&namzdf        !   vertical physics                                     (default: NO selection) 
     190!----------------------------------------------------------------------- 
     191   !                       ! type of vertical closure 
     192   ln_zdfcst   = .true.       !  constant mixing 
     193   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     194   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     195   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     196   ! 
     197   !                       ! convection 
     198   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     199   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     200   ! 
     201   !                       ! time-stepping 
     202   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     203   ! 
     204   !                       ! coefficients 
     205   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     206   rn_avt0     =   0.         !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     207   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     208   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    200209/ 
    201210!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/OVERFLOW/EXP00/namelist_sco_FCT2_flux_ubs_cfg

    r7640 r7954  
    6666   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    6767                           !                              = 2 : nonlinear friction 
    68 / 
    69 !----------------------------------------------------------------------- 
    70 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    71 !----------------------------------------------------------------------- 
    72 / 
    73 !----------------------------------------------------------------------- 
    74 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    75 !----------------------------------------------------------------------- 
    7668/ 
    7769!----------------------------------------------------------------------- 
     
    191183   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    192184/ 
    193 !----------------------------------------------------------------------- 
    194 &namzdf        !   vertical physics 
    195 !----------------------------------------------------------------------- 
    196    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    197    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    198    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    199    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     185!!====================================================================== 
     186!!                     vertical physics namelists                     !! 
     187!!====================================================================== 
     188!----------------------------------------------------------------------- 
     189&namzdf        !   vertical physics                                     (default: NO selection) 
     190!----------------------------------------------------------------------- 
     191   !                       ! type of vertical closure 
     192   ln_zdfcst   = .false.      !  constant mixing 
     193   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     194   ln_zdftke   = .true.       !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     195   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     196   ! 
     197   !                       ! convection 
     198   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     199   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     200   ! 
     201   !                       ! time-stepping 
     202   ln_zdfexp   = .false.   ! split-explicit (T) or implicit (F) scheme 
     203   ! 
     204   !                       ! coefficients 
     205   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     206   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     207   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     208   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    200209/ 
    201210!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/OVERFLOW/EXP00/namelist_zps_FCT2_flux_ubs_cfg

    r7640 r7954  
    6666   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    6767                           !                              = 2 : nonlinear friction 
    68 / 
    69 !----------------------------------------------------------------------- 
    70 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    71 !----------------------------------------------------------------------- 
    72 / 
    73 !----------------------------------------------------------------------- 
    74 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    75 !----------------------------------------------------------------------- 
    7668/ 
    7769!----------------------------------------------------------------------- 
     
    191183   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    192184/ 
    193 !----------------------------------------------------------------------- 
    194 &namzdf        !   vertical physics 
    195 !----------------------------------------------------------------------- 
    196    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    197    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    198    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    199    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     185!!====================================================================== 
     186!!                     vertical physics namelists                     !! 
     187!!====================================================================== 
     188!----------------------------------------------------------------------- 
     189&namzdf        !   vertical physics                                     (default: NO selection) 
     190!----------------------------------------------------------------------- 
     191   !                       ! type of vertical closure 
     192   ln_zdfcst   = .true.       !  constant mixing 
     193   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     194   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     195   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     196   ! 
     197   !                       ! convection 
     198   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     199   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     200   ! 
     201   !                       ! time-stepping 
     202   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     203   ! 
     204   !                       ! coefficients 
     205   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     206   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     207   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     208   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    200209/ 
    201210!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/OVERFLOW/EXP00/namelist_zps_FCT4_flux_ubs_cfg

    r7640 r7954  
    6666   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    6767                           !                              = 2 : nonlinear friction 
    68 / 
    69 !----------------------------------------------------------------------- 
    70 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    71 !----------------------------------------------------------------------- 
    72 / 
    73 !----------------------------------------------------------------------- 
    74 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    75 !----------------------------------------------------------------------- 
    7668/ 
    7769!----------------------------------------------------------------------- 
     
    191183   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    192184/ 
    193 !----------------------------------------------------------------------- 
    194 &namzdf        !   vertical physics 
    195 !----------------------------------------------------------------------- 
    196    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    197    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    198    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    199    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     185!!====================================================================== 
     186!!                     vertical physics namelists                     !! 
     187!!====================================================================== 
     188!----------------------------------------------------------------------- 
     189&namzdf        !   vertical physics                                     (default: NO selection) 
     190!----------------------------------------------------------------------- 
     191   !                       ! type of vertical closure 
     192   ln_zdfcst   = .true.       !  constant mixing 
     193   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     194   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     195   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     196   ! 
     197   !                       ! convection 
     198   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     199   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     200   ! 
     201   !                       ! time-stepping 
     202   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     203   ! 
     204   !                       ! coefficients 
     205   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     206   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     207   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     208   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    200209/ 
    201210!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/OVERFLOW/EXP00/namelist_zps_FCT4_vect_een_cfg

    r7640 r7954  
    6666   nn_bfr      =    0      !  type of bottom friction :   = 0 : free slip,  = 1 : linear friction 
    6767                           !                              = 2 : nonlinear friction 
    68 / 
    69 !----------------------------------------------------------------------- 
    70 &nambbc        !   bottom temperature boundary condition                (default: NO) 
    71 !----------------------------------------------------------------------- 
    72 / 
    73 !----------------------------------------------------------------------- 
    74 &nambbl        !   bottom boundary layer scheme                         ("key_trabbl") 
    75 !----------------------------------------------------------------------- 
    7668/ 
    7769!----------------------------------------------------------------------- 
     
    191183   rn_bhm_0      =     1.e+12  !  horizontal bilaplacian eddy viscosity [m4/s] 
    192184/ 
    193 !----------------------------------------------------------------------- 
    194 &namzdf        !   vertical physics 
    195 !----------------------------------------------------------------------- 
    196    rn_avm0     =   1.e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    197    rn_avt0     =   0.     !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    198    ln_zdfevd   = .false.  !  enhanced vertical diffusion (evd) 
    199    ln_zdfnpc   = .false.  !  Non-Penetrative Convective algorithm 
     185!!====================================================================== 
     186!!                     vertical physics namelists                     !! 
     187!!====================================================================== 
     188!----------------------------------------------------------------------- 
     189&namzdf        !   vertical physics                                     (default: NO selection) 
     190!----------------------------------------------------------------------- 
     191   !                       ! type of vertical closure 
     192   ln_zdfcst   = .true.       !  constant mixing 
     193   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     194   ln_zdftke   = .false.      !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     195   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     196   ! 
     197   !                       ! convection 
     198   ln_zdfevd   = .false.      !  enhanced vertical diffusion 
     199   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     200   ! 
     201   !                       ! time-stepping 
     202   ln_zdfexp   = .false.      ! split-explicit (T) or implicit (F) scheme 
     203   ! 
     204   !                       ! coefficients 
     205   rn_avm0     =   1.e-4      !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     206   rn_avt0     =   0.e0       !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     207   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     208   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
    200209/ 
    201210!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/OVERFLOW/cpp_OVERFLOW.fcm

    r7423 r7954  
    1 bld::tool::fppkeys key_zdfcst key_mpp_mpi key_iomput key_nosignedzero 
     1bld::tool::fppkeys   key_mpp_mpi key_iomput key_nosignedzero 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/1_namelist_cfg

    r7821 r7954  
    3333/ 
    3434!----------------------------------------------------------------------- 
    35 &namcrs        !   Grid coarsening for dynamics output and/or 
    36                !   passive tracer coarsened online simulations 
     35&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    3736!----------------------------------------------------------------------- 
    3837/ 
     
    141140&namtra_ldfeiv !   eddy induced velocity param. 
    142141!---------------------------------------------------------------------------------- 
    143    ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
    144    ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
    145    rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    146    nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
    147    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    148    !                                !   =  0           constant  
    149    !                                !   = 10 F(k)      =ldf_c1d  
    150    !                                !   = 20 F(i,j)    =ldf_c2d  
    151    !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    152    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    153142/ 
    154143!----------------------------------------------------------------------- 
     
    209198/ 
    210199!----------------------------------------------------------------------- 
    211 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    212 !----------------------------------------------------------------------- 
    213 / 
    214 !----------------------------------------------------------------------- 
    215 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    216 !----------------------------------------------------------------------- 
    217 / 
    218 !----------------------------------------------------------------------- 
    219 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    220 !----------------------------------------------------------------------- 
    221 / 
    222 !----------------------------------------------------------------------- 
    223200&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
    224201!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/namelist_cfg

    r7822 r7954  
    3333/ 
    3434!----------------------------------------------------------------------- 
    35 &namcrs        !   Grid coarsening for dynamics output and/or 
    36                !   passive tracer coarsened online simulations 
     35&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    3736!----------------------------------------------------------------------- 
    3837/ 
     
    210209/ 
    211210!----------------------------------------------------------------------- 
    212 &namzdf        !   vertical physics 
    213 !----------------------------------------------------------------------- 
    214 / 
    215 !----------------------------------------------------------------------- 
    216 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    217 !----------------------------------------------------------------------- 
    218 / 
    219 !----------------------------------------------------------------------- 
    220 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    221 !----------------------------------------------------------------------- 
    222 / 
    223 !----------------------------------------------------------------------- 
    224 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    225 !----------------------------------------------------------------------- 
    226 / 
    227 !----------------------------------------------------------------------- 
    228 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     211&namzdf        !   vertical physics                                     (default: NO selection) 
     212!----------------------------------------------------------------------- 
     213   !                       ! type of vertical closure 
     214   ln_zdfcst   = .true.       !  constant mixing 
     215/ 
     216!----------------------------------------------------------------------- 
     217&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     218!----------------------------------------------------------------------- 
     219/ 
     220!----------------------------------------------------------------------- 
     221&namzdf_ddm    !   double diffusive mixing parameterization             (ln_zdfddm =T) 
     222!----------------------------------------------------------------------- 
     223/ 
     224!----------------------------------------------------------------------- 
     225&namzdf_tmx    !   tidal mixing parameterization                        (ln_zdftmx =T) 
     226!----------------------------------------------------------------------- 
     227/ 
     228!----------------------------------------------------------------------- 
     229&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    229230!----------------------------------------------------------------------- 
    230231/ 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/WAD/EXP00/namelist_cfg

    r7609 r7954  
    3232/ 
    3333!----------------------------------------------------------------------- 
    34 &namzgr        !   vertical coordinate 
    35 !----------------------------------------------------------------------- 
    36    ln_zps      = .false.   !  z-coordinate - partial steps 
    37    ln_sco      = .true.    !  s-coordinate 
    38 / 
    39 !----------------------------------------------------------------------- 
    4034&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    4135!----------------------------------------------------------------------- 
     
    4640/ 
    4741!----------------------------------------------------------------------- 
    48 &namcrs        !   Grid coarsening for dynamics output and/or 
    49                !   passive tracer coarsened online simulations 
     42&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T) 
    5043!----------------------------------------------------------------------- 
    5144/ 
     
    6255   nn_fsbc     = 1         !  frequency of surface boundary condition computation 
    6356   !                       !     (also = the frequency of sea-ice model call) 
    64    ln_usr      = .true.    !  analytical formulation                    (T => fill namsbc_ana ) 
    65    ln_blk      = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
     57   ln_usr      = .true.    !  analytical formulation                    (T => check usrdef_sbc) 
     58   ln_blk      = .false.   !  CORE bulk formulation                     (T => fill namsbc_blk ) 
    6659   nn_ice      = 0         !  =0 no ice boundary condition   , 
    67    ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf) 
    68    ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
     60   ln_rnf      = .false.   !  runoffs                                   (T => fill namsbc_rnf ) 
     61   ln_ssr      = .false.   !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr ) 
    6962   nn_fwb      = 0         !  FreshWater Budget: =0 unchecked 
    70 / 
    71 !----------------------------------------------------------------------- 
    72 &namsbc_ana    !   analytical surface boundary condition 
    73 !----------------------------------------------------------------------- 
    74    nn_tau000   =   100     !  gently increase the stress over the first ntau_rst time-steps 
    75    rn_utau0    =   0.0e0   !  uniform value for the i-stress 
    76 / 
    77 !----------------------------------------------------------------------- 
    78 &namsbc_flx    !   surface boundary condition : flux formulation 
    79 !----------------------------------------------------------------------- 
    80 / 
    81 !----------------------------------------------------------------------- 
    82 &namsbc_clio   !   namsbc_clio  CLIO bulk formulae 
    83 !----------------------------------------------------------------------- 
    84 / 
    85 !----------------------------------------------------------------------- 
    86 &namsbc_core   !   namsbc_core  CORE bulk formulae 
    87 !----------------------------------------------------------------------- 
    88 / 
    89 !----------------------------------------------------------------------- 
    90 &namsbc_mfs   !   namsbc_mfs  MFS bulk formulae 
    91 !----------------------------------------------------------------------- 
    9263/ 
    9364!----------------------------------------------------------------------- 
     
    350321   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
    351322/ 
    352 !----------------------------------------------------------------------- 
    353 &namzdf        !   vertical physics 
    354 !----------------------------------------------------------------------- 
    355    nn_evdm     =    1      !  evd apply on tracer (=0) or on tracer and momentum (=1) 
    356 / 
    357 !----------------------------------------------------------------------- 
    358 &namzdf_ric    !   richardson number dependent vertical diffusion       ("key_zdfric" ) 
    359 !----------------------------------------------------------------------- 
    360 / 
    361 !----------------------------------------------------------------------- 
    362 &namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
     323!!====================================================================== 
     324!!                     vertical physics namelists                     !! 
     325!!====================================================================== 
     326!!    namzdf        vertical physics 
     327!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T) 
     328!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T) 
     329!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T) 
     330!!    namzdf_tmx    tidal mixing parameterization                       (ln_zdftmx=T) 
     331!!====================================================================== 
     332!----------------------------------------------------------------------- 
     333&namzdf        !   vertical physics                                     (default: NO selection) 
     334!----------------------------------------------------------------------- 
     335   !                       ! type of vertical closure 
     336   ln_zdfcst   = .false.      !  constant mixing 
     337   ln_zdfric   = .false.      !  local Richardson dependent formulation (T =>   fill namzdf_ric) 
     338   ln_zdftke   = .true.       !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke) 
     339   ln_zdfgls   = .false.      !  Generic Length Scale closure           (T =>   fill namzdf_gls) 
     340   ! 
     341   !                       ! convection 
     342   ln_zdfevd   = .true.       !  enhanced vertical diffusion 
     343      nn_evdm     =    1         ! apply on tracer (=0) or on tracer and momentum (=1) 
     344      rn_evd      =  100.        ! mixing coefficient [m2/s] 
     345   ln_zdfnpc   = .false.      !  Non-Penetrative Convective algorithm 
     346      nn_npc      =    1         ! frequency of application of npc 
     347      nn_npcp     =  365         ! npc control print frequency 
     348   ! 
     349   ln_zdfddm   = .false.   ! double diffusive mixing 
     350      rn_avts  =    1.e-4     !  maximum avs (vertical mixing on salinity) 
     351      rn_hsbfr =    1.6       !  heat/salt buoyancy flux ratio 
     352   ! 
     353   ln_zdftmx   = .false.   ! tidal mixing parameterization              (T =>   fill namzdf_tmx) 
     354   ! 
     355   ln_zdfqiao  = .false.   ! surface wave-induced mixing (Qiao et al. 2010) (T =>   ln_wave=ln_sdw=T. & fill namsbc_wave) 
     356   ! 
     357   !                       ! time-stepping 
     358   ln_zdfexp   = .false.   ! split-explicit (T) or implicit (F) scheme 
     359   ! 
     360   !                       ! coefficients 
     361   rn_avm0     =   1.2e-4     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F) 
     362   rn_avt0     =   1.2e-5     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F) 
     363   nn_avb      =    0         !  profile for background avt & avm (=1) or not (=0) 
     364   nn_havtb    =    0         !  horizontal shape for avtb (=1) or not (=0) 
     365/ 
     366!----------------------------------------------------------------------- 
     367&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
    363368!----------------------------------------------------------------------- 
    364369   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    365370/ 
    366 !----------------------------------------------------------------------- 
    367 &namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
    368 !----------------------------------------------------------------------- 
    369 / 
    370 !----------------------------------------------------------------------- 
    371 &namzdf_ddm    !   double diffusive mixing parameterization             ("key_zdfddm") 
    372 !----------------------------------------------------------------------- 
    373 / 
    374 !----------------------------------------------------------------------- 
    375 &namzdf_tmx    !   tidal mixing parameterization                        ("key_zdftmx") 
    376 !----------------------------------------------------------------------- 
    377    ln_tmx_itf  = .false.   !  ITF specific parameterisation 
    378 / 
     371 
     372!!====================================================================== 
     373!!                  ***  Miscellaneous namelists  *** 
     374!!====================================================================== 
    379375!----------------------------------------------------------------------- 
    380376&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     
    411407!!gm 
    412408!----------------------------------------------------------------------- 
    413 &namflo       !   float parameters                                      ("key_float") 
    414 !----------------------------------------------------------------------- 
    415 / 
    416 !----------------------------------------------------------------------- 
    417 &namptr       !   Poleward Transport Diagnostic 
    418 !----------------------------------------------------------------------- 
    419 / 
    420 !----------------------------------------------------------------------- 
    421409&namhsb       !  Heat and salt budgets 
    422410!----------------------------------------------------------------------- 
     
    430418/ 
    431419!----------------------------------------------------------------------- 
    432 &namobs       !  observation usage switch                               ('key_diaobs') 
     420&namobs       !  observation usage switch 
    433421!----------------------------------------------------------------------- 
    434422/ 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/TEST_CASES/WAD/cpp_WAD.fcm

    r7645 r7954  
    1  bld::tool::fppkeys key_zdftke key_iomput key_mpp_mpi key_nosignedzero 
     1 bld::tool::fppkeys   key_iomput key_mpp_mpi key_nosignedzero 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/CONFIG/cfg.txt

    r7953 r7954  
    55ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC 
    66ORCA2_OFF_TRC OPA_SRC OFF_SRC TOP_SRC 
    7 GYRE_PISCES OPA_SRC TOP_SRC 
    87ORCA2_LIM3_PISCES OPA_SRC LIM_SRC_3 TOP_SRC NST_SRC 
    98GYRE_PISCES_XIOS OPA_SRC TOP_SRC 
     9GYRE_PISCES OPA_SRC TOP_SRC 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r7646 r7954  
    167167      CALL zdf_mxl( kt )                                                   ! In any case, we need mxl 
    168168      ! 
    169       hmld(:,:)         = sf_dyn(jf_mld)%fnow(:,:,1) * tmask(:,:,1)    ! mixed layer depht 
    170       avt(:,:,:)        = sf_dyn(jf_avt)%fnow(:,:,:) * tmask(:,:,:)    ! vertical diffusive coefficient  
    171       ! 
    172 #if defined key_trabbl && ! defined key_c1d 
    173       ahu_bbl(:,:)      = sf_dyn(jf_ubl)%fnow(:,:,1) * umask(:,:,1)    ! bbl diffusive coef 
    174       ahv_bbl(:,:)      = sf_dyn(jf_vbl)%fnow(:,:,1) * vmask(:,:,1) 
    175 #endif 
     169      hmld(:,:)       = sf_dyn(jf_mld)%fnow(:,:,1) * tmask(:,:,1)    ! mixed layer depht 
     170      avt(:,:,:)      = sf_dyn(jf_avt)%fnow(:,:,:) * tmask(:,:,:)    ! vertical diffusive coefficient  
     171      ! 
     172      IF( ln_trabbl .AND. .NOT.lk_c1d ) THEN       ! diffusive Bottom boundary layer param  
     173         ahu_bbl(:,:) = sf_dyn(jf_ubl)%fnow(:,:,1) * umask(:,:,1)    ! bbl diffusive coef 
     174         ahv_bbl(:,:) = sf_dyn(jf_vbl)%fnow(:,:,1) * vmask(:,:,1) 
     175      ENDIF 
    176176      ! 
    177177      ! 
     
    275275      ENDIF 
    276276      ! 
    277       IF( lk_trabbl ) THEN 
     277      IF( ln_trabbl ) THEN 
    278278                 jf_ubl  = jfld + 1    ;         jf_vbl  = jfld + 2     ;      jfld = jf_vbl 
    279279           slf_d(jf_ubl) = sn_ubl      ;   slf_d(jf_vbl) = sn_vbl 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90

    r7953 r7954  
    315315 
    316316                            CALL tra_qsr_init   ! penetrative solar radiation qsr 
    317       IF( lk_trabbl     )   CALL tra_bbl_init   ! advective (and/or diffusive) bottom boundary layer scheme 
     317      IF( ln_trabbl     )   CALL tra_bbl_init   ! advective (and/or diffusive) bottom boundary layer scheme 
    318318 
    319319                            CALL trc_nam_run    ! Needed to get restart parameters for passive tracers 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

    r7753 r7954  
    1313   !!             -   ! 2010-11  (G. Madec) add mbk. arrays associated to the deepest ocean level 
    1414   !!             -   ! 2013-04  (F. Roquet, G. Madec)  use of eosbn2 instead of local hard coded alpha and beta 
     15   !!            4.0  ! 2017-04  (G. Madec)  ln_trabbl namelist variable instead of a CPP key 
    1516   !!---------------------------------------------------------------------- 
    16 #if   defined key_trabbl 
    17    !!---------------------------------------------------------------------- 
    18    !!   'key_trabbl'   or                             bottom boundary layer 
     17 
    1918   !!---------------------------------------------------------------------- 
    2019   !!   tra_bbl_alloc : allocate trabbl arrays 
     
    4948   PUBLIC   bbl           !  routine called by trcbbl.F90 and dtadyn.F90 
    5049 
    51    LOGICAL, PUBLIC, PARAMETER ::   lk_trabbl = .TRUE.    !: bottom boundary layer flag 
    52  
    5350   !                                !!* Namelist nambbl * 
     51   LOGICAL , PUBLIC ::   ln_trabbl   !: bottom boundary layer flag 
    5452   INTEGER , PUBLIC ::   nn_bbl_ldf  !: =1   : diffusive bbl or not (=0) 
    5553   INTEGER , PUBLIC ::   nn_bbl_adv  !: =1/2 : advective bbl or not (=0) 
     
    8280      !!                ***  FUNCTION tra_bbl_alloc  *** 
    8381      !!---------------------------------------------------------------------- 
    84       ALLOCATE( utr_bbl  (jpi,jpj) , ahu_bbl  (jpi,jpj) , mbku_d  (jpi,jpj) , mgrhu(jpi,jpj) ,     & 
    85          &      vtr_bbl  (jpi,jpj) , ahv_bbl  (jpi,jpj) , mbkv_d  (jpi,jpj) , mgrhv(jpi,jpj) ,     & 
    86          &      ahu_bbl_0(jpi,jpj) , ahv_bbl_0(jpi,jpj) ,                                          & 
    87          &      e3u_bbl_0(jpi,jpj) , e3v_bbl_0(jpi,jpj) ,                                      STAT=tra_bbl_alloc ) 
     82      ALLOCATE( utr_bbl  (jpi,jpj) , ahu_bbl  (jpi,jpj) , mbku_d(jpi,jpj) , mgrhu(jpi,jpj) ,     & 
     83         &      vtr_bbl  (jpi,jpj) , ahv_bbl  (jpi,jpj) , mbkv_d(jpi,jpj) , mgrhv(jpi,jpj) ,     & 
     84         &      ahu_bbl_0(jpi,jpj) , ahv_bbl_0(jpi,jpj) ,                                        & 
     85         &      e3u_bbl_0(jpi,jpj) , e3v_bbl_0(jpi,jpj) ,                                    STAT=tra_bbl_alloc ) 
    8886         ! 
    8987      IF( lk_mpp            )   CALL mpp_sum ( tra_bbl_alloc ) 
     
    111109      IF( nn_timing == 1 )  CALL timing_start( 'tra_bbl') 
    112110      ! 
    113       IF( l_trdtra )   THEN                         !* Save the input trends 
     111      IF( l_trdtra )   THEN                         !* Save the T-S input trends 
    114112         CALL wrk_alloc( jpi, jpj, jpk, ztrdt, ztrds ) 
    115113         ztrdt(:,:,:) = tsa(:,:,:,jp_tem) 
     
    301299            ! 
    302300         END DO 
    303          !                                                       ! =========== 
    304       END DO                                                     ! end tracer 
    305       !                                                          ! =========== 
     301         !                                                  ! =========== 
     302      END DO                                                ! end tracer 
     303      !                                                     ! =========== 
     304      ! 
    306305      IF( nn_timing == 1 )  CALL timing_stop( 'tra_bbl_adv') 
    307306      ! 
     
    498497      INTEGER ::   ios                  !   -      - 
    499498      REAL(wp), POINTER, DIMENSION(:,:) :: zmbk 
    500       ! 
    501       NAMELIST/nambbl/ nn_bbl_ldf, nn_bbl_adv, rn_ahtbbl, rn_gambbl 
     499      !! 
     500      NAMELIST/nambbl/ ln_trabbl, nn_bbl_ldf, nn_bbl_adv, rn_ahtbbl, rn_gambbl 
    502501      !!---------------------------------------------------------------------- 
    503502      ! 
     
    519518         WRITE(numout,*) 'tra_bbl_init : bottom boundary layer initialisation' 
    520519         WRITE(numout,*) '~~~~~~~~~~~~' 
    521          WRITE(numout,*) '   Namelist nambbl : set bbl parameters' 
    522          WRITE(numout,*) '      diffusive bbl (=1)   or not (=0)    nn_bbl_ldf = ', nn_bbl_ldf 
    523          WRITE(numout,*) '      advective bbl (=1/2) or not (=0)    nn_bbl_adv = ', nn_bbl_adv 
    524          WRITE(numout,*) '      diffusive bbl coefficient           rn_ahtbbl  = ', rn_ahtbbl, ' m2/s' 
    525          WRITE(numout,*) '      advective bbl coefficient           rn_gambbl  = ', rn_gambbl, ' s' 
    526       ENDIF 
    527  
     520         WRITE(numout,*) '       Namelist nambbl : set bbl parameters' 
     521         WRITE(numout,*) '          bottom boundary layer flag          ln_trabbl  = ', ln_trabbl 
     522      ENDIF 
     523      IF( .NOT.ln_trabbl )   RETURN 
     524      ! 
     525      IF(lwp) THEN 
     526         WRITE(numout,*) '          diffusive bbl (=1)   or not (=0)    nn_bbl_ldf = ', nn_bbl_ldf 
     527         WRITE(numout,*) '          advective bbl (=1/2) or not (=0)    nn_bbl_adv = ', nn_bbl_adv 
     528         WRITE(numout,*) '          diffusive bbl coefficient           rn_ahtbbl  = ', rn_ahtbbl, ' m2/s' 
     529         WRITE(numout,*) '          advective bbl coefficient           rn_gambbl  = ', rn_gambbl, ' s' 
     530      ENDIF 
     531      ! 
    528532      !                              ! allocate trabbl arrays 
    529533      IF( tra_bbl_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'tra_bbl_init : unable to allocate arrays' ) 
    530  
     534      ! 
    531535      IF( nn_bbl_adv == 1 )    WRITE(numout,*) '       * Advective BBL using upper velocity' 
    532536      IF( nn_bbl_adv == 2 )    WRITE(numout,*) '       * Advective BBL using velocity = F( delta rho)' 
    533  
     537      ! 
    534538      !                             !* vertical index of  "deep" bottom u- and v-points 
    535539      DO jj = 1, jpjm1                    ! (the "shelf" bottom k-indices are mbku and mbkv) 
     
    544548      zmbk(:,:) = REAL( mbkv_d(:,:), wp )   ;   CALL lbc_lnk(zmbk,'V',1.)   ;   mbkv_d(:,:) = MAX( INT( zmbk(:,:) ), 1 ) 
    545549      CALL wrk_dealloc( jpi, jpj, zmbk ) 
    546  
     550      ! 
    547551      !                                 !* sign of grad(H) at u- and v-points 
    548552      mgrhu(jpi,:) = 0   ;   mgrhu(:,jpj) = 0   ;   mgrhv(jpi,:) = 0   ;   mgrhv(:,jpj) = 0 
     
    565569      ahu_bbl_0(:,:) = rn_ahtbbl * e2_e1u(:,:) * e3u_bbl_0(:,:) * umask(:,:,1) 
    566570      ahv_bbl_0(:,:) = rn_ahtbbl * e1_e2v(:,:) * e3v_bbl_0(:,:) * vmask(:,:,1) 
    567  
    568571      ! 
    569572      IF( nn_timing == 1 )  CALL timing_stop( 'tra_bbl_init') 
    570573      ! 
    571574   END SUBROUTINE tra_bbl_init 
    572  
    573 #else 
    574    !!---------------------------------------------------------------------- 
    575    !!   Dummy module :                      No bottom boundary layer scheme 
    576    !!---------------------------------------------------------------------- 
    577    LOGICAL, PUBLIC, PARAMETER ::   lk_trabbl = .FALSE.   !: bbl flag 
    578 CONTAINS 
    579    SUBROUTINE tra_bbl_init               ! Dummy routine 
    580    END SUBROUTINE tra_bbl_init 
    581    SUBROUTINE tra_bbl( kt )              ! Dummy routine 
    582       WRITE(*,*) 'tra_bbl: You should not have seen this print! error?', kt 
    583    END SUBROUTINE tra_bbl 
    584 #endif 
    585575 
    586576   !!====================================================================== 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf.F90

    r7953 r7954  
    4343#  include "vectopt_loop_substitute.h90" 
    4444   !!---------------------------------------------------------------------- 
    45    !! NEMO/OPA 3.7 , NEMO Consortium (2015) 
     45   !! NEMO/OPA 4.0 , NEMO Consortium (2017) 
    4646   !! $Id$ 
    4747   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
     
    6464      ! 
    6565      IF( neuler == 0 .AND. kt == nit000 ) THEN     ! at nit000 
    66          r2dt =  rdt                          ! = rdt (restarting with Euler time stepping) 
     66         r2dt =  rdt                                     ! = rdt (restarting with Euler time stepping) 
    6767      ELSEIF( kt <= nit000 + 1) THEN                ! at nit000 or nit000+1 
    68          r2dt = 2. * rdt                      ! = 2 rdt (leapfrog) 
     68         r2dt = 2. * rdt                                 ! = 2 rdt (leapfrog) 
    6969      ENDIF 
    7070      ! 
     
    120120      !!          an implicit scheme is required. 
    121121      !!---------------------------------------------------------------------- 
    122       USE zdftke 
    123       USE zdfgls 
    124       !!---------------------------------------------------------------------- 
    125122      ! 
    126123      ! Choice from ln_zdfexp (namzdf namelist variable read in zdfphy module) 
     
    129126      ENDIF 
    130127      ! 
    131       ! Force implicit schemes 
    132       IF( ln_zdftke .OR. ln_zdfgls   )   nzdf = 1   ! TKE, or GLS physics 
     128      ! Force implicit schemes when absolutely needed 
     129      IF( .NOT.ln_zdfcst             )   nzdf = 1   ! turbulent closure schemes 
    133130      IF( ln_traldf_iso              )   nzdf = 1   ! iso-neutral lateral physics 
    134131      IF( ln_traldf_hor .AND. ln_sco )   nzdf = 1   ! horizontal lateral physics in s-coordinate 
    135       IF( ln_zdfexp .AND. nzdf == 1 )   CALL ctl_stop( 'tra_zdf : If using the rotation of lateral mixing operator',   & 
    136             &                         ' GLS or TKE scheme, the implicit scheme is required, set ln_zdfexp = .false.' ) 
    137             ! 
     132      ! 
    138133      IF(lwp) THEN 
    139134         WRITE(numout,*) 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r7953 r7954  
    467467                            CALL tra_qsr_init      ! penetrative solar radiation qsr 
    468468                            CALL tra_bbc_init      ! bottom heat flux 
    469       IF( lk_trabbl     )   CALL tra_bbl_init      ! advective (and/or diffusive) bottom boundary layer scheme 
     469      IF( ln_trabbl     )   CALL tra_bbl_init      ! advective (and/or diffusive) bottom boundary layer scheme 
    470470                            CALL tra_dmp_init      ! internal tracer damping 
    471471                            CALL tra_adv_init      ! horizontal & vertical advection 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7953 r7954  
    232232      IF( ln_traqsr  )   CALL tra_qsr       ( kstp )  ! penetrative solar radiation qsr 
    233233      IF( ln_trabbc  )   CALL tra_bbc       ( kstp )  ! bottom heat flux 
    234       IF( lk_trabbl  )   CALL tra_bbl       ( kstp )  ! advective (and/or diffusive) bottom boundary layer scheme 
     234      IF( ln_trabbl  )   CALL tra_bbl       ( kstp )  ! advective (and/or diffusive) bottom boundary layer scheme 
    235235      IF( ln_tradmp  )   CALL tra_dmp       ( kstp )  ! internal damping trends 
    236236      IF( ln_bdy     )   CALL bdy_tra_dmp   ( kstp )  ! bdy damping trends 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/TOP_SRC/TRP/trcbbl.F90

    r7753 r7954  
    55   !!                                  layer scheme 
    66   !!====================================================================== 
    7    !!============================================================================== 
    87   !! History :  OPA  !  1996-06  (L. Mortier)  Original code 
    98   !!            8.0  !  1997-11  (G. Madec)    Optimization 
     
    1312   !!             -   !  2010-04  (G. Madec)  Campin & Goosse advective bbl  
    1413   !!             -   !  2010-06  (C. Ethe, G. Madec)  merge TRA-TRC 
     14   !!            4.0  !  2017-04  (G. Madec)  ln_trabbl namelist variable instead of a CPP key 
    1515   !!---------------------------------------------------------------------- 
    16 #if  defined key_top &&  defined key_trabbl  
     16#if  defined key_top 
    1717   !!---------------------------------------------------------------------- 
    18    !!   'key_trabbl                      diffusive or/and adevective bottom boundary layer 
     18   !!   'key_top'                                                TOP models 
    1919   !!---------------------------------------------------------------------- 
    20    !!    trc_bbl       : update the tracer trends due to the bottom boundary layer (advective and/or diffusive) 
     20   !!    trc_bbl      : update the tracer trends due to the bottom boundary layer (advective and/or diffusive) 
    2121   !!---------------------------------------------------------------------- 
    22    USE oce_trc             ! ocean dynamics and active tracers variables 
    23    USE trc                 ! ocean passive tracers variables 
    24    USE trabbl              !  
    25    USE prtctl_trc          ! Print control for debbuging 
    26    USE trd_oce 
    27    USE trdtra 
     22   USE oce_trc        ! ocean dynamics and active tracers variables 
     23   USE trc            ! ocean passive tracers variables 
     24   USE trd_oce        ! trends: ocean variables 
     25   USE trdtra         ! tracer trends 
     26   USE trabbl         ! bottom boundary layer  
     27   USE prtctl_trc     ! Print control for debbuging 
    2828 
    29    PUBLIC   trc_bbl       !  routine called by step.F90 
     29   PUBLIC   trc_bbl   !  routine called by trctrp.F90 
    3030 
    3131   !!---------------------------------------------------------------------- 
    32    !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
     32   !! NEMO/TOP 4.0 , NEMO Consortium (2017) 
    3333   !! $Id$  
    3434   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3535   !!---------------------------------------------------------------------- 
    36  
    3736CONTAINS 
    38  
    3937 
    4038   SUBROUTINE trc_bbl( kt ) 
     
    9896   END SUBROUTINE trc_bbl 
    9997 
    100 #else 
    101    !!---------------------------------------------------------------------- 
    102    !!   Dummy module :                      No bottom boundary layer scheme 
    103    !!---------------------------------------------------------------------- 
    104 CONTAINS 
    105    SUBROUTINE trc_bbl( kt )              ! Empty routine 
    106       WRITE(*,*) 'tra_bbl: You should not have seen this print! error?', kt 
    107    END SUBROUTINE trc_bbl 
    10898#endif 
    10999 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90

    r7646 r7954  
    1515   USE oce_trc         ! ocean dynamics and active tracers variables 
    1616   USE trc             ! ocean passive tracers variables  
    17    USE trabbl          ! bottom boundary layer               (trc_bbl routine) 
    1817   USE trcbbl          ! bottom boundary layer               (trc_bbl routine) 
    1918   USE trcdmp          ! internal damping                    (trc_dmp routine) 
     
    6362         ! 
    6463                                CALL trc_sbc    ( kt )      ! surface boundary condition 
    65          IF( lk_trabbl )        CALL trc_bbl    ( kt )      ! advective (and/or diffusive) bottom boundary layer scheme 
     64         IF( ln_trabbl )        CALL trc_bbl    ( kt )      ! advective (and/or diffusive) bottom boundary layer scheme 
    6665         IF( ln_trcdmp )        CALL trc_dmp    ( kt )      ! internal damping trends 
    6766         IF( ln_bdy )           CALL trc_bdy_dmp( kt )      ! BDY damping trends 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/TOP_SRC/trc.F90

    r7953 r7954  
    1717   PUBLIC   trc_alloc   ! called by nemogcm.F90 
    1818 
    19    !! parameters for the control of passive tracers 
    20    !! ---------------------------------------------    
    21    INTEGER, PUBLIC                                                 ::   numnat_ref = -1   !: logical unit for the reference passive tracer namelist_top_ref 
    22    INTEGER, PUBLIC                                                 ::   numnat_cfg = -1   !: logical unit for the reference passive tracer namelist_top_cfg 
    23    INTEGER, PUBLIC                                                 ::   numont     = -1   !: logical unit for the reference passive tracer namelist output output.namelist.top 
    24    INTEGER, PUBLIC                                                 ::   numtrc_ref = -1   !: logical unit for the reference passive tracer namelist_top_ref 
    25    INTEGER, PUBLIC                                                 ::   numtrc_cfg = -1   !: logical unit for the reference passive tracer namelist_top_cfg 
    26    INTEGER, PUBLIC                                                 ::   numonr     = -1   !: logical unit for the reference passive tracer namelist output output.namelist.top 
    27    INTEGER, PUBLIC                                                 ::   numstr        !: logical unit for tracer statistics 
    28    INTEGER, PUBLIC                                                 ::   numrtr        !: logical unit for trc restart (read ) 
    29    INTEGER, PUBLIC                                                 ::   numrtw        !: logical unit for trc restart ( write ) 
     19   !                                     !!- logical units of passive tracers 
     20   INTEGER, PUBLIC ::   numnat_ref = -1   !: reference passive tracer namelist_top_ref 
     21   INTEGER, PUBLIC ::   numnat_cfg = -1   !: reference passive tracer namelist_top_cfg 
     22   INTEGER, PUBLIC ::   numont     = -1   !: reference passive tracer namelist output output.namelist.top 
     23   INTEGER, PUBLIC ::   numtrc_ref = -1   !: reference passive tracer namelist_top_ref 
     24   INTEGER, PUBLIC ::   numtrc_cfg = -1   !: reference passive tracer namelist_top_cfg 
     25   INTEGER, PUBLIC ::   numonr     = -1   !: reference passive tracer namelist output output.namelist.top 
     26   INTEGER, PUBLIC ::   numstr            !: tracer statistics 
     27   INTEGER, PUBLIC ::   numrtr            !: trc restart (read ) 
     28   INTEGER, PUBLIC ::   numrtw            !: trc restart ( write ) 
    3029 
    3130   !! passive tracers fields (before,now,after) 
    3231   !! -------------------------------------------------- 
    33    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)               ::  trai           !: initial total tracer 
    34    REAL(wp), PUBLIC                                                ::  areatot        !: total volume  
    35    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  )         ::  cvol           !: volume correction -degrad option-  
    36    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:)         ::  trn            !: tracer concentration for now time step 
    37    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:)         ::  tra            !: tracer concentration for next time step 
    38    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:)         ::  trb            !: tracer concentration for before time step 
    39    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  )         ::  sbc_trc_b      !: Before sbc fluxes for tracers 
    40    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  )         ::  sbc_trc        !: Now sbc fluxes for tracers 
     32   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)       ::  trai           !: initial total tracer 
     33   REAL(wp), PUBLIC                                        ::  areatot        !: total volume  
     34   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::  cvol           !: volume correction -degrad option-  
     35   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  trn            !: tracer concentration for now time step 
     36   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  tra            !: tracer concentration for next time step 
     37   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  trb            !: tracer concentration for before time step 
     38   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::  sbc_trc_b      !: Before sbc fluxes for tracers 
     39   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::  sbc_trc        !: Now sbc fluxes for tracers 
    4140 
    42    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  )         ::  trc_i          !: prescribed tracer concentration in sea ice for SBC 
    43    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  )         ::  trc_o          !: prescribed tracer concentration in ocean for SBC 
    44    INTEGER             , PUBLIC                                    ::  nn_ice_tr      !: handling of sea ice tracers 
     41   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::  trc_i          !: prescribed tracer concentration in sea ice for SBC 
     42   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::  trc_o          !: prescribed tracer concentration in ocean for SBC 
     43   INTEGER             , PUBLIC                            ::  nn_ice_tr      !: handling of sea ice tracers 
    4544 
    4645   !! interpolated gradient 
    4746   !!--------------------------------------------------   
    48    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)           ::  gtru           !: hor. gradient at u-points at bottom ocean level 
    49    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)           ::  gtrv           !: hor. gradient at v-points at bottom ocean level 
    50    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)           ::  gtrui          !: hor. gradient at u-points at top    ocean level 
    51    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)           ::  gtrvi          !: hor. gradient at v-points at top    ocean level 
    52    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)             ::  qsr_mean        !: daily mean qsr 
     47   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  gtru           !: hor. gradient at u-points at bottom ocean level 
     48   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  gtrv           !: hor. gradient at v-points at bottom ocean level 
     49   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  gtrui          !: hor. gradient at u-points at top    ocean level 
     50   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  gtrvi          !: hor. gradient at v-points at top    ocean level 
     51   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  qsr_mean        !: daily mean qsr 
    5352    
    5453   !! passive tracers  (input and output) 
    5554   !! ------------------------------------------   
    56    LOGICAL             , PUBLIC                                    ::  ln_rsttr       !: boolean term for restart i/o for passive tracers (namelist) 
    57    LOGICAL             , PUBLIC                                    ::  lrst_trc       !: logical to control the trc restart write 
    58    INTEGER             , PUBLIC                                    ::  nn_writetrc    !: time step frequency for concentration outputs (namelist) 
    59    INTEGER             , PUBLIC                                    ::  nutwrs         !: output FILE for passive tracers restart 
    60    INTEGER             , PUBLIC                                    ::  nutrst         !: logical unit for restart FILE for passive tracers 
    61    INTEGER             , PUBLIC                                    ::  nn_rsttr       !: control of the time step ( 0 or 1 ) for pass. tr. 
    62    CHARACTER(len = 80) , PUBLIC                                    ::  cn_trcrst_in   !: suffix of pass. tracer restart name (input) 
    63    CHARACTER(len = 256), PUBLIC                                    ::  cn_trcrst_indir  !: restart input directory 
    64    CHARACTER(len = 80) , PUBLIC                                    ::  cn_trcrst_out  !: suffix of pass. tracer restart name (output) 
    65    CHARACTER(len = 256), PUBLIC                                    ::  cn_trcrst_outdir  !: restart output directory 
    66    REAL(wp)            , PUBLIC                                    ::  rdttrc         !: passive tracer time step 
    67    REAL(wp)            , PUBLIC                                    ::  r2dttrc        !: = 2*rdttrc except at nit000 (=rdttrc) if neuler=0 
    68    LOGICAL             , PUBLIC                                    ::  ln_top_euler   !: boolean term for euler integration  
    69    LOGICAL             , PUBLIC                                    ::  ln_trcdta      !: Read inputs data from files 
    70    LOGICAL             , PUBLIC                                    ::  ln_trcdmp      !: internal damping flag 
    71    LOGICAL             , PUBLIC                                    ::  ln_trcdmp_clo  !: internal damping flag on closed seas 
    72    INTEGER             , PUBLIC                                    ::  nittrc000      !: first time step of passive tracers model 
    73    LOGICAL             , PUBLIC                                    ::  l_trcdm2dc     !: Diurnal cycle for TOP 
     55   LOGICAL             , PUBLIC ::   ln_rsttr           !: boolean term for restart i/o for passive tracers (namelist) 
     56   LOGICAL             , PUBLIC ::   lrst_trc           !: logical to control the trc restart write 
     57   INTEGER             , PUBLIC ::   nn_writetrc        !: time step frequency for concentration outputs (namelist) 
     58   INTEGER             , PUBLIC ::   nutwrs             !: output FILE for passive tracers restart 
     59   INTEGER             , PUBLIC ::   nutrst             !: logical unit for restart FILE for passive tracers 
     60   INTEGER             , PUBLIC ::   nn_rsttr           !: control of the time step ( 0 or 1 ) for pass. tr. 
     61   CHARACTER(len = 80) , PUBLIC ::   cn_trcrst_in       !: suffix of pass. tracer restart name (input) 
     62   CHARACTER(len = 256), PUBLIC ::   cn_trcrst_indir    !: restart input directory 
     63   CHARACTER(len = 80) , PUBLIC ::   cn_trcrst_out      !: suffix of pass. tracer restart name (output) 
     64   CHARACTER(len = 256), PUBLIC ::   cn_trcrst_outdir   !: restart output directory 
     65   REAL(wp)            , PUBLIC ::   rdttrc             !: passive tracer time step 
     66   REAL(wp)            , PUBLIC ::   r2dttrc            !: = 2*rdttrc except at nit000 (=rdttrc) if neuler=0 
     67   LOGICAL             , PUBLIC ::   ln_top_euler       !: boolean term for euler integration  
     68   LOGICAL             , PUBLIC ::   ln_trcdta          !: Read inputs data from files 
     69   LOGICAL             , PUBLIC ::   ln_trcdmp          !: internal damping flag 
     70   LOGICAL             , PUBLIC ::   ln_trcdmp_clo      !: internal damping flag on closed seas 
     71   INTEGER             , PUBLIC ::   nittrc000          !: first time step of passive tracers model 
     72   LOGICAL             , PUBLIC ::   l_trcdm2dc         !: Diurnal cycle for TOP 
    7473 
    7574   !! Information for the ice module for tracers 
     
    8079         CHARACTER(len=2) :: ctrc_o     ! choice of ocean trc cc 
    8180   END TYPE 
    82  
    83    REAL(wp)        , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)  :: trc_ice_ratio      ! ice-ocean tracer ratio 
    84    REAL(wp)        , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)  :: trc_ice_prescr     ! prescribed ice trc cc 
    85    CHARACTER(len=2), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)  :: cn_trc_o ! choice of ocean tracer cc 
     81   ! 
     82   REAL(wp)        , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   trc_ice_ratio    !: ice-ocean tracer ratio 
     83   REAL(wp)        , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   trc_ice_prescr   !: prescribed ice trc cc 
     84   CHARACTER(len=2), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   cn_trc_o         !: choice of ocean tracer cc 
    8685 
    8786 
    8887   !! information for outputs 
    8988   !! -------------------------------------------------- 
    90    TYPE, PUBLIC :: PTRACER                                                            !: Passive tracer type 
     89   TYPE, PUBLIC :: PTRACER                                                         !: Passive tracer type 
    9190       CHARACTER(len = 20)  :: clsname  !: short name 
    9291       CHARACTER(len = 80)  :: cllname  !: long name 
     
    9796       LOGICAL              :: llobc   !: read in a file or not 
    9897   END TYPE PTRACER 
    99  
    100    CHARACTER(len = 20), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ctrcnm         !: tracer name  
    101    CHARACTER(len = 80), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ctrcln         !: trccer field long name 
    102    CHARACTER(len = 20), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ctrcun         !: tracer unit 
    103  
    104    TYPE, PUBLIC :: DIAG                                                               !: passive trcacer ddditional diagnostic type 
     98   ! 
     99   CHARACTER(len = 20), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ctrcnm     !: tracer name  
     100   CHARACTER(len = 80), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ctrcln     !: trccer field long name 
     101   CHARACTER(len = 20), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ctrcun     !: tracer unit 
     102   ! 
     103   TYPE, PUBLIC :: DIAG                                                           !: passive trcacer ddditional diagnostic type 
    105104      CHARACTER(len = 20)  :: sname    !: short name 
    106105      CHARACTER(len = 80)  :: lname    !: long name 
    107106      CHARACTER(len = 20)  :: units    !: unit 
    108107   END TYPE DIAG 
    109  
     108   ! 
    110109   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:)         ::  trc3d          !: 3D diagnostics for tracers 
    111110   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:)         ::  trc2d          !: 2D diagnostics for tracers 
     
    113112   !! information for inputs 
    114113   !! -------------------------------------------------- 
    115    LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ln_trc_ini     !: Initialisation from data input file 
    116    LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ln_trc_obc     !: Use open boundary condition data 
    117    LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ln_trc_sbc     !: Use surface boundary condition data 
    118    LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::  ln_trc_cbc     !: Use coastal boundary condition data 
    119    LOGICAL , PUBLIC                                     ::  ln_rnf_ctl     !: remove runoff dilution on tracers 
    120    REAL(wp), PUBLIC                                     ::  rn_bc_time     !: Time scaling factor for SBC and CBC data (seconds in a day) 
     114   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::   ln_trc_ini     !: Initialisation from data input file 
     115   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::   ln_trc_obc     !: Use open boundary condition data 
     116   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::   ln_trc_sbc     !: Use surface boundary condition data 
     117   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)    ::   ln_trc_cbc     !: Use coastal boundary condition data 
     118   LOGICAL , PUBLIC                                     ::   ln_rnf_ctl     !: remove runoff dilution on tracers 
     119   REAL(wp), PUBLIC                                     ::   rn_bc_time     !: Time scaling factor for SBC and CBC data (seconds in a day) 
    121120 
    122  
    123    !! variables to average over physics over passive tracer sub-steps. 
    124    !! ---------------------------------------------------------------- 
    125    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  un_tm       !: i-horizontal velocity average     [m/s] 
    126    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  vn_tm       !: j-horizontal velocity average     [m/s] 
    127    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  tsn_tm      !: t/s average     [m/s] 
    128    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  avs_tm      !: vertical diffusivity coeff. at  w-point   [m2/s] 
    129    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  rhop_tm     !:  
    130 #if defined key_trabbl 
    131    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  ahu_bbl_tm  !: u-, w-points 
    132    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  ahv_bbl_tm  !: j-direction slope at u-, w-points 
    133    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  utr_bbl_tm  !: j-direction slope at u-, w-points 
    134    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  vtr_bbl_tm  !: j-direction slope at u-, w-points 
    135 #endif 
    136    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  sshn_tm     !: average ssh for the now step [m] 
    137    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  sshb_hold   !:hold sshb from the beginning of each sub-stepping[m]   
    138  
    139    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  rnf_tm     !: river runoff 
    140    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  h_rnf_tm   !: depth in metres to the bottom of the relevant grid box 
    141    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  hmld_tm    !: mixed layer depth average [m] 
    142    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  fr_i_tm    !: average ice fraction     [m/s] 
    143    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  emp_tm     !: freshwater budget: volume flux [Kg/m2/s] 
    144    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  fmmflx_tm  !: freshwater budget: freezing/melting [Kg/m2/s] 
    145    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  emp_b_hold !: hold emp from the beginning of each sub-stepping[m]   
    146    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  qsr_tm     !: solar radiation average [m] 
    147    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  wndm_tm    !: 10m wind average [m] 
    148    ! 
    149  
    150    ! Temporary physical arrays for sub_stepping 
    151    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  tsn_temp 
    152    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  un_temp,vn_temp,wn_temp     !: hold current values of avt, un, vn, wn 
    153    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  avs_temp, rhop_temp     !: hold current values of avt, un, vn, wn 
    154    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  sshn_temp, sshb_temp, ssha_temp, rnf_temp,h_rnf_temp 
    155    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  hdivn_temp, rotn_temp 
    156    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  hdivb_temp, rotb_temp 
    157    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  hmld_temp, qsr_temp, fr_i_temp,wndm_temp 
    158    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  emp_temp, fmmflx_temp, emp_b_temp 
    159    ! 
    160 #if defined key_trabbl 
    161    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  ahu_bbl_temp, ahv_bbl_temp, utr_bbl_temp, vtr_bbl_temp !: hold current values  
    162 #endif 
    163121   ! 
    164122   ! 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

    r7812 r7954  
    88#if defined key_top 
    99   !!---------------------------------------------------------------------- 
    10    !!   trc_stp      : passive tracer system time-stepping 
    11    !!---------------------------------------------------------------------- 
    12    USE oce_trc          ! ocean dynamics and active tracers variables 
     10   !!   trc_stp       : passive tracer system time-stepping 
     11   !!---------------------------------------------------------------------- 
     12   USE oce_trc        ! ocean dynamics and active tracers variables 
    1313   USE sbc_oce 
    1414   USE trc 
    15    USE trctrp           ! passive tracers transport 
    16    USE trcsms           ! passive tracers sources and sinks 
     15   USE trctrp         ! passive tracers transport 
     16   USE trcsms         ! passive tracers sources and sinks 
    1717   USE trcwri 
    1818   USE trcrst 
     19   USE trcsub         ! 
    1920   USE trdtrc_oce 
    2021   USE trdmxl_trc 
    21    USE prtctl_trc       ! Print control for debbuging 
    22    USE iom 
    23    USE in_out_manager 
    24    USE trcsub 
     22   ! 
     23   USE prtctl_trc     ! Print control for debbuging 
     24   USE iom            ! 
     25   USE in_out_manager ! 
    2526 
    2627   IMPLICIT NONE 
     
    2930   PUBLIC   trc_stp    ! called by step 
    3031 
    31    REAL(wp), DIMENSION(:,:,:), SAVE, ALLOCATABLE ::   qsr_arr ! save qsr during TOP time-step 
    32    REAL(wp) :: rdt_sampl 
    33    INTEGER  :: nb_rec_per_day, ktdcy 
    34    REAL(wp) :: rsecfst, rseclast 
    35    LOGICAL  :: llnew 
     32   LOGICAL  ::   llnew                   ! ??? 
     33   REAL(wp) ::   rdt_sampl               ! ??? 
     34   INTEGER  ::   nb_rec_per_day, ktdcy   ! ??? 
     35   REAL(wp) ::   rsecfst, rseclast       ! ??? 
     36   REAL(wp), DIMENSION(:,:,:), SAVE, ALLOCATABLE ::   qsr_arr   ! save qsr during TOP time-step 
    3637 
    3738   !!---------------------------------------------------------------------- 
     
    4647      !!                     ***  ROUTINE trc_stp  *** 
    4748      !!                       
    48       !! ** Purpose : Time loop of opa for passive tracer 
     49      !! ** Purpose :   Time loop of opa for passive tracer 
    4950      !!  
    50       !! ** Method  :  
    51       !!              Compute the passive tracers trends  
    52       !!              Update the passive tracers 
     51      !! ** Method  :   Compute the passive tracers trends  
     52      !!                Update the passive tracers 
    5353      !!------------------------------------------------------------------- 
    54       INTEGER, INTENT( in ) ::  kt      ! ocean time-step index 
    55       INTEGER               ::  jk, jn  ! dummy loop indices 
    56       REAL(wp)              ::  ztrai 
    57       CHARACTER (len=25)    ::  charout  
    58  
     54      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
     55      ! 
     56      INTEGER ::   jk, jn   ! dummy loop indices 
     57      REAL(wp)::   ztrai    ! local scalar 
     58      CHARACTER (len=25) ::   charout   ! 
    5959      !!------------------------------------------------------------------- 
    6060      ! 
     
    115115      ! 
    116116   END SUBROUTINE trc_stp 
     117 
    117118 
    118119   SUBROUTINE trc_mean_qsr( kt ) 
     
    128129      !!              In coupled mode, the sampling is done at every coupling frequency  
    129130      !!---------------------------------------------------------------------- 
    130       INTEGER, INTENT(in) ::   kt 
    131       INTEGER  :: jn 
    132       REAL(wp) :: zkt, zrec 
    133       CHARACTER(len=1)               ::   cl1                      ! 1 character 
    134       CHARACTER(len=2)               ::   cl2                      ! 2 characters 
    135  
     131      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
     132      ! 
     133      INTEGER  ::   jn   ! dummy loop indices 
     134      REAL(wp) ::   zkt, zrec     ! local scalars 
     135      CHARACTER(len=1) ::   cl1   ! 1 character 
     136      CHARACTER(len=2) ::   cl2   ! 2 characters 
     137      !!---------------------------------------------------------------------- 
     138      ! 
    136139      IF( kt == nittrc000 ) THEN 
    137140         IF( ln_cpl )  THEN   
     
    143146         ENDIF 
    144147         ! 
    145          IF( lwp ) THEN 
     148         IF(lwp) THEN 
    146149            WRITE(numout,*)  
    147150            WRITE(numout,*) ' Sampling frequency dt = ', rdt_sampl, 's','   Number of sampling per day  nrec = ', nb_rec_per_day 
     
    171174                    CALL iom_get( numrtr, jpdom_autoglo, 'qsr_arr_'//cl2, qsr_arr(:,:,jn) )   !  A mean of qsr 
    172175                  ENDIF 
    173               ENDDO 
     176              END DO 
    174177            ELSE 
    175178               DO jn = 1, nb_rec_per_day 
     
    184187            DO jn = 1, nb_rec_per_day 
    185188               qsr_arr(:,:,jn) = qsr_mean(:,:) 
    186             ENDDO 
     189            END DO 
    187190         ENDIF 
    188191         ! 
     
    220223               CALL iom_rstput( kt, nitrst, numrtw, 'qsr_arr_'//cl2, qsr_arr(:,:,jn) ) 
    221224             ENDIF 
    222          ENDDO 
     225         END DO 
    223226         CALL iom_rstput( kt, nitrst, numrtw, 'qsr_mean', qsr_mean(:,:) ) 
    224227      ENDIF 
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/TOP_SRC/trcsub.F90

    r7953 r7954  
    22   !!====================================================================== 
    33   !!                       ***  MODULE trcsubstp  *** 
    4    !!TOP :   Averages physics variables for TOP substepping.  
     4   !!   TOP :   Averages physics variables for TOP substepping.  
    55   !!====================================================================== 
    66   !! History :  1.0  !  2011-10  (K. Edwards)  Original 
     
    88#if defined key_top 
    99   !!---------------------------------------------------------------------- 
    10    !!   trc_sub    : passive tracer system sub-stepping  
     10   !!   trc_sub       : passive tracer system sub-stepping  
    1111   !!---------------------------------------------------------------------- 
    12    USE oce_trc          ! ocean dynamics and active tracers variables 
     12   USE oce_trc        ! ocean dynamics and active tracers variables 
    1313   USE trc 
    14    USE prtctl_trc       ! Print control for debbuging 
    15    USE iom 
    16    USE in_out_manager 
    17    USE lbclnk 
    18    USE trabbl 
     14   USE trabbl         ! bottom boundary layer 
    1915   USE zdf_oce 
    2016   USE domvvl 
    21    USE divhor          ! horizontal divergence            (div_hor routine) 
    22    USE sbcrnf    , ONLY: h_rnf, nk_rnf    ! River runoff 
    23    USE bdy_oce   , ONLY: ln_bdy, bdytmask ! BDY 
     17   USE divhor         ! horizontal divergence 
     18   USE sbcrnf   , ONLY: h_rnf, nk_rnf    ! River runoff 
     19   USE bdy_oce  , ONLY: ln_bdy, bdytmask ! BDY 
     20   ! 
     21   USE prtctl_trc     ! Print control for debbuging 
     22   USE in_out_manager !  
     23   USE iom 
     24   USE lbclnk 
    2425#if defined key_agrif 
    2526   USE agrif_opa_update 
     
    2930   IMPLICIT NONE 
    3031 
    31    PUBLIC   trc_sub_stp      ! called by trc_stp 
    32    PUBLIC   trc_sub_ini      ! called by trc_ini to initialize substepping arrays. 
    33    PUBLIC   trc_sub_reset    ! called by trc_stp to reset physics variables 
    34    PUBLIC   trc_sub_ssh      ! called by trc_stp to reset physics variables 
    35  
    36    REAL(wp)  :: r1_ndttrc     !    1 /  nn_dttrc  
    37    REAL(wp)  :: r1_ndttrcp1   !    1 / (nn_dttrc+1)  
    38  
    39    !                                                       !* iso-neutral slopes (if l_ldfslp=T) 
    40    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  uslp_temp, vslp_temp, wslpi_temp, wslpj_temp   !: hold current values  
    41    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  uslp_tm  , vslp_tm  , wslpi_tm  , wslpj_tm     !: time mean  
     32   PUBLIC   trc_sub_stp     ! called by trc_stp 
     33   PUBLIC   trc_sub_ini     ! called by trc_ini to initialize substepping arrays. 
     34   PUBLIC   trc_sub_reset   ! called by trc_stp to reset physics variables 
     35   PUBLIC   trc_sub_ssh     ! called by trc_stp to reset physics variables 
     36 
     37   REAL(wp) :: r1_ndttrc     ! = 1 /  nn_dttrc  
     38   REAL(wp) :: r1_ndttrcp1   ! = 1 / (nn_dttrc+1)  
     39 
     40 
     41   !! averaged and temporary saved variables  (needed when a larger passive tracer time-step is used) 
     42   !! ---------------------------------------------------------------- 
     43   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::    un_tm ,   un_temp   !: i-horizontal velocity average     [m/s] 
     44   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::    vn_tm ,   vn_temp   !: j-horizontal velocity average     [m/s] 
     45   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::              wn_temp   !: hold current values of avt, un, vn, wn 
     46   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::   tsn_tm ,  tsn_temp   !: t/s average     [m/s] 
     47   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   avs_tm ,  avs_temp   !: vertical diffusivity coeff. at  w-point   [m2/s] 
     48   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  rhop_tm , rhop_temp   !:  
     49   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  sshn_tm , sshn_temp   !: average ssh for the now step [m] 
     50 
     51   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::    rnf_tm ,    rnf_temp   !: river runoff 
     52   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  h_rnf_tm ,  h_rnf_temp   !: depth in metres to the bottom of the relevant grid box 
     53   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::   hmld_tm ,   hmld_temp   !: mixed layer depth average [m] 
     54   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::   fr_i_tm ,   fr_i_temp   !: average ice fraction     [m/s] 
     55   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::    emp_tm ,    emp_temp   !: freshwater budget: volume flux [Kg/m2/s] 
     56   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     :: fmmflx_tm , fmmflx_temp   !: freshwater budget: freezing/melting [Kg/m2/s] 
     57   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     :: emp_b_hold,  emp_b_temp   !: hold emp from the beginning of each sub-stepping[m]   
     58   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::    qsr_tm ,    qsr_temp   !: solar radiation average [m] 
     59   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::   wndm_tm ,   wndm_temp   !: 10m wind average [m] 
     60   ! 
     61   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  sshb_hold   !:hold sshb from the beginning of each sub-stepping[m]   
     62   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::   sshb_temp, ssha_temp 
     63   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  hdivn_temp, rotn_temp 
     64   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  hdivb_temp, rotb_temp 
     65   ! 
     66   !                                                    !!- bottom boundary layer param (ln_trabbl=T) 
     67   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  ahu_bbl_tm, ahu_bbl_temp  ! BBL diffusive i-coef. 
     68   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  ahv_bbl_tm, ahv_bbl_temp  ! BBL diffusive j-coef. 
     69   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  utr_bbl_tm, utr_bbl_temp  ! BBL u-advection 
     70   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)     ::  vtr_bbl_tm, vtr_bbl_temp  ! BBL v-advection 
     71 
     72   !                                                      !!- iso-neutral slopes (if l_ldfslp=T) 
     73   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   uslp_temp, vslp_temp, wslpi_temp, wslpj_temp   !: hold current values  
     74   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   uslp_tm  , vslp_tm  , wslpi_tm  , wslpj_tm     !: time mean  
     75 
    4276 
    4377   !!---------------------------------------------------------------------- 
    44    !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
     78   !! NEMO/TOP 4.0 , NEMO Consortium (2017) 
    4579   !! $Id$  
    4680   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     
    5791      !!              on TOP steps, calculate averages. 
    5892      !!------------------------------------------------------------------- 
    59       INTEGER, INTENT( in ) ::  kt        ! ocean time-step index 
    60       INTEGER               ::  ji,jj,jk  ! dummy loop indices 
    61       REAL(wp)              ::  z1_ne3t, z1_ne3u, z1_ne3v, z1_ne3w 
     93      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
     94      ! 
     95      INTEGER ::   ji, jj, jk   ! dummy loop indices 
     96      REAL(wp)::   z1_ne3t, z1_ne3u, z1_ne3v, z1_ne3w   ! local scalars 
    6297      !!------------------------------------------------------------------- 
    6398      ! 
     
    74109           r1_ndttrc        = 1._wp / REAL( nn_dttrc    , wp )  
    75110           r1_ndttrcp1      = 1._wp / REAL( nn_dttrc + 1, wp ) 
    76            ! 
    77111      ENDIF   
    78112 
    79        IF( MOD( kt , nn_dttrc ) /= 0 ) THEN 
    80           ! 
    81           un_tm   (:,:,:)        = un_tm   (:,:,:)        + un   (:,:,:)        * e3u_n(:,:,:)  
    82           vn_tm   (:,:,:)        = vn_tm   (:,:,:)        + vn   (:,:,:)        * e3v_n(:,:,:)  
    83           tsn_tm  (:,:,:,jp_tem) = tsn_tm  (:,:,:,jp_tem) + tsn  (:,:,:,jp_tem) * e3t_n(:,:,:)   
    84           tsn_tm  (:,:,:,jp_sal) = tsn_tm  (:,:,:,jp_sal) + tsn  (:,:,:,jp_sal) * e3t_n(:,:,:)   
    85           rhop_tm (:,:,:)        = rhop_tm (:,:,:)        + rhop (:,:,:)        * e3t_n(:,:,:)   
    86           avs_tm  (:,:,:)        = avs_tm  (:,:,:)        + avs  (:,:,:)        * e3w_n(:,:,:)   
     113      IF( MOD( kt , nn_dttrc ) /= 0 ) THEN 
     114         ! 
     115         un_tm   (:,:,:)        = un_tm   (:,:,:)        + un   (:,:,:)        * e3u_n(:,:,:)  
     116         vn_tm   (:,:,:)        = vn_tm   (:,:,:)        + vn   (:,:,:)        * e3v_n(:,:,:)  
     117         tsn_tm  (:,:,:,jp_tem) = tsn_tm  (:,:,:,jp_tem) + tsn  (:,:,:,jp_tem) * e3t_n(:,:,:)   
     118         tsn_tm  (:,:,:,jp_sal) = tsn_tm  (:,:,:,jp_sal) + tsn  (:,:,:,jp_sal) * e3t_n(:,:,:)   
     119         rhop_tm (:,:,:)        = rhop_tm (:,:,:)        + rhop (:,:,:)        * e3t_n(:,:,:)   
     120         avs_tm  (:,:,:)        = avs_tm  (:,:,:)        + avs  (:,:,:)        * e3w_n(:,:,:)   
    87121         IF( l_ldfslp ) THEN 
    88122            uslp_tm (:,:,:)      = uslp_tm (:,:,:)        + uslp (:,:,:) 
     
    91125            wslpj_tm(:,:,:)      = wslpj_tm(:,:,:)        + wslpj(:,:,:) 
    92126         ENDIF 
    93 # if defined key_trabbl 
    94           IF( nn_bbl_ldf == 1 ) THEN 
    95              ahu_bbl_tm(:,:)     = ahu_bbl_tm(:,:)        + ahu_bbl(:,:)  
    96              ahv_bbl_tm(:,:)     = ahv_bbl_tm(:,:)        + ahv_bbl(:,:)  
    97           ENDIF 
    98           IF( nn_bbl_adv == 1 ) THEN 
    99              utr_bbl_tm(:,:)     = utr_bbl_tm(:,:)        + utr_bbl(:,:)  
    100              vtr_bbl_tm(:,:)     = vtr_bbl_tm(:,:)        + vtr_bbl(:,:)  
    101           ENDIF 
    102 # endif 
    103           ! 
    104           sshn_tm  (:,:)         = sshn_tm  (:,:)         + sshn  (:,:)  
    105           rnf_tm   (:,:)         = rnf_tm   (:,:)         + rnf   (:,:)  
    106           h_rnf_tm (:,:)         = h_rnf_tm (:,:)         + h_rnf (:,:)  
    107           hmld_tm  (:,:)         = hmld_tm  (:,:)         + hmld  (:,:) 
    108           fr_i_tm  (:,:)         = fr_i_tm  (:,:)         + fr_i  (:,:) 
    109           emp_tm   (:,:)         = emp_tm   (:,:)         + emp   (:,:)  
    110           fmmflx_tm(:,:)         = fmmflx_tm(:,:)         + fmmflx(:,:) 
    111           qsr_tm   (:,:)         = qsr_tm   (:,:)         + qsr   (:,:) 
    112           wndm_tm  (:,:)         = wndm_tm  (:,:)         + wndm  (:,:) 
    113  
     127         IF( ln_trabbl ) THEN 
     128            IF( nn_bbl_ldf == 1 ) THEN 
     129               ahu_bbl_tm(:,:)     = ahu_bbl_tm(:,:)        + ahu_bbl(:,:)  
     130               ahv_bbl_tm(:,:)     = ahv_bbl_tm(:,:)        + ahv_bbl(:,:)  
     131            ENDIF 
     132            IF( nn_bbl_adv == 1 ) THEN 
     133               utr_bbl_tm(:,:)     = utr_bbl_tm(:,:)        + utr_bbl(:,:)  
     134               vtr_bbl_tm(:,:)     = vtr_bbl_tm(:,:)        + vtr_bbl(:,:)  
     135            ENDIF 
     136         ENDIF  
     137         ! 
     138         sshn_tm  (:,:)         = sshn_tm  (:,:)         + sshn  (:,:)  
     139         rnf_tm   (:,:)         = rnf_tm   (:,:)         + rnf   (:,:)  
     140         h_rnf_tm (:,:)         = h_rnf_tm (:,:)         + h_rnf (:,:)  
     141         hmld_tm  (:,:)         = hmld_tm  (:,:)         + hmld  (:,:) 
     142         fr_i_tm  (:,:)         = fr_i_tm  (:,:)         + fr_i  (:,:) 
     143         emp_tm   (:,:)         = emp_tm   (:,:)         + emp   (:,:)  
     144         fmmflx_tm(:,:)         = fmmflx_tm(:,:)         + fmmflx(:,:) 
     145         qsr_tm   (:,:)         = qsr_tm   (:,:)         + qsr   (:,:) 
     146         wndm_tm  (:,:)         = wndm_tm  (:,:)         + wndm  (:,:) 
     147         ! 
    114148      ELSE                           !  It is time to substep  
    115          !   1. set temporary arrays to hold physics variables 
     149         !   1. set temporary arrays to hold physics/dynamical variables 
    116150         un_temp    (:,:,:)      = un    (:,:,:) 
    117151         vn_temp    (:,:,:)      = vn    (:,:,:) 
     
    124158            vslp_temp  (:,:,:)   = vslp  (:,:,:)   ;   wslpj_temp (:,:,:)   = wslpj (:,:,:) 
    125159         ENDIF 
    126 # if defined key_trabbl 
    127           IF( nn_bbl_ldf == 1 ) THEN 
    128              ahu_bbl_temp(:,:)   = ahu_bbl(:,:)   
    129              ahv_bbl_temp(:,:)   = ahv_bbl(:,:)  
    130           ENDIF 
    131           IF( nn_bbl_adv == 1 ) THEN 
    132              utr_bbl_temp(:,:)   = utr_bbl(:,:)  
    133              vtr_bbl_temp(:,:)   = vtr_bbl(:,:)  
    134           ENDIF 
    135 # endif 
     160         IF( ln_trabbl ) THEN 
     161            IF( nn_bbl_ldf == 1 ) THEN 
     162               ahu_bbl_temp(:,:)   = ahu_bbl(:,:)   
     163               ahv_bbl_temp(:,:)   = ahv_bbl(:,:)  
     164            ENDIF 
     165            IF( nn_bbl_adv == 1 ) THEN 
     166               utr_bbl_temp(:,:)   = utr_bbl(:,:)  
     167               vtr_bbl_temp(:,:)   = vtr_bbl(:,:)  
     168            ENDIF 
     169         ENDIF  
    136170         sshn_temp  (:,:)        = sshn  (:,:) 
    137171         sshb_temp  (:,:)        = sshb  (:,:) 
     
    162196            wslpj_tm (:,:,:)     = wslpj_tm(:,:,:)        + wslpj(:,:,:)  
    163197         ENDIF 
    164 # if defined key_trabbl 
    165           IF( nn_bbl_ldf == 1 ) THEN 
    166              ahu_bbl_tm(:,:)     = ahu_bbl_tm(:,:)        + ahu_bbl(:,:)  
    167              ahv_bbl_tm(:,:)     = ahv_bbl_tm(:,:)        + ahv_bbl(:,:)  
    168           ENDIF 
    169           IF( nn_bbl_adv == 1 ) THEN 
    170              utr_bbl_tm(:,:)     = utr_bbl_tm(:,:)        + utr_bbl(:,:)  
    171              vtr_bbl_tm(:,:)     = vtr_bbl_tm(:,:)        + vtr_bbl(:,:)  
    172           ENDIF 
    173 # endif 
     198         IF( ln_trabbl ) THEN 
     199            IF( nn_bbl_ldf == 1 ) THEN 
     200               ahu_bbl_tm(:,:)     = ahu_bbl_tm(:,:)        + ahu_bbl(:,:)  
     201               ahv_bbl_tm(:,:)     = ahv_bbl_tm(:,:)        + ahv_bbl(:,:)  
     202            ENDIF 
     203            IF( nn_bbl_adv == 1 ) THEN 
     204               utr_bbl_tm(:,:)     = utr_bbl_tm(:,:)        + utr_bbl(:,:)  
     205               vtr_bbl_tm(:,:)     = vtr_bbl_tm(:,:)        + vtr_bbl(:,:)  
     206            ENDIF 
     207         ENDIF  
    174208         sshn_tm  (:,:)          = sshn_tm    (:,:)       + sshn  (:,:)  
    175209         rnf_tm   (:,:)          = rnf_tm     (:,:)       + rnf   (:,:)  
     
    195229            fmmflx(:,:)          = fmmflx_tm  (:,:) * r1_ndttrc  
    196230            fr_i  (:,:)          = fr_i_tm    (:,:) * r1_ndttrc 
    197 # if defined key_trabbl 
    198             IF( nn_bbl_ldf == 1 ) THEN 
    199                ahu_bbl(:,:)      = ahu_bbl_tm (:,:) * r1_ndttrc   
    200                ahv_bbl(:,:)      = ahv_bbl_tm (:,:) * r1_ndttrc  
    201             ENDIF 
    202             IF( nn_bbl_adv == 1 ) THEN 
    203                utr_bbl(:,:)      = utr_bbl_tm (:,:) * r1_ndttrc   
    204                vtr_bbl(:,:)      = vtr_bbl_tm (:,:) * r1_ndttrc  
    205             ENDIF 
    206 # endif 
     231            IF( ln_trabbl ) THEN 
     232               IF( nn_bbl_ldf == 1 ) THEN 
     233                  ahu_bbl(:,:)      = ahu_bbl_tm (:,:) * r1_ndttrc   
     234                  ahv_bbl(:,:)      = ahv_bbl_tm (:,:) * r1_ndttrc  
     235               ENDIF 
     236               IF( nn_bbl_adv == 1 ) THEN 
     237                  utr_bbl(:,:)      = utr_bbl_tm (:,:) * r1_ndttrc   
     238                  vtr_bbl(:,:)      = vtr_bbl_tm (:,:) * r1_ndttrc  
     239               ENDIF 
     240            ENDIF 
    207241         ELSE 
    208242            wndm  (:,:)          = wndm_tm    (:,:) * r1_ndttrcp1  
     
    211245            fmmflx(:,:)          = fmmflx_tm  (:,:) * r1_ndttrcp1  
    212246            fr_i  (:,:)          = fr_i_tm    (:,:) * r1_ndttrcp1  
    213 # if defined key_trabbl 
    214             IF( nn_bbl_ldf == 1 ) THEN 
    215                ahu_bbl(:,:)      = ahu_bbl_tm (:,:) * r1_ndttrcp1   
    216                ahv_bbl(:,:)      = ahv_bbl_tm (:,:) * r1_ndttrcp1  
    217             ENDIF 
    218             IF( nn_bbl_adv == 1 ) THEN 
    219                utr_bbl(:,:)      = utr_bbl_tm (:,:) * r1_ndttrcp1   
    220                vtr_bbl(:,:)      = vtr_bbl_tm (:,:) * r1_ndttrcp1  
    221             ENDIF 
    222 # endif 
     247            IF( ln_trabbl ) THEN 
     248               IF( nn_bbl_ldf == 1 ) THEN 
     249                  ahu_bbl(:,:)      = ahu_bbl_tm (:,:) * r1_ndttrcp1   
     250                  ahv_bbl(:,:)      = ahv_bbl_tm (:,:) * r1_ndttrcp1  
     251               ENDIF 
     252               IF( nn_bbl_adv == 1 ) THEN 
     253                  utr_bbl(:,:)      = utr_bbl_tm (:,:) * r1_ndttrcp1   
     254                  vtr_bbl(:,:)      = vtr_bbl_tm (:,:) * r1_ndttrcp1  
     255               ENDIF 
     256            ENDIF 
    223257         ENDIF 
    224258         ! 
     
    298332 
    299333      ! Physics variables that are set after initialization: 
    300       fr_i_tm(:,:) = 0._wp 
    301       emp_tm (:,:) = 0._wp 
     334      fr_i_tm  (:,:) = 0._wp 
     335      emp_tm   (:,:) = 0._wp 
    302336      fmmflx_tm(:,:)  = 0._wp 
    303       qsr_tm (:,:) = 0._wp 
    304       wndm_tm(:,:) = 0._wp 
    305 # if defined key_trabbl 
    306       IF( nn_bbl_ldf == 1 ) THEN 
    307          ahu_bbl_tm(:,:) = 0._wp 
    308          ahv_bbl_tm(:,:) = 0._wp 
    309       ENDIF 
    310       IF( nn_bbl_adv == 1 ) THEN 
    311          utr_bbl_tm(:,:) = 0._wp 
    312          vtr_bbl_tm(:,:) = 0._wp 
    313       ENDIF 
    314 # endif 
     337      qsr_tm   (:,:) = 0._wp 
     338      wndm_tm  (:,:) = 0._wp 
     339      IF( ln_trabbl ) THEN 
     340         IF( nn_bbl_ldf == 1 ) THEN 
     341            ahu_bbl_tm(:,:) = 0._wp 
     342            ahv_bbl_tm(:,:) = 0._wp 
     343         ENDIF 
     344         IF( nn_bbl_adv == 1 ) THEN 
     345            utr_bbl_tm(:,:) = 0._wp 
     346            vtr_bbl_tm(:,:) = 0._wp 
     347         ENDIF 
     348      ENDIF 
    315349      ! 
    316350      IF( nn_timing == 1 )  CALL timing_stop('trc_sub_ini') 
     
    359393      qsr   (:,:)     =  qsr_temp   (:,:) 
    360394      wndm  (:,:)     =  wndm_temp  (:,:) 
    361 # if defined key_trabbl 
    362       IF( nn_bbl_ldf == 1 ) THEN 
    363          ahu_bbl(:,:) = ahu_bbl_temp(:,:)  
    364          ahv_bbl(:,:) = ahv_bbl_temp(:,:)  
    365       ENDIF 
    366       IF( nn_bbl_adv == 1 ) THEN 
    367          utr_bbl(:,:) = utr_bbl_temp(:,:)  
    368          vtr_bbl(:,:) = vtr_bbl_temp(:,:)  
    369       ENDIF 
    370 # endif 
     395      IF( ln_trabbl ) THEN 
     396         IF( nn_bbl_ldf == 1 ) THEN 
     397            ahu_bbl(:,:) = ahu_bbl_temp(:,:)  
     398            ahv_bbl(:,:) = ahv_bbl_temp(:,:)  
     399         ENDIF 
     400         IF( nn_bbl_adv == 1 ) THEN 
     401            utr_bbl(:,:) = utr_bbl_temp(:,:)  
     402            vtr_bbl(:,:) = vtr_bbl_temp(:,:)  
     403         ENDIF 
     404      ENDIF 
    371405      ! 
    372406      hdivn (:,:,:)   =  hdivn_temp (:,:,:) 
     
    397431      qsr_tm     (:,:) = qsr   (:,:) 
    398432      wndm_tm    (:,:) = wndm  (:,:) 
    399 # if defined key_trabbl 
    400       IF( nn_bbl_ldf == 1 ) THEN 
    401          ahu_bbl_tm(:,:) = ahu_bbl(:,:)  
    402          ahv_bbl_tm(:,:) = ahv_bbl(:,:)  
    403       ENDIF 
    404       IF( nn_bbl_adv == 1 ) THEN 
    405          utr_bbl_tm(:,:) = utr_bbl(:,:)  
    406          vtr_bbl_tm(:,:) = vtr_bbl(:,:)  
    407       ENDIF 
    408 # endif 
     433      IF( ln_trabbl ) THEN 
     434         IF( nn_bbl_ldf == 1 ) THEN 
     435            ahu_bbl_tm(:,:) = ahu_bbl(:,:)  
     436            ahv_bbl_tm(:,:) = ahv_bbl(:,:)  
     437         ENDIF 
     438         IF( nn_bbl_adv == 1 ) THEN 
     439            utr_bbl_tm(:,:) = utr_bbl(:,:)  
     440            vtr_bbl_tm(:,:) = vtr_bbl(:,:)  
     441         ENDIF 
     442      ENDIF 
    409443      ! 
    410444      ! 
     
    509543      !!------------------------------------------------------------------- 
    510544      USE lib_mpp, ONLY: ctl_warn 
    511       INTEGER ::  ierr 
    512       !!------------------------------------------------------------------- 
    513       ! 
    514       ALLOCATE( un_temp(jpi,jpj,jpk)        ,  vn_temp(jpi,jpj,jpk)  ,   & 
    515          &      wn_temp(jpi,jpj,jpk)        ,                            & 
    516          &      rhop_temp(jpi,jpj,jpk)      ,  rhop_tm(jpi,jpj,jpk) ,   & 
    517          &      sshn_temp(jpi,jpj)          ,  sshb_temp(jpi,jpj) ,      & 
    518          &      ssha_temp(jpi,jpj)          ,                           & 
    519 #if defined key_trabbl 
    520          &      ahu_bbl_temp(jpi,jpj)       ,  ahv_bbl_temp(jpi,jpj),    & 
    521          &      utr_bbl_temp(jpi,jpj)       ,  vtr_bbl_temp(jpi,jpj),    & 
    522 #endif 
    523          &      rnf_temp(jpi,jpj)           ,  h_rnf_temp(jpi,jpj) ,     & 
    524          &      tsn_temp(jpi,jpj,jpk,2)     ,  emp_b_temp(jpi,jpj),      & 
    525          &      emp_temp(jpi,jpj)           ,  fmmflx_temp(jpi,jpj),     & 
    526          &      hmld_temp(jpi,jpj)          ,  qsr_temp(jpi,jpj) ,       & 
    527          &      fr_i_temp(jpi,jpj)          ,  fr_i_tm(jpi,jpj) ,        & 
    528          &      wndm_temp(jpi,jpj)          ,  wndm_tm(jpi,jpj) ,        & 
    529          &      avs_tm(jpi,jpj,jpk)         ,  avs_temp(jpi,jpj,jpk) ,   & 
    530          &      hdivn_temp(jpi,jpj,jpk)     ,  hdivb_temp(jpi,jpj,jpk),  & 
    531          &      un_tm(jpi,jpj,jpk)          ,  vn_tm(jpi,jpj,jpk)  ,     & 
    532          &      sshn_tm(jpi,jpj)            ,  sshb_hold(jpi,jpj) ,      & 
    533          &      tsn_tm(jpi,jpj,jpk,2)       ,                            & 
    534          &      emp_tm(jpi,jpj)             ,  fmmflx_tm(jpi,jpj)  ,     & 
    535          &      emp_b_hold(jpi,jpj)         ,                            & 
    536          &      hmld_tm(jpi,jpj)            ,  qsr_tm(jpi,jpj) ,         & 
    537 #if defined key_trabbl 
    538          &      ahu_bbl_tm(jpi,jpj)         ,  ahv_bbl_tm(jpi,jpj),      & 
    539          &      utr_bbl_tm(jpi,jpj)         ,  vtr_bbl_tm(jpi,jpj),      & 
    540 #endif 
    541          &      rnf_tm(jpi,jpj)             ,  h_rnf_tm(jpi,jpj) , STAT=trc_sub_alloc )   
     545      INTEGER ::  ierr(3) 
     546      !!------------------------------------------------------------------- 
     547      ! 
     548      ierr(:) = 0 
     549      ! 
     550      ALLOCATE( un_temp(jpi,jpj,jpk)      ,  vn_temp(jpi,jpj,jpk)   ,     & 
     551         &      wn_temp(jpi,jpj,jpk)      ,                               & 
     552         &      rhop_temp(jpi,jpj,jpk)    ,  rhop_tm(jpi,jpj,jpk)   ,     & 
     553         &      sshn_temp(jpi,jpj)        ,  sshb_temp(jpi,jpj)     ,     & 
     554         &      ssha_temp(jpi,jpj)        ,                               & 
     555         &      rnf_temp(jpi,jpj)         ,  h_rnf_temp(jpi,jpj)    ,     & 
     556         &      tsn_temp(jpi,jpj,jpk,2)   ,  emp_b_temp(jpi,jpj)    ,     & 
     557         &      emp_temp(jpi,jpj)         ,  fmmflx_temp(jpi,jpj)   ,     & 
     558         &      hmld_temp(jpi,jpj)        ,  qsr_temp(jpi,jpj)      ,     & 
     559         &      fr_i_temp(jpi,jpj)        ,  fr_i_tm(jpi,jpj)       ,     & 
     560         &      wndm_temp(jpi,jpj)        ,  wndm_tm(jpi,jpj)       ,     & 
     561         &      avs_tm(jpi,jpj,jpk)       ,  avs_temp(jpi,jpj,jpk)  ,     & 
     562         &      hdivn_temp(jpi,jpj,jpk)   ,  hdivb_temp(jpi,jpj,jpk),     & 
     563         &      un_tm(jpi,jpj,jpk)        ,  vn_tm(jpi,jpj,jpk)     ,     & 
     564         &      sshn_tm(jpi,jpj)          ,  sshb_hold(jpi,jpj)     ,     & 
     565         &      tsn_tm(jpi,jpj,jpk,2)     ,                               & 
     566         &      emp_tm(jpi,jpj)           ,  fmmflx_tm(jpi,jpj)     ,     & 
     567         &      emp_b_hold(jpi,jpj)       ,                               & 
     568         &      hmld_tm(jpi,jpj)          ,  qsr_tm(jpi,jpj)        ,     & 
     569         &      rnf_tm(jpi,jpj)           ,  h_rnf_tm(jpi,jpj)      , STAT=ierr(1) )   
     570      ! 
     571      IF( l_ldfslp ) THEN 
     572         ALLOCATE( uslp_temp(jpi,jpj,jpk) ,  wslpi_temp(jpi,jpj,jpk),     & 
     573            &      vslp_temp(jpi,jpj,jpk) ,  wslpj_temp(jpi,jpj,jpk),     & 
     574            &      uslp_tm  (jpi,jpj,jpk) ,  wslpi_tm  (jpi,jpj,jpk),     & 
     575            &      vslp_tm  (jpi,jpj,jpk) ,  wslpj_tm  (jpi,jpj,jpk), STAT=ierr(2) ) 
     576      ENDIF 
     577      IF( ln_trabbl ) THEN 
     578         ALLOCATE( ahu_bbl_temp(jpi,jpj)  , utr_bbl_temp(jpi,jpj)   ,     & 
     579            &      ahv_bbl_temp(jpi,jpj)  , vtr_bbl_temp(jpi,jpj)   ,     & 
     580            &      ahu_bbl_tm  (jpi,jpj)  , utr_bbl_tm  (jpi,jpj)   ,     & 
     581            &      ahv_bbl_tm  (jpi,jpj)  , vtr_bbl_tm  (jpi,jpj)   , STAT=ierr(3) )  
     582      ENDIF 
     583      ! 
     584      trc_sub_alloc = MAXVAL( ierr ) 
    542585      ! 
    543586      IF( trc_sub_alloc /= 0 )   CALL ctl_warn('trc_sub_alloc: failed to allocate arrays') 
    544       ! 
    545       IF( l_ldfslp ) THEN 
    546          ALLOCATE( uslp_temp(jpi,jpj,jpk)   ,  wslpi_temp(jpi,jpj,jpk),      & 
    547             &      vslp_temp(jpi,jpj,jpk)   ,  wslpj_temp(jpi,jpj,jpk),      & 
    548             &      uslp_tm  (jpi,jpj,jpk)   ,  wslpi_tm  (jpi,jpj,jpk),      & 
    549             &      vslp_tm  (jpi,jpj,jpk)   ,  wslpj_tm  (jpi,jpj,jpk),  STAT=trc_sub_alloc ) 
    550       ENDIF 
    551       ! 
    552       IF( trc_sub_alloc /= 0 )   CALL ctl_warn('trc_sub_alloc: failed to allocate ldf_slp arrays') 
    553587      ! 
    554588   END FUNCTION trc_sub_alloc 
Note: See TracChangeset for help on using the changeset viewer.