Dear all, I am using GMT5 and I am trying to get on the same figure LAB depth and surface topography. But I did not understand what I misunderstood or miss to get this figure. Any comments or suggestions.
Cheers
/bin/rm tmp.ps tmp.cpt
gmt makecpt -Crainbow -T40/300/50 -Z >tmp.cpt
Make a basemap of central and southern Africa:
gmt psbasemap -R-1.5/50.5/-35.5/11.5 -JM14 -Ba2f1 -P -X2.0 -Y15 -K >tmp.ps
Convert the xyz file into GRD file using xyz2grd
awk ‘{print $1, $2, $3}’ LAB.xyz | gmt xyz2grd -R -I0.5 -Gtmp.grd
#gmt grdimage tmp.grd -JM -Ctmp.cpt -E50 -K -O >>tmp.ps
Resample tmp.grd from coarse (0.5 deg.) into finer resolution (0.1 degree) using grdsample.
The output file is tmph.grd.
gmt grdsample tmp.grd -Gtmph.grd -I0.1
Cut out the portion of elevation data from a global database
The output is topo0.grd
gmt grdcut etopo1_bedrock_africa.grd -Gtopo0.grd -R
Resample it into the same interval as tmph.grd
grdsample topo0.grd -Gtopoh.grd -I0.1
Produce the gradient file. This time I usd -Nt to make a more dramatic shade
gmt grdgradient topoh.grd -Gtopoh.int -A45 -Nt
Produce the image using the resampled grd file and tmp.int
gmt grdimage tmph.grd -Itopoh.int -Ctmp.cpt -R -J -O -K >>tmp.ps
Overlay with coastlines and political boudnaries
gmt pscoast -R -JM -B2f2 -A1000 -W1 -N1 -D -P -V -K -O >> tmp.ps
Plot the scale bar at the bottom of the plot
gmt psscale -Ctmp.cpt -Dx8c/1c+w6c/0.5c+jTC+h -Bxaf+l"LAB " -By+lkm -K -O -Y-2.5 >> tmp.ps