docs: def Field_to_xarray(field):
      '''Convert a climlab.Field object to xarray.DataArray'''
      dom = field.domain
      dims = []; dimlist = []; coords = {};
      for axname in dom.axes:
          dimlist.append(axname)
          try:
              assert field.interfaces[dom.axis_index[axname]]
              bounds_name = axname + '_bounds'
              dims.append(bounds_name)
              coords[bounds_name] = dom.axes[axname].bounds
          except:
              dims.append(axname)
              coords[axname] = dom.axes[axname].points
      #  Might need to reorder the data
      da = DataArray(field.transpose([dom.axis_index[name] for name in dimlist]),
                        dims=dims, coords=coords)
      for name in dims:
          try:
              da[name].attrs['units'] = dom.axes[name].units
          except:
              pass
      return da