#!/bin/sh

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

  MDL="${HOME}/hysplit4"
  MAP="${MDL}/graphics/arlmap"
  OUT="."
  MET="${HOME}/Tutorial/volcano"
  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 "'${MD}L/bdyfiles/'" >>ASCDATA.CFG
  fi
  echo "### $0 ###"

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

  syr=10
  smo=04
  sda=16
  shr=00 
          
  run=12
  ztop=32000.0
  data="apr1420.bin"

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

  echo "$syr $smo $sda $shr    " >CONTROL
  echo "2                      ">>CONTROL
  echo "63.63 -19.62  100.0    ">>CONTROL
  echo "63.63 -19.62 6000.0    ">>CONTROL
  echo "$run                   ">>CONTROL
  echo "0                      ">>CONTROL
  echo "$ztop                  ">>CONTROL
  echo "1                      ">>CONTROL
  echo "$MET/                  ">>CONTROL
  echo "$data                  ">>CONTROL
  echo "2                      ">>CONTROL
  echo "SO2                    ">>CONTROL
  echo "208333.3               ">>CONTROL
  echo "12.0                   ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "SO4                    ">>CONTROL
  echo "0.0                    ">>CONTROL
  echo "0.0                    ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "1                      ">>CONTROL
  echo "45.0 0.0               ">>CONTROL
  echo "0.25 0.25              ">>CONTROL
  echo "90.0 360.0             ">>CONTROL
  echo "$OUT/                  ">>CONTROL
  echo "volcano.bin            ">>CONTROL
  echo "1                      ">>CONTROL
  echo "10000                  ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 00 00 00 00         ">>CONTROL
  echo "00 12 00               ">>CONTROL
  echo "2                      ">>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 "1.0 1.0 1.0            ">>CONTROL
  echo "0.001 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 "ichem=2,            " >>SETUP.CFG
  echo "maxdim=2,           " >>SETUP.CFG
  echo "numpar=5000,        " >>SETUP.CFG
  echo "maxpar=50000,       " >>SETUP.CFG
  echo "/                   " >>SETUP.CFG

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

  for R in 05 10 20; do

     rm -f volcano.bin
     echo "1 2 0.$R 1.5"    >CHEMRATE.TXT
     ${MDL}/exec/hycs_std  

     echo "'TITLE&','### $0 ### SO2 $R %/h &'"  >LABELS.CFG  
     echo "'UNITS&','ppb&'"                    >>LABELS.CFG
     echo "'VOLUM&','&'"                       >>LABELS.CFG
     ${MDL}/exec/concplot -ivolcano.bin -oplot_so2_$R -j${MAP} -x0.453E+09 -s1 -z80 -k2 -c4 -v10+5+2+1+0.5+0.2+0.1

     echo "'TITLE&','### $0 ### SO4 $R %/h &'"  >LABELS.CFG
     echo "'UNITS&','ug&'"                     >>LABELS.CFG
     echo "'VOLUM&','/m3&'"                    >>LABELS.CFG
     ${MDL}/exec/concplot -ivolcano.bin -oplot_so4_$R -j${MAP} -x1.0E+09 -s2 -z80 -k2 -c4 -v10+5+2+1+0.5+0.2+0.1

  done