Changeset 13627 for NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3
- Timestamp:
- 2020-10-16T16:26:17+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/DYN/dynvor.F90
r13621 r13627 647 647 648 648 649 650 649 SUBROUTINE vor_eeT( kt, Kmm, kvor, pu, pv, pu_rhs, pv_rhs ) 651 650 !!---------------------------------------------------------------------- … … 665 664 !! References : Arakawa and Lamb 1980, Mon. Wea. Rev., 109, 18-36 666 665 !!---------------------------------------------------------------------- 667 INTEGER , INTENT(in ) :: kt ! ocean time-step index666 INTEGER , INTENT(in ) :: kt ! ocean time-step index 668 667 INTEGER , INTENT(in ) :: Kmm ! ocean time level index 669 INTEGER , INTENT(in ) :: kvor ! total, planetary, relative, or metric670 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pu, pv ! now velocities671 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pu_rhs, pv_rhs 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 672 671 ! 673 672 INTEGER :: ji, jj, jk ! dummy loop indices … … 682 681 IF( kt == nit000 ) THEN 683 682 IF(lwp) WRITE(numout,*) 684 IF(lwp) WRITE(numout,*) 'dyn:vor_ee n: vorticity term: energy and enstrophy conserving scheme'683 IF(lwp) WRITE(numout,*) 'dyn:vor_eeT : vorticity term: energy and enstrophy conserving scheme' 685 684 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 686 685 ENDIF … … 745 744 zwx(:,:) = e2u(:,:) * e3u(:,:,jk,Kmm) * pu(:,:,jk) 746 745 zwy(:,:) = e1v(:,:) * e3v(:,:,jk,Kmm) * pv(:,:,jk) 747 746 ! 748 747 ! !== compute and add the vorticity term trend =! 749 748 DO_2D( 0, 1, 0, 1 ) … … 875 874 CASE( np_ENE ) ; WRITE(numout,*) ' ==>>> energy conserving scheme (Coriolis at F-points) (ENE)' 876 875 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') 877 877 CASE( np_EET ) ; WRITE(numout,*) ' ==>>> energy conserving scheme (EEN scheme using e3t) (EET)' 878 878 CASE( np_EEN ) ; WRITE(numout,*) ' ==>>> energy and enstrophy conserving scheme (EEN)'
Note: See TracChangeset
for help on using the changeset viewer.