- Timestamp:
- 2016-11-25T16:41:40+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r7171 r7340 39 39 USE wrk_nemo ! Memory Allocation 40 40 USE timing ! Timing 41 USE sbcwave , ONLY: usd3dt, vsd3dt,wsd3d41 USE sbcwave ! Stokes velocities 42 42 43 43 IMPLICIT NONE … … 162 162 REAL(wp), POINTER, DIMENSION(:,: ) :: z2d 163 163 REAL(wp), POINTER, DIMENSION(:,:,:) :: z3d, zhdiv 164 REAL(wp) :: dsshnu, dsshnv164 REAL(wp) :: zdsshu, zdsshv 165 165 !!---------------------------------------------------------------------- 166 166 ! … … 209 209 END DO 210 210 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 228 224 ENDIF 229 225 !
Note: See TracChangeset
for help on using the changeset viewer.