Back to grdout examples

#! /bin/csh
gmtset VERBOSE TRUE LABEL_FONT_SIZE 12 ANOT_FONT_SIZE 12
set output = Togrd.ps
set B = a10f5/a1g.5SWne
set B2 = a10f5/a1g0.5swne
set R = 166400/177000/2.5/6.5
set R = 0/63/2.5/6.5
#-------------------------------------------------
# grd2xyz togrd.grd=1 >! togrd.xyz
# xyz2grd togrd.xyz -I1/0.012 -Gtogrd2.grd -R0/360/3/6
#-------------------------------------------------
# xyz2grd interpretes x values as deg thus 0-360 only
#-------------------------------------------------

set cpt = gray.cpt
#makecpt -Cgray -T-20/20/20 -Z -V -I >! $cpt

psbasemap -Ba10f5:"Distance (km)":/a1f0.5:"Two-Way Time (s)":SWne -R$R -JX10/-3 -K -X0.75 -Y2 -Vl >! Togrd.ps
grdimage  -B$B2 togrd.grd=1 -JX10/3 -R -C$cpt -O -K >> Togrd.ps
psxy -R -JX -O -N -G0 -Sv0.05c/0.24c/0.12c -K <>Togrd.ps
3.0  3     -50   0.25
31   3.0    30    0.3
10   3.95   -90  0.25
28   3.90   -90  0.25
48   3.8   -90  0.25
15   5.2   90 0.25  
35   5.2   90 0.25 
55   5.3   90 0.25 
END

pstext -R -JX -O -N -W255 -G255 -K <>Togrd.ps
1    3.15    4  0    0    1 seafloor multiple    seafloor   multiple side
30   2.95   4  0    0    3 side swipe   side swipe si
10   4.05   3  0    0    2 MohoMohoMoho
28   4.05   3  0    0    2 MohoMohoMoho
48   3.95   3  0    0    2 MohoMohoMoho
15   5.05   3  0    0    2 base of 2Abase of 2Abase of 2A
35   5.05  3  0    0    2  base of 2Abase of 2Abase of 2A
55   5.05   3  0    0    2  base of 2Abase of 2Abase of 2A
END


pstext -R -JX -O -N <>Togrd.ps
32.5 6.8   14  0    0    2 Galapagos Spreading Center (G-PRIME): ridge flank MCS profile
1    3.1   10  0    0    1 seafloor multiple
30   2.9   10  0    0    3 side swipe
10   4.0   10  0    0    2 Moho
28   4.0   10  0    0    2 Moho
48   3.9   10  0    0    2 Moho
15   5.0   10  0    0    2 base of 2A
35   5.0   10  0    0    2  base of 2A
55   5.0   10  0    0    2  base of 2A
END

imagetool Togrd.ps