$in = "ruapehu_subset.dat"; 
$out = "ruapehu_subset.eps";
$loc = "ruapehu_loc.xyz";
     
#ruapehu_subset.dat: N = 27	<2730400/2744200>	<6210700/6240500>	<0.002/3087.3>                                                         
$west = 2720000;
$east = 2740000;
$south = 6210000;
$north = 6230000;

#system "surface $in -Gaso.grd -R$west/$east/$south/$north -I500 -V";
#system "makecpt -Chaxby -T-60/0/2 -Z -V > aso.cpt";
#system "grdimage aso.grd -R -X1i -Y1i -JX6i -E150 -Caso.cpt -K -V > $out";
#system "grdcontour aso.grd -Caso.cpt -JX6i -R -A- -W.5p,0 -K -O >> $out";
#system "psscale -D6.8i/3i/3i/.1i -Caso.cpt -B20:'':/:mGal: -O -K -V >> $out";
system "gmt psxy $in -R$west/$east/$south/$north -JX6i -W1p -Sc0.1i -G255/0/0 -B4000g2000 -V -K > $out";
system "gmt pstext $in -R -JX6i -D0.1i/0.1i -F+f10p -V -K -O >> $out";
system "gmt psxy $loc -R -JX6i -W1p -St0.1i -G0 -V -O >> $out";