diff --git a/xwavelet/wavelet.py b/xwavelet/wavelet.py index ee345a6..8a1790c 100644 --- a/xwavelet/wavelet.py +++ b/xwavelet/wavelet.py @@ -154,11 +154,6 @@ def dwvlt(da, s, spacing_tol=1e-3, dim=None, xo=50e3, a=1.0, ntheta=16, wtype="m else: raise NotImplementedError("Only the Morlet wavelet is implemented for now.") - new_coords = dict() - for d in dim: - new_coords[d] = da[d] - wavelet = xr.DataArray(wavelet.data, dims=dim, coords=new_coords) - dawt = (da * np.conj(wavelet)).sum(dim, skipna=True) * np.prod(delta_x) / s dawt = dawt.drop_vars(sdim) dawt[sdim] = xo * s