code: def xrmerge(das, accept_new=True):
      
      da = das[0]
      for new_da in das[1:]:
          # Expand both to have same dimensions, padding with NaN
          da, new_da = xr.align(da, new_da, join='outer')
          # Fill NaNs one way or the other re. accept_new
          da = new_da.fillna(da) if accept_new else da.fillna(new_da)
      return da
docs: Merges xarrays with different dimension sets
      Parameters
      ----------
      das : list of data_arrays
  
      accept_new
  
      Returns
      -------
      da : an xarray that is the merge of das
  
      References
      ----------
      Thanks to @jcmgray https://github.com/pydata/xarray/issues/742#issue-130753818
  
      In the future, we may not need this as xarray may provide the merge for us.