issue_comments: 575605008

This ```python def quantile(self, q, dim=None, interpolation="linear", keep_attrs=None): from xarray.core.computation import apply_ufunc from xarray.core import utils

scalar = utils.is_scalar(q) 
q = np.atleast_1d(np.asarray(q, dtype=np.float64)) 
if dim is None: 
    dim = self.dims

if utils.is_scalar(dim): 
    dim = [dim]

def _wrapper(npa, **kwargs): 
    # move quantile axis to end. required for apply_ufunc 
    return np.moveaxis(np.nanpercentile(npa, **kwargs), 0, -1)

axis = np.arange(-1, -1 * len(dim) - 1, -1) 
result = apply_ufunc( 
    _wrapper, 
    self, 
    input_core_dims=[dim], 
    exclude_dims=set(dim), 
    output_core_dims=[["quantile"]], 
    output_dtypes=[np.float64], 
    output_sizes={"quantile": len(q)}, 
    dask="parallelized", 
    kwargs={"q": q * 100, "axis": axis, "interpolation": interpolation}, 
)

# for backward compatibility 
result = result.transpose("quantile", ...) 
if scalar: 
    result = result.squeeze("quantile") 
if keep_attrs: 
    result.attrs = self._attrs

return result

xr.Variable.quantile = quantile `` seems to makequantile` work, but it might be a bad idea to monkeypatch. Thoughts, @dcherian?