@echo off
setLocal EnableDelayedExpansion

cd ..
set WRK=%CD%

set DIR=c:
set PGM=%DIR%\hysplit4
cd %PGM%\working

IF EXIST ASCDATA.CFG DEL ASCDATA.CFG
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

IF NOT EXIST PARDUMP (
   echo file PARDUMP is missing
   echo Go back and run the cust_pard batch file!
   pause
)

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

SET CEN=20
SET SYR=10
SET SMO=04
SET SDA=15
SET SHR=00

SET LAT=63.63
SET LON=-19.62

SET HST=jwocky.gsfc.nasa.gov
SET DAT=L3_aersl_omi_%CEN%%SYR%%SMO%%SDA%.txt

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

echo user anonymous john.smith@gmail.com   >ftp.txt
echo cd /pub/omi/data/aerosol/Y%CEN%%SYR% >>ftp.txt
echo binary                               >>ftp.txt  
echo prompt                               >>ftp.txt
echo get %DAT% ga%SYR%%SMO%%SDA%.txt      >>ftp.txt
echo bye                                  >>ftp.txt

ftp -n %HST% <ftp.txt

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

ECHO 'TITLE^&','### %0: shift at 1000 km ### ^&' >LABELS.CFG
%PGM%\exec\parshift -iPARDUMP -oPARINIT -r10.0:1000.0:%LAT%:%LON%

IF EXIST parxplot.ps DEL parxplot.ps
%PGM%\exec\parxplot -iPARINIT -j%PGM%\graphics\arlmap -k1 -z80 
parxplot.ps
PAUSE

%PGM%\exec\autoview -i0.5 -m0.5 -cvolcano.bin -pPARINIT

ECHO 'TITLE^&','### %0: blended shift at 1000 km ### ^&' >LABELS.CFG
%PGM%\exec\parshift -iPARDUMP -oPARINIT -b -r10.0:1000.0:%LAT%:%LON%

IF EXIST parxplot.ps DEL parxplot.ps
%PGM%\exec\parxplot -iPARINIT -j%PGM%\graphics\arlmap -k1 -z80 
parxplot.ps
PAUSE

%PGM%\exec\autoview -i0.5 -m0.5 -cvolcano.bin -pPARINIT