if exists(sdilesfolder, filename):
sdiles.update({v.name: Dataset(sdilesfolder+filename, mode='a',clobber = False)})
else:
d = [d for d in basenc.variables.values() if d.name in v.dimensions and d.name not in ['lat','lon']]
sdiles.update({v.name:self.initSdile(dile, v, d, sdilesfolder, filename)})
return sdiles
if __name__ == '__main__':
kg = KellyGrinder()
kc = KellyConverter()
timer = Chrono()
baselink = "http://users.rcc.uchicago.edu/~davidkelly999/agmerra.origgrid.2deg.tile"
basefilelink = baselink + "/0004/clim_0004_0047.tile.nc4"
sdilesfolder = "results/sdiles/"
# basefile setup
print "downloading base file..."
timer.start()
basefile = kg.downloadFile(basefilelink, "results/base/")
basenc = ncOpen(basefile)
basefolder = "results/"
timer.stop()
y += 1 x += 1 return iter_len def onClose(self): self.client.close() if __name__ == '__main__': ng = NetcdfGeometry() df = DileFactory() timer = Chrono() path = "/sdiles/ubuntu/sdiles/sdile_tasmax_2_0_1.nc" fname = pathLeaf(path,False) print "computing md5 for ", fname, "..." timer.start() md5 = getMD5(path) timer.stop() print "md5 computed in: ", timer.formatted() rgrp = ncOpen(path, mode='r') bb = ng.getBoundingBox(rgrp['lat'],rgrp['lon']) zoom = ng.getZoomLevel(rgrp['lat'], rgrp['lon'])
