mscross / pysplit

A package for HYSPLIT air parcel trajectory analysis.
BSD 3-Clause "New" or "Revised" License
154 stars 82 forks source link

NAM 12km Meteorology #40

Open wyowind opened 6 years ago

wyowind commented 6 years ago

Hi Mellissa,

I have been using the web-based version of HYSPLIT to generate many trajectories with the NAM 12 km (pressure) meteorology. Because this process is onerous, I was wondering if your trajectory_generator function supports these NAM files, which are produced daily at 6 hr intervals. Following your instructions for the GDAS files, I get an error, presumably because the filenames are formatted differently. Any help is appreciated!

Thanks, Tim

mscross commented 6 years ago

Hi Tim,

Thanks for your interest in PySPLIT! PySPLIT's trajectory generator theoretically supports any ARL-packed meteorology file, BUT the filenames have to be formatted appropriately in order for PySPLIT to find them, and the bulk generator, which goes month by month, is right now (coming updates will extend the PySPLIT trajectory generation functionality) optimized for weekly, semi-monthly, and monthly files in the way it searches for meteorology. That means that with daily/sub-daily files, depending on how comprehensive your coverage is per month, you might run into an issue where PySPLIT finds a large number of matching meteorology files and supplying all of those to HYSPLIT makes the latter error. Short answer, right now it's a MAYBE.

Mellissa

wyowind commented 6 years ago

Hi Mellissa,

Thank you for the prompt reply! I look forward to future updates.

Tim

brendano257 commented 6 years ago

Hey Tim,

I had a similar issue to yours, where I needed to use higher resolution met data (3Km HRRR 6-hour files in my case), and also had to call it based on the times of specific air pollution events, so bulk generation was unnecessary. My fork of Mellissa's Pysplit has "generate_singletraj" which lets you pass a list of met files for each run instead of having generate_bulktraj parse the files for you. This means you can use any ARL-packaged met files so long as you parse the names yourself. I can also send you test code parsing and using the HRRR files if you'd like.

Cheers, Brendan