Ignore:
Timestamp:
2017-04-23T09:30:41+02:00 (3 years ago)
Author:
gm
Message:

#1880 (HPC-09): add zdfphy (the ZDF manager) + remove all key_…

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7931 r7953  
    7474      !!              -8- Outputs and diagnostics 
    7575      !!---------------------------------------------------------------------- 
    76       INTEGER ::   ji,jj,jk ! dummy loop indice 
    77       INTEGER ::   indic    ! error indicator if < 0 
    78       INTEGER ::   kcall    ! optional integer argument (dom_vvl_sf_nxt) 
     76      INTEGER ::   ji, jj, jk  ! dummy loop indice 
     77      INTEGER ::   indic        ! error indicator if < 0 
     78      INTEGER ::   kcall        ! optional integer argument (dom_vvl_sf_nxt) 
    7979      !! --------------------------------------------------------------------- 
    8080#if defined key_agrif 
     
    125125                         CALL bn2    ( tsn, rab_n, rn2  ) ! now    Brunt-Vaisala frequency 
    126126 
    127       ! 
    128127      !  VERTICAL PHYSICS 
    129                          CALL zdf_bfr( kstp )         ! bottom friction (if quadratic) 
    130       !                                               ! Vertical eddy viscosity and diffusivity coefficients 
    131       IF( lk_zdfric  )   CALL zdf_ric ( kstp )             ! Richardson number dependent Kz 
    132       IF( lk_zdftke  )   CALL zdf_tke ( kstp )             ! TKE closure scheme for Kz 
    133       IF( lk_zdfgls  )   CALL zdf_gls ( kstp )             ! GLS closure scheme for Kz 
    134       IF( ln_zdfqiao )   CALL zdf_qiao( kstp )             ! Qiao vertical mixing  
    135       ! 
    136       IF( lk_zdfcst  ) THEN                                ! Constant Kz (reset avt, avm[uv] to the background value) 
    137          avt (:,:,:) = rn_avt0 * wmask (:,:,:) 
    138          avm (:,:,:) = rn_avm0 * wmask (:,:,:) 
    139          avmu(:,:,:) = rn_avm0 * wumask(:,:,:) 
    140          avmv(:,:,:) = rn_avm0 * wvmask(:,:,:) 
    141       ENDIF 
    142       ! 
    143       IF( ln_rnf_mouth ) THEN                         ! increase diffusivity at rivers mouths 
    144          DO jk = 2, nkrnf   ;   avt(:,:,jk) = avt(:,:,jk) + 2._wp * rn_avt_rnf * rnfmsk(:,:) * tmask(:,:,jk)   ;   END DO 
    145       ENDIF 
    146       ! 
    147       IF( ln_zdfevd  )   CALL zdf_evd( kstp )         ! enhanced vertical eddy diffusivity 
    148       ! 
    149       IF( ln_zdfddm  ) THEN                           ! double diffusive mixing 
    150                          CALL zdf_ddm( kstp ) 
    151       ELSE                                            ! avs=avt 
    152          DO jk = 2, jpkm1   ;   avs(:,:,jk) = avt(:,:,jk)   ;   END DO 
    153       ENDIF 
    154       ! 
    155       IF( lk_zdftmx  )   CALL zdf_tmx( kstp )         ! tidal vertical mixing 
    156  
    157                          CALL zdf_mxl( kstp )         ! mixed layer depth 
    158  
    159                                                       ! write TKE or GLS information in the restart file 
    160       IF( lrst_oce .AND. lk_zdftke )   CALL tke_rst( kstp, 'WRITE' ) 
    161       IF( lrst_oce .AND. lk_zdfgls )   CALL gls_rst( kstp, 'WRITE' ) 
    162       ! 
     128                         CALL zdf_phy( kstp )         ! vertical physics update (bfr, avt, avs, avm + MLD) 
     129 
    163130      !  LATERAL  PHYSICS 
    164131      ! 
Note: See TracChangeset for help on using the changeset viewer.