search

sharppy / SHARPpy

Sounding/Hodograph Analysis and Research Program in Python
https://sharppy.github.io/SHARPpy/index.html
Other
216 stars 112 forks source link

Issue When Opening Newer RAOBs #117

Closed jpan84 closed 7 years ago

jpan84 commented 7 years ago

I take FSL format files from esrl.noaa.gov/raobs/ and convert to SHARPpy's tabular form using my own program. Older data have no issues. OTOH, often, when I try to open RAOB data from after 2011 or so, I get the following error:

"index -1 is out of bounds for axis 0 with size 0"

87111112.txt

12103003.txt

georgephillips1 commented 7 years ago

Don't know if this is the problem, but the surface level at 961 mb shown in the FSL output is not included in your text file in the 2012 sounding data.

jpan84 commented 7 years ago

@georgephillips1 you're right. I deleted that because it caused height out of order error. Should I rearrange it instead of deleting?

georgephillips1 commented 7 years ago

@jpan84 that's what I did. I had to do a sort to get the right format. I was working with FSL output file from your 2012 case, and there are other issues with a direct conversion that you will likely find. One is that there are two cases with the same pressure level/height duplicated. SHARPpy doesn't like that. The other is that the FSL output has 360 degrees for a due north wind, whereas SHARPpy expects 0 degrees. I opened a ticket on the 2nd issue and maybe should on the first as well. I can code workarounds for both in my conversion program, but that doesn't help the next person writing a conversion program.

jpan84 commented 7 years ago

@georgephillips1 Oh yes, I realized the 360.00 issue and have just used find and replace to deal with that. As for sorting instead of deleting, it worked. Thanks for your help.