[cdftools] 146/228: JMM +AL+TCD : Bug fix in cdfcensus : reset cumulating variable to 0 between each time frame

[Alastair McKinstry]

This is an automated email from the git hooks/post-receive script.

mckinstry pushed a commit to branch master
in repository cdftools.

commit 0139c4401efa863cfd4ddd5e9a7c12df4c34f994
Author: molines <molines at 1055176f-818a-41d9-83e1-73fbe5b947c5>
Date:   Thu Apr 19 13:57:19 2012 +0000

    JMM +AL+TCD : Bug fix in cdfcensus : reset cumulating variable to 0 between each time frame
    
    
    git-svn-id: http://servforge.legi.grenoble-inp.fr/svn/CDFTOOLS/trunk@594 1055176f-818a-41d9-83e1-73fbe5b947c5
---
 cdfcensus.f90 | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/cdfcensus.f90 b/cdfcensus.f90
index e8a2d1a..45218c4 100644
--- a/cdfcensus.f90
+++ b/cdfcensus.f90
@@ -89,7 +89,6 @@ PROGRAM cdfcensus
   !!----------------------------------------------------------------------
   CALL ReadCdfNames()
 
-  dvoltotal=0.d0
 
   narg = iargc()
   IF ( narg == 0 ) THEN
@@ -191,11 +190,11 @@ PROGRAM cdfcensus
      CASE ( '-srange' )
         CALL getarg(ijarg,cldum) ; READ(cldum,*) zsmin ; ijarg = ijarg+1
         CALL getarg(ijarg,cldum) ; READ(cldum,*) zsmax ; ijarg = ijarg+1
-        CALL getarg(ijarg,cldum) ; READ(cldum,*) zds  ; ijarg = ijarg+1
+        CALL getarg(ijarg,cldum) ; READ(cldum,*) zds   ; ijarg = ijarg+1
      CASE ( '-trange' )
         CALL getarg(ijarg,cldum) ; READ(cldum,*) ztmin ; ijarg = ijarg+1
         CALL getarg(ijarg,cldum) ; READ(cldum,*) ztmax ; ijarg = ijarg+1
-        CALL getarg(ijarg,cldum) ; READ(cldum,*) zdt  ; ijarg = ijarg+1
+        CALL getarg(ijarg,cldum) ; READ(cldum,*) zdt   ; ijarg = ijarg+1
      CASE ( '-full' )
         lfull = .TRUE.
      CASE DEFAULT
@@ -222,7 +221,6 @@ PROGRAM cdfcensus
   ALLOCATE ( dcensus (ns,nt), ddump(ns,nt) )
   ALLOCATE ( rsigma0(ns,nt), rsigma2(ns,nt), rsigma4(ns,nt) )
   ALLOCATE ( zsx (ns,nt), zty(ns,nt), rdumdep(1), tim(npt))
-  dcensus(:,:)=0.d0
 
   ! fill up rsigma0 array with theoretical density
   DO ji=1,ns
@@ -270,6 +268,9 @@ PROGRAM cdfcensus
   ENDIF
 
   DO jt = 1, npt
+     ! reset cumulating variables to 0
+     dcensus(:,:) = 0.d0
+     dvoltotal    = 0.d0
      ! Enter main loop
      DO jk=ik1,ik2
         zt(:,:) = getvar(cf_tfil, cn_votemper, jk, npiglo, npjglo, ktime = jt)

-- 
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-science/packages/cdftools.git