Skip to content

van der Grinten projection (-JV) fails in non -Rg domains #7282

@KristofKoch

Description

@KristofKoch

Description of the problem

When plotting a global map with the van der Grinten projection (-JV) everything behaves as expected, when using it to plot quadrants or hemispheres things go wild: Instead of quadrants or hemispheres I get modern art.

Full script that generated the error

gmt begin van_der_Grinten_quadrants_hemispheres png
  gmt subplot begin 2x4 -Fs8c/16c -M8c/1c -T"van der Grinten quadrant & hemisphere maps"
    while read minX maxX minY maxY title
    do
    gmt subplot set
      gmt coast -R${minX}/${maxX}/${minY}/${maxY} -JV0/8c -Bg -Dc -Glightgray -Scornsilk -A10000 -Wthinnest -B+t"${title}"
    done << EOF
-180   0   0 90 NW
   0 180   0 90 NE
   0 180 -90  0 SE
-180   0 -90  0 SW
-180   0 -90 90 W
-180 180   0 90 N
   0 180 -90 90 E
-180 180 -90  0 S
EOF
  gmt subplot end
gmt end show

If you prefer a one-liner use the following for a plot of the NW quadrant:

gmt coast -R-180/0/0/90 -JV0/10c -Bg -Dc -Glightgray -Scornsilk -A10000 -Wthinnest -view van_der_Grinten_NW_quadrant

Full error message

No error messages generated

Actual outcome

van_der_Grinten_quadrants_hemispheres

Expected outcome

Four quadrant maps in the first row and four hemisphere maps in the second row in above plot.

System information

  • Operating system: MacOS 13.2.1
  • GMT version (gmt --version): 6.5.0_51b8cba_2023.02.25 also fails with 6.4.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