########################################################################
##
## SCRIPT: upaair.sh            NCO/PMB/SDM
##
## Syntax:  upaair $1 $2 $3
##          $1 = model hour (00 06 12 18)
##          $2 = station ID
##          $3 = radius of data box
##
## CHANGE LOG:
##
##   11/2008   MALONE      Changed data retrieval directory to gpfs
##                         sdmedit area
##   04/2010   MALONE      Changed data retrieval directory from gpfs
##                         sdmedit area to /nco/ops/sdm[2]
##
########################################################################
set -xa

envir=prod

## ERROR CHECK
## Exit if 3 valid arguments are not entered.
## $1 = model hour... 00 06 12 18
## $2 = aircraft ID... found in column 1 of ListairZZ/LISTACARSZZ
## $3 = radius of data box
[ $# -ne 3 ] && { echo "$0: <hour> <id> <len> "; exit 99; }

## DEFINE THE WORKING DIRECTORY.
WKDIR=/stmp/d$$
mkdir $WKDIR
cd $WKDIR

## GATHER AIRCRAFT DATA
ls -t /nco/ops/sdm/Listair$1 /nco/ops/sdm2/Listair$1 > fileair
ls -t /nco/ops/sdm/LISTACARS$1 /nco/ops/sdm2/LISTACARS$1 > fileacr
Listair=`head -1 fileair`
Listacr=`head -1 fileacr`
## echo " RUN $1 " >> $WKDIR/airout

## RUN UPAAIR PROGRAM TO PRESENT DATA
cp /nw${envir}/exec/sdm_upaair $WKDIR

HOUR=$1
ID=$2
LEN=$3
IFILE=1
cat <<EOH > cardec1
${HOUR} ${ID} ${LEN} ${IFILE}
EOH

IFILE=2
cat <<EOH > cardec2
${HOUR} ${ID} ${LEN} ${IFILE}
EOH

cat /nwprod/dictionaries/sonde.land.tbl > $WKDIR/dict

ln -sf "$WKDIR/dict"                             "fort.15"
ln -sf "$WKDIR/crun"                             "fort.17"
ln -sf "$Listair"                                "fort.20"
ln -sf "$WKDIR/header1"                          "fort.49"
ln -sf "$WKDIR/selairc"                          "fort.50"

$WKDIR/sdm_upaair <cardec1 1>upaairc.out 2>upaairc.err

cat          $WKDIR/header1 >> $WKDIR/airout
sort +4n -5  $WKDIR/selairc >> $WKDIR/airout

ln -sf "$WKDIR/dict"                             "fort.15"
ln -sf "$WKDIR/crun"                             "fort.17"
ln -sf "$Listacr"                                "fort.20"
ln -sf "$WKDIR/header2"                          "fort.49"
ln -sf "$WKDIR/selacar"                          "fort.50"

$WKDIR/sdm_upaair <cardec2 1>upaacar.out 2>upaacar.err

cat          $WKDIR/header2 >> $WKDIR/airout
sort +4n -5  $WKDIR/selacar >> $WKDIR/airout

more $WKDIR/airout

## END OF SCRIPT