@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 2                       >>CONTROL
echo RNUC                    >>CONTROL
echo 1.0                     >>CONTROL
echo 51.0                    >>CONTROL
echo 00 00 00 00 00          >>CONTROL
echo NGAS                    >>CONTROL
echo 1.0                     >>CONTROL
echo 51.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_unit.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 2                       >>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 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 maxdim = 2,          >>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 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 ---------------------------------------------------------------------

%PGM%\exec\con2rem -ifdnpp_unit.bin -ofdnpp_dose.bin -a..\files\activity.txt -d1

%PGM%\exec\concsum -ifdnpp_dose.bin -ofdnpp_total.bin -l -pDOSE

%PGM%\exec\conavgpd -ifdnpp_total.bin -ofdnpp_sums.bin -a11031409 -b11031612 -r1


REM -------- Dose Time Series at Tokai-mura ---------

echo Tokai 36.4356 140.6025 >latlon.txt
%PGM%\exec\con2stn  -ifdnpp_total.bin -oTokai.txt -slatlon.txt

ECHO 'TITLE^&','Total Dose^&' >LABELS.CFG
ECHO 'UNITS^&','R^&'         >>LABELS.CFG
ECHO 'VOLUM^&','^&'          >>LABELS.CFG
%PGM%\exec\timeplot -iTokai.txt -p2 -y
IF EXIST dose_time.ps DEL dose_time.ps
RENAME timeplot.ps dose_time.ps

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

ECHO 'TITLE^&','Total Dose^&' >LABELS.CFG
ECHO 'UNITS^&','R^&'         >>LABELS.CFG
ECHO 'VOLUM^&','^&'          >>LABELS.CFG
%PGM%\exec\concplot -ifdnpp_sums.bin -odose_plot -j%PGM%\graphics\arlmap -z95 -h37.0:140.0 -g0:500 -k2 -c4 -v1.0+0.3+0.1+0.03+0.01+0.003+0.001