set +x
#########################################################
# Run Hysplit over predefined locations
#########################################################

cd $DATA
set -x

city=$1

loc=`echo $city | awk -F\; '{print $1}'`
lat=`echo $city | awk -F\; '{print $2}'`
lon=`echo $city | awk -F\; '{print $3}'`
base=`echo $city | awk -F\; '{print $4}'`
top=`echo $city | awk -F\; '{print $5}'`

if [ "$base" = "" ] ; then
  base=10
fi
if [ "$top" = "" ] ; then
  top=50
fi

export DATA=$DATA/$loc
mkdir $DATA
cp tmp.predefined $DATA
cd $DATA
ln -s ../NAMA NAMA
ln -s ../NAMF NAMF
# Run setup to initialize working directory and utility scripts
sh $utilscript/setup.sh

sed s/#LOC#/$loc/g tmp.predefined > tmp1
sed s/#LAT#/$lat/g tmp1 > tmp2
sed s/#LON#/$lon/g tmp2 > tmp1
sed s/#BASE#/$base/g tmp1 > tmp2
sed s/#TOP#/$top/g tmp2 > hysplit.ini

cat hysplit.ini | while read tmp
do
  let k=$k+1
  case $k in
    1)  export title=$tmp;;
    2)  export meteo=$tmp;;
    3)  export site=$tmp;;
    4)  export olat=$tmp;;
    5)  export olon=$tmp;;
    6)  export nfile=$tmp;;
    7)  export dir1=$tmp;;
    8)  export met1=$tmp;;
    9)  export dir2=$tmp;;
    10) export met2=$tmp;;
    11) export year=$tmp;;
    12) export month=$tmp;;
    13) export day=$tmp;;
    14) export hour=$tmp;;
    15) export run=$tmp;;
    16) export dur=$tmp;;
    17) export rate=$tmp;;
    18) export rbot=$tmp;;
    19) export rtop=$tmp;;
    20) export poll=$tmp;;
    21) export syr=$tmp;;
    22) export smo=$tmp;;
    23) export sda=$tmp;;
    24) export ctop=$tmp;;
  esac
done

export met1=NAMA
export met2=NAMF
export nfile=2

#############################################################
# Run Trajectory and Concentration Models (HYSPLIT)
#############################################################
$USHhysplit/hysplit_trajconc_1.sh