Changeset 5983 for branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

Timestamp:

2015-12-02T17:12:45+01:00 (8 years ago)

Author:

emanuelaclementi

Message:

ticket #1643 Aligned branch INGV1-WAVE at trunk revision 5936

File:

• branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 (modified) (4 diffs)

branches/2015/dev_r5936_INGV1_WAVE/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 
    !!----------------------------------------------------------------------
 
@@ -38 +39 @@
    USE wrk_nemo        ! Memory Allocation
    USE timing          ! Timing
+   USE sbcwave,  ONLY: usd2dt, vsd2dt,wsd3d
 
    IMPLICIT NONE
@@ -160 +162 @@
       REAL(wp), POINTER, DIMENSION(:,:  ) ::  z2d
       REAL(wp), POINTER, DIMENSION(:,:,:) ::  z3d, zhdiv
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:  ) ::  sshnu, sshnv, dsshnu, dsshnv
       !!----------------------------------------------------------------------
       !
@@ -206 +209 @@
          END DO
       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,jj  ,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) + usd2dt(:,:) * dsshnu(:,:)     &
+                 &           + vsd2dt(:,:) * dsshnv(:,:)      &
+                 &           - ( wsd3d (:,:,1) ) * tmask(:,:,1)
+      ENDIF
+      !
 #if defined key_bdy
       IF( lk_bdy ) THEN