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

	
15

	
16

	
17

	
18

	
19

	
20

	
21

	
22

	
23

	
            #!/usr/bin/env bash
# Testing gmt grdfilter if openmp is used.

ps=openmp.ps

if ! [[ ${HAVE_OPENMP} =~ TRUE|ON ]]; then
  # Since a PS-producing test is judge by comparison to original PS
  # and since no PS is produced without OpenMP, we simply duplicate
  # the original so they both exist and are identical.
  echo "OpenMP not available - just duplicating PS to pass test"
  cp "${GMT_SRCDIR:-.}"/$ps .
  exit 0
fi
FILT=g			# Gaussian filter
INC=1			# 1x1 degree output
D=1000			# 1000 km filter width

# Run gmt grdfilter as specified
gmt grdfilter -D4 -F${FILT}$D -I$INC @earth_relief_10m -Gt.nc -fg
gmt makecpt -Cglobe > t.cpt
gmt grdimage t.nc -JQ0/7i -Ba -BWSne+t"$D km Gaussian filter" -Ct.cpt -P -K -Xc -Y1.5i > $ps
gmt psscale -Ct.cpt -D3.5i/-0.5i+w6i/0.1i+h+jTC -O -K -Bxa -By+l"m" >> $ps
gmt grdimage @earth_relief_10m -JQ0/7i -Ba -BWSne+t"Original data" -Ct.cpt -O -Y4.75i >> $ps