Changeset 7422 for branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

Timestamp:

2016-12-01T18:17:41+01:00 (7 years ago)

Author:

gm

Message:

#1805 dev_INGV_UKMO_2016: improve Stokes drift (including dynspg_ts , Stokes-Coriolis force , and GLS surface roughness

File:

• branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90 (modified) (1 diff)

branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

@@ -107 +107 @@
       !                                         !==  effective transport  ==!
       IF( ln_wave .AND. ln_sdw )  THEN
-         DO jk = 1, jpkm1
-            zun(:,:,jk) = e2u(:,:) * e3u_n(:,:,jk) *      &
-                        &  ( un(:,:,jk) + usd3d(:,:,jk) )                       ! eulerian transport + Stokes Drift
-            zvn(:,:,jk) = e1v(:,:) * e3v_n(:,:,jk) *      &
-                        &  ( vn(:,:,jk) + vsd3d(:,:,jk) )
-            zwn(:,:,jk) = e1e2t(:,:) *                    &
-                        &  ( wn(:,:,jk) + wsd3d(:,:,jk) )
+         DO jk = 1, jpkm1                                                       ! eulerian transport + Stokes Drift
+            zun(:,:,jk) = e2u  (:,:) * e3u_n(:,:,jk) * ( un(:,:,jk) + usd(:,:,jk) )
+            zvn(:,:,jk) = e1v  (:,:) * e3v_n(:,:,jk) * ( vn(:,:,jk) + vsd(:,:,jk) )
+            zwn(:,:,jk) = e1e2t(:,:)                 * ( wn(:,:,jk) + wsd(:,:,jk) )
          END DO
       ELSE