@@ -60,14 +60,14 @@ ds_fields = parcels.tutorial.open_dataset("CopernicusMarine_data_for_Argo_tutori
6060ds_fields.load() # load the dataset into memory
6161
6262# Create an idealised wind field and add it to the dataset
63- tdim, ydim, xdim = ( len(ds_fields.time),len(ds_fields. latitude), len(ds_fields.longitude) )
63+ ydim, xdim = len(ds_fields.latitude), len(ds_fields.longitude)
6464ds_fields["UWind"] = xr.DataArray(
65- data=0.5 * np.ones((tdim, ydim, xdim)) * np.sin(ds_fields.latitude.values - ds_fields.latitude.values.mean())[None, :, None],
66- coords=[ds_fields.time, ds_fields. latitude, ds_fields.longitude])
65+ data=0.5 * np.ones((ydim, xdim)) * np.sin(ds_fields.latitude.values - ds_fields.latitude.values.mean())[:, None],
66+ coords=[ds_fields.latitude, ds_fields.longitude])
6767
6868ds_fields["VWind"] = xr.DataArray(
69- data=np.zeros((tdim, ydim, xdim)),
70- coords=[ds_fields.time, ds_fields. latitude, ds_fields.longitude])
69+ data=np.zeros((ydim, xdim)),
70+ coords=[ds_fields.latitude, ds_fields.longitude])
7171
7272fields = {
7373 "U": ds_fields["uo"],
0 commit comments