(28) Mixing UTM and geographic data sets

Like the [(22) World-wide seismicity of last month] case, the GMT.jl interface offers a much handier way of achieving what otherwise had to be the translation of Example 28

using GMT
resetGMT() # hide

G = gmtread("@Kilauea.utm.nc");
C = makecpt(cmap="copper", range=(0,1500));
imshow(G, cmap=C, shade=true, frame=(axes="WS", annot=true),
       coast=(shore=true, ocean=:lightblue, frame=(axes="EN", annot=true, grid=true)), show=true)