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 7334 for branches/UKMO/r5936_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/step.F90 – NEMO

Ignore:
Timestamp:
2016-11-25T10:51:35+01:00 (7 years ago)
Author:
jcastill
Message:

Changes after branch review before merging

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/r5936_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7279 r7334  
    133133      IF( lk_zdftke  )   CALL zdf_tke( kstp )            ! TKE closure scheme for Kz 
    134134      IF( lk_zdfgls  )   CALL zdf_gls( kstp )            ! GLS closure scheme for Kz 
    135       IF( ln_zdfqiao )   THEN 
    136          !Activated Qiao enhanced turbulence but neither ln_wave or ln_sdw are activated  
    137          IF ( .NOT.( ln_wave .AND. ln_sdw ) )   THEN  
    138             CALL ctl_stop ( 'Ask for wave Qiao enhanced turbulence but ln_wave and ln_sdw have to be activated')  
    139          ELSE  
    140             CALL zdf_qiao(kstp )                             ! Qiao vertical mixing  
    141             DO jk = 1, jpkm1 
    142                DO jj = 1, jpj 
    143                   DO ji = 1, jpi 
    144                      avmu(ji,jj,jk) = (avmu(ji,jj,jk) + QBvu(ji,jj,jk)) * umask(ji,jj,jk) 
    145                      avmv(ji,jj,jk) = (avmv(ji,jj,jk) + QBvv(ji,jj,jk)) * vmask(ji,jj,jk) 
    146                      avt( ji,jj,jk) = (avt( ji,jj,jk) + QBv(ji,jj,jk))  * tmask(ji,jj,jk) 
    147                   END DO 
    148                END DO 
    149             END DO 
    150          ENDIF 
    151       ENDIF 
     135      IF( ln_zdfqiao )   CALL zdf_qiao( kstp )             ! Qiao vertical mixing  
    152136      ! 
    153137      IF( lk_zdfcst  ) THEN                              ! Constant Kz (reset avt, avm[uv] to the background value) 
Note: See TracChangeset for help on using the changeset viewer.