#!/bin/csh -f
#
set YMD=`date +%Y%m%d`
limit stacksize 150000
source /usr/nawips/GEMPAK/Gemenviron

###############################################
# This program retrieves the 2-meter, 850mb, ##
# and the 1000-970mb layer temperatures in   ##
# Kelvin from the ruc 211 model.             ##
# from the ruc211 model.                     ##
###############################################
cd $NAWIPS/../radar_testing/usptype
set RUC=`ls -1t $HDS/ruc/*ruc211.gem | head -1`
gdinfo << EOF > out_pdlytmpk.txt
restore info_pdly.nts
gdfile = $RUC
r

e
EOF
gdinfo << EOF > out_2mtmpk.txt
restore info_2m.nts
gdfile = $RUC
r

e
EOF
gdinfo << EOF > out_850tmpk.txt
restore info_850.nts
gdfile = $RUC
r

e
EOF
set M2=`grep -c TMPK out_2mtmpk.txt`
set PDLY=`grep -c TMPK out_pdlytmpk.txt`
set PRES=`grep -c TMPK out_850tmpk.txt`
set answer=`echo "${M2}_${PDLY}_${PRES}"`
if ($answer == "1_1_1") then
cp $RUC ruc.gem
else
echo NO
echo $answer
endif

###########################################
# This program grids all CONUS radars to ##
# create the US and regional composites. ##
###########################################
cd $NAWIPS/../radar_testing
if (-e $NAWIPS/../radar_testing/grid/rradus.gem) rm $NAWIPS/../radar_testing/grid/rradus.gem
## set radarType to either n0r or n0q
set radarType="$1"
echo "radarType: $1"
gdradr << EOF
restore ${radarType}_gdradr.nts
l
r

e
EOF

cp $NAWIPS/../radar_testing/grid/${YMD}_radr.gem $NAWIPS/../radar_testing/grid/rradus.gem

gddiag << EOFb
restore gddiag_${radarType}.nts
gdoutf=grid/rradus.gem
gdfile=grid/rradus.gem
l
r

e
EOFb
gpend

## Create clutter removed precip grid
set RUC=`ls -1t $NAWIPS/../radar_testing/usptype/ruc.gem | head -1`
echo $RUC
set RADC=`ls -1t $NAWIPS/../radar_testing/grid/*us*.gem | head -1`
echo $RADC
#
########################################################################
# create a new grid file of the same projection as the radar grid
# ----------------------------------------------------------------
echo ""
echo "Running GEMPAK prog gdcfil"
echo ""
$GEMEXE/gdcfil << EOF
gdoutf = usptype.gem
cpyfil = ${RADC}
maxgrd = 600
anlyss =
l
r

e
EOF