| 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!!!
|
|
|
|
| Re: Partial summation of a data set! [message #51369 is a reply to message #51366] |
Thu, 16 November 2006 03:24  |
enod
Messages: 41
Registered: November 2004
|
Member |
|
|
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
