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 5983 for branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/step.F90 – NEMO

Ignore:
Timestamp:
2015-12-02T17:12:45+01:00 (8 years ago)
Author:
emanuelaclementi
Message:

ticket #1643 Aligned branch INGV1-WAVE at trunk revision 5936

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/step.F90

    r5930 r5983  
    2626   !!            3.6  !  2012-07  (J. Simeon, G. Madec. C. Ethe)  Online coarsening of outputs 
    2727   !!            3.6  !  2014-04  (F. Roquet, G. Madec) New equations of state 
     28   !!            3.6  !  2014-10  (E. Clementi, P. Oddo) Add Qiao vertical mixing in case of waves 
    2829   !!            3.7  !  2014-10  (G. Madec)  LDF simplication  
    2930   !!             -   !  2014-12  (G. Madec) remove KPP scheme 
     
    7576      !!              -8- Outputs and diagnostics 
    7677      !!---------------------------------------------------------------------- 
    77       INTEGER ::   jk      ! dummy loop indice 
     78      INTEGER ::   ji,jj,jk ! dummy loop indice 
    7879      INTEGER ::   indic    ! error indicator if < 0 
    7980      INTEGER ::   kcall    ! optional integer argument (dom_vvl_sf_nxt) 
     
    132133      IF( lk_zdftke  )   CALL zdf_tke( kstp )            ! TKE closure scheme for Kz 
    133134      IF( lk_zdfgls  )   CALL zdf_gls( kstp )            ! GLS closure scheme for Kz 
     135      IF( ln_zdfqiao )   THEN 
     136        CALL zdf_qiao(kstp )                             ! Qiao vertical mixing  
     137         DO jk = 1, jpkm1 
     138          DO jj = 1, jpj 
     139            DO ji = 1, jpi 
     140               avmu(ji,jj,jk) = (avmu(ji,jj,jk) + QBvu(ji,jj,jk)) * umask(ji,jj,jk) 
     141               avmv(ji,jj,jk) = (avmv(ji,jj,jk) + QBvv(ji,jj,jk)) * vmask(ji,jj,jk) 
     142               avt( ji,jj,jk) = (avt( ji,jj,jk) + QBv(ji,jj,jk))  * tmask(ji,jj,jk) 
     143            END DO 
     144          END DO 
     145         END DO 
     146      ENDIF 
     147      ! 
    134148      IF( lk_zdfcst  ) THEN                              ! Constant Kz (reset avt, avm[uv] to the background value) 
    135149         avt (:,:,:) = rn_avt0 * wmask (:,:,:) 
     
    218232                         CALL dyn_vor       ( kstp )  ! vorticity term including Coriolis 
    219233                         CALL dyn_ldf       ( kstp )  ! lateral mixing 
     234      IF( ln_stcor    )  CALL dyn_stcor     ( kstp )  ! Stokes-Coriolis forcing 
    220235                         CALL dyn_hpg       ( kstp )  ! horizontal gradient of Hydrostatic pressure 
    221236                         CALL dyn_spg       ( kstp )  ! surface pressure gradient 
Note: See TracChangeset for help on using the changeset viewer.