#!/bin/csh -f
if(x${1} != x) then
gpmap << SELECT1
garea = $1
e
SELECT1
else
gpmap << SELECT2
garea = usnps
proj = nps
map = 6
title = 5/-1/Select Map Area
dev = xw|SELECT_AREA
clear = y
\$mapfil = hipowo.gsf
r
cursor garea
e
SELECT2
gpend
endif
if(x${2} != x) then
set DEVICE=$2
else
echo -n "Select an output device (xw,nc,gf) : "
set DEVICE=$<
endif
set EXT=""
if($DEVICE == nc) then
set EXT=".nc"
endif
if($DEVICE == gf) then
set EXT=".gif"
endif
# This is the directory of our lightning data
set NLDN=$GEMDATA/nldn
# set the current year, month day
set YMD=`date -u '+%Y%m%d'`
# get the day of week, month abbrev., day and 4 digit year for title
set TITLE = `date -u '+%A %B %d, %Y'`
# we'll start at hour 0 and loop through hour 23
@ HOUR = 0
@ STOP_HOUR = 23
@ COLOR=30
while($HOUR <= $STOP_HOUR)
if($HOUR < 10) then
set SFHOUR=0${HOUR}
else
set SFHOUR=${HOUR}
endif
set FILES=`ls $NLDN/${YMD}${SFHOUR}*_nldn.gem`
if($#FILES > 0) then
foreach FILE ($FILES)
set CLEAR = yes
echo "plotting $FILE"
set TIMES=`sfctime $FILE`
foreach TIME ($TIMES)
sfmap << EOF
map = 1
dev = ${DEVICE}|lightning_${YMD}${EXT}
proj = nps
area = dset
filter = no
clear = $CLEAR
text = 1.2/23/1/hw
title = 5/-1/Lightning Strikes $TITLE ${HOUR}Z
panel = .02;.02;.97;.97/1/1/4
sffile = $FILE
sfparm = sgnl;mark:15:4:1
color = 0;${COLOR}
dattim = $TIME
r
e
EOF
set CLEAR = no
end # foreach TIME
end # foreach file
endif
@ HOUR = $HOUR + 1
@ COLOR = $COLOR - 1
if($DEVICE == "gf") gpend
end
exit(0)