1#!/usr/bin/env bash 2# GMT EXAMPLE 42 3# 4# Purpose: Illustrate Antarctica and stereographic projection 5# GMT modules: makecpt, grdimage, coast, legend, colorbar, set, text, plot 6# Unix progs: [curl grdconvert] 7# 8gmt begin ex42 9 gmt set PROJ_ELLIPSOID WGS-84 10 # Data obtained via website and converted to netCDF thus: 11 # curl http://www.antarctica.ac.uk//bas_research/data/access/bedmap/download/bedelev.asc.gz 12 # gunzip bedelev.asc.gz 13 # grdconvert bedelev.asc BEDMAP_elevation.nc=ns -V 14 gmt makecpt -Cearth -T-7000/4000 15 gmt grdimage @BEDMAP_elevation.nc -Jx1:60000000 -Q 16 gmt coast -R-180/180/-90/-60 -Js0/-90/-71/1:60000000 -Bafg -Di -W0.25p 17 gmt colorbar -DJRM+w6.5c/0.5c+o1c/0+mc -F+p+i -Bxa1000+lELEVATION -By+lm 18 # GSHHG 19 gmt coast -Glightblue -Sroyalblue2 -X5c -Y12c 20 gmt coast -Glightbrown -A+ag -Bafg 21 gmt legend -DjLM+w3.5c+jRM+o1c/0 -F+p+i <<- EOF 22 H 12 Legend 23 D 0.25c 1p 24 S 0.4c s 0.3c blue 0.25p 0.75c Ocean 25 S 0.4c s 0.3c lightblue 0.25p 0.75c Ice front 26 S 0.4c s 0.3c lightbrown 0.25p 0.75c Grounding line 27 EOF 28 # Fancy line 29 gmt plot -R0/19/0/25 -Jx1c -B0 -W2p -X-6c -Y-13.5c <<- EOF 30 0 14 31 6.5 14 32 13 11.5 33 19 11.5 34 EOF 35 gmt text -F+f18p+jBL -Dj8p/0 <<- EOF 36 0 13 BEDMAP 37 0 24 GSHHG 38 EOF 39gmt end show 40 41