zhx215
commented
6 years ago Hi,
I am glad to answer your question... "fpath" should be replaced by the directory of trajectory files (generated by HYSPLIT GUI); for example: trajgroup = pysplit.make_trajectorygroup(r'C:/trajectories/colgate/feb')
here "feb" refers all February trajectories (I use PYSPLIT codes to generate trajectories and this program will rename all files to certain criteria specified). Ideally, you can add other information to "group" your trajectories, such as "1000" to read only trajectories with starting-point altitude of 1000m AGL. Here are the code I use to run the example codes (to plot the trajectories on map):
import pysplit
%matplotlib inline
trajgroup = pysplit.make_trajectorygroup(r'C:/trajectories/colgate/*feb*')
mapcorners = [-150, 15, -40, 75]
standard_pm = None
bmap_params = pysplit.MapDesign(mapcorners, standard_pm)
bmap = bmap_params.make_basemap()
color_dict = {500.0 : 'blue',
1000.0 : 'orange',
1500.0 : 'red'}
for traj in trajgroup:
altitude0 = traj.data.geometry.apply(lambda p: p.z)[0]
traj.trajcolor = color_dict[altitude0]
for traj in trajgroup[::4]:
bmap.plot(*traj.path.xy, c=traj.trajcolor, latlon=True, zorder=20)
#add two lines of codes to plot coast and country lines; the zorder matters here.
bmap.drawcoastlines(linewidth=1.0, color='black', zorder=20)
bmap.drawcountries(linewidth=0.5, color='grey', zorder=20)I recommend you learn how to generate trajectories using pysplit code interacting with hysplit program. There is an example code available already. I also recommend you to learn some python basics and the use of jupyter notebook... in the code files, there are some explanatory notes that can help you understand the structure of pysplit....
Best, Ryan
heatherguy
mscross
Hello Mellissa,
This is my first time using HYSPLIT and pysplit. I am trying to create some trajectory plots from HYSPLIT trajectories that I have run previously using the HYSPLIT Windows GUI. I started by trying to follow your example basic_plotting_example.py.
When I try to run the command: trajgroup = pysplit.make_trajectorygroup(fpath)
It returns the following error: OutOfBoundsDatetime: Out of bounds nanosecond timestamp: 4016-08-31 23:00:00
I am not sure where this error comes from or how I can trace it back to its source. It is possible that it could be an error in how I initialized the HYSPLIT trajectory in the first place. For instance the datetime for this trajectory should be 2016-08-31 23:00:00 not 4016. I can't find any information in the HYSPLIT documentation about this.
Have you seen this problem before? Do you have any ideas for how I might fix it in python without having to re-run all my trajectories? I am not sure if this is the right place to post this question but I hope you might point me in the right direction.
I am using python 3.6.
Thank you!
Heather