docs: def loadnetcdf(filename, copy=True):
      """Load a dataarray from a NetCDF file.
  
      Args:
          filename (str): Filename (*.nc).
          copy (bool): If True, dataarray is copied in memory. Default is True.
  
      Returns:
          dataarray (xarray.DataArray): Loaded dataarray.
      """
      filename = str(Path(filename).expanduser())
  
      if copy:
          dataarray = xr.open_dataarray(filename).copy()
      else:
          dataarray = xr.open_dataarray(filename, chunks={})
  
      if dataarray.name is None:
          dataarray.name = filename.rstrip('.nc')
  
      for key, val in dataarray.coords.items():
          if val.dtype.kind == 'S':
              dataarray[key] = val.astype('U')
          elif val.dtype == np.int32:
              dataarray[key] = val.astype('i8')
  
      return dataarray