---------------------------------------------------------------------------
KeyboardInterrupt Traceback (most recent call last)
<ipython-input-24-5b6192ad99b7> in <module>()
1 ## Step 2 - Geometric Adjustments
2 # We make the OC data apply to spatial resolution of SST, because the latter is lower.
----> 3 chla_da_low = ect_adjust_geometry(master=sst_da, slave=chla_da)
4 chla_da
/home/marcoz/Projects/BC/CCI-Tools/sandbox/notebooks/Marco/use_case_9.py in ect_adjust_geometry(master, slave)
36
37 # Help! This is soooo slow... few minutes on Norman's PC
---> 38 temp_da = slave.groupby('time').apply(resample)
39 temp_lon = scipy.ndimage.zoom(slave.lon, [lon_factor])
40 temp_lat = scipy.ndimage.zoom(slave.lat, [lat_factor])
/home/marcoz/Projects/Anaconda/miniconda3/lib/python3.5/site-packages/xarray/core/groupby.py in apply(self, func, shortcut, **kwargs)
381 grouped = self._iter_grouped()
382 applied = (maybe_wrap_array(arr, func(arr, **kwargs)) for arr in grouped)
--> 383 combined = self._concat(applied, shortcut=shortcut)
384 result = self._maybe_restore_empty_groups(combined)
385 return result
/home/marcoz/Projects/Anaconda/miniconda3/lib/python3.5/site-packages/xarray/core/groupby.py in _concat(self, applied, shortcut)
393 combined = self._concat_shortcut(applied, concat_dim, positions)
394 else:
--> 395 combined = concat(applied, concat_dim, positions=positions)
396
397 if isinstance(combined, type(self.obj)):
/home/marcoz/Projects/Anaconda/miniconda3/lib/python3.5/site-packages/xarray/core/combine.py in concat(objs, dim, data_vars, coords, compat, positions, indexers, mode, concat_over)
112 raise TypeError('can only concatenate xarray Dataset and DataArray '
113 'objects, got %s' % type(first_obj))
--> 114 return f(objs, dim, data_vars, coords, compat, positions)
115
116
/home/marcoz/Projects/Anaconda/miniconda3/lib/python3.5/site-packages/xarray/core/combine.py in _dataarray_concat(arrays, dim, data_vars, coords, compat, positions)
281 def _dataarray_concat(arrays, dim, data_vars, coords, compat,
282 positions):
--> 283 arrays = list(arrays)
284
285 if data_vars != 'all':
/home/marcoz/Projects/Anaconda/miniconda3/lib/python3.5/site-packages/xarray/core/groupby.py in <genexpr>(.0)
380 else:
381 grouped = self._iter_grouped()
--> 382 applied = (maybe_wrap_array(arr, func(arr, **kwargs)) for arr in grouped)
383 combined = self._concat(applied, shortcut=shortcut)
384 result = self._maybe_restore_empty_groups(combined)
/home/marcoz/Projects/BC/CCI-Tools/sandbox/notebooks/Marco/use_case_9.py in resample(x)
32
33 def resample(x):
---> 34 y = scipy.ndimage.zoom(x, [lat_factor, lon_factor])
35 return xr.DataArray(y)
36
/home/marcoz/Projects/Anaconda/miniconda3/lib/python3.5/site-packages/scipy/ndimage/interpolation.py in zoom(input, zoom, output, order, mode, cval, prefilter)
560 zoom = numpy.asarray(zoom, dtype=numpy.float64)
561 zoom = numpy.ascontiguousarray(zoom)
--> 562 _nd_image.zoom_shift(filtered, zoom, None, output, order, mode, cval)
563 return return_value
564
KeyboardInterrupt: