#!/bin/ksh
#
# Get ALAMO in situ observations
#
# Set error trapping on
#
set -euvx
#
# Document command line options
#
usage="Usage:\n$0 -o <output directory> -i <input directory> -r <execution directory> -b <initial date in yyyymm format> -e <final date in yyyymm format> -D <dir on disk> -B <bin dir> -S <data set number> -I <start day of input data> -E <end day of input data>"
#
# Set defaults for command line options
#
outdir=unknown
rundir=unknown
inpdir=unknown
dtadir=unknown
inidate=unknown
enddate=unknown
bindir=unknown
dataset=1
typeset -Z2 dataset
errflg=0
#
# Get command line options
#
while getopts b:e:f:i:o:r:D:B:S: option
do
  case $option in
    b) inidate=$OPTARG;;
    e) enddate=$OPTARG;;
    i) inpdir=$OPTARG;;
    o) outdir=$OPTARG;;
    r) rundir=$OPTARG;;
    B) bindir=$OPTARG;;
    D) dtadir=$OPTARG;;
    S) dataset=$OPTARG;;
    \?) errflg=1;;
  esac
done
#
# Check for errors
#
if [[ $errflg -ne 0 || $rundir = unknown || $outdir = unknown || $inpdir = unknown || $inidate = unknown || $enddate = unknown ]] ; then
  print $usage
  exit 1
fi
#
# Check that directories exists.
#
CHECK_DIR $outdir
CHECK_DIR $rundir
#
cd $rundir
#
# Check if file exists locally. If not, fetch from ECFS
CP_ARCH $inpdir/fb_alamo.nc .

cat >namsel.in<<EOF
&namsel
  ninidate=$inidate
  nenddate=$enddate
/
EOF

$bindir/fbsel fb_alamo.nc $outdir/profiles_${inidate}_${enddate}_${dataset}.nc