[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: 20011115: NEXRAD reflectivity mosaic availability



David,

The land/water coloring is done in gpmap using the BND parameter.
Attacthed is the script (NEXRIII.csh) which uses gdplot2 and gpmap.

The individual NIDS station plot is created from a cgi I created that
calls: radmap_sw -s input_nids_file output_gif_file "${TITLE}"

The radmap_sw program must have the map and station tables
in the working directory (eg it doesn't search for them).
If you need more information on radmap, see the
$GEMPAK/source/contrib/tdl/radmap directory.

The files that must exist in the radmap working directory are:
drwxrwxrwx 9219/2000      0 Sep 18 10:57 2001 .radmap_work/
-rw-r--r-- 9219/2000  82485 Mar 20 13:58 2001 .radmap_work/cities.rad
-rw-r--r-- 9219/2000  35353 Mar 20 13:58 2001 .radmap_work/geog.tbl
-rw-r--r-- 9219/2000 1756160 Mar 20 13:58 2001 .radmap_work/hicnus.nws
-rw-r--r-- 9219/2000 1600512 Mar 20 13:58 2001 .radmap_work/hipona.nws
-rw-r--r-- 9219/2000 750592 Mar 20 13:58 2001 .radmap_work/hiponw.cia
-rw-r--r-- 9219/2000 289792 Mar 20 13:58 2001 .radmap_work/lakes.ncp
-rw-r--r-- 9219/2000  64512 Mar 20 13:58 2001 .radmap_work/loconh.gsf
-rw-r--r-- 9219/2000 221184 Mar 20 13:58 2001 .radmap_work/mepowo.gsf
-rw-r--r-- 9219/2000   6831 Mar 20 13:58 2001 .radmap_work/nidprd.tbl

You can find those files wither in the directory mentioned above,
or the standard $GEMTBL directories.

The tdl program is just a hardcoded version of gpmap which they created to
eliminate some of the time it takes to search for the files. Just makes
getting the output to the web faster.

If you need more info, I'll be back after the 28th. Kind of in a hurry now.

Steve Chiswell
Unidata User SUpport


On Sat, 17 Nov 2001, David Ovens wrote:

> Steve,
> 
> These radar mosaics are fantastic.  Is there any chance you could send
> me a copy of the scripts that generate the basic US mosaic image and
> any of the individual station plots.  I see the c-shell script for
> generating the gdradr file, and I've been able to generate my own, but
> now I am interested in the land/water mask on the former and the
> station names, counties, color table, and color-bar labels in the
> latter.  The former I think I can do, but the latter looks like
> something other than gdplot2.  Is it nmap?  If so, how is this gif
> generated in batch mode?
> 
> Thanks,
> 
> David 
> -- 
> 
> David Ovens           e-mail: address@hidden
> (206) 685-8108          plan: Real-time MM5 forecasting for Pacific Northwest
> Research Meteorologist
> Dept of Atmospheric Sciences, Box 351640
> University of Washington 
> Seattle, WA  98195
> 
> 
> Steve Chiswell wrote:
> > 
> > 
> > As a service to the community, 2 NEXRAD reflectivity mosaic products
> > are available on the FNEXRAD data feed. One is a 6km national composite.
> > The other is a 1km regional floater. These are the products used to
> > create the example images shown at:
> > http://motherlode.ucar.edu/unidata/images/nids/gdradr.html
> > 
> > Millersville University has accepted the job for selecting the floater 
> > region
> > which is being created for the IDD.
> > 
> > The radar composites on FNEXRAD are distributed in standard GRIB format
> > so that users can access the data locally using a variety of data
> > display tools.
> > 
> > The composites are created using the GEMPAK program GDRADR which is 
> > described
> > in more detail at:
> > http://www.unidata.ucar.edu/packages/gempak/examples/gdradr/
> > Examples of other uses of the products are also shown on the above web page.
> > 
> > Any site receiving NEXRAD data can create their own mosaics locally
> > using GDRADR. The products on the FNEXRAD feed provide an alternative
> > for sites that do not want to receive all the NEXRAD radar data,
> > or may only want to create their own local region composites.
> > 
> > LDM pqact.conf examples, and further information on decoding and displaying
> > the products can be found here:
> > http://www.unidata.ucar.edu/packages/gempak/examples/gdradr/index.html#DECODING
> > 
> > 
> > Steve Chiswell
> > Unidata User Support
> > 
> > 
> 
> 
> 
#!/bin/csh -f

source /home/gempak/NAWIPS/Gemenviron
set WEBDIR=/opt/apache/htdocs/motherlode/unidata/images/nids

cd /tmp
set WORK=.nexrwork.$$
if(! -e $WORK) mkdir $WORK

cd $WORK


gpmap << EOF_BND
 MAP      = 32
 GAREA    = 22.0;-122.0;47.0;-65.0
 PROJ     = lcc/25;-103;60/3;3;0;0
 \$mapfil = hipowo.cia
 SATFIL   =  
 RADFIL   =  
 LATLON   = 0
 PANEL    = 0
 TITLE    = 0
 TEXT     = 1.3/23/1/1/s/l/hw
 CLEAR    = y
 DEVICE   = gif|radar_mosaic.gif|900;700
 LUTFIL   =  
 STNPLT   =  
 VGFILE   =  
 AFOSFL   =  
 AWPSFL   =  
 LINE     =  
 WATCH    =  
 WARN     =  
 HRCN     =  
 ISIG     =  
 LTNG     =  
 ATCF     =  
 AIRM     =  
 NCON     =  
 SVRL     =  
 BND      = bg/24+world_bnds/8+lakes/24
 ATCO     =
 r

 e
EOF_BND


gdplot2 << EOF
device   = gif|radar_mosaic.gif|900;700
gdfile   = nexr
\$mapfil = hipowo.cia
proj     = lcc/25;-103;60/3;3;0;0
garea    = 22.0;-122.0;47.0;-65.0
map      = 0
clear    = n
gdattim  = LAST
text     = 1/1/1
TITLE    = 1/-1/~ Unidata 4km composite
GDPFUN   = n0r
GLEVEL   = 0
GVCORD   = none
TYPE     = F
!FINT     = 5;10;15;20;25;30;35;40;45;50;55;60;65;70;75
!FLINE    = 0;26-24;21-23;20;18;17;15-7
FINT     = 20;25;30;35;40;45;50;55;60;65;70;75
FLINE    = 0;21-23;20;18;17;15-7
!CLRBAR   = 1|1/1/1
CLRBAR   = 1/V/LL/0.001;0.15/0.82;0.010/1|1/1/1
PANEL    = 0
SKIP     = 0
SCALE    = 0
CONTUR   = 0
HILO     =  
LATLON   = 0
STNPLT   =  
SATFIL   =  
RADFIL   = 
LUTFIL   =  
STREAM   =  
POSN     = 0
COLORS   = 1
MARKER   = 0
GRDLBL   = 0
FILTER   = YES
r

e
EOF

gpend

mv radar_mosaic.gif $WEBDIR
rm *.nts 
cd ..
rmdir $WORK

cd $WEBDIR
@ COUNT = 0
if(-e mosaicloop${COUNT}) rm -f mosaicloop${COUNT}
while ($COUNT < 8)
   @ NEXT = $COUNT + 1
   mv mosaicloop${NEXT} mosaicloop${COUNT}
   @ COUNT = $NEXT
end
cp radar_mosaic.gif mosaicloop8