| Re: Partial summation of a data set! [message #51366] |
Thu, 16 November 2006 04:22  |
kim20026
Messages: 54
Registered: November 2006
|
Member |
|
|
Thank you, Tian!!! You are right...!!!
tianyf@gmail.com wrote:
> The problem may lies in this:
> total(values[chk])
>
> Try:
> total(values[*,chk])
>
>
>
> Regards,
> Tian
>
>
> DirtyHarry wrote:
>> One more time, thanks for your invaluable suggestions for my previous
>> post. I have another data file with 50 columns * 8767 rows. This is a
>> part of it.
>> Day month year data1 data2 data3 ... data47
>> 1 1 1981 ...
>> .
>> .
>> .
>> 12 30 2005
>>
>> What I am trying to do now is to add up the data in 28 to 35 columns in
>> each year. So what I have done so far is...
>>
>> pro calout
>>
>> ;CONTROL VARIABLES
>> rootpath="D:\IDL_Documents\Temp\" ;working directory
>> ncol= 50
>> nrow= 8767
>> h=''
>> data=fltarr(ncol,nrow-1)
>>
>> close, 1
>> openr, 1, rootpath + 'f1_grow_stratum.daily' ; This is my data file.
>> readf, 1, h
>> readf, 1, data
>> close, 1
>>
>> years=reform(data[2,*]) ;3rd column is year.
>> values=reform(data[28:35,*])
>>
>> for i=1981, 2005 do begin
>> chk=where(years eq i, n_chk)
>> if n_chk eq 0 then continue
>> print, 'Year :', i, ' -> Total = ', total(values[chk])
>> endfor
>> end
>> --------------------------
>> It seemed that there was no problem. No stop, no errors. However,
>> results were totally different form what I expected. Any suggestions or
>> comments? Thanks, in advance!!!
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
