def __init__(self, variables=None, axes=None) :

if axes == None :

self.axes = Axes()

else :

self.axes = axes

if variables == None :

self.variables = {}

else :

self.variables = variables

def __getattr__(self, attributeName) :

if 'variables' in self.__dict__ :

if attributeName in self.variables :

return self.variables[attributeName]

if 'axes' in self.__dict__ :

return self.axes[attributeName]

raise AttributeError

def __getitem__(self, item) :

return getattr(self, item)

def close(self) :

pass

def write(self, filePath, compress = True) :

from netCDF4 import Dataset

# making sure the variable's axes are given

# to the file is the user's responsability

with Dataset(filePath, 'w') as output :

for axisName, axis in self.axes.items() :

axisName = Axes.ncStandardize(axisName)

output.createDimension(axisName, len(axis))

if axisName == 'time' :

output.createVariable('time', int, ('time',))

output.variables['time'].units = 'seconds since 1970-1-1'

from datetime import datetime

epoch = datetime(1970, 1, 1)

output.variables['time'][:] = [

(instant-epoch).total_seconds()

for instant in axis.data]

else :

# this is a real axis with information worth recording

if axis.units != 'indices' or len(axis) - 1 != axis[-1] :

output.createVariable(

axisName,

type(np.asscalar(axis.data.ravel()[0])),

(axisName,))

output.variables[axisName][:] = axis.data

if axis.units != None :

output.variables[axisName].units = axis.units

# TODO metadata...

for variableName, variable in self.variables.items() :

#variable.censor_nans(ratio = 1./3)

variable.censor_nans(ratio = 1)

if variableName == '~' :

variableName = 'unknown'

output.createVariable(

variableName,

type(np.asscalar(variable.data.ravel()[0])),

tuple(

[Axes.ncStandardize(axisName) for axisName in

variable.axes.keys()]),

zlib = compress)

if 'units' in variable.metadata :

output.variables[variableName].units = variable.units