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

Timestamp:

2016-11-28T11:26:34+01:00 (8 years ago)

Author:

emanuelaclementi

Message:

#1643 removal of surface wn modification due to Stokes Drift in sshsshwzv.F90 and minor modifications in development branch 2015/dev_r5936_INGV1_WAVE

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 
    !!----------------------------------------------------------------------
 
@@ -39 +38 @@
    USE wrk_nemo        ! Memory Allocation
    USE timing          ! Timing
-   USE sbcwave         ! Stokes velocities
 
    IMPLICIT NONE
@@ -162 +160 @@
       REAL(wp), POINTER, DIMENSION(:,:  ) ::  z2d
       REAL(wp), POINTER, DIMENSION(:,:,:) ::  z3d, zhdiv
-      REAL(wp) ::  zdsshu, zdsshv
       !!----------------------------------------------------------------------
       !
@@ -210 +207 @@
       ENDIF
 
-      IF( ln_wave .AND. ln_sdw ) THEN
-         ! Compute d(ssh)/dx  and d(ssh)/dy  
-         ! Compute the surface vertical velocity accounting for the Stokes Drift
-         DO jj = 1 , jpjm1
-            DO ji = 1 , fs_jpim1
-              zdsshu = ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj)
-              zdsshv = ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj)
-              wn(ji,jj,1) = wn(ji,jj,1) + ( usd3dt(ji,jj,1) * zdsshu   &
-                 &                        + vsd3dt(ji,jj,1) * zdsshv   &
-                 &        - wsd3d (ji,jj,1) ) * tmask(ji,jj,1)
-            END DO
-         END DO
-      ENDIF
-      !
 #if defined key_bdy
       IF( lk_bdy ) THEN