Selective summation problem

I've got a nested loop that is selecting which points to add to a summation.
However, it is not actually doing to addition, and I can't figure out what
I'm doing wrong. Any help would be appreciated!

Here is my loop, followed by the beginning of the output:
do t=0,1800-1
    print("Averaging month "+t)
    Tcnt = new(16,integer,0)
    Tsum = new(16,float,0.)
    Pcnt = new(16,integer,0)
    Psum = new(16,float,0.)
    do e=0,15
        Tcnt(e) = num(eco(:,:).eq.(e+1))
        Pcnt(e) = num(eco(:,:).eq.(e+1))
    end do
    print("Tcnt="+Tcnt)
    print("Pcnt="+Pcnt)
    do i=0,222-1
        do j=0,462-1
            if ((.not.ismissing(eco(i,j))).and.(eco(i,j).gt.0)) then
                e = eco(i,j) - 1
                print("e="+e+" t="+temps(t,i,j)+" pr="+precip(t,i,j))
                Tsum(e) = Tsum(e) + temps(t,i,j)
                Psum(e) = Psum(e) + precip(t,i,j)
                print("Tsum(e)="+Tsum(e))
                print("Psum(e)="+Psum(e))
            end if
        end do
    end do
    ;print("Tsum="+Tsum)
    ;print("Psum="+Psum)
    ts(t,:) = Tsum/Tcnt
    pr(t,:) = Psum/Pcnt
end do

 Copyright (C) 1995-2009 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 5.1.1
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
(0) Averaging month 0
(0) Tcnt=1023
(1) Tcnt=760
(2) Tcnt=1564
(3) Tcnt=4155
(4) Tcnt=2927
(5) Tcnt=1811
(6) Tcnt=1194
(7) Tcnt=1898
(8) Tcnt=324
(9) Tcnt=753
(10) Tcnt=166
(11) Tcnt=408
(12) Tcnt=1250
(13) Tcnt=1428
(14) Tcnt=828
(15) Tcnt=1148
(0) Pcnt=1023
(1) Pcnt=760
(2) Pcnt=1564
(3) Pcnt=4155
(4) Pcnt=2927
(5) Pcnt=1811
(6) Pcnt=1194
(7) Pcnt=1898
(8) Pcnt=324
(9) Pcnt=753
(10) Pcnt=166
(11) Pcnt=408
(12) Pcnt=1250
(13) Pcnt=1428
(14) Pcnt=828
(15) Pcnt=1148
(0) e=15 t=12.09 pr=0.0921936
(0) Tsum(e)=0
(0) Psum(e)=0
(0) e=15 t=11.89 pr=0.0980645
(0) Tsum(e)=0
(0) Psum(e)=0
(0) e=15 t=12 pr=0.102774
(0) Tsum(e)=0
(0) Psum(e)=0
(0) e=15 t=10.42 pr=0.0926452
(0) Tsum(e)=0
(0) Psum(e)=0
.
.
.

Thanks!
Bridget

-- 
Bridget Thrasher, PhD
Postdoctoral Researcher
Climate Central
www.climatecentral.org

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri, 2 Oct 2009 11:24:12 -0700