set +x
#########################################################
# Run Hysplit over predefined locations
#########################################################
set -x
cd $DATA
DATA_JOB=$DATA

for city in $1
do
   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_JOB/$loc
   mkdir $DATA
   cd $DATA

   cp ../tmp.predefined $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

   k=0
   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

   #############################################################
   # When Source ID starts with letter Z, force continuous release
   #############################################################
   C1=`echo $site | cut -c1-1`
   if [ "${C1}" = "Z" ]
   then
      export dur=$run
   fi

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

done