maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git


  
1

	
2

	
3

	
4

	
5

	
6

	
7

	
8

	
9

	
10

	
11

	
12

	
13

	
14

	
            #!/usr/bin/env bash
# GMT_KNOWN_FAILURE
# Demonstrate issue #536 with labeling of periodic contours
# Turns out non-issue due to different contour lengths and -Gd
# Using -Gl shows annotation machinery works fine.
ps=periodic.ps
gmt grdmath -Rg -I1 32 22 SBAZ = t.nc
# Plot as non-periodic.  This plots labels.  The pile of blue contours around
# 180 occurs because the grid goes rapidly from -180 to +180 here and -Z+p is not used
# This is expected and correct.
gmt grdcontour t.nc -JG0/0/5.0i -A20+u@. -Bgaf -K -Wa1,blue -Gl-90/0/0/0,0/0/90/0 -P -Xc -Y0.3i > $ps
# Plot as periodic with -Zp.  Transition from -180 to +180 goes well for contouring
# but we have no contour labels anymore. This is the main issue here.
gmt grdcontour t.nc -J -A20+u@. -Wa1,blue -Bgaf -Gl-90/0/0/0,0/0/90/0 -Z+p -O -Y5.2i >> $ps