#!/bin/bash cd SBAS #gmt makecpt -Celevation -T-100/400/10 -Z -D > new.cpt #gmt makecpt -Cterra -T200/500/10 -Z -D > new.cpt gmt makecpt -Cpolar -T-10/10/1 -Z -D > new.cpt # Variables format="jpg" file="vel_ll_nan2.grd" output="velocity3" gmt begin $output $format gmt set FONT_ANNOT 10p,Times-Roman gmt set FONT_LABEL 10p,Times-Roman gmt grdimage $file -JM15c -I+d -B -BWSne -Cnew.cpt gmt coast -B -W1p,black #-Scornflowerblue #-Gbrown #gmt colorbar -DjTL+w5c+o.7c/1c+h+ml+e -Bxaf+l"LOS Displacement" -By+l"(rad)" #-Cnew.cpt gmt colorbar -DJTC+w5c+o.7c/1c+h+ml -Bxaf+l"Mean velocity (mm/yr)" -By -Cnew.cpt gmt end show