Changeset 1551

Timestamp:

2009-07-29T12:12:02+02:00 (15 years ago)

Author:

smasson

Message:

bugfix in the new 300m heat contain, see ticket:342

File:

• trunk/NEMO/OPA_SRC/DIA/diahth.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/DIA/diahth.F90

@@ -161 +161 @@
       ! deepest layer
       zthick(:,:) = 300. - zthick(:,:)   !   remaining thickness to reach 300m
-      htc3(:,:) = htc3(:,:) + tn(:,:,ilevel+1) * MIN( fse3t(:,:,ilevel+1), zthick(:,:) ) * tmask(:,:,ilevel+1)
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            htc3(ji,jj) = htc3(ji,jj) + tn(ji,jj,ilevel+1) * MIN( fse3t(ji,jj,ilevel+1), zthick(ji,jj) ) * tmask(ji,jj,ilevel+1)
+         END DO
+      END DO
       ! from temperature to heat contain
       zcoef = rau0 * rcp