Re: 20040309: range rings

Stonie,

When displaying radar or satellite data, you must use clear=yes
to display the image (eg in GPMAP), then in GPANOT you
will set CLEAR=no and overlay the rings on the image
produced in GPMAP.

As I also mentioned, if you want to expedite this process,
you can set the DISPLAY in GPANOT to a vgf file (remembering to
run gpend when finished!). Then, you can use the vgf file in GPMAP
in the future so that you only need to run 1 program instead of 2.

Steve Chiswell
Unidata User Support




On Tue, 2004-03-09 at 14:23, Stonie Cooper wrote:
> O.K. - I don't mean to be a party crasher, but I have been following this 
> thread with interest.  I can't get range rings to display per directions, so 
> let's look at the test:
> 
> Here is the gpmap:
> # gpmap
> MAP      = 1//1 + 2//4 + 26/// + 5///
>  GAREA    = 32.0;-104.2;35.8;-98.7
>  PROJ     = RAD
>  SATFIL   =
>  RADFIL   = /data/ldm/gempak/nexrad/NIDS/LBB/N0R/N0R_20040309_2047
>  LATLON   = 0
>  PANEL    = 0
>  TITLE    = 1/-3
>  TEXT     = .8/2/1/hw
>  CLEAR    = no
>  DEVICE   = xw
>  LUTFIL   = DEFAULT
>  STNPLT   = 17/.7|1|nexrad.tbl
>  VGFILE   =
>  AFOSFL   =
>  AWPSFL   =
>  LINE     = 3
>  WATCH    =
>  WARN     =
>  HRCN     =
>  ISIG     =
>  LTNG     =
>  ATCF     =
>  AIRM     =
>  NCON     =
>  CSIG     =
>  SVRL     =
>  BND      =
>  ATCO     =
>  TCMG     =
>  QSCT     =
>  WSTM     =
>  WOU      =
>  WCN      =
>  IMBAR    = 1
> r
> 
> exit
> Note that I have set clear=no, as I assume the gpanot will be an overlay.
> 
> Then:
> # gpanot
> GDFILE   =
>  SATFIL   =
>  RADFIL   = /data/ldm/gempak/nexrad/NIDS/LBB/N0R/N0R_20040309_2047
>  PROJ     = RAD
>  GAREA    = 32.0;-104.2;35.8;-98.7
>  PANEL    = 0
>  DEVICE   = xw
>  CLEAR    = no
>  SHAPE    = regpoly
>  INFO     = 1/1.0/360/360
>  LOCI     = #39.78;-104.53
>  LINE     = 3
>  CTYPE    = C
> r
> 
> And I get nothing.  What am I missing?  I'm using xw instead of gif so I can 
> actually watch the process . . . 
> 
> 
> 
> 
> On Tuesday 09 March 2004 21:00, Steve Chiswell wrote:
> > Bruce,
> >
> > As in my message, the paremeters LOCI, SHAPE and INFO
> >
> > are fore the GPANOT program:
> > > >then use use GPANOT
> > > > to draw navigated rings, eg with your PROJ=rad and RADFIL specifying
> > > > your radar file:
> > > >
> > > >  SHAPE    = regpoly
> > > >  INFO     = 1/1.0/360/360
> > > >  LOCI     = #39.78;-104.53
> >
> > Steve Chiswell
> >
> > On Tue, 2004-03-09 at 13:25, Haynie, Bruce wrote:
> > > Oops!  Thanks for pointing that out.  Made the correction, but still get
> > > the same result.
> > >
> > > Bruce
> > >
> > > -----Original Message-----
> > > From: Jeff Weber [mailto:address@hidden]
> > > Sent: Tuesday, March 09, 2004 20:18
> > > To: Haynie, Bruce
> > > Cc: address@hidden
> > > Subject: RE: 20040309: range rings
> > >
> > > Bruce,
> > >
> > > Not sure if this is an "exact" dup of your script, but if so:
> > >
> > > LOCI        = #33.65;-10180
> > >
> > > seems to be missing a decimal point...:)
> > >
> > > Jeff
> > > ---------------------------------------------------------------------
> > > Jeff Weber                                    address@hidden        :
> > > Unidata Program Center                        PH:303-497-8676        :
> > > University Corp for Atmospheric Research      3300 Mitchell Ln       :
> > > http://www.unidata.ucar.edu/staff/jweber      Boulder,Co 80307-3000  :
> > > ---------------------------------------------------------------------
> > >
> > > On Tue, 9 Mar 2004, Haynie, Bruce wrote:
> > > > Guess I'm missing something here.  Part of my script follows:
> > > >
> > > > gpmap << EOF
> > > >  MAP    = 1//1 + 2//4 + 26/// + 5///
> > > >  GAREA  = 32.0;-104.2;35.8;-98.7
> > > >  PROJ   = RAD
> > > >  SATFIL =
> > > >  RADFIL = $RADDIR/$FILE
> > > >   SHAPE = regpoly
> > > >   INFO  = 1/1.0/360/360
> > > >   LOCI  = #33.65;-10180
> > > >  LATLON = 0
> > > >  PANEL  = 0
> > > >  TITLE  = 1/-3
> > > >  TEXT   = .8/2/1/hw
> > > >  CLEAR  = y
> > > >  DEVICE = gif|lbb_rr_${IMGT}.gif|800;600
> > > >  LUTFIL = DEFAULT
> > > >  IMBAR  = 1
> > > >  STNPLT = 17/.7|1|$STATIONS
> > > >  VGFILE =
> > > >  AWPSFL =
> > > >  LINE   = 3
> > > > \$MAPFIL = hicnus.nws + hipowo.gsf + louhus.nws + loisus.nws r
> > > >
> > > > e
> > > > EOF
> > > >
> > > > gpend
> > > >
> > > > Tried running the above script, but no range rings were added.
> > > >
> > > > Bruce
> > > >
> > > > -----Original Message-----
> > > > From: Steve Chiswell [mailto:address@hidden]
> > > > Sent: Tuesday, March 09, 2004 19:14
> > > > To: Haynie, Bruce
> > > > Cc: address@hidden; GEMPAK support
> > > > Subject: 20040309: range rings
> > > >
> > > > Bruce,
> > > >
> > > > One possibility is to plot your radar with GPMAP,
> > > >
> > > >
> > > > The LOCI point would be the location of the radar. The radius of "1"
> > > > in the INFO line will draw a ring at 1.0 degrees radius from the LOCI
> > > > point (KFTG's lat/lon in this example).
> > > >
> > > > To use repeatedly in your script, you could create a VGF file with the
> > > >
> > > > rings drawn at the desired radii, and use that in the GPMAP parameter
> > > > VGFIL.
> > > >
> > > > Steve Chiswell
> > > > Unidata User Support
> > > >
> > > > On Tue, 2004-03-09 at 11:27, Haynie, Bruce wrote:
> > > > > Howdy!
> > > > >
> > > > > I'd like to display range rings on my GEMPAK-generated radar images.
> > > > > A quick search through the online literature yielded unfavorable
> > > > > results.  Is there a handy method that will allow me to display
> > > > > range rings using csh scripts??
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Bruce
> > > > > ------------------
> > > > > Bruce Haynie                    eml: address@hidden
> > > > > Research Associate              phn: 806-885-4644 x224
> > > > > West TX MesoNet         fax: 806-885-4468
> > > > > http://www.mesonet.ttu.edu