Contour maps
# Purpose: Make two contour maps based on the data in the file osu91a1f_16.nc
# GMT progs: gmtset, grdcontour, basemap, coast
gmtset(MAP_GRID_CROSS_SIZE_PRIMARY=0, FONT_ANNOT_PRIMARY=10, PS_CHAR_ENCODING="Standard+")
basemap(limits=(0,16.5,0,19.0), frame=0, figscale=1)
coast!(limits=:global360, proj=(name=:Hammer, center=0), frame=(ticks=30,),
resolution=:crude, land=:lightbrown, water=:lightblue, figsize=15,
x_off=0.64, y_off=0.5)
grdcontour!("@osu91a1f_16.nc", cont=10, annot=(int=50,labels=(font=7,)), labels=(dist=10,),
range=(-1000,-1), pen=((contour=1,pen="thinnest,-"), (annot=1, pen="thin,-")),
ticks=(gap=(0.25,0.05),))
grdcontour!("@osu91a1f_16.nc", cont=10, annot=(int=50,labels=(font=7,)), labels=(dist=10,),
range=(-1,1000), ticks=(gap=(0.25,0.05),))
coast!(limts=:d, proj=:Hammer, frame=(title="Low Order Geoid", ticks=30), res=:crude,
land=:lightbrown, water=:lightblue, y_off=8.5, figsize=15)
grdcontour!("@osu91a1f_16.nc", cont=10, annot=(int=50,labels=(font=7,)), labels=(dist=10,),
range=(-1000,-1), pen=((contour=1,pen="thinnest,-"), (annot=1, pen="thin,-")),
ticks=(gap=(0.25,0.05),labels=""))
grdcontour!("@osu91a1f_16.nc", cont=10, annot=(int=50,labels=(font=7,)), labels=(dist=10,),
range=(-1,100), ticks=(gap=(0.25,0.05),labels=""), show=true)
See also GMT ex01