station map

Mon Jun 1 23:05:17 EDT 2009

Dear M.Mahakur,

Thank you, The picture looks very good.
The  radii influence value I set is 8, and it displayed normally.
Thank you for your time, you are so kind hearted.

Best wishes.
2009/6/1 M.Mahakur <mmahakur at tropmet.res.in>

> Hello,
>
> Yes, for station data the default display domain is global. The small lines
> are due to
> exitance undef grid boxes. That is due to no. of stations not well
> distributed or less.
> You can try starting with 2-3 large scan radii, and many small radii later
> to have the
> fine structures, check if plot is ok after that. you can use
> smth9(oacres(....)) also
> for a smooth plot (using the 9-point smoother)..
>
> With Regards,
>
> --
> M. Mahakur
> Forecasting Research Division
> INDIAN INSTITUTE OF TROPICAL METEOROLOGY
> (Ministry of Earth Sciences, Govt. of India)
> Dr. Homi Bhabha Road, NCL Post, Pashan
> PUNE - 411 008, INDIA
>
>
> ---------- Original Message -----------
> From: sh zhou <sh.maeam at GMAIL.COM>
> To: GRADSUSR at LIST.CINECA.IT
> Sent: Mon, 1 Jun 2009 22:57:50 +0800
> Subject: Re: station map
>
> > Dear, M. Mahakur,
> >
> > Thank you for your time and replay.
> > Yes, as you said, the resolution I set is too high, and I made a mistake
> at
> > the same time when set the area location, that the stations' location are
> > out of " background grid".
> > The radii of influence should be set between "0" and "50", my Grads
> version
> > is 1.8sl11, maybe in other versions are different values.
> > Now, I had set background grid data like this:
> > ###
> > xdef  350 linear  104.61  0.0450
> > ydef  250 linear   21.98  0.0450
> > ###
> > then the display seems much better now, but there are so many small lines
> > appeared, why does those small lines appear?
> > In this web page " http://www.iges.org/grads/gadoc/users.html" ( my
> manual
> > ), I haven't found the reason, If you know, could you tell me?
> > By the way, in my first step, when display station data, the background
> is
> > the whole word, is that default display?
> >
> > Thank you for your time.
> >
> > Sincerely.
> >
> > 2009/6/1 M.Mahakur <mmahakur at tropmet.res.in>
> >
> > > Hello,
> > >
> > > How many stations do you have to inerpolate to such a hi resolution map
> > > (0.05 x 0.01)?
> > > Try with several scan radii, starting from a bigger for example (I
> guess
> > > there is limit
> > > to nos. of radii and minimum value of radii, see if 0.01 is allowed,
> > > although I have
> > > written!!):
> > >
> > > qq=oacres(r,rain.2, 5,4,2,1,0.5,0.3
> > > ,0.1,0.08,0.06,0.05,0.04,0.03,0.02,0.01)
> > >
> > > Default radii are 10,7,4,2,1
> > >
> > > http://www.iges.org/grads/gadoc/gradfuncoacres.html
> > >
> > > --
> > > M. Mahakur
> > >
> > > ---------- Original Message -----------
> > > From: sh zhou <sh.maeam at GMAIL.COM>
> > > To: GRADSUSR at LIST.CINECA.IT
> > > Sent: Mon, 1 Jun 2009 20:49:19 +0800
> > > Subject: station map
> > >
> > > > Dear All,
> > > >
> > > > I want to draw a station map,  firstly, I created a binary file, and
> > > wrote a
> > > > "ctl" file as following:
> > > > ###
> > > > dset st_data.dat
> > > > title station data
> > > > dtype station
> > > > stnmap station_data.map
> > > > undef -999.9
> > > > tdef 1 linear 24May2008 7dy
> > > > vars 1
> > > > rain 0 99 rain data
> > > > endvars
> > > > ###
> > > > but the picture's background is the whole world( there is a picture
> in
> > > > attachment), and stations values displayed in a small area,  when I
> added
> > > > "xdef,yedf,zdef" like this:
> > > >
> > > > ###
> > > > ...
> > > > undef -999.9
> > > > xdef 100 linear 100.0 0.05
> > > > ydef 100 linear 20.0 0.01
> > > > zdef 1 linear 1000.0
> > > > tdef 1 linear 24May2008 7dy
> > > > ...
> > > > ###
> > > > however, the background didn't change, Does "xdef" and "ydef"'s range
> > >  set
> > > > by default as following?
> > > > ###
> > > > ga-> open station.ctl
> > > > Scanning description file:  station.ctl
> > > > Data file st_data.dat is open as file 1
> > > > LON set to 0 360
> > > > LAT set to -90 90
> > > > LEV set to 500 500
> > > > Time values set: 2008:5:24:0 2008:5:24:0
> > > > ga-> q file
> > > > ###
> > > >
> > > > Secondly, I created a background grid file, and I used it to display
> > > station
> > > > values, but I got a error when used command " oacres", the following
> > > steps
> > > > is what I did.
> > > >
> > > > ###
> > > > ga-> open grid.ctl
> > > > Scanning description file:  grid.ctl
> > > > Data file grid.dat is open as file 1
> > > > LON set to 99 106.45
> > > > LAT set to 19 20.49
> > > > LEV set to 1000 1000
> > > > Time values set: 2008:5:24:0 2008:5:24:0
> > > > ga-> open station.ctl
> > > > Scanning description file:  station.ctl
> > > > Data file st_data.dat is open as file 2
> > > > ga-> define qq=oacres(r,rain.2)
> > > > Warning from OACRES:  Less than two stations
> > > >     Grid will be all missing values
> > > > Define memory allocation size = 90000 bytes
> > > > ga-> d qq
> > > > Cannot contour grid - all undefined values
> > > > ###
> > > > according to the error message:"Warning from OACRES:  Less than two
> > > stations"
> > > > , and it occupied so much memory  , it seems there are some mistakes
> in
> > > the
> > > > first step, am I right? if yes, could you tell me what should I do ?
> and
> > > is
> > > > there any problem in my  steps?
> > > >
> > > > Thank you.
> > > >
> > > > Sincerely.
> > > >
> > > > --
> > > >
> > > > Shenghui                Student
> > > >
> > > > HuNan University .  China
> > > ------- End of Original Message -------
> > >
> >
> > --
> >
> > Shenghui                 Student
> > Major:                    Physics
> > HuNan University         China
> ------- End of Original Message -------
>
>

--

Shenghui                 Student
Major:                    Physics
HuNan University         China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090602/9a59d710/attachment.html 