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