HYSPLIT Forum: hysplitbbs.arl.noaa.gov

I'm running a simulation from Eastern Europe and I want to visualize the output over Europe. To do this I've set up the control file like so:

#!/bin/sh
now=$(date)
echo "Current date and time: $now"
#-------------------------------------------------------------

WEB=""
if [ -f /usr/bin/firefox ];then WEB="/usr/bin";fi
if [ -f /usr/local/bin/firefox ];then WEB="/usr/local/bin";fi
MDL="${HOME}/Dispersion//hysplit.v5.2.0_CentOS7.9.2009"
OUT="${MDL}/outfiles"
MET="${HOME}/Dispersion/"
cd $OUT

if [ ! -f ASCDATA.CFG ]; then
echo "-90.0 -180.0" >ASCDATA.CFG
echo "1.0 1.0" >>ASCDATA.CFG
echo "180 360" >>ASCDATA.CFG
echo "2" >>ASCDATA.CFG
echo "0.2" >>ASCDATA.CFG
echo "'$MDL/bdyfiles/'" >>ASCDATA.CFG
fi
echo "### $0 ###"

#--------------------------------------------------------------

syr=$(date "+%y")
smo=$(date "+%m")
sda=$(date "+%d")
shr=00

olat=50.00
olon=26.00
lvl1=10.0

run=240
ztop=10000.0
data=$file

#----------------------------------------------------------

echo "$syr $smo $sda $shr " >CONTROL
echo "1 ">>CONTROL
echo "$olat $olon $lvl1 ">>CONTROL
echo "$run ">>CONTROL
echo "0 ">>CONTROL
echo "$ztop ">>CONTROL
echo "1 ">>CONTROL
echo "$MET/ ">>CONTROL
echo "$data ">>CONTROL
echo "1 ">>CONTROL
echo "C137 ">>CONTROL
echo "1.0E+16 ">>CONTROL
echo "1.0 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "1 ">>CONTROL
echo "0.0 0.0 ">>CONTROL
echo "0.05 0.05 ">>CONTROL
echo "30.0 30.0 ">>CONTROL
echo "$OUT/ ">>CONTROL
echo "cdump ">>CONTROL
echo "1 ">>CONTROL
echo "1000 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "00 12 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 "10960.0 ">>CONTROL
echo "0.0 ">>CONTROL

#----------------------------------------------------------

rm -f cdump
rm -f SETUP.CFG

${MDL}/exec/hycs_std

echo "'TITLE&','10-Day Dispersion&'" >LABELS.CFG
${MDL}/exec/concplot +g1 -icdump -c1 -z40 -h50:25 -j${MDL}/graphics/arlmap
${MDL}/exec/splitsvg -iconcplot.html -otemp.svg

/usr/bin/convert -trim +repage -density 150 -loop 100 -delay 50 -background white \
-dispose Background F??-temp.svg concplot.gif

#UNCOMMENT IF YOU WANT TO DISPLAY IN FIREFOX
#if [[ "$OSTYPE" == "darwin"* ]]; then
# open concplot.html
#else
# ${WEB}/firefox concplot.gif &
# ${WEB}/firefox concplot.html &
#fi



Which produces the following output image:

But even though the plume only travels over Europe, the output shows the entire global map.

I have been playing with the settings and can adjust the parameters in the control file like so:

#!/bin/sh
now=$(date)
echo "Current date and time: $now"
#-------------------------------------------------------------

WEB=""
if [ -f /usr/bin/firefox ];then WEB="/usr/bin";fi
if [ -f /usr/local/bin/firefox ];then WEB="/usr/local/bin";fi
MDL="${HOME}/hysplit.v5.2.0_CentOS7.9.2009"
OUT="${MDL}/outfiles"
MET="${HOME}/Dispersion/data"
cd $OUT

if [ ! -f ASCDATA.CFG ]; then
echo "-90.0 -180.0" >ASCDATA.CFG
echo "1.0 1.0" >>ASCDATA.CFG
echo "180 360" >>ASCDATA.CFG
echo "2" >>ASCDATA.CFG
echo "0.2" >>ASCDATA.CFG
echo "'$MDL/bdyfiles/'" >>ASCDATA.CFG
fi
echo "### $0 ###"

#--------------------------------------------------------------

syr=$(date "+%y")
smo=$(date "+%m")
sda=$(date "+%d")
shr=00

olat=50.30
olon=26.64
lvl1=10.0

run=240
ztop=10000.0
data=$file

#----------------------------------------------------------

echo "$syr $smo $sda $shr " >CONTROL
echo "1 ">>CONTROL
echo "$olat $olon $lvl1 ">>CONTROL
echo "$run ">>CONTROL
echo "0 ">>CONTROL
echo "$ztop ">>CONTROL
echo "1 ">>CONTROL
echo "$MET/ ">>CONTROL
echo "$data ">>CONTROL
echo "1 ">>CONTROL
echo "C137 ">>CONTROL
echo "1.0E+16 ">>CONTROL
echo "1.0 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "1 ">>CONTROL
echo "15.0 50.0 ">>CONTROL
echo "0.05 0.05 ">>CONTROL
echo "30.0 50.0 ">>CONTROL
echo "$OUT/ ">>CONTROL
echo "cdump ">>CONTROL
echo "1 ">>CONTROL
echo "1000 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "00 12 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 "10960.0 ">>CONTROL
echo "0.0 ">>CONTROL

#----------------------------------------------------------

rm -f cdump
rm -f SETUP.CFG

${MDL}/exec/hycs_std

echo "'TITLE&','10-Day Dispersion &'" >LABELS.CFG
${MDL}/exec/concplot +g1 -icdump -c1 -z40 -h50:25 -j${MDL}/graphics/arlmap
${MDL}/exec/splitsvg -iconcplot.html -otemp.svg

/usr/bin/convert -trim +repage -density 150 -loop 100 -delay 50 -background white \
-dispose Background F??-temp.svg concplot.gif

#UNCOMMENT IF YOU WANT TO DISPLAY IN FIREFOX
#if [[ "$OSTYPE" == "darwin"* ]]; then
# open concplot.html
#else
# ${WEB}/firefox concplot.gif &
# ${WEB}/firefox concplot.html &
#fi


And it shows only the desired area (Europe) - but by adjusting the grid span like this the plume is cut off around some locations in northern and western Europe, which is not what I want...:

Is there something I am missing here? I have adjusted the zoom from 0 to 90 and it has had no real effect on the scaling of the map, when the grid span covers the area I want, the output map ends up being the entire northern hemisphere.

It's really difficult to distinguish individual countries on the northern hemisphere map. Is it possible to limit the view to just western Europe while having a larger extent for the plume, no matter where the plume goes?

Try keeping the grid center at 15.0, 50.0 but increasing the span to 80.

change this
echo "15.0 50.0 ">>CONTROL
echo "0.05 0.05 ">>CONTROL
echo "30.0 50.0 ">>CONTROL

to
echo "15.0 50.0 ">>CONTROL
echo "0.05 0.05 ">>CONTROL
echo "30.0 80.0 ">>CONTROL

Oh dear I have made an error in my copying of the control file - the should be

echo "50.0 15.0 ">>CONTROL

for the center of grid, not the other way around.

When I run it like this as you've said, with:
echo "50.0 15.0 ">>CONTROL
echo "0.05 0.05 ">>CONTROL
echo "30.0 80.0 ">>CONTROL

I still get a map of the entire northern hemisphere

You might want to use the -g option in addition to other options you are already using for the CONCPLOT program. For example, with -g0:250 -h50:25, the plot will be spatially limited to an approximate radius of 250 km from lat=50 deg and lon=25 deg.

Oh thank you so much - I think this has worked!

Thank you.