[gradsusr] [BULK] Re: Masking Domain in Grads

Mon Mar 7 22:04:18 EST 2011

On Mon, Mar 7, 2011 at 9:28 PM, Zablone Owiti <zablonowiti at yahoo.com> wrote:

> Hi Yaqiang,
> Thanks for the response, but the masking (Blanking) we meant here cannot be
> handled by the function you  suggested.
> 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.
> The question is how to generate a blanking file in grads for a given
> boundary (polygon).
>
>

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.

   Arlindo

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

-- 
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110307/a17c0258/attachment-0003.html 