$in1 = "r-1000.dat"; 
$in2 = "r-100.dat"; 
$in3 = "r-10.dat"; 
$in4 = "r-1.dat"; 
$in5 = "r-0.1.dat"; 
$in6 = "r-0.01.dat"; 
$in7 = "r-trace.dat"; 
$in8 = "r.txt"; 
$out = "ruapehu_dotmap.eps";
$loc = "ruapehu_loc.xyz";
     
#ruapehu_subset.dat: N = 27	<2730400/2744200>	<6210700/6240500>	<0.002/3087.3>                                                         
$west = 2710000;
$east = 2870000;
$south = 6210000;
$north = 6370000;

#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 $in7 -R$west/$east/$south/$north -JX4i -W0.75p -Sc0.05i -G0 -Bag20000wesn -V -K -P > $out";
system "gmt psxy $in6 -R -JX -W0.75p -Sc0.075i -G0/114/178 -V -K -O -P >> $out";
system "gmt psxy $in5 -R -JX -W1p -Sc0.1i -G230/159/0 -V -K -O -P >> $out";
system "gmt psxy $in4 -R -JX -W1p -Sc0.125i -G204/121/167 -V -K -O -P >> $out";
system "gmt psxy $in3 -R -JX -W1p -Sc0.15i -G0/158/115 -V -K -O -P >> $out";
system "gmt psxy $in2 -R -JX -W1p -Sc0.175i -G213/94/0 -V -K -O -P >> $out";
#system "gmt psxy $in1 -R -JX6i -W1p -Sc0.1i -G0 -V -K -O >> $out";
system "gmt psxy $loc -R -JX -W1p -St0.3i -V -O -K -P >> $out";
system "gmt pstext $in8 -R -JX -D0.25i/0i -F+f10p -V -O -P >> $out";
system "gmt psconvert $out -A+m0.5i/0.5i -Tg";