Changeset 10094 for NEMO/branches/UKMO/dev_r10037_dynvor_EEUV/src/OCE/DYN/dynspg_ts.F90

Timestamp:

2018-09-06T15:24:19+02:00 (6 years ago)

Author:

davestorkey

Message:

UKMO dev_r10037_dynvor_EEUV branch : add new dynvor scheme.

File:

• NEMO/branches/UKMO/dev_r10037_dynvor_EEUV/src/OCE/DYN/dynspg_ts.F90 (modified) (5 diffs)

NEMO/branches/UKMO/dev_r10037_dynvor_EEUV/src/OCE/DYN/dynspg_ts.F90

@@ -104 +104 @@
       ALLOCATE( wgtbtp1(3*nn_baro), wgtbtp2(3*nn_baro), zwz(jpi,jpj), STAT=ierr(1) )
       !
-      IF( ln_dynvor_een .OR. ln_dynvor_eeT )   &
+      IF( ln_dynvor_een .OR. ln_dynvor_eeT .OR. ln_dynvor_eeUV)   &
          &     ALLOCATE( ftnw(jpi,jpj) , ftne(jpi,jpj) , & 
          &               ftsw(jpi,jpj) , ftse(jpi,jpj) , STAT=ierr(2) )
@@ -156 +156 @@
       REAL(wp) ::   za0, za1, za2, za3              !   -      -
       REAL(wp) ::   zmdi, zztmp            , z1_ht  !   -      -
+      REAL(wp) ::   zmask_ne,zmask_nw,zmask_se,zmask_sw
+      REAL(wp) ::   z1_huv_ne,z1_huv_nw,z1_huv_se,z1_huv_sw
       REAL(wp), DIMENSION(jpi,jpj) :: zsshp2_e, zhf
       REAL(wp), DIMENSION(jpi,jpj) :: zwx, zu_trd, zu_frc, zssh_frc
@@ -283 +285 @@
                   ftse(ji,jj) = ( ff_f(ji  ,jj  ) + ff_f(ji  ,jj-1) + ff_f(ji-1,jj-1) ) * z1_ht
                   ftsw(ji,jj) = ( ff_f(ji  ,jj-1) + ff_f(ji-1,jj-1) + ff_f(ji-1,jj  ) ) * z1_ht
+               END DO
+            END DO
+         CASE( np_EEUV )                  != EEN scheme using e3u and e3v (energy conserving scheme)
+            ftne(1,:) = 0._wp ; ftnw(1,:) = 0._wp ; ftse(1,:) = 0._wp ; ftsw(1,:) = 0._wp
+            DO jj = 2, jpj
+               DO ji = 2, jpi
+                  zmask_ne = MAX(ssumask(ji  ,jj),ssvmask(ji,jj  ))
+                  zmask_nw = MAX(ssumask(ji-1,jj),ssvmask(ji,jj  ))
+                  zmask_se = MAX(ssumask(ji  ,jj),ssvmask(ji,jj-1))
+                  zmask_sw = MAX(ssumask(ji-1,jj),ssvmask(ji,jj-1))
+                  z1_huv_ne = 2._wp * zmask_ne / ( hu_n(ji  ,jj) + hv_n(ji,jj  ) + 1._wp - zmask_ne )
+                  z1_huv_nw = 2._wp * zmask_nw / ( hu_n(ji-1,jj) + hv_n(ji,jj  ) + 1._wp - zmask_nw )
+                  z1_huv_se = 2._wp * zmask_se / ( hu_n(ji  ,jj) + hv_n(ji,jj-1) + 1._wp - zmask_se )
+                  z1_huv_sw = 2._wp * zmask_sw / ( hu_n(ji-1,jj) + hv_n(ji,jj-1) + 1._wp - zmask_sw )
+                  ftne(ji,jj) = ( ff_f(ji-1,jj  ) + ff_f(ji  ,jj  ) + ff_f(ji  ,jj-1) ) * z1_huv_ne
+                  ftnw(ji,jj) = ( ff_f(ji-1,jj-1) + ff_f(ji-1,jj  ) + ff_f(ji  ,jj  ) ) * z1_huv_nw
+                  ftse(ji,jj) = ( ff_f(ji  ,jj  ) + ff_f(ji  ,jj-1) + ff_f(ji-1,jj-1) ) * z1_huv_se
+                  ftsw(ji,jj) = ( ff_f(ji  ,jj-1) + ff_f(ji-1,jj-1) + ff_f(ji-1,jj  ) ) * z1_huv_sw
                END DO
             END DO
@@ -426 +446 @@
          END DO
          !
-      CASE( np_EET , np_EEN )      ! energy & enstrophy scheme (using e3t or e3f)         
+      CASE( np_EET , np_EEN,  np_EEUV )      ! energy & enstrophy scheme (using e3t or e3f or e3u and e3v)         
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -1035 +1055 @@
             END DO
             !
-         CASE( np_EET , np_EEN )   ! energy & enstrophy scheme (using e3t or e3f)
+         CASE( np_EET , np_EEN,  np_EEUV )   ! energy & enstrophy scheme (using e3t or e3f or e3u and e3v)
             DO jj = 2, jpjm1
                DO ji = fs_2, fs_jpim1   ! vector opt.