Changeset 13627


Ignore:
Timestamp:
2020-10-16T16:26:17+02:00 (6 months ago)
Author:
techene
Message:

#2385 cosmetics and warning for ENT uncompatibility with vector form

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/DYN/dynvor.F90

    r13621 r13627  
    647647 
    648648 
    649  
    650649   SUBROUTINE vor_eeT( kt, Kmm, kvor, pu, pv, pu_rhs, pv_rhs ) 
    651650      !!---------------------------------------------------------------------- 
     
    665664      !! References : Arakawa and Lamb 1980, Mon. Wea. Rev., 109, 18-36 
    666665      !!---------------------------------------------------------------------- 
    667       INTEGER                         , INTENT(in   ) ::   kt          ! ocean time-step index 
     666      INTEGER                         , INTENT(in   ) ::   kt               ! ocean time-step index 
    668667      INTEGER                         , INTENT(in   ) ::   Kmm              ! ocean time level index 
    669       INTEGER                         , INTENT(in   ) ::   kvor        ! total, planetary, relative, or metric 
    670       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pu, pv    ! now velocities 
    671       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pu_rhs, pv_rhs    ! total v-trend 
     668      INTEGER                         , INTENT(in   ) ::   kvor             ! total, planetary, relative, or metric 
     669      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pu, pv           ! now velocities 
     670      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pu_rhs, pv_rhs   ! total v-trend 
    672671      ! 
    673672      INTEGER  ::   ji, jj, jk     ! dummy loop indices 
     
    682681      IF( kt == nit000 ) THEN 
    683682         IF(lwp) WRITE(numout,*) 
    684          IF(lwp) WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 
     683         IF(lwp) WRITE(numout,*) 'dyn:vor_eeT : vorticity term: energy and enstrophy conserving scheme' 
    685684         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
    686685      ENDIF 
     
    745744         zwx(:,:) = e2u(:,:) * e3u(:,:,jk,Kmm) * pu(:,:,jk) 
    746745         zwy(:,:) = e1v(:,:) * e3v(:,:,jk,Kmm) * pv(:,:,jk) 
    747  
     746         ! 
    748747         !                                   !==  compute and add the vorticity term trend  =! 
    749748         DO_2D( 0, 1, 0, 1 ) 
     
    875874         CASE( np_ENE )   ;   WRITE(numout,*) '   ==>>>   energy conserving scheme (Coriolis at F-points) (ENE)' 
    876875         CASE( np_ENT )   ;   WRITE(numout,*) '   ==>>>   energy conserving scheme (Coriolis at T-points) (ENT)' 
     876                              IF( ln_dynadv_vec )   CALL ctl_warn('dyn_vor_init: ENT scheme may not work in vector form') 
    877877         CASE( np_EET )   ;   WRITE(numout,*) '   ==>>>   energy conserving scheme (EEN scheme using e3t) (EET)' 
    878878         CASE( np_EEN )   ;   WRITE(numout,*) '   ==>>>   energy and enstrophy conserving scheme (EEN)' 
Note: See TracChangeset for help on using the changeset viewer.