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)