@echo off
setLocal EnableDelayedExpansion

cd ..
set WRK=%CD%
cd %WRK%\work

set DIR=c:
set PGM=%DIR%\hysplit4

echo -90.0   -180.0  lat/lon of lower left corner      >ASCDATA.CFG
echo 1.0     1.0     lat/lon spacing in degrees       >>ASCDATA.CFG
echo 180     360     lat/lon number of data points    >>ASCDATA.CFG
echo 2               default land use category        >>ASCDATA.CFG
echo 0.2             default roughness length (m)     >>ASCDATA.CFG
echo '%PGM%\bdyfiles\'  directory of files            >>ASCDATA.CFG
 
REM ----------------------------------------------------------

  echo 10 04 16 00              >CONTROL
  echo 2                       >>CONTROL
  echo 63.63 -19.62  100.0     >>CONTROL
  echo 63.63 -19.62 6000.0     >>CONTROL
  echo 12                      >>CONTROL
  echo 0                       >>CONTROL
  echo 32000.0                 >>CONTROL
  echo 1                       >>CONTROL
  echo ../../Tutorial/volcano/ >>CONTROL
  echo apr1420.bin             >>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 .\                      >>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

REM ----------------------------------------------------------

for %%R in (05 10 20) do (

   echo 1 2 0.%%R 1.5 >CHEMRATE.TXT

   %PGM%\exec\hycs_std  

   echo 'TITLE^&','Sulfur Dioxide %%R pcnt/h ^&' >LABELS.CFG
   echo 'UNITS^&','ppb^&'            >>LABELS.CFG
   echo 'VOLUM^&','^&'               >>LABELS.CFG
   %PGM%\exec\concplot -ivolcano.bin -oplot_so2_%%R -j%PGM%\graphics\arlmap -x0.453E+09 -s1 -z80 -k2 -c4 -v10+5+2+1+0.5+0.2+0.1

   echo 'TITLE^&','Sulfate  %%R pcnt/h ^&'        >LABELS.CFG
   echo 'UNITS^&','ug^&'             >>LABELS.CFG
   echo 'VOLUM^&','/m3^&'            >>LABELS.CFG
   %PGM%\exec\concplot -ivolcano.bin -oplot_so4_%%R -j%PGM%\graphics\arlmap -x1.0E+09 -s2 -z80 -k2 -c4 -v10+5+2+1+0.5+0.2+0.1

)