Skip to content

grdselect -Ai unable to report intersection when one grid is global? #8894

@mwoot

Description

@mwoot

Description

grdselect -Ai returns an unexpected error when trying to find the intersection between a global grid and another grid. Issue comes from my forum post.

Full script that generated the error

gmt grdcut @earth_geoid_01d_p -Rg -Gglobal.grd
gmt grdcut global.grd -R0/1/0/1 -Gsubset.grd
gmt grdselect -Ai global.grd subset.grd

Full error message

grdselect [WARNING]: No common area for all grids is possible.

Expected outcome

-R0/1/0/1

Running the intersection between two non-global grids works as expected:

gmt grdcut global.grd -R0/15/0/15 -Gbigger_subset.grd
gmt grdselect -Ai subset.grd bigger_subset.grd

System information

  • Operating system: Ubuntu 22.04.5 LTS
  • GMT version (gmt --version): 6.6.0

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions