<br><br><div class="gmail_quote">On Mon, Mar 7, 2011 at 9:28 PM, Zablone Owiti <span dir="ltr"><<a href="mailto:zablonowiti@yahoo.com">zablonowiti@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="color:#000;font-size:14pt"><div><span>Hi Yaqiang,</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:garamond,new york,times,serif;background-color:transparent;font-style:normal"><span>Thanks for the response, but the masking (Blanking) we meant here cannot be handled by the function you suggested.</span></div>
<div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:garamond,new york,times,serif;background-color:transparent;font-style:normal"><span>The objective is to have the contours or shading restricted within a country or state boundary and NOT to mask out certain values, land or ocean areas. <br>
</span></div><div style="color:rgb(0, 0, 0);font-size:18.6667px;font-family:garamond,new york,times,serif;background-color:transparent;font-style:normal"><span>The question is how to generate a
blanking file in grads for a given boundary (polygon).<br></span></div><div> </div></div></div></blockquote><div><br></div><div><br></div><div>Hmm. I've never done this but here how would go about it. Select a rectangular domain and draw contours as usual. Next, get a shapefile with polygons describing your states. Next, draw filled polygons (say, color white) for each state that you want to mask out --- that is, exclude your state of interest. This is not a very elegant solution, but might do the trick. Read the shapefile documentation carefully. YMMV.</div>
<div><br></div><div> Arlindo</div><div><br></div><div><br></div><div><br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="color:#000;font-size:14pt">
<div><div><br></div></div><div><br><blockquote style="border-left:2px solid rgb(16, 16, 255);margin-left:5px;padding-left:5px"><div style="font-family:garamond,new york,times,serif;font-size:14pt"><div style="font-family:times new roman,new york,times,serif;font-size:12pt">
<font face="Arial" size="2"><hr size="1"><b><span style="font-weight:bold">From:</span></b> Yaqiang Wang <<a href="mailto:yaqiang.wang@gmail.com" target="_blank">yaqiang.wang@gmail.com</a>><div class="im"><br><b><span style="font-weight:bold">To:</span></b> Zablone Owiti <<a href="mailto:zablonowiti@yahoo.com" target="_blank">zablonowiti@yahoo.com</a>>; GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
</div><b><span style="font-weight:bold">Sent:</span></b> Monday, 7 March 2011 9:55 PM<br><b><span style="font-weight:bold">Subject:</span></b> Re: [gradsusr] [BULK] Re: Masking Domain in Grads<br></font><div><div></div><div class="h5">
<br> maskout function could be used to solve your problem.<br><a href="http://www.iges.org/grads/gadoc/gradfuncmaskout.html" target="_blank">http://www.iges.org/grads/gadoc/gradfuncmaskout.html</a><br><br>On Mon, Mar 7, 2011 at 9:21 PM, Zablone Owiti <<a href="mailto:zablonowiti@yahoo.com" target="_blank">zablonowiti@yahoo.com</a>> wrote:<br>
> Hi Sujata,<br>><br>> I have not been able to mask-out or blank a particular boundary (country or<br>> state) so that the contour or shading does not spill out of the boundary<br>> country boundary.<br>> I had earlier posted this query to Grads Users forum but did not get a<br>
> response. I hope this time there will be a response from the Grads Gurus.<br>><br>> ===================================================<br>> ZABLONE OWITI<br>><br>> From: Sujata Pattanayak <<a href="mailto:sujata05@gmail.com" target="_blank">sujata05@gmail.com</a>><br>
> To: Zablone Owiti
<<a href="mailto:zablonowiti@yahoo.com" target="_blank">zablonowiti@yahoo.com</a>><br>> Cc:<br>> Sent: Monday, 7 March 2011 8:22 PM<br>> Subject: Re: [gradsusr] [BULK] Re: opening text file data<br>><br>
> Hi Zablone,<br>> I am attaching my model domain. I want to plot certain model variables<br>> in a particular state (i.e ORISSA) as mentioned in the plot. Can I plot on a<br>> particular state (like Orissa) and mask the other area.<br>
><br>> regards,<br>> Sujata<br>><br>> On Mon, Mar 7, 2011 at 5:14 PM, Sujata Pattanayak <<a href="mailto:sujata05@gmail.com" target="_blank">sujata05@gmail.com</a>><br>> wrote:<br>><br>> Hi Zablone,<br>
> Thank you very much for your kind help.<br>><br>> Regards,<br>> Sujata<br>><br>> On Mon, Mar 7, 2011 at 4:32 PM, Zablone Owiti <<a href="mailto:zablonowiti@yahoo.com" target="_blank">zablonowiti@yahoo.com</a>> wrote:<br>
><br>> Hi Sujata,<br>> You can use the shapefiles (admin98) in Grads to generate State boundaries<br>> after you display something.<br>> The command is:<br>>> draw shp admin98<br>> This command is available with Grads V. 2.0.a8 or later. You can get more<br>
> details here: <a href="http://opengrads.org/doc/udxt/shape/shape.html" target="_blank">http://opengrads.org/doc/udxt/shape/shape.html</a><br>><br>> Hope this will help you<br>><br>><br>> ===================================================<br>
> ZABLONE OWITI<br>> GRADUATE STUDENT<br>> Nanjing University of Information, Science and Technology<br>> College of Atmospheric Science<br>> Add: 219 Ning Liu Rd, Nanjing, Jiangsu, 21004, P.R. China<br>> Tel: +86-25-58731402<br>
> Fax: +86-25-58731456<br>> Mob.
15077895632<br>> Website: <a href="http://www.nuist.edu.cn" target="_blank">www.nuist.edu.cn</a><br>> ====================================================<br>> "The real art of conversation is not only to say the right thing at the<br>
> right time, but also to leave unsaid the wrong thing at the tempting<br>> moment."<br>><br>><br>> From: Sujata Pattanayak <<a href="mailto:sujata05@gmail.com" target="_blank">sujata05@gmail.com</a>><br>
> To: Zablone Owiti <<a href="mailto:zablonowiti@yahoo.com" target="_blank">zablonowiti@yahoo.com</a>>; GrADS Users Forum<br>> <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>
> Sent: Monday, 7 March 2011 1:13 PM<br>> Subject: Re: [gradsusr] [BULK] Re: opening text file data<br>><br>> Hi Zablone,<br>> How do you get the state boundaries with GrADS? When I am plotting<br>> through grads, it
just display the country boundary.<br>> Can you please help me in this regard??<br>><br>> Regards,<br>> Sujata<br>><br>> On Wed, Feb 9, 2011 at 10:39 AM, Zablone Owiti <<a href="mailto:zablonowiti@yahoo.com" target="_blank">zablonowiti@yahoo.com</a>><br>
> wrote:<br>><br>> Dear Ramaphane,<br>> As Muhammad has indicated, the script runs. BUT your Objective is to plot<br>> the station Locations on the map, with station names. If that is the<br>> case, you need to open a file (eg. open model.ctl) to set the environment<br>
> (so that you don't get the error massage of No scaling environment) and<br>> display a variable. The trick I normally use so that I get a blank map is to<br>> set cmin a value higher than the maximum value of the variable (there could<br>
> be a better way of doing this, I hope somebody will point it out). For<br>> example I set cmin 10 before diplaying
pr variable in model.ctl.<br>> You also need to zoom in your domain using set lat and set lon.<br>> See attached map(city.png) I have done with the script.<br>> Hope this will help.<br>><br>> ===================================================<br>
> ZABLONE OWITI<br>> GRADUATE STUDENT<br>> Nanjing University of Information, Science and Technology<br>> College of International Education<br>> Add: 219 Ning Liu Rd, Nanjing, Jiangsu, 21004, P.R. China<br>
> Tel: +86-25-58731402<br>> Fax: +86-25-58731456<br>> Mob. 15077895632<br>> Website: <a href="http://www.nuist.edu.cn" target="_blank">www.nuist.edu.cn</a><br>> ====================================================<br>
> "The real art of conversation is not only to say the right thing at the<br>> right time, but also to leave unsaid the wrong thing at the tempting<br>> moment."<br>><br>><br>> From: Muhammad Rahiz <<a href="mailto:muhammad.rahiz@ouce.ox.ac.uk" target="_blank">muhammad.rahiz@ouce.ox.ac.uk</a>><br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>> Cc:<br>> Sent: Tuesday, 8 February 2011 11:23 PM<br>> Subject: Re: [gradsusr] [BULK] Re: opening text file data<br>
><br>> I am able to run the script. See results below<br>><br>> lat is -18.38 lon is 21.85 and the city is Shakawe<br>> start is 1<br>> lat is -17.82 lon is 25.15 and the city is Kasane<br>> start is 2<br>
> lat is -19.98 lon is 23.42 and the city is Maun<br>> start is 3<br>> lat is -18.55 lon is 25.63 and the city is Pandamatenga<br>> start is 4<br>> lat is -20.52 lon is 26.10 and the city is Sua_Pan<br>> start is 5<br>
> lat is -21.42 lon is 25.62 and the city is Letlhakane<br>> start is 6<br>> lat is -21.17 lon is 27.52 and the city is Francistown<br>> start
is 7<br>> lat is -21.97 lon is 27.83 and the city is Selibe_Phikwe<br>> start is 8<br>> lat is -23.12 lon is 26.83 and the city is Mahalapye<br>> start is 9<br>> lat is -24.67 lon is 25.92 and the city is SSKA<br>
> start is 10<br>> lat is -24.58 lon is 24.80 and the city is Jwaneng<br>> start is 11<br>> lat is -25.47 lon is 25.43 and the city is Good_Hope<br>> start is 12<br>> lat is -26.22 lon is 22.40 and the city is Tsabong<br>
> start is 13<br>> lat is -24.02 lon is 21.88 and the city is Tshane<br>> start is 14<br>> lat is -25.27 lon is 23.25 and the city is Werda<br>> start is 15<br>> lat is -21.68 lon is 21.63 and the city is Gantsi<br>
> start is 16<br>> lat is -22.38 lon is 26.72 and the city is Serowe<br>> start is 17<br>> lat is -24.38 lon is 25.50 and the city is Molepolole<br>> start is 18<br>> lat is -25.75 lon is 28.20 and the city is Pretoria<br>
> start is 19<br>>
lat is -33.92 lon is 18.45 and the city is Cape_Town<br>> start is 20<br>> lat is -22.57 lon is 17.10 and the city is Windhoek<br>> start is 21<br>> lat is -17.83 lon is 31.05 and the city is Harare<br>> start is 22<br>
> lat is -29.32 lon is 27.48 and the city is Maseru<br>> start is 23<br>> lat is -26.33 lon is 31.13 and the city is Manzini<br>> start is 24<br>> lat is -25.97 lon is 32.58 and the city is Maputo<br>> start is 25<br>
> last city is 24<br>> No scaling environment<br>><br>> This means x is environment and y is<br>> DRAW error: Syntax is DRAW MARK marktype x y size<br>> Non-numeric args to numeric operation<br>> Error occurred on line 53<br>
> In file <a href="http://city_plot.gs" target="_blank">city_plot.gs</a><br>><br>><br>> --<br>> Muhammad Rahiz<br>> Researcher & DPhil Candidate (Climate Systems & Policy)<br>> School of Geography & the Environment<br>
> University of
Oxford<br>><br>> On Tue, 8 Feb 2011, Galebonwe Ramaphane wrote:<br>><br>>> Muhammad,<br>>><br>>> I'm also running GrADS on linux. I have gone over the code several times.<br>>> It's one of those annoying problems. I'll attach the sample script and the<br>
>> text file for your inspection.<br>>><br>>> Galebonwe Ramaphane<br>>> Department of Meteorological Services<br>>> P.O. BOX 10100<br>>> Gaborone<br>>> Botswana<br>>> Tel: +267 3612200<br>
>> Mobile: +267 74435261<br>>><br>>><br>>> -----Original Message-----<br>>> From: <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>]<br>
>> On Behalf Of Muhammad Rahiz<br>>> Sent:
Tuesday, February 08, 2011 10:36 AM<br>>> To: GrADS Users Forum<br>>> Subject: [BULK] Re: [gradsusr] opening text file data<br>>> Importance: Low<br>>><br>>> Galebonwe,<br>>><br>>> I'm using GrADS on linux and as far as I know, if your text file is in the<br>
>> same directory as your script, you don't have to put in the full path.<br>>> I'm not sure if it's different on Windows although I doubt so. I presume<br>>> you've checked your code?<br>
>><br>>> Muhammad Rahiz<br>>> Researcher & DPhil Candidate (Climate Systems & Policy) School of<br>>> Geography & the Environment University of Oxford<br>>><br>>> On Tue, 8 Feb 2011, Galebonwe Ramaphane wrote:<br>
>><br>>>> I have a similar but minor problem. My calling script does not see the<br>>>> text file. I tried putting the text file in the same directory as the<br>>>>
calling script. I still failed to read the text file. Do I include the full<br>>>> path of the text file?<br>>>><br>>>> Galebonwe Ramaphane<br>>>> Department of Meteorological Services<br>
>>> P.O. BOX 10100<br>>>> Gaborone<br>>>> Botswana<br>>>> Tel: +267 3612200<br>>>> Mobile: +267 74435261<br>>>><br>>>> -----Original Message-----<br>>>> From: <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a><br>
>>> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Muhammad Rahiz<br>>>> Sent: Monday, February 07, 2011 6:32 PM<br>>>> To: GrADS Users Forum<br>
>>> Subject: Re: [gradsusr] opening text file data<br>>>><br>>>> James,
you can try the following.<br>>>><br>>>> *****<br>>>><br>>>> 'sdfopen <a href="http://ifile.nc" target="_blank">ifile.nc</a>'<br>>>><br>>>> * ttot = number of lines in text file<br>
>>> ttot = 755<br>>>> tx = 1<br>>>> while(tx<=ttot)<br>>>> fname ='7hA01.out.txt'<br>>>> rec = read(fname)<br>>>> IO = sublin(rec,1)<br>>>> if(IO = 1)<br>
>>> say fname' > File open error'<br>>>> 'quit'<br>>>> endif<br>>>> if(IO = 0)<br>>>> t1 = subwrd(rec,2)<br>>>> t2 = subwrd(rec,3)<br>
>>> endif<br>>>> say ''t1' 't2''<br>>>><br>>>> tx = tx +1<br>>>> endwhile<br>>>><br>>>> *****<br>>>><br>>>> --<br>>>> Muhammad Rahiz<br>
>>> Researcher & DPhil Candidate
(Climate Systems & Policy) School of<br>>>> Geography & the Environment University of Oxford<br>>>><br>>>> On Mon, 7 Feb 2011, James Ciarlo` wrote:<br>>>><br>>>>> Hi,<br>
>>>> I have several text files such as the one I attached and I need to<br>>>>> open them which GrADS but I cannot understand how to do so. Can<br>>>>> someone help me please?<br>>>>><br>
>>>> James<br>>>>><br>>>>><br>>>> _______________________________________________<br>>>> gradsusr mailing list<br>>>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
>>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>>>><br>>>>
_______________________________________________<br>>>> gradsusr mailing list<br>>>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>>>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
>>><br>>><br>><br>><br>><br>> _______________________________________________<br>> gradsusr mailing list<br>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>><br>><br>><br>><br>><br>><br>><br>><br>><br>><br>><br>> _______________________________________________<br>
> gradsusr mailing list<br>> <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>><br><br><br><br>-- <br>*************************************************<br>Dr. Ya-Qiang Wang<br>Centre for Atmosphere Watch and Services (CAWAS), CMA<br>Chinese Academy of Meteorological Sciences (CAMS)<br>46, Zhong-Guan-Cun South Avenue<br>
Beijing, 100081<br>China<br><br><a href="mailto:yaqiang.wang@gmail.com" target="_blank">yaqiang.wang@gmail.com</a><br>**************************************************<br><br><br></div></div></div></div></blockquote></div>
</div><br></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>