[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: 20040309: range rings
- Subject: Re: 20040309: range rings
- Date: 09 Mar 2004 14:57:38 -0700
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