Page 1 of 1

Particle script

Posted: July 21st, 2015, 8:03 am
by poprostudaria
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).

Re: Particle script

Posted: July 22nd, 2015, 11:37 am
by ariel.stein
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

Re: Particle script

Posted: July 30th, 2015, 7:26 am
by poprostudaria
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?

Re: Particle script

Posted: July 30th, 2015, 9:00 am
by ariel.stein
Are you setting up the model to write the particle position file (PARDUMP) in the SETUP.CFG?

Re: Particle script

Posted: July 30th, 2015, 9:26 am
by poprostudaria
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,
/

Re: Particle script

Posted: July 30th, 2015, 11:29 am
by ariel.stein
You have to set ndump=1 and ncycl=1 if you want an hourly particle position file

Re: Particle script

Posted: July 30th, 2015, 11:53 am
by poprostudaria
Thank you very much! Now it works!