DEM'den kesit alma
From GMT Türkiye Wiki
N40E029C.grd isimli srtm grd'den kesit alma. Bu grd'nin elde edilmesi için bakınız.
project -C29.8/40.3 -E29.8/40.5 -G0.000855555 > ! profil.xy grdtrack profil.xy -GN40E029C.grd > ! profile.xydz set r = `minmax -C profile.xydz | awk ' {print $1,$2,$3,$4,$5,$6,$7-100,$8+100}'` #plot x axis as latitute awk '{print $2,$4}' profile.xydz| psxy -R$r[3]/$r[4]/$r[7]/$r[8] -JX10/2 -W -B0.1/250 > ! x.ps xv -rotate -90 x.ps #plot x axis as distance in degree awk '{print $3,$4}' profile.xydz| psxy -R$r[5]/$r[6]/$r[7]/$r[8] -JX10/2 -W -B0.05f.025/200f100 > ! x.ps xv -rotate -90 x.ps #plot x axis as distance in degree fill it with gray echo 0 0 0 0 > ! tmp echo 180 180 180 0 > ! tmp2 cat tmp profile.xydz tmp2 > ! tmp3 awk '{print $3,$4}' tmp3 | psxy -R$r[5]/$r[6]/$r[7]/$r[8] -JX10/2 -W -B0.05f.025/200f100 -G120 > ! x.ps xv -rotate -90 x.ps set name = marmara.ps set R = "-R29/30/40/41 -JM6" psbasemap $R -B0.2 -K -P> ! $name grdimage $R -Cgray.cpt N40E029C.grd -IN40E029CI.grd -O -K>> $name pscoast -Df -S120 -W -O $R -Ia/2/blue -K >> $name psxy profile.xydz $R -O -K -W8/red>> $name psbasemap -Lf29.8/40.1/29.8/40.1/20 $R -O -B.2 >> $name xv $name &