def _make_dataset(varname, seed=None, **var_kws):
rs = np.random.RandomState(seed)
_dims = {'time': 10, 'x': 5, 'y': 5}
_dim_keys = ('time', 'x', 'y')
ds = xr.Dataset()
ds['time'] = ('time', pd.date_range('2000-01-01', periods=_dims['time']))
ds['x'] = np.linspace(0, 10, _dims['x'])
ds['y'] = np.linspace(0, 10, _dims['y'])
data = rs.normal(size=tuple(_dims[d] for d in _dim_keys))
ds[varname] = (_dim_keys, data)
ds.coords['numbers'] = ('time',
np.array(range(_dims['time']), dtype='int64'))
return ds
评论列表
文章目录