Particle script

Post any questions or comments regarding the LINUX version of HYSPLIT. This includes the model execution, GUIs, results, or graphics. Be sure to mention the version of LINUX you are using.
Post Reply
poprostudaria
Posts: 13
Joined: May 26th, 2015, 7:30 am
Registered HYSPLIT User: No

Particle script

Post by poprostudaria » July 21st, 2015, 8:03 am

Hi,
at the begining I would like to highlight, that I have been using HYSPLIT for one month, so I'm sorry if my question will be too easy for you. And sorry for my english too, I hope you'll understand me.
I need to run model for episode (which lasted several days) for Linux. I would like to show it as particles. And my question is: can I find somewhere script whereby I can do it? I tried to do it with concbyday.sh file (examples/scripts/unix) but it doesn't work as I wish (mainly - I didn't receive any plot).

ariel.stein
Posts: 623
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes
Contact:

Re: Particle script

Post by ariel.stein » July 22nd, 2015, 11:37 am

Please take a look at the following link for a description of the use of the particle file:
http://ready.arl.noaa.gov/documents/Tut ... _pard.html

Also, in the following link you will find a script that uses the particle plot routines:
http://ready.arl.noaa.gov/documents/Tut ... st_pard.sh

poprostudaria
Posts: 13
Joined: May 26th, 2015, 7:30 am
Registered HYSPLIT User: No

Re: Particle script

Post by poprostudaria » July 30th, 2015, 7:26 am

Hi one more time, thank you very much for your answer. I tried to do it in that way:
#!/bin/sh
# set default directory structure if not passed through

MDL="home/HYSPLIT"
OUT="/home/conc"
MET="/home/meteo"
echo test $OUT
# set model simulation variables

poll="IND"
syr=05
smo=09
shr=10
olat=51.00
olon=17.05
olvl=10.0
run=840
ztop=10000.0
met1="ARLDATA2005.bin"

for sda in 07
do
echo "$syr $smo $sda $shr " >CONTROL
echo "1 ">>CONTROL
echo "$olat $olon $olvl ">>CONTROL
echo "$run ">>CONTROL
echo "0 ">>CONTROL
echo "$ztop ">>CONTROL
echo "1 ">>CONTROL
echo "$MET/ ">>CONTROL
echo "$met1 ">>CONTROL
echo "1 ">>CONTROL
echo "$poll ">>CONTROL
echo "10.0 ">>CONTROL
echo "0.1 ">>CONTROL
echo "$syr $smo $sda $shr 00 ">>CONTROL
echo "1 ">>CONTROL
echo "30.0 70.0 ">>CONTROL
echo "0.50 0.50 ">>CONTROL
echo "20.0 20.0 ">>CONTROL
echo "$OUT/ ">>CONTROL
echo "cdump ">>CONTROL
echo "1 ">>CONTROL
echo "10 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "00 00 00 00 00 ">>CONTROL
echo "00 24 00 ">>CONTROL
echo "1 ">>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 "5.27 ">>CONTROL
echo "0.0 ">>CONTROL

if [ -f cdump ];then rm cdump; fi
${MDL}/exec/hycs_std
mv cdump ${poll}${smo}${sda}
mv CONC.CFG SETUP.CFG

${MDL}/exec/parhplot -a1
${MDL}/exec/parxplot -iPARDUMP -j${MDL}/graphics/arlmap -oparxplot
${MDL}/exec/parvplot -iPARDUMP -oparvplot.ps

echo $run
run=`expr $run - 24`
echo $run
done "

And I received a massage "PARDUMP file not found in local directory". What's wrong with my script?

ariel.stein
Posts: 623
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes
Contact:

Re: Particle script

Post by ariel.stein » July 30th, 2015, 9:00 am

Are you setting up the model to write the particle position file (PARDUMP) in the SETUP.CFG?

poprostudaria
Posts: 13
Joined: May 26th, 2015, 7:30 am
Registered HYSPLIT User: No

Re: Particle script

Post by poprostudaria » July 30th, 2015, 9:26 am

My setup.cfg looks like this:
&SETUP
tratio = 0.75,
delt = 0.0,
initd = 0,
kpuff = 0,
khmax = 9999,
khinp = 0,
numpar = 2500,
maxpar = 10000,
nbptyp = 1,
qcycle = 0.0,
efile = '',
k10m = 1,
kdef = 0,
krand = 2,
kzmix = 0,
kbls = 1,
kblt = 0,
isot = -99,
vscale = 200.0,
vscales = 200.0,
vscaleu = 200.0,
hscale = 10800.0,
capemin = -1.0,
tvmix = 1.00,
tkerd = 0.18,
tkern = 0.18,
kmix0 = 250,
kmixd = 0,
ninit = 1,
ndump = 0,
ncycl = 0,
pinbc = 'PARINBC',
pinpf = 'PARINIT',
poutf = 'PARDUMP',
messg = 'MESSAGE',
vdist = 'VMSDIST',
mgmin = 10,
conage = 48,
kmsl = 0,
kwet = 1,
ichem = 0,
cpack = 1,
cmass = 0,
kspl = 1,
krnd = 6,
frhmax = 3.00,
splitf = 1.00,
frhs = 1.00,
frvs = 0.01,
frts = 0.10,
dxf = 1.00,
dyf = 1.00,
dzf = 0.01,
/

ariel.stein
Posts: 623
Joined: November 7th, 2012, 3:14 pm
Registered HYSPLIT User: Yes
Contact:

Re: Particle script

Post by ariel.stein » July 30th, 2015, 11:29 am

You have to set ndump=1 and ncycl=1 if you want an hourly particle position file

poprostudaria
Posts: 13
Joined: May 26th, 2015, 7:30 am
Registered HYSPLIT User: No

Re: Particle script

Post by poprostudaria » July 30th, 2015, 11:53 am

Thank you very much! Now it works!

Post Reply