#!/bin/sh

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

  MDL="${HOME}/hysplit4"
  MAP="${MDL}/graphics/arlmap"
  MET="${HOME}/Tutorial/smoke"

  OUT="."
  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 ###"

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

  syr=07
  smo=04
  sda=16
  shr=00
     
  olat=31.15
  olon=-82.40
  lvl1=0.0
        
  run=96
  ztop=10000.0
  data=EDAS16to20.bin

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

  echo "$syr $smo $sda $shr    " >CONTROL
  echo "1                      ">>CONTROL
  echo "$olat $olon $lvl1      ">>CONTROL
  echo "$run                   ">>CONTROL
  echo "0                      ">>CONTROL
  echo "$ztop                  ">>CONTROL
  echo "1                      ">>CONTROL
  echo "$MET/                  ">>CONTROL
  echo "$data                  ">>CONTROL
  echo "1                      ">>CONTROL
  echo "PM25                   ">>CONTROL
  echo "1.0                    ">>CONTROL
  echo "72.0                   ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "1                      ">>CONTROL
  echo "0.0 0.0                ">>CONTROL
  echo "0.1 0.1                ">>CONTROL
  echo "20.0 20.0              ">>CONTROL
  echo "$OUT/                  ">>CONTROL
  echo "smoke.bin              ">>CONTROL
  echo "1                      ">>CONTROL
  echo "100                    ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 01 00               ">>CONTROL
  echo "1                      ">>CONTROL
  echo "0.0 0.0 0.0            ">>CONTROL
  echo "0.0 0.0 0.0 0.0 0.0    ">>CONTROL
  echo "0.0 0.0 0.0            ">>CONTROL
  echo "0.0                    ">>CONTROL
  echo "0.0                    ">>CONTROL

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

  echo "&SETUP "                        >SETUP.CFG
  echo "initd = 3, "                   >>SETUP.CFG
  echo "khmax = 24, "                  >>SETUP.CFG
  echo "numpar = 25000, "              >>SETUP.CFG
  echo "maxpar = 50000, "              >>SETUP.CFG
  echo "efile = 'EMITIMES',"           >>SETUP.CFG
  echo "/ "                            >>SETUP.CFG

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

  for H in 1 2 3 4; do

      if [ $H -eq 1 ];then heat="0.215E+11"; fi
      if [ $H -eq 2 ];then heat="0.430E+11"; fi
      if [ $H -eq 3 ];then heat="0.860E+11"; fi
      if [ $H -eq 4 ];then heat="1.720E+11"; fi

      echo "YYYY MM DD HH DURATION RECORDS"                                    >EMITIMES
      echo "YYYY MM DD HH MM DURATION LAT LON  HGT  RATE  AREA   HEAT"        >>EMITIMES
      echo "2007 04 16 00 0096 0001"                                          >>EMITIMES
      echo "2007 04 16 18 00 7200 31.150 -82.400 0. 0.24E+15 0.20E+07 $heat"  >>EMITIMES

      rm -f smoke.bin

      ${MDL}/exec/hycs_std  

      grep "mixd,rise" MESSAGE >rise_$H.txt

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

      echo "'TITLE&','### $H ###&'"  >LABELS.CFG
      ${MDL}/exec/concplot -ismoke.bin -j${MDL}/graphics/arlmap -c4 -k1 -z90 -uug -n40:40 -v500+200+100+50+20+10
      rm -f plot_$H.ps
      mv concplot.ps plot_$H.ps

      echo Mayo 30.2556 -81.4533 >samplers.txt
      ${MDL}/exec/con2stn -ismoke.bin -ocon2stn.txt -r1 -ssamplers.txt

      ${MDL}/exec/timeplot -icon2stn.txt -s../../Tutorial/smoke/mayo_pm25.txt
      rm -f time_$H.ps 
      mv timeplot.ps time_$H.ps

 done