@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 11 03 14 09              >CONTROL0
echo 2                       >>CONTROL0
echo 37.4206 141.0329 1.0    >>CONTROL0
echo 37.4206 141.0329 100.0  >>CONTROL0
echo 192                     >>CONTROL0
echo 0                       >>CONTROL0
echo 10000.0                 >>CONTROL0
echo 1                       >>CONTROL0
echo ../../Tutorial/japan/   >>CONTROL0
echo gdas11-22.bin           >>CONTROL0
echo 1                       >>CONTROL0
echo CPAR                    >>CONTROL0
echo 0.0                     >>CONTROL0
echo 0.0                     >>CONTROL0
echo 00 00 00 00 00          >>CONTROL0
echo 2                       >>CONTROL0
echo 38.0 140.0              >>CONTROL0
echo 0.25 0.25               >>CONTROL0
echo 20.0 30.0               >>CONTROL0
echo ./                      >>CONTROL0
echo deposit.bin             >>CONTROL0
echo 1                       >>CONTROL0
echo 0                       >>CONTROL0
echo 00 00 00 00 00          >>CONTROL0
echo 00 00 00 00 00          >>CONTROL0
echo 00 192 00               >>CONTROL0
echo 38.0 140.0              >>CONTROL0
echo 0.25 0.25               >>CONTROL0
echo 20.0 30.0               >>CONTROL0
echo ./                      >>CONTROL0
echo airconc.bin             >>CONTROL0
echo 1                       >>CONTROL0
echo 100                     >>CONTROL0
echo 00 00 00 00 00          >>CONTROL0
echo 00 00 00 00 00          >>CONTROL0
echo 00 03 00                >>CONTROL0
echo 1                       >>CONTROL0
echo 1.0 1.0 1.0             >>CONTROL0
echo 0.001 0.0 0.0 0.0 0.0   >>CONTROL0
echo 0.0 8.0E-05 8.0E-05     >>CONTROL0
echo 0.0                     >>CONTROL0

echo ^&SETUP                       >SETUP.CFG
echo delt = 15.0,                 >>SETUP.CFG
echo initd = 0,                   >>SETUP.CFG
echo khmax = 12,                  >>SETUP.CFG
echo numpar = -2000,              >>SETUP.CFG
echo maxpar = 1000000,            >>SETUP.CFG
echo efile = '../files/EMITIMES', >>SETUP.CFG
echo /                            >>SETUP.CFG

REM -------- Calculation with NO resuspension ----------------------------------------------------------------

COPY CONTROL0    CONTROL
ECHO 0.0       >>CONTROL

%PGM%\exec\hycs_std

echo 'TITLE^&','Without Resuspension^&'  >LABELS.CFG 
%PGM%\exec\concplot.exe -ideposit.bin -oemit_surf1.ps -j%PGM%\graphics\arlmap -c4 -y0.001 -z80 -ukBq -g0:1000 -h37.0:144.0 -v100+50+20+10+5+2+1
%PGM%\exec\concplot.exe -iairconc.bin -oemit_surf2.ps -j%PGM%\graphics\arlmap -c4 -x1.0 -z80 -uBq -g0:1000 -h37.0:144.0 -v100+10+1+0.1+0.01+0.001

%PGM%\exec\c2datem  -iairconc.bin -oTokai.txt -m..\files\JAEA_C137_all.txt -c1000.0
ECHO 'TITLE^&','Cs-137 time series at Tokai-mura^&' >LABELS.CFG
ECHO 'UNITS^&',' mBq^&'                            >>LABELS.CFG
%PGM%\exec\timeplot -iTokai.txt -s..\files\JAEA_C137_all.txt -p2 -y
RENAME timeplot.ps emit_surf3.ps

REM -------- Calculation WITH resuspension ----------------------------------------------------------------

COPY CONTROL0    CONTROL
echo 1.0E-06   >>CONTROL

%PGM%\exec\hycs_std

ECHO 'TITLE^&','With Resuspension^&'  >LABELS.CFG    
%PGM%\exec\concplot.exe -ideposit.bin -oemit_surf4.ps -j%PGM%\graphics\arlmap -c4 -y0.001 -z80 -ukBq -g0:1000 -h37.0:144.0 -v100+50+20+10+5+2+1
%PGM%\exec\concplot.exe -iairconc.bin -oemit_surf5.ps -j%PGM%\graphics\arlmap -c4 -x1.0 -z80 -uBq -g0:1000 -h37.0:144.0 -v100+10+1+0.1+0.01+0.001

%PGM%\exec\c2datem  -iairconc.bin -oTokai.txt -m..\files\JAEA_C137_all.txt -c1000.0
ECHO 'TITLE^&','Cs-137 time series at Tokai-mura^&' >LABELS.CFG
ECHO 'UNITS^&',' mBq^&'                            >>LABELS.CFG
%PGM%\exec\timeplot -iTokai.txt -s..\files\JAEA_C137_all.txt -p2 -y
RENAME timeplot.ps emit_surf6.ps