##############################################################################
##
## SCRIPT: sdmair.sh                NCO/PMB/SDM
##
## Syntax:  sdmair $1 $2 $3 $4
##          $1 = model hour (00 06 12 18)
##          $2 = Latitude
##          $3 = Longitude
##          $4 = Length (width of data box)
##
## CHANGE LOG:
##
##   11/2008   MALONE      Changed data retrieval directory to gpfs sdmedit
##                         area
##
##   04/2010   MALONE      Changed data retrieval directory to /nco/ops/sdm[2]
##
##############################################################################
set -xa
envir=prod

## ERROR CHECK.
## Exit if 4 valid arguments are not entered.
## $1 = model hour... 00 06 12 18
## $2 = Latitude
## $3 = Longitude
## $4 = Length
[ $# -ne 4 ] && { echo "$0: <hour> <lat> <lon> <len> "; exit 99; }

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

## RETRIEVE DATA FROM OPERATIONAL DIRECTORY.
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`

## SETUP AND RUN SDM_SDMAIR_NEW EXECUTABLE ON THE LISTAIR AND LISTACR FILES
## CREATED ABOVE.
cp /nw${envir}/exec/sdm_sdmair_new $WKDIR

HOUR=$1
LAT=$2
LON=$3
LEN=$4
IFILE=1
cat <<EOH > cardec1
${HOUR} ${LAT} ${LON} ${LEN} ${IFILE}
EOH

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

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_sdmair_new <cardec1 1>sdmairc.out 2>sdmairc.err

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

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_sdmair_new <cardec2 1>sdmacar.out 2>sdmacar.err

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

more $WKDIR/airout

## END OF SCRIPT