Alan D. Snow

You could try re-projecting to the South Pole Stereographic projection and then clipping.

I investigated and the latest version on pypi is 0.3.4 from 2018 with a vendored older version of pybind11. There are more recent releases: https://github.com/r-barnes/richdem/releases If you install from source...

> And how to install the more recent versions in the releases The readme has the instructions

To build richdem, I downloaded the assets from the releases here: https://github.com/r-barnes/richdem/releases Then compiled richdem and then installed the python wrapper.

In addition to the instructions, [this](https://github.com/conda-forge/richdem-feedstock/blob/e2d439f262abbe1054a2241ec7228ecc079d7700/recipe/build.sh#L18-L28) is a helpful reference. If compiling the code is problematic, I suggest using [conda](https://conda.io/en/latest/) with the [conda-forge](https://conda-forge.org/) channel to install richdem. ```bash conda config...

Note: You have different versions of GDAL on each system. Please make sure that you have the same versions of GDAL and other libraries when making the comparison.

I would like to add my vote `+1` to this idea.

It would be nice to have this as an extension to xarray. That way you could call something like ```python with xr.open_dataset('/path/to/file.nc') as xd: t = xd.wrf.getvar('temp') proj4 = xd.wrf.proj4...

Looks like it would not be bad to try to merge with salem. See: https://github.com/fmaussion/salem/issues/82

Should GDAL default to XSCALE=1 & YSCALE=1? Or, is there a better solution? This is what has been done for datacube (https://github.com/opendatacube/datacube-core/issues/1456).