Changeset 6673

Timestamp:

2016-06-07T17:07:47+02:00 (8 years ago)

Author:

timgraham

Message:

Bugfix for heat transport from vT

File:

• branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaprod.F90 (modified) (3 diffs)

branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaprod.F90

@@ -112 +112 @@
             DO jj = 2, jpjm1
                DO ji = fs_2, fs_jpim1   ! vector opt.
-                  z3d(ji,jj,jk) = vn(ji,jj,jk) * fse3v(ji,jj,jk) * e1v(ji,jj)
+                  z3d(ji,jj,jk) = z3d(ji,jj,jk) * fse3v(ji,jj,jk) * e1v(ji,jj)
                END DO
             END DO
@@ -148 +148 @@
       ENDIF
 
-      IF( iom_use("vs") ) THEN
+      IF( iom_use("vs") .OR. iom_use("sopst_vs") ) THEN
          z3d(:,:,:) = 0.e0 
          DO jk = 1, jpkm1
@@ -158 +158 @@
          END DO
          CALL iom_put( "vs", z3d )                  ! product of salinity and meridional velocity at V points
+         DO jk = 1, jpkm1
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  z3d(ji,jj,jk) = z3d(ji,jj,jk) * fse3v(ji,jj,jk) * e1v(ji,jj)
+               END DO
+            END DO
+         END DO
+         CALL dia_ptr_ohst_components( jp_sal, 'vts', z3d)
       ENDIF