Register
Login
Resources
Docs Blog Datasets Glossary Case Studies Tutorials & Webinars
Product
Data Engine LLMs Platform Enterprise
Pricing Explore
Connect to our Discord channel

oblique.sh 630 B

You have to be logged in to leave a comment. Sign In
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
  1. #!/usr/bin/env bash
  2. #
  3. # Tests gmt project to make oblique lines
  4. ps=oblique.ps
  5. # Since GMT4 gmt project has no option for small circle we must do it differently
  6. gmt math -T0/360/1 45 = t.txt
  7. cat << EOF | gmt project -T30/60 -C0/0 -Fpq > p.txt
  8. 0 90
  9. 0 0
  10. EOF
  11. ppole=`$AWK '{if (NR == 1) printf "%s/%s\n", $1, $2}' p.txt`
  12. centr=`$AWK '{if (NR == 2) printf "%s/%s\n", $1, $2}' p.txt`
  13. gmt pscoast -Rg -JG30/50/7i -P -K -Glightgray -Dc -Bg > $ps
  14. echo 30 50 | gmt psxy -R -J -O -K -Sa0.2i -Gred -W0.25p >> $ps
  15. echo 0 0 | gmt psxy -R -J -O -K -Sc0.1i -Gblack >> $ps
  16. gmt project t.txt -T$ppole -C$centr -Fpq | gmt psxy -R -J -O -W3p >> $ps
Tip!

Press p or to see the previous file or, n or to see the next file

Comments

Loading...