Changeset 13167


Ignore:
Timestamp:
2020-06-26T18:38:12+02:00 (6 weeks ago)
Author:
techene
Message:

#2385 : pass sette tests without key_qco, oops should be true this time

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/sshwzv.F90

    r13164 r13167  
    307307      ! 
    308308      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    309       REAL(wp)             ::   zCu, zcff, z1_e3t                     ! local scalars 
     309      REAL(wp)             ::   zCu, zcff, z1_e3t, zdt                ! local scalars 
    310310      REAL(wp) , PARAMETER ::   Cu_min = 0.15_wp                      ! local parameters 
    311311      REAL(wp) , PARAMETER ::   Cu_max = 0.30_wp                      ! local parameters 
     
    324324      ! 
    325325      ! Calculate Courant numbers 
     326      zdt = 2._wp * rn_Dt                            ! 2*rn_Dt and not rDt (for restartability) 
    326327      IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 
    327328         DO_3D_00_00( 1, jpkm1 ) 
    328329            z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 
    329             Cu_adv(ji,jj,jk) =   rDt *                                                         & 
     330            Cu_adv(ji,jj,jk) =   zdt *                                                         & 
    330331               &  ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) )            & 
    331332               &  + ( MAX( e2u(ji  ,jj) * e3u(ji  ,jj,jk,Kmm)                                  & 
     
    344345         DO_3D_00_00( 1, jpkm1 ) 
    345346            z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 
    346             Cu_adv(ji,jj,jk) =   rDt *                                                         & 
     347            Cu_adv(ji,jj,jk) =   zdt *                                                      & 
    347348               &  ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) )         & 
    348349               &  + ( MAX( e2u(ji  ,jj)*e3u(ji  ,jj,jk,Kmm)*uu(ji  ,jj,jk,Kmm), 0._wp ) -   & 
Note: See TracChangeset for help on using the changeset viewer.