$in = "simple.out"; 
$out = "simple.ps";

     
#simple.out: N = 100	<-20/25>	<-20/25>	<-150/361>
                                                       
$west = -25;
$east = 30;
$south = -25;
$north = 30;

#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 -B10:Easting:/10:Northing:WsNe -V -K -P > $out";
system "gmt pstext $in -R -JX6i -D0.1i/0.1i -F+f10p -V -O -P >> $out";
system "gmt psconvert $out -A+m1i -Tg";
