Station data to grid data

[Siraj Ul Islam]

Dear All ,

I am also trying to convert my station data into gridded data.
For this I have the stations lat/lon ,time , and variable information but
all this is in txt file. For using any  functions which deals with this kind
of conversion in Grads like
oacres (*gexpr, sexpr)*
oabin (*gexpr, sexpr *)

I first need to make this .txt file readable for Grads. Is any one have
information about this kind of problem


On 11/23/05, Qiong Zheng <qiongzh at gmail.com> wrote:
>
> Hi,
> I have used "oacres", it works very well. One thing you can try is, open
> the grid file (cal.controla.1dom.nc) firstly, then open the station file (
> a.ctl) (don't forget to change the OACRES command!).
> Good Luck
> Qiong
>
>
>  On 11/21/05, Hui Du <huidu at ucdavis.edu> wrote:
> >
> > Hello, grads users
> >   I am trying to convert the station data to grid data using the
> > function 'oacres'.
> >   I created a station file for one variable.
> >   I also successfully got the  map file using 'stnmap'
> >   Then I opened two files to try to get the grid data, the following are
> > my command and error information.
> > *************************************************************************************
> >
> >  ga-> open a.ctl
> > Scanning description file:  a.ctl
> > Data file controla.dat is open as file 1
> > LON set to 0 360
> > LAT set to -90 90
> > LEV set to 500 500
> > Time values set: 1995:8:1:0 1995:8:1:0
> > ga-> sdfopen ../cal.controla.1dom.nc
> > Scanning self-describing file:  ../cal.controla.1dom.nc
> > SDF file ../cal.controla.1dom.nc is open as file 2
> > ga-> set lon 230 250
> > LON set to 230 250
> > ga-> set lat 30 45
> > LAT set to 30 45
> > ga-> set t 100
> > Time values set: 1995:8:25:18 1995:8:25:18
> > ga-> d oacres(t2.2,t2)
> > Segmentation fault (core dumped)
> > *******************************************************************
> > then grads broke, I have no idea what made this happen.
> >
> > For the a.clt, its script is the following
> > ***********************************************
> > dset controla.dat
> > title t2
> > dtype station
> > stnmap t2.map
> > undef -999
> > tdef 1697 linear 00z1aug1995 6hr
> > vars 1
> > t2 0 99 t2 data
> > endvars
> > *****************************************************
> > I don't think it is a good idea to convert my station data because with
> > only one surface variable t2, the size of station data is about 580M.
> >
> > Has anyone have any program to do this outside grads? Or can you tell me
> > where the source code of the function 'oacres' is
> >
> > Thanks a lot
> >
> >
>


--
Siraj Ul Islam
1: Research Fellow (Climate Section)
Global Change Impact Studies Centre
First Floor Saudi Pak tower
Blue Area Islamabad Pakistan
Ph :   92-51-9219785
Fax : 92-51-9219787

2: Mphil Scholar
Computational Physics and Simulation Group
Department of Physics
Quaid e Azam University
Islamabad.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20051123/7c7ef4ea/attachment.html