code: def max_n_day_precipitation_amount(pr, window=1, freq='YS'):
      r
  
      # rolling sum of the values
      arr = pr.rolling(time=window, center=False).sum()
      out = arr.resample(time=freq).max(dim='time', keep_attrs=True)
  
      out.attrs['units'] = pr.units
      # Adjust values and units to make sure they are daily
      return utils.pint_multiply(out, 1 * units.day, 'mm')
docs: r"""Highest precipitation amount cumulated over a n-day moving window.
  
      Calculate the n-day rolling sum of the original daily total precipitation series
      and determine the maximum value over each period.
  
      Parameters
      ----------
      da : xarray.DataArray
        Daily precipitation values [Kg m-2 s-1] or [mm]
      window : int
        Window size in days.
      freq : str, optional
        Resampling frequency : default 'YS' (yearly)
  
      Returns
      -------
      xarray.DataArray
        The highest cumulated n-day precipitation value at the given time frequency.
  
      Examples
      --------
      The following would compute for each grid cell of file `pr.day.nc` the highest 5-day total precipitation
      at an annual frequency:
  
      >>> da = xr.open_dataset('pr.day.nc').pr
      >>> window = 5
      >>> output = max_n_day_precipitation_amount(da, window, freq="YS")