set -x

#  Kitzmiller/Samplatsky  Jul/Dec 2000   IBM AIX
#                         May 2001  tgsv19

#  Collects most recent set of surface obs in METAR form.
#  Places codes for precip obs in the 460x360 FMH-definition RCM grid.

#  Command-line argument: number of hours prior to current time
#  for which to seek METARs, such as 01 for 1 hour, 02 for 2 hours

#  Input file: /com/hourly/prod/hourly.{yyyymmdd}/rawmeter.{hh}
#  Output file: metarrcm.dat

#  Output codes in metarrcm.dat:  7 - no ob in the box
#  0 = ob, no precip, 1 - rain, 2 - snow

#     Locate current ASP raw METAR file.  Put appropriate 
#     date-time into control file

rm -f fline metarlist.ctl

dset="/com/hourly/prod/hourly.${PDY}/rawmetar.${cyc}"
cp $dset METARS
line="${PDY} ${cyc}55"
echo $line > fline
cat $PARMradar/radar_metarlist.ctl0 fline > metarlist.ctl
cp $PARMradar/radar_sfc_stns.txt sfc_stns.txt

echo "METAR dataset: " $dset
echo "Date-time: " $line

#   Interpolate precip status obs to 10-km radar grid

export pgm=radar_metarlist
. prep_step

$EXECradar/radar_metarlist
err=$? ; export err; err_chk