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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7931 for branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfini.F90 – NEMO

Ignore:
Timestamp:
2017-04-19T12:15:07+02:00 (7 years ago)
Author:
gm
Message:

#1880 (HPC-09): remove key_zdfddm + phasing with last changes of HPC08 branch

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfini.F90

    r7646 r7931  
    5151      INTEGER ::   ioptio, ios       ! local integers 
    5252      !! 
    53       NAMELIST/namzdf/ rn_avm0, rn_avt0, nn_avb, nn_havtb, ln_zdfexp, nn_zdfexp,  & 
    54          &        ln_zdfevd, nn_evdm, rn_avevd, ln_zdfnpc, nn_npc, nn_npcp,       & 
    55          &        ln_zdfqiao 
     53      NAMELIST/namzdf/ ln_zdfcst, ln_zdfric, ln_zdftke, ln_zdfgls,   &     ! type of closure scheme 
     54         &             ln_zdfevd, nn_evdm, rn_evd ,                  &     ! convection : evd 
     55         &             ln_zdfnpc, nn_npc , nn_npcp,                  &     ! convection : npc 
     56         &             ln_zdfddm, rn_avts, rn_hsbfr,                 &     ! double diffusion 
     57         &             ln_zdftmx,                                    &     ! tidal mixing 
     58         &             ln_zdfqiao,                                   &     ! surface wave-induced mixing 
     59         &             ln_zdfexp, nn_zdfexp,                         &     ! time-stepping 
     60         &             rn_avm0, rn_avt0, nn_avb, nn_havtb                  ! coefficients 
     61 
     62 
     63!!org      NAMELIST/namzdf/ rn_avm0, rn_avt0, nn_avb, nn_havtb, ln_zdfexp, nn_zdfexp,  & 
     64!!org         &        ln_zdfevd, nn_evdm, rn_avevd, ln_zdfnpc, nn_npc, nn_npcp,       & 
     65!!org         &        ln_zdfqiao 
    5666      !!---------------------------------------------------------------------- 
    5767 
    … …  
    7080         WRITE(numout,*) '~~~~~~~~' 
    7181         WRITE(numout,*) '   Namelist namzdf : set vertical mixing mixing parameters' 
     82         WRITE(numout,*) '      vertical closure scheme' 
     83         WRITE(numout,*) '         constant vertical mixing coefficient    ln_zdfcst = ', ln_zdfcst 
     84         WRITE(numout,*) '         constant vertical mixing coefficient    ln_zdfric = ', ln_zdfric 
     85         WRITE(numout,*) '         constant vertical mixing coefficient    ln_zdftke = ', ln_zdftke 
     86         WRITE(numout,*) '         constant vertical mixing coefficient    ln_zdfgls = ', ln_zdfgls 
     87         WRITE(numout,*) '      convection: ' 
     88         WRITE(numout,*) '         enhanced vertical diffusion             ln_zdfevd = ', ln_zdfevd 
     89         WRITE(numout,*) '            applied on momentum (=1/0)             nn_evdm = ', nn_evdm 
     90         WRITE(numout,*) '            vertical coefficient for evd           rn_evd  = ', rn_evd 
     91         WRITE(numout,*) '         non-penetrative convection (npc)        ln_zdfnpc = ', ln_zdfnpc 
     92         WRITE(numout,*) '            npc call  frequency                    nn_npc  = ', nn_npc 
     93         WRITE(numout,*) '            npc print frequency                    nn_npcp = ', nn_npcp 
     94         WRITE(numout,*) '      double diffusive mixing                    ln_zdfddm = ', ln_zdfddm 
     95         WRITE(numout,*) '         maximum avs for dd mixing                 rn_avts = ', rn_avts 
     96         WRITE(numout,*) '         heat/salt buoyancy flux ratio             rn_hsbfr= ', rn_hsbfr 
     97         WRITE(numout,*) '      surface wave-induced mixing                ln_zdfqiao= ', ln_zdfqiao                                          ! surface wave induced mixing 
     98         WRITE(numout,*) '      tidal mixing                               ln_zdftmx = ', ln_zdftmx 
     99         WRITE(numout,*) '      time splitting / backward scheme           ln_zdfexp = ', ln_zdfexp 
     100         WRITE(numout,*) '         number of sub-time step (ln_zdfexp=T)   nn_zdfexp = ', nn_zdfexp 
     101         WRITE(numout,*) '      coefficients : ' 
    72102         WRITE(numout,*) '      vertical eddy viscosity             rn_avm0   = ', rn_avm0 
    73103         WRITE(numout,*) '      vertical eddy diffusivity           rn_avt0   = ', rn_avt0 
    74104         WRITE(numout,*) '      constant background or profile      nn_avb    = ', nn_avb 
    75105         WRITE(numout,*) '      horizontal variation for avtb       nn_havtb  = ', nn_havtb 
    76          WRITE(numout,*) '      time splitting / backward scheme    ln_zdfexp = ', ln_zdfexp 
    77          WRITE(numout,*) '      number of time step                 nn_zdfexp = ', nn_zdfexp 
    78          WRITE(numout,*) '      enhanced vertical diffusion         ln_zdfevd = ', ln_zdfevd 
    79          WRITE(numout,*) '         applied on momentum (=1/0)       nn_evdm   = ', nn_evdm 
    80          WRITE(numout,*) '      vertical coefficient for evd        rn_avevd  = ', rn_avevd 
    81          WRITE(numout,*) '      non-penetrative convection (npc)    ln_zdfnpc = ', ln_zdfnpc 
    82          WRITE(numout,*) '      npc call  frequency                 nn_npc    = ', nn_npc 
    83          WRITE(numout,*) '      npc print frequency                 nn_npcp   = ', nn_npcp 
    84          WRITE(numout,*) '      Qiao formulation flag               ln_zdfqiao=', ln_zdfqiao 
    85       ENDIF 
     106      ENDIF 
     107 
     108      IF(ln_zdfddm) THEN                    ! double diffusive mixing' 
     109         ALLOCATE( avs(jpi,jpj,jpk) ) 
     110         avs(:,:,:) = rn_avt0 * wmask(:,:,:)  
     111      ENDIF 
     112 
    86113 
    87114      !                          !* Parameter & logical controls 
Note: See TracChangeset for help on using the changeset viewer.