@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

echo 07 04 16 00              >CONTROL
echo 1                       >>CONTROL
echo 31.15 -82.40 100.0      >>CONTROL
echo 96                      >>CONTROL
echo 0                       >>CONTROL
echo 10000.0                 >>CONTROL
echo 1                       >>CONTROL
echo ../../Tutorial/smoke/   >>CONTROL
echo EDAS16to20.bin          >>CONTROL
echo 1                       >>CONTROL
echo PM25                    >>CONTROL
echo 1.0                     >>CONTROL
echo 72.0                    >>CONTROL
echo 00 00 00 18 00          >>CONTROL
echo 1                       >>CONTROL
echo 0.0 0.0                 >>CONTROL
echo 0.1 0.1                 >>CONTROL
echo 20.0 20.0               >>CONTROL    
echo ./                      >>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 == 1 SET heat=0.215E+11
IF %%H == 2 SET heat=0.430E+11
IF %%H == 3 SET heat=0.860E+11
IF %%H == 4 SET heat=1.720E+11

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

%PGM%\exec\hycs_std

FINDSTR /c:"mixd,rise" MESSAGE >rise_%%H.txt

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

echo 'TITLE^&','### %%H ###^&'  >LABELS.CFG
%PGM%\exec\concplot -ismoke.bin -j%PGM%\graphics\arlmap -c4 -k1 -z90 -uug -n40:40 -v500+200+100+50+20+10
IF EXIST plot_%%H.ps DEL plot_%%H.ps
RENAME concplot.ps plot_%%H.ps

echo Mayo 30.2556 -81.4533 >samplers.txt
%PGM%\exec\con2stn -ismoke.bin -ocon2stn.txt -r1 -ssamplers.txt

%PGM%\exec\timeplot -icon2stn.txt -s..\..\Tutorial\smoke\mayo_pm25.txt
IF EXIST time_%%H.ps DEL time_%%H.ps 
RENAME timeplot.ps time_%%H.ps 

)