REM draw seismic spatial distribution map                for windows
psbasemap -R75/160/-10/50 -JM7.8i -BNEWSa10:xx:/a10:yy: -K -G >p.ps
pscoast -R -JM -S100/149/237 -G240/248/255 -W1 -Di -N1/1P -N2/0.5P -V -O -K>>p.ps 
gawk "{print $1,$2,0.02*$3}" p.txt |psxy -R -JM -O -K -Sci -G255/255/0 -W0.45p>>p.ps    
echo 118.25 25.06 > pratt.d
gawk "{print $1, $2, 0, 1000, 1000}" pratt.d |psxy -R -J -O -K -SE -Wthinnest/red >>p.ps
echo 118.25 25.06 |psxy -R -JM -O -K -Sa0.25i -Gred -W0.45p>>p.ps

 

 

REM draw topography map

makecpt -Crainbow -T-10000/10000/400 -V -Z > sc.cpt

grdgradient etopo2.grd -Ne0.8 -A100 -M -Gsc_i.grd 
grdimage -JM6i -R65/160/-10/50 -B10/10 -Csc.cpt etopo2.grd -Isc_i.grd -K -P -V>topo.ps