Concentration plot showing the whole Northern Hemisphere when I only want Northern Europe

General questions and postings pertaining to the use of HYSPLIT regardless of the platform. For platform specific questions, use the HYSPLIT Platform forums.
Post Reply
padraigflattery
Posts: 5
Joined: January 27th, 2022, 7:47 am
Registered HYSPLIT User: Yes

Concentration plot showing the whole Northern Hemisphere when I only want Northern Europe

Post by padraigflattery »

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: 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...: Image

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?
alicec
Posts: 416
Joined: February 8th, 2016, 12:56 pm
Registered HYSPLIT User: Yes

Re: Concentration plot showing the whole Northern Hemisphere when I only want Northern Europe

Post by alicec »

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
padraigflattery
Posts: 5
Joined: January 27th, 2022, 7:47 am
Registered HYSPLIT User: Yes

Re: Concentration plot showing the whole Northern Hemisphere when I only want Northern Europe

Post by padraigflattery »

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 :(
sonny.zinn
Posts: 354
Joined: May 8th, 2019, 1:31 pm
Registered HYSPLIT User: Yes

Re: Concentration plot showing the whole Northern Hemisphere when I only want Northern Europe

Post by sonny.zinn »

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.
padraigflattery
Posts: 5
Joined: January 27th, 2022, 7:47 am
Registered HYSPLIT User: Yes

Re: Concentration plot showing the whole Northern Hemisphere when I only want Northern Europe

Post by padraigflattery »

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

Thank you.
Post Reply

Return to “Users”