#!/bin/sh

#-------------------------------------------------------------

  WEB=""
  if [ -f /usr/bin/firefox ];then WEB="/usr/bin";fi
  if [ -f /usr/local/bin/firefox ];then WEB="/usr/local/bin";fi
  MDL=${1:-"${HOME}/hysplit"}
  OUT=${2:-"${MDL}/working"}
  DSP=${3:-"YES"}
  TTR=${TTR:-"${HOME}/Tutorial"}
  MET="${TTR}/captex"
  cd $OUT
  echo "### ${0##*/} ###"

#-------------------------------------------------------------

  rm -f INFILE
  for run in 001 002 003 004 005; do       
      if [ ! -f hysplit2.$run ];then
         echo "Missing HYSPLIT file: hysplit2.$run"
         echo "Rerun previous section: ens_data.sh"
         exit
      else
         echo "hysplit2.$run" >>INFILE
      fi
  done

#----------------------------------------------------------

  outf="temp1"
  let num=0
  cat INFILE | while read fname; do
  let num=$num+1
  if [ $num == 1 ];then 
     inp1=$fname
  else
     inp2=$fname
     ${MDL}/exec/concadd -i$inp1 -b$inp2 -o$outf -p0
     if [ "$outf" == "temp1" ];then
        inp1="temp1"
        outf="temp2"
     else
        inp1="temp2"
        outf="temp1"
     fi
  fi
  done

  ${MDL}/exec/conlight -i$outf -ohysplit2.006 -m0.20
  rm -f INFILE temp1 temp2

#----------------------------------------------------------

  rm -f sumstat.txt
  for run in 001 002 003 004 005 006; do       
      ${MDL}/exec/c2datem -ihysplit2.$run -odatem.$run -c1.0E+12 -m${MET}/captex2_meas.txt -xi -s
      ${MDL}/exec/statmain -t0 -rdatem.$run -d${MET}/captex2_meas.txt -l10.0 -o1 -ssumstat.txt
  done
  cat sumstat.txt