Edad del fondo oceánico

GMT progs : psscale, psxy

#!/bin/bash

PS=grafo.ps
AGEGRID=globalage_1.5.grd
AGECPT=age.cpt
AGEFILE=isochrons.xy
SCALING=0.10c
PLATES=nuvel_plates_mod
WESN=-80/20/-60/60

psbasemap -R$WESN -Jm$SCALING -Ba20/a20SneW -K -V >$PS
grdimage $AGEGRID -Jm -C$AGECPT -R$WESN -B -K -O -V >>$PS
pscoast -R -Jm -B -C10 -Dl -G0 -K -O -V >>$PS
#psxy $AGEFILE -Jm  -R  -: -M -W0.2p/0/0/0 -O -K -V >>$PS
psxy $PLATES -Jm -R -: -M -H3 -W1.5p/0/0/0ta -O -V -K >>$PS
psscale -C$AGECPT -D12/7.5/15/0.5 -L  -N300 -V -O >>$PS

Explicación
Este script grafica la edad del suelo oceánico. En lugar de la topografía global se usa el archivo globalage_1.5.grd que contiene la edad del suelo oceánico en millones de años asignando el valor correspondiente a cada célula de la grilla (x,y). El archivo nuvel_plates_mod contiene los bordes de las principales placas tectónicas.