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).
Particle script
-
- Posts: 19
- Joined: May 26th, 2015, 7:30 am
- Registered HYSPLIT User: No
-
- Posts: 660
- Joined: November 7th, 2012, 3:14 pm
- Registered HYSPLIT User: Yes
Re: Particle script
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
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
-
- Posts: 19
- Joined: May 26th, 2015, 7:30 am
- Registered HYSPLIT User: No
Re: Particle script
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?
#!/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?
-
- Posts: 660
- Joined: November 7th, 2012, 3:14 pm
- Registered HYSPLIT User: Yes
Re: Particle script
Are you setting up the model to write the particle position file (PARDUMP) in the SETUP.CFG?
-
- Posts: 19
- Joined: May 26th, 2015, 7:30 am
- Registered HYSPLIT User: No
Re: Particle script
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,
/
&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,
/
-
- Posts: 660
- Joined: November 7th, 2012, 3:14 pm
- Registered HYSPLIT User: Yes
Re: Particle script
You have to set ndump=1 and ncycl=1 if you want an hourly particle position file
-
- Posts: 19
- Joined: May 26th, 2015, 7:30 am
- Registered HYSPLIT User: No
Re: Particle script
Thank you very much! Now it works!