code: def daily_temperature_range_variability(tasmax, tasmin, freq="YS"):
      r
  
      vdtr = abs((tasmax - tasmin).diff(dim='time'))
      out = vdtr.resample(time=freq).mean(dim='time')
      out.attrs['units'] = tasmax.units
      return out
docs: r"""Mean absolute day-to-day variation in daily temperature range.
  
      Mean absolute day-to-day variation in daily temperature range.
  
      Parameters
      ----------
      tasmax : xarray.DataArray
        Maximum daily temperature values [℃] or [K]
      tasmin : xarray.DataArray
        Minimum daily temperature values [℃] or [K]
      freq : str, optional
        Resampling frequency
  
      Returns
      -------
      xarray.DataArray
        The average day-to-day variation in daily temperature range for the given time period.
  
      Notes
      -----
      Let :math:`TX_{ij}` and :math:`TN_{ij}` be the daily maximum and minimum temperature at
      day :math:`i` of period :math:`j`. Then calculated is the absolute day-to-day differences in
      period :math:`j` is:
  
      .. math::
  
         vDTR_j = \frac{ \sum_{i=2}^{I} |(TX_{ij}-TN_{ij})-(TX_{i-1,j}-TN_{i-1,j})| }{I}