In [34]:
import netCDF4
from pyudl.cdmr import Dataset
from pyudl.tds import TDSCatalog

tstvar = "Geopotential_height_isobaric"

In [55]:
cat = TDSCatalog('http://localhost:8080/thredds/catalog.xml?dataset=HRRR_nc3')
url = cat.datasets.values()[0].accessUrls['CdmRemote']

In [56]:
%%timeit -n 1 -r 1
ds = Dataset(url)
#lcc_info = ds.variables['LambertConformal_Projection']
temp = ds.variables[tstvar]
#x = ds.variables['x'][:] * 1000.
#y = ds.variables['y'][:] * 1000.
tempData = temp[0, :]


1 loops, best of 1: 4.88 s per loop

In [57]:
%%timeit -n 1 -r 1
ds = netCDF4.Dataset(cat.datasets.values()[0].accessUrls['OPENDAP'])
#lcc_info = ds.variables['LambertConformal_Projection']
temp = ds.variables[tstvar]
#x = ds.variables['x'][:] * 1000.
#y = ds.variables['y'][:] * 1000.
tempData = temp[1, :]


1 loops, best of 1: 8.9 s per loop

In [58]:
#%timeit -n 1 -r 1 ds = Dataset(url)

In [59]:
#%timeit -n 1 -r 1 ds = netCDF4.Dataset(cat.datasets.values()[0].accessUrls['OPENDAP'])

In [60]:
%%timeit -n 1 -r 1
ds = netCDF4.Dataset("/Users/lesserwhirls/.unidata/tds/content/thredds/public/testdata/HRRR.nc")
#lcc_info = ds.variables['LambertConformal_Projection']
temp = ds.variables[tstvar]
#x = ds.variables['x'][:] * 1000.
#y = ds.variables['y'][:] * 1000.
tempData = temp[2, :]


1 loops, best of 1: 922 ms per loop


netCDF-4



In [61]:
cat = TDSCatalog('http://localhost:8080/thredds/catalog.xml?dataset=HRRR_nc4')
url = cat.datasets.values()[0].accessUrls['CdmRemote']

In [62]:
%%timeit -n 1 -r 1
ds = Dataset(url)
#lcc_info = ds.variables['LambertConformal_Projection']
temp = ds.variables[tstvar]
#x = ds.variables['x'][:] * 1000.
#y = ds.variables['y'][:] * 1000.
tempData = temp[0, :]


1 loops, best of 1: 7.39 s per loop

In [63]:
%%timeit -n 1 -r 1
ds = netCDF4.Dataset(cat.datasets.values()[0].accessUrls['OPENDAP'])
#lcc_info = ds.variables['LambertConformal_Projection']
temp = ds.variables[tstvar]
#x = ds.variables['x'][:] * 1000.
#y = ds.variables['y'][:] * 1000.
tempData = temp[1, :]


1 loops, best of 1: 10.8 s per loop

In [64]:
%%timeit -n 1 -r 1
ds = netCDF4.Dataset("/Users/lesserwhirls/.unidata/tds/content/thredds/public/testdata/HRRR.nc4")
#lcc_info = ds.variables['LambertConformal_Projection']
temp = ds.variables[tstvar]
#x = ds.variables['x'][:] * 1000.
#y = ds.variables['y'][:] * 1000.
tempData = temp[2, :]


1 loops, best of 1: 3.28 s per loop