New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7340 for branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90 – NEMO

Ignore:
Timestamp:
2016-11-25T16:41:40+01:00 (7 years ago)
Author:
emanuelaclementi
Message:

#1643 Correction after review in development branch 2015/dev_r5936_INGV1_WAVE

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r7171 r7340  
    3939   USE wrk_nemo        ! Memory Allocation 
    4040   USE timing          ! Timing 
    41    USE sbcwave,  ONLY:  usd3dt, vsd3dt,wsd3d 
     41   USE sbcwave         ! Stokes velocities 
    4242 
    4343   IMPLICIT NONE 
     
    162162      REAL(wp), POINTER, DIMENSION(:,:  ) ::  z2d 
    163163      REAL(wp), POINTER, DIMENSION(:,:,:) ::  z3d, zhdiv 
    164       REAL(wp) ::  dsshnu, dsshnv 
     164      REAL(wp) ::  zdsshu, zdsshv 
    165165      !!---------------------------------------------------------------------- 
    166166      ! 
     
    209209         END DO 
    210210      ENDIF 
    211 ! 
    212 !     In case ln_wave and ln_sdw, the surface vertical velocity is modified 
    213 !     accounting for Sokes Drift velocity 
    214 ! 
    215       IF ( ln_wave .AND. ln_sdw )  THEN 
    216       ! Compute d(ssh)/dx  and d(ssh)/dy   
    217       ! Compute the surface vertical velocity accounting for the Stokes Drift 
    218       !--------------------------------- 
    219        DO jj = 2 , jpjm1 
    220         DO ji = 2 , jpim1 
    221           dsshnu = ( ssha(ji+1,jj) - ssha(ji,jj) ) / e1u(ji,jj) 
    222           dsshnv = ( ssha(ji,jj+1) - ssha(ji,jj) ) / e2v(ji,jj) 
    223           wn(ji,jj,1) = wn(ji,jj,1) +( usd3dt(ji,jj,1) * dsshnu     & 
    224                 &     + vsd3dt(ji,jj,1) * dsshnv                    & 
    225                 &     - ( wsd3d (ji,jj,1) )) * tmask(ji,jj,1) 
    226         ENDDO 
    227        ENDDO 
     211 
     212      IF( ln_wave .AND. ln_sdw ) THEN 
     213         ! Compute d(ssh)/dx  and d(ssh)/dy   
     214         ! Compute the surface vertical velocity accounting for the Stokes Drift 
     215         DO jj = 1 , jpjm1 
     216            DO ji = 1 , fs_jpim1 
     217              zdsshu = ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj) 
     218              zdsshv = ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj) 
     219              wn(ji,jj,1) = wn(ji,jj,1) + ( usd3dt(ji,jj,1) * zdsshu   & 
     220                 &                        + vsd3dt(ji,jj,1) * zdsshv   & 
     221                 &        - wsd3d (ji,jj,1) ) * tmask(ji,jj,1) 
     222            END DO 
     223         END DO 
    228224      ENDIF 
    229225      ! 
Note: See TracChangeset for help on using the changeset viewer.