pygmt.grdpaste: Join two grids along their common edge

[Chuan1937]

pygmt.grdpaste: Join two grids along their common edge

Preview: https://pygmt-dev--4399.org.readthedocs.build/en/4399/api/generated/pygmt.grdpaste.html#pygmt.grdpaste

[seisman]

Also need to add it to doc/api/index.rst.

[Chuan1937]

Could you also please check the min and max values in the result?
Must use temporary files for xarray.DataArray inputs because virtualfile_in returns all zeros when handling grdpaste.

[seisman]

Could you also please check the min and max values in the result? Must use temporary files for xarray.DataArray inputs because virtualfile_in returns all zeros when handling grdpaste.

I think this points to an upstream GMT bug, which should be fixed.

[Chuan1937]

Due to a bug fix in the upstream GMT GenericMappingTools/gmt#8901, I was able to run test_grdpaste successfully after recompiling GMT locally. The test failure here likely also requires the latest GMT.

@seisman

[seisman]

Could you please mark these tests as xfail for GMT older than 6.6 and also document that this method need GMT dev version

[seisman]

There are 4 cases: file in file out, file in xarray out, xarray in file out, xarray in xarray out. At least the first case works with gmt 6.6, right?

[Chuan1937]

There are 4 cases: file in file out, file in xarray out, xarray in file out, xarray in xarray out. At least the first case works with gmt 6.6, right?

Yes, But test_grdpaste only include case 3, 4. Do I need to add case1,2?

[seisman]

There are 4 cases: file in file out, file in xarray out, xarray in file out, xarray in xarray out. At least the first case works with gmt 6.6, right?

Yes, But test_grdpaste only include case 3, 4. Do I need to add case1,2?

Yes please