Dear, M. Mahakur,<br><br>Thank you for your time and replay.<br>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&#39; location are out of &quot; background grid&quot;. <br>
The radii of influence should be set between &quot;0&quot; and &quot;50&quot;, my Grads version is 1.8sl11, maybe in other versions are different values.<br>Now, I had set background grid data like this:<br>###<br>xdef  350 linear  104.61  0.0450<br>
ydef  250 linear   21.98  0.0450<br>###<br>then the display seems much better now, but there are so many small lines appeared, why does those small lines appear? <br>In this web page &quot; <a href="http://www.iges.org/grads/gadoc/users.html">http://www.iges.org/grads/gadoc/users.html</a>&quot; ( my manual ), I haven&#39;t found the reason, If you know, could you tell me?<br>
By the way, in my first step, when display station data, the background is the whole word, is that default display?<br><br>Thank you for your time.<br><br>Sincerely.<br><br><div class="gmail_quote">2009/6/1 M.Mahakur <span dir="ltr">&lt;<a href="mailto:mmahakur@tropmet.res.in">mmahakur@tropmet.res.in</a>&gt;</span><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hello,<br>
<br>
How many stations do you have to inerpolate to such a hi resolution map (0.05 x 0.01)?<br>
Try with several scan radii, starting from a bigger for example (I guess there is limit<br>
to nos. of radii and minimum value of radii, see if 0.01 is allowed, although I have<br>
written!!):<br>
<br>
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)<br>
<br>
Default radii are 10,7,4,2,1<br>
<br>
<a href="http://www.iges.org/grads/gadoc/gradfuncoacres.html" target="_blank">http://www.iges.org/grads/gadoc/gradfuncoacres.html</a><br>
<br>
--<br>
M. Mahakur<br>
<div><div></div><div class="h5"><br>
---------- Original Message -----------<br>
From: sh zhou &lt;<a href="mailto:sh.maeam@GMAIL.COM">sh.maeam@GMAIL.COM</a>&gt;<br>
To: <a href="mailto:GRADSUSR@LIST.CINECA.IT">GRADSUSR@LIST.CINECA.IT</a><br>
Sent: Mon, 1 Jun 2009 20:49:19 +0800<br>
Subject: station map<br>
<br>
&gt; Dear All,<br>
&gt;<br>
&gt; I want to draw a station map,  firstly, I created a binary file, and wrote a<br>
&gt; &quot;ctl&quot; file as following:<br>
&gt; ###<br>
&gt; dset st_data.dat<br>
&gt; title station data<br>
&gt; dtype station<br>
&gt; stnmap station_data.map<br>
&gt; undef -999.9<br>
&gt; tdef 1 linear 24May2008 7dy<br>
&gt; vars 1<br>
&gt; rain 0 99 rain data<br>
&gt; endvars<br>
&gt; ###<br>
&gt; but the picture&#39;s background is the whole world( there is a picture in<br>
&gt; attachment), and stations values displayed in a small area,  when I added<br>
&gt; &quot;xdef,yedf,zdef&quot; like this:<br>
&gt;<br>
&gt; ###<br>
&gt; ...<br>
&gt; undef -999.9<br>
&gt; xdef 100 linear 100.0 0.05<br>
&gt; ydef 100 linear 20.0 0.01<br>
&gt; zdef 1 linear 1000.0<br>
&gt; tdef 1 linear 24May2008 7dy<br>
&gt; ...<br>
&gt; ###<br>
&gt; however, the background didn&#39;t change, Does &quot;xdef&quot; and &quot;ydef&quot;&#39;s range  set<br>
&gt; by default as following?<br>
&gt; ###<br>
&gt; ga-&gt; open station.ctl<br>
&gt; Scanning description file:  station.ctl<br>
&gt; Data file st_data.dat is open as file 1<br>
&gt; LON set to 0 360<br>
&gt; LAT set to -90 90<br>
&gt; LEV set to 500 500<br>
&gt; Time values set: 2008:5:24:0 2008:5:24:0<br>
&gt; ga-&gt; q file<br>
&gt; ###<br>
&gt;<br>
&gt; Secondly, I created a background grid file, and I used it to display station<br>
&gt; values, but I got a error when used command &quot; oacres&quot;, the following steps<br>
&gt; is what I did.<br>
&gt;<br>
&gt; ###<br>
&gt; ga-&gt; open grid.ctl<br>
&gt; Scanning description file:  grid.ctl<br>
&gt; Data file grid.dat is open as file 1<br>
&gt; LON set to 99 106.45<br>
&gt; LAT set to 19 20.49<br>
&gt; LEV set to 1000 1000<br>
&gt; Time values set: 2008:5:24:0 2008:5:24:0<br>
&gt; ga-&gt; open station.ctl<br>
&gt; Scanning description file:  station.ctl<br>
&gt; Data file st_data.dat is open as file 2<br>
&gt; ga-&gt; define qq=oacres(r,rain.2)<br>
&gt; Warning from OACRES:  Less than two stations<br>
&gt;     Grid will be all missing values<br>
&gt; Define memory allocation size = 90000 bytes<br>
&gt; ga-&gt; d qq<br>
&gt; Cannot contour grid - all undefined values<br>
&gt; ###<br>
&gt; according to the error message:&quot;Warning from OACRES:  Less than two stations&quot;<br>
&gt; , and it occupied so much memory  , it seems there are some mistakes in the<br>
&gt; first step, am I right? if yes, could you tell me what should I do ? and is<br>
&gt; there any problem in my  steps?<br>
&gt;<br>
&gt; Thank you.<br>
&gt;<br>
&gt; Sincerely.<br>
&gt;<br>
&gt; --<br>
&gt;<br>
&gt; Shenghui                Student<br>
&gt;<br>
&gt; HuNan University .  China<br>
</div></div>------- End of Original Message -------<br>
</blockquote></div><br><br clear="all"><br>-- <br><br>Shenghui                 Student<br>Major:                    Physics<br>HuNan University         China<br>