from osgeo import gdal
import numpy as np
hrrrPath = 'https://noaa-hrrr-bdp-pds.s3.amazonaws.com/hrrr.20201117/conus/hrrr.t00z.wrfsfcf00.grib2'
ds = gdal.Open('/vsicurl/' + hrrrPath)