HRRR : Failed to get CRS from crs variable LambertConfor mal_Projection #38
Description
import psyplot.project
import matplotlib.pyplot as plt
n1 = 'GFS/Global_onedeg'
# n1 = 'HRRR/CONUS_2p5km' # XXX psy_maps/plotters.py:309: UserWarning: Failed to get CRS from crs variable LambertConformal_Projection
#n1 = 'NAM/Polar_90km' # XXX psy_maps/plotters.py:309: UserWarning: Failed to get CRS from crs variable PolarStereographic_Projection
nwp0 = f'https://thredds.ucar.edu/thredds/dodsC/grib/NCEP/{n1}/Best'
map0 = psyplot.project.plot.mapplot(nwp0,
name = 'Geopotential_height_isobaric',
clabel='{desc}', title = n1 + ' %FT%H:%M'
)
map0.update(z=500, method='nearest') # XXX KeyError: 'Unknown formatoption keyword lev! Possible similiar frasings are levels.'
map0.update(t = '2021-10-02', method='nearest') # XXX KeyError: 'Unknown formatoption keyword time! Possible similiar frasings are title.'
plt.show()
n0 = 'u-component_of_wind_isobaric'
vectors = psyplot.project.plot.mapvector(nwp0, name=[[n0, n0.replace('u-', 'v-')]])
plt.show()Problem description
First off : cool software !
However, I've noticed some technical difficulties :
- The first plot (w/ GFS) works fine, but if you uncomment the "n1 = 'HRRR/CONUS_2p5km'", it gives a warning :
psy_maps/plotters.py:309: UserWarning: Failed to get CRS from crs variable LambertConformal_Projection
... and then kind of hangs.
-
The plot.mapvector call dies (after some time) w/ RuntimeError: NetCDF: Access failure
-
The calls to select time & level don't seem to work. RFE : some way to put the "level" info in the title.
-
From the command line :
psyplot https://thredds.ucar.edu/thredds/dodsC/grib/NCEP/GFS/Global_onedeg/Best -n Pressure_reduced_to_MSL_msl
... seems like it should work, but I get : AttributeError: 'NoneType' object has no attribute 'variables'
FYI, this is on Fedora 34.