Abstracting time series from Grided data

[Sachin Ghude]

Hi!

You can try following,


'set x 1'
'set y 1'
'set t 1 57'
'set gxout print'
'set prnopts %g 1 1'
'd tloop(aave(maskout(no2,no2-100),lon=72,lon=72,lat=23,lat=23))'
say result
write ('c:timeseries.txt',result)

This will write to file in txt formate.

to print on screen you can try following

'set x 1'
'set y 1'
'set t 1 57'
'set prnopts %g 1 1'
'd tloop(aave(maskout(no2,no2-100),lon=72,lon=72,lat=23,lat=23))'


bye

Sachin D Ghude
Indian Institute of Tropical Meteorology
Pashan, Pune 411008, INDIA

---------- Original Message -----------
From: Dhaval Prajapati <dhaval.prajapati at ICCSIR.ORG>
To: GRADSUSR at LIST.CINECA.IT
Sent: Fri, 5 Jun 2009 06:07:49 +0000
Subject: Abstracting time series from Grided data

>
>  Dear All,
>
> I want to abstract time series from the grided data in  Grads.
>
> Suppose for lat 23 and long 72
> I want to have time series values in Ascii for this position
> if I give one time at a time
> like set t 1
> or set t 2
> I am getting the value for this time
> But If I use
> set t 1 57
> d rf
> I am getting a time series plote
> But I want all the values to be written on the screen how to do it??????
>
> Thanks
> Dhaval
> Indian Centre for Climate and Societal Impacts Research (ICCSIR)
> Ahmedabad Education Society compound,
> Opp. H.L. Commerce College,
> Navrangpura,
> Ahmedabad 380 009,
> Gujarat,
> India
> Contact No. : +91(079)40045462
>
>
------- End of Original Message -------