@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 ^&SETUP               >SETUP.CFG
echo maxdim = 2,          >>SETUP.CFG
echo delt = 5.0,          >>SETUP.CFG
echo initd = 0,           >>SETUP.CFG
echo khmax = 9999,        >>SETUP.CFG
echo numpar = -2500,      >>SETUP.CFG
echo maxpar = 10000,      >>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

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

set RUN=72
for %%D in (14 15 16) do (
for %%H in (00 03 06 09 12 15 18 21) do (

echo 11 03 %%D %%H            >CONTROL
echo 2                       >>CONTROL
echo 37.4206 141.0329 1.0    >>CONTROL
echo 37.4206 141.0329 100.0  >>CONTROL
echo !RUN!                   >>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 3.0                     >>CONTROL
echo 00 00 00 00 00          >>CONTROL
echo NGAS                    >>CONTROL
echo 1.0                     >>CONTROL
echo 3.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 TG_03%%D%%H             >>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

%PGM%\exec\hycs_std

set /a RUN=!RUN!-3

REM ---- terminate day hour loop
)
)

rem ---------------------------------------------

%PGM%\exec\condecay -1:1:11025.8:C137 -3:1:8.02330:I131 -4:1:754.020:C134 -5:1:0.86958:I133 -6:1:1.10167:B140 -7:1:1.44980:L140 -8:1:3.23000:T132 -9:1:2.18000:A110 -10:1:3.0200:NB95 -14:2:5.2474:X133 +e..\files\cfactors_all +oCG_ +t031106

for %%D in (14 15 16) do (
for %%H in (00 03 06 09 12 15 18 21) do (
  %PGM%\exec\con2rem -iCG_03%%D%%H -oDG_03%%D%%H -s1 -t0 -d1 -a..\files\activity_unit.txt
  %PGM%\exec\concsum -iDG_03%%D%%H -oSG_03%%D%%H -l -pDOSE
  echo Finished: %%D %%H
)
)

dir /b SG_??????.bin >merg_list.txt
%PGM%\exec\conmerge -imerg_list.txt -ofdnpp_total.bin

%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 Rate 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 -umR/h -k2 -c4 -v1.0+0.3+0.1+0.03+0.01+0.003+0.001