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

	
24

	
25

	
26

	
27

	
28

	
29

	
30

	
31

	
32

	
33

	
34

	
35

	
36

	
37

	
38

	
39

	
            #!/usr/bin/env bash
# Show what (x,y,p,q,r,s) in project are for man page
gmt begin project_setup ps
	a=$(gmt math -Q 4 2.5 DIV ATAN R2D =)
	X=(`echo 0 0 | gmt project -C0/-1 -E2.5/3 -N`)
	gmt plot -JX4i/0 -R-3.5/4/-2.7/2.6 -Glightgray -B -Bx+l"@%7%x@%% or @%7%r@%%" -By+l"@%7%y@%% or @%7%s@%%" <<- EOF
	-1.5	-1.0625
	0	-2
	2	1.2
	0.5	2.1375
	EOF
	gmt plot -Sc0.3c -Gorange -Bg10 <<- EOF
	0	-1
	2	2.2
	EOF
	gmt plot -Sv16p+e+h1+s -W2p -Gblack <<- EOF
	0	-1	2	2.2
	0	-1	-2.5	0.5625
	EOF
	echo 0 0 | gmt plot -Sc0.3c -Gred
	# Get coordinates of the (0,q) point as well so we can dash the line
	x=$(gmt math -Q ${X[3]} $a SIND MUL NEG =)
	y=$(gmt math -Q ${X[3]} $a COSD MUL 1 SUB =)
	echo ${X[4]} ${X[5]} | gmt plot -Sc0.2c -Gblue
	gmt text -F+f12p,Times-Italic+a+j -Dj0.15c <<- EOF
	0		-1			0		TL @%7%C@%%
	2		2.2			0		BR @%7%E@%%
	1.9		1.9			${a}	TC	p
	-2.3	0.4			${a}	RM	q
	2		1.2			${a}	TC	L@-max@-
	0		-2			${a}	TC	L@-min@-
	0		-2			${a}	RB	W@-min@-
	-1.5	-1.0625		${a}	RM	W@-max@-
	0.45	0.8			-16		TC	@~a@~
	EOF
	echo 0 0 ${X[4]} ${X[5]} | gmt plot -Sv12p+s -W0.25p,red,-
	echo 0 0 ${x} ${y} | gmt plot -Sv12p+s -W0.25p,red,-
	echo 0 -1 1i ${a} 90 | gmt plot -Sm9p+b -W0.5p -Gblack
gmt end show