def read_base(filename): f=mygis.Dataset(filename,mode="r") output=Bunch() variables=Bunch() for k in f.variables.keys(): curvar=f.variables[k] if (k!="pr") and (k!="tasmin") and (k!="tasmax") and (k!="TREFMN") and (k!="TREFMX") and (k!="tas"): data=Bunch(data=curvar[:]) else: data=Bunch() data.atts=copy.deepcopy(curvar.__dict__) # data.dims=copy.deepcopy(curvar.dimensions) data.dims=curvar.dimensions data.dtype=curvar.dtype variables[k]=data output.variables=variables output.global_attributes=copy.deepcopy(f.__dict__) # output.dimensions=copy.deepcopy(f.dimensions) output.dimensions=f.dimensions output.ncfile=f # f.close() return output