#!/bin/sh

MDL="${HOME}/hysplit4"
MET="${HOME}/Tutorial/captex"
OUT="${MDL}/working"
cd $OUT

if [ ! -f ASCDATA.CFG ]; then
   echo "-90.0  -180.0"     >ASCDATA.CFG
   echo "1.0     1.0"      >>ASCDATA.CFG
   echo "180     360"      >>ASCDATA.CFG
   echo "2"                >>ASCDATA.CFG
   echo "0.2"              >>ASCDATA.CFG
   echo "'$MDL/bdyfiles/'" >>ASCDATA.CFG
fi
echo "### $0 ###"

echo "83 09 25 17         "  >CONTROL
echo "1                   " >>CONTROL
echo "39.90 -84.22 750.0  " >>CONTROL
echo "68                  " >>CONTROL
echo "0                   " >>CONTROL
echo "10000.0             " >>CONTROL
echo "1                   " >>CONTROL
echo "$MET/               " >>CONTROL
echo "captex2_narr.bin    " >>CONTROL
echo "./                  " >>CONTROL
echo "tdump               " >>CONTROL

rm -f SETUP.CFG
rm -f tdump*       
rm -f CONTROL.???

${MDL}/exec/dat2cntl -i${MET}/captex2_meas.txt

let run=1
while [ $run -lt 1000 ];do
   RUN=`printf %3.3d $run`
   if [ -f CONTROL.$RUN ];then ${MDL}/exec/hyts_std $RUN; fi
   rm -f CONTROL.$RUN MESSAGE.$RUN
   let run=$run+1
done

ls tdump*  >INFILE
${MDL}/exec/trajfreq -ftfreq.bin -g0.5 -iINFILE -r0 -s0:99999 

echo "'TITLE&','### $0 ### Traj frequency&'"  >LABELS.CFG
echo "'MAPID&',' Values &'"             >>LABELS.CFG
echo "'UNITS&',' %&'"                   >>LABELS.CFG
echo "'VOLUM&',' &'"                    >>LABELS.CFG

${MDL}/exec/concplot -m0 -k1 -z80 -j${MDL}/graphics/arlmap -itfreq.bin
gs concplot.ps

# rm -f LABELS.CFG INFILE
# rm -f tfreq.bin tdump?????_????????
# rm -f WARNING.???