pysplit
pysplit copied to clipboard
Ensemble mode in back trajectory
Will it be possible to add ensemble mode in the code, I was going through one of the forum of HYSPLIT and the following is suggested in one of the posts using bash script,
#!/bin/sh
MDL="${HOME}/hysplit4"
MET="${HOME}/Tutorial/captex"
OUT="."
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
echo "83 09 25 17 " >CONTROL
echo "3 " >>CONTROL
echo "38.00 -84.00 750.0 " >>CONTROL
echo "44.00 -74.00 750.0 " >>CONTROL
echo "39.00 -83.00 750.0 " >>CONTROL
echo "6 " >>CONTROL
echo "0 " >>CONTROL
echo "10000.0 " >>CONTROL
echo "1 " >>CONTROL
echo "$MET/ " >>CONTROL
echo "captex2_narr.bin " >>CONTROL
echo "./ " >>CONTROL
echo "fdump " >>CONTROL
rm -f SETUP.CFG
${MDL}/exec/latlon
${MDL}/exec/hyts_std
${MDL}/exec/trajplot -v0 -z80 -j${MDL}/graphics/arlmap -ifdump
cat trajplot.ps >>results.ps
rm -f fdump trajplot.ps
# gs results.ps
Apparently, a file named latlon creates the latitude and longitude for the ensemble and then 'hyts_std' is run for each of the coordinates. It will be great if its possible to include something called in bulk_trajectory.py so that user can choose the mode he wants.
Thanks, Nirav
Thanks for the question, Nirav!
I'm planning to add all the remaining trajectory modes- Trajectory matrices (like in your script), meteorological grid offset ensembles, and multiple trajectories in both space and time. These may be separate functions in order to keep complexity down, but we'll see. I have recently moved house again, and within the next week or two when my new place and home office are functional I will start a series of updates, which will be sketched out in a "Coming Soon" section of the README.
Cheers, Mellissa
Has a method for computing ensembles been added? thanks