Changeset 7099 for branches/UKMO/r5518_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

Timestamp:

2016-10-25T18:15:30+02:00 (7 years ago)

Author:

jcastill

Message:

Changes as in r7078 of the original ING branch: branches/2015/dev_r5936_INGV1_WAVE

File:

• branches/UKMO/r5518_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 (modified) (4 diffs)

branches/UKMO/r5518_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

@@ -9 +9 @@
    !!             -   !  2010-09  (D.Storkey and E.O'Dea) bug fixes for BDY module
    !!            3.3  !  2011-10  (M. Leclair) split former ssh_wzv routine and remove all vvl related work
+   !!            3.6  !  2014-10  (E. Clementi, P. Oddo) add wave contribution to surface vertical velocity
    !!----------------------------------------------------------------------
 
@@ -39 +40 @@
    USE wrk_nemo        ! Memory Allocation
    USE timing          ! Timing
+   USE sbcwave,  ONLY:  zusd2dt, zvsd2dt,wsd3d
 
    IMPLICIT NONE
@@ -164 +166 @@
       REAL(wp), POINTER, DIMENSION(:,:  ) ::  z2d
       REAL(wp), POINTER, DIMENSION(:,:,:) ::  z3d, zhdiv
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:  ) ::  sshnu, sshnv, dsshnu, dsshnv
       !
       INTEGER             ::   ji, jj, jk   ! dummy loop indices
@@ -216 +219 @@
       ENDIF
 
+! 
+!     In case ln_wave and ln_sdw the surface vertical velocity is modified 
+!     accounting for Sokes Drift velocity 
+      IF ( ln_wave .AND. ln_sdw )  THEN 
+        ALLOCATE(sshnu(jpi,jpj),sshnv(jpi,jpj),dsshnu(jpi,jpj),dsshnv(jpi,jpj)) 
+        sshnu (:,:) = 0._wp 
+        sshnv (:,:) = 0._wp 
+        dsshnu(:,:) = 0._wp 
+        dsshnv(:,:) = 0._wp 
+        ! sshn interpolated on U-V grid 
+        !--------------------------------- 
+        DO jj = 1 , jpjm1 
+          DO ji = 1 , jpim1 
+            sshnu(ji,jj) =  0.5 * ( 2. - umask(ji,jj,1) ) * & 
+                         &        ( sshn(ji  ,jj) * tmask(ji  ,jj,1)      & 
+                         &        + sshn(ji+1,jj) * tmask(ji+1,jj,1) ) 
+            sshnv(ji,jj) =  0.5 * ( 2. - vmask(ji,jj,1) ) * & 
+                         &        ( sshn(ji,jj  ) * tmask(ji,jji ,1)      & 
+                         &        + sshn(ji,jj+1) * tmask(ji,jj+1,1) ) 
+          ENDDO 
+        ENDDO 
+        ! Compute d(ssh)/dx  and d(ssh)/dy   
+        !--------------------------------- 
+        DO jj = 1 , jpjm1 
+          DO ji = 1 , jpim1 
+            dsshnu(ji,jj) = ( sshnu(ji+1,jj) - sshnu(ji,jj) ) / e1u(ji,jj) 
+            dsshnv(ji,jj) = ( sshnv(ji,jj+1) - sshnv(ji,jj) ) / e2v(ji,jj) 
+          ENDDO 
+        ENDDO 
+        ! Compute the surface vertical velocity accounting for the Stokes Drift  
+        !--------------------------------------------------------------------- 
+        wn(:,:,1) = wn(:,:,1) + zusd2dt(:,:) * dsshnu(:,:)     & 
+                  &           + zvsd2dt(:,:) * dsshnv(:,:)      & 
+                  &           - ( wsd3d (:,:,1) ) * tmask(:,:,1) 
+      ENDIF 
+      ! 
 #if defined key_bdy
       IF (lk_bdy) THEN