From 3004b1b42707fa0b1d36f69770e89bb2e0201602 Mon Sep 17 00:00:00 2001
From: Javier Vegas-Regidor <javier.vegas@bsc.es>
Date: Thu, 12 Jan 2017 15:51:18 +0100
Subject: [PATCH] Fixed error in heatc

---
 src/cdfheatc.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/cdfheatc.f90 b/src/cdfheatc.f90
index 2171be4..815a9fb 100644
--- a/src/cdfheatc.f90
+++ b/src/cdfheatc.f90
@@ -168,13 +168,13 @@ PROGRAM cdfheatc
      PRINT * ,'TIME : ', tim(jt)
      IF (mxloption /= 0) rmxldep(:,:) = getvar(cf_tfil, cn_somxl010, 1, npiglo, npjglo, ktime=jt)
 
-     DO jk = ikmin,ikmax
+     DO jk = 1,nvpk
         ik = jk + ikmin -1
         ! Get velocities v at ik
         temp( :,:)   = getvar(cf_tfil, cn_votemper, ik, npiglo, npjglo, kimin=iimin, kjmin=ijmin, ktime=jt)
         IF ( jt == 1 ) THEN
             PRINT *, 'Read mask'
-            tmask(:,:,jk)   = getvar(cn_fmsk, 'tmask', ik, npiglo, npjglo, kimin=iimin, kjmin=ijmin)
+            tmask(:,:,jk)   = getvar(cn_fmsk, 'tmask', ik, npiglo, npjglo, kimin=iimin, kjmin=ijmin, ktime=jt)
 
             ! get e3t at level ik ( ps...)
             PRINT *, 'Load e3t'