#!/bin/bash paleta="/opt/gmt/share/cpt/GMT_relief.cpt" grilla="/vox/201001/linux/topografia/ETOPO1_Bed_g_gdal.grd" ps="cile.ps" zona="-74/-71/-40/-30" grdcut $grilla -R$zona -G$ps.grd grdgradient $ps.grd -A50 -G$ps.grad -Nt grdview $ps.grd -JM2.5i -C$paleta -Qs -JZ1i -I$ps.grad -E230/45 -Wc -B1a2 -K > $ps rm *.grad *.grd
Explicación
Generamos una perspectiva 3D del borde continental de Chile, donde se aprecia la fosa de subducción y la topografá continental.
Formamos una sub-grilla de la zona a grafica con grdcut, luego añadimos iluminación con un acimut de 50 grados con el
comando grdgradient para finalmente graficar la perspectiva con el comando grdview donde los parámetros
funcionan del mismo modo como el ejemplo 2.