def init_files(self,period,extent): # if self._init: v = mt.Variable(self.nc_var,'mm') cs = mt.gen_coordset(period,extent) mvar = mt.MappedVariable(v,cs,self.dtype) #in_dtype) self.fm = SplitFileManager(self.path,mvar) clobber = self._file_mode == 'w' self.fm.create_files(self.schema,False,clobber,chunksize=self.chunksizes,ncparams=self.ncparams)#,chunksize=(1,16,16)) #None) #(64,32,32)) #(1,16,16)) self._open = True
def __setitem__(self, k, extent): group = self.f.create_group(k) g = DatasetManager(group) lats = mt.Coordinates(mt.latitude, extent.latitudes.to_degrees()) lons = mt.Coordinates(mt.longitude, extent.longitudes.to_degrees()) cs = mt.CoordinateSet([lats, lons]) v = mt.Variable( 'area', 'm2', dict(long_name='Metres squared covered by cell', cell_size=extent.cell_size.to_degrees(), lat_orient=extent.parent_ref.lat_orient)) mapped_var = mt.MappedVariable(v, cs, np.float64) ncv = g.create_variable(mapped_var) ncv[...] = extent.areas.data
def init_files(self, time_coords): #period,extent): # if self._init: v = mt.Variable(self.nc_var, 'mm') mvar = mt.MappedVariable(v, (time_coords, self.coords[0], self.coords[1]), np.float32) #in_dtype) self.fm = self.file_manager_type(self.path, mvar) clobber = self._file_mode == 'w' self.fm.create_files(False, clobber, chunksize=(1, 16, 16)) #None) #(64,32,32)) #(1,16,16)) self._open = True
def setup_var_coords(): global m_tvar global extent global test_path georef = geo.GeoReference((0,0),1000,1000,0.05) extent = extents.Extent(georef).ioffset[0:10,0:10] period = dt.dates('dec 2000 - jan 25 2001') tvar = mt.Variable('test_var','mm') m_tvar = mt.MappedVariable(tvar,mt.gen_coordset(period,extent),np.float32) test_path = os.path.join(os.path.dirname(__file__),'file_tests') shutil.rmtree(test_path,True)