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'])