def from_netcdf(cls, path):
"""
Construct a fuel moisture model from data stored in a netCDF file.
:param path: the path to the netCDF4 file
"""
import netCDF4
print "reading from netCDF file", path
d = netCDF4.Dataset(path)
ncfmc = d.variables['FMC_GC'][:,:,:]
print "FuelMoistureModel.from_netcdf: reading FMC_GC shape",ncfmc.shape
d0, d1, k = ncfmc.shape
P = d.variables['FMC_COV'][:,:,:,:]
Tk = np.array([1.0, 10.0, 100.0]) * 3600
fm = FuelMoistureModel(ncfmc[:,:,:k-2], Tk)
fm.m_ext[:,:,k-2:] = ncfmc[:,:,k-2:]
fm.P[:,:,:,:] = P
return fm
评论列表
文章目录