Changeset 2468 for trunk/NEMO/OPA_SRC

Timestamp:

2010-12-08T15:42:16+01:00 (13 years ago)

Author:

cetlod

Message:

v3.2:bug correction on diffusive heat flux for AR5 diag. see ticket #771

File:

• trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/TRA/traldf_iso.F90

@@ -114 +114 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdit, zdjt, ztfw     ! 3D workspace
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zdis, zdjs, zsfw     !  "      "
-# if defined key_diaar5
-# endif  
       !!----------------------------------------------------------------------
 
@@ -253 +251 @@
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
-               z2d(ji,jj) = z2d(ji,jj) + zztmp * zftu(ji,jj,jk) * ( tn(ji,jj,jk) + tn(ji+1,jj,jk) ) * e1u(ji,jj) * fse3u(ji,jj,jk) 
-            END DO
-         END DO
-      END DO
+               z2d(ji,jj) = z2d(ji,jj) + zftu(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+      z2d(:,:) =  z2d(:,:) * zztmp
       CALL lbc_lnk( z2d, 'U', -1. )
       CALL iom_put( "udiff_heattr", z2d )                  ! heat transport in i-direction
@@ -263 +262 @@
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
-               z2d(ji,jj) = z2d(ji,jj) + zztmp * zftv(ji,jj,jk) * ( tn(ji,jj,jk) + tn(ji,jj+1,jk) ) * e2v(ji,jj) * fse3v(ji,jj,jk) 
-            END DO
-         END DO
-      END DO
+               z2d(ji,jj) = z2d(ji,jj) + zftv(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+      z2d(:,:) =  z2d(:,:) * zztmp
       CALL lbc_lnk( z2d, 'V', -1. )
       CALL iom_put( "vdiff_heattr", z2d )                  !  heat transport in i-direction