@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              >CONTROL
echo 2                       >>CONTROL
echo 37.4206 141.0329 1.0    >>CONTROL
echo 37.4206 141.0329 100.0  >>CONTROL
echo 51                      >>CONTROL
echo 0                       >>CONTROL
echo 10000.0                 >>CONTROL
echo 1                       >>CONTROL
echo ../files/               >>CONTROL
echo wrf11031400.bin         >>CONTROL
echo 1                       >>CONTROL
echo CPAR                    >>CONTROL
echo 0.0                     >>CONTROL
echo 0.0                     >>CONTROL
echo 00 00 00 00 00          >>CONTROL
echo 1                       >>CONTROL
echo 38.0 140.0              >>CONTROL
echo 0.05 0.05               >>CONTROL
echo 20.0 30.0               >>CONTROL    
echo ./                      >>CONTROL
echo fdnpp.bin               >>CONTROL
echo 2                       >>CONTROL
echo 0 100                   >>CONTROL
echo 00 00 00 00 00          >>CONTROL
echo 00 00 00 00 00          >>CONTROL
echo 00 03 00                >>CONTROL
echo 1                       >>CONTROL
echo 1.0 1.0 1.0             >>CONTROL
echo 0.001 0.0 0.0 0.0 0.0   >>CONTROL
echo 0.0 8.0E-05 8.0E-05     >>CONTROL
echo 0.0                     >>CONTROL
echo 0.0                     >>CONTROL

echo ^&SETUP                       >SETUP.CFG
echo delt = 5.0,                  >>SETUP.CFG
echo initd = 0,                   >>SETUP.CFG
echo khmax = 24,                  >>SETUP.CFG
echo numpar = -2500,              >>SETUP.CFG
echo maxpar = 300000,             >>SETUP.CFG
echo efile = '../files/EMITIMES', >>SETUP.CFG
echo krand = 1,                   >>SETUP.CFG
echo vscales = 5.0,               >>SETUP.CFG
echo vscaleu = 200.0,             >>SETUP.CFG
echo kmix0 = 150,                 >>SETUP.CFG
echo rhb = 100,                   >>SETUP.CFG
echo rht = 80,                    >>SETUP.CFG
echo /                            >>SETUP.CFG

%PGM%\exec\hycs_std

REM -------- Time Series at JAEA ---------

%PGM%\exec\c2datem  -ifdnpp.bin -oTokai.txt -m..\files\JAEA_C137.txt -c1000.0 -z2
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.txt -p2 -r2 -y
RENAME timeplot.ps plot_time.ps

REM -------- Deposition Plot -------------

ECHO 'TITLE^&','Cs-137 deposition ^&' >LABELS.CFG
ECHO 'UNITS^&',' kBq^&'  >>LABELS.CFG
%PGM%\exec\concplot -ifdnpp.bin -oplot_dept -j%PGM%\graphics\arlmap -z95 -h37.0:140.0 -g0:500 -y0.001 -ukBq -b0 -t0 -k2 -c4 -v1000+500+200+100+50+20+10 -r3

REM -------- Simple Dose Plot -------------

ECHO 'TITLE^&','Cs-137 dose rate^&' >LABELS.CFG
ECHO 'UNITS^&',' uR^&' >>LABELS.CFG
ECHO 'VOLUM^&','/h^&'  >>LABELS.CFG
%PGM%\exec\concplot -ifdnpp.bin -oplot_dose -j%PGM%\graphics\arlmap -z95 -h37.0:140.0 -g0:500 -k2 -x3.340E-05 -y1.08E-06 -c4 -v0.1+0.05+0.02+0.01+0.005+0.002+0.001+0.0005+0.0002+0.0001 -r2