Re: Selective summation problem

That did it. Thanks!

-Bridget

On Fri, Oct 2, 2009 at 1:20 PM, Rick Brownrigg <brownrig@ucar.edu> wrote:

> Hi Bridget,
> I believe a statement such as:
>
> Tsum(e) = Tsum(e) + temps(t,i,j)
>
> is not "working" because Tsum(e) on the right-hand side of the expression
> contains missing values, which in turn resulted from the third parameter of
> the creation of Tsum:
>
> Tsum = new(16,float,0.)
>
> Note that the optional third param to new() is not an initializer as much
> as it is a specifier for what constitutes missing_values.
>
> One way to get what you're after might be to initialize your summation
> arrays as:
>
> Tsum=new(16,float)
> Tsum=0.
>
> Hope this helps....
>
> Rick
>
> On Oct 2, 2009, at 12:24 PM, Bridget Thrasher wrote:
>
> 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
>
>
>

-- 
Bridget Thrasher, PhD
Postdoctoral Researcher
Climate Central
www.climatecentral.org
(650) 376-7130

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