#!/bin/csh ###THis is based on my Pamirs plots, but is designed to do the ###Altyn Tagh Fault ###This one is designed to overlay the seismicity on the topo ###Hopefully this will do the focalmechs as well #Here are some variables to save some trouble set region = '-R74/130/19/45' set filename = 'instrumentalhistoricseismicityandtopo' set projection = '-JM8' ###These need to be done only the first time (or whenever you change the #range of the grd. ##First make the grd from the gtopo30 data set: #here is the location: #http://edcwww.cr.usgs.gov/landdaac/gtopo30/gtopo30.html #note that you need to be sure that the grdraster.info file is #set up correctly. On alai it is in /disk2/opt/GMT3.1/lib/grdraster.info #these are the western grids #echo 'making the southern grd from the raster' #grdraster 1 -R70/100/19/40 -G/tmp/atfS.grd -I2.5m -V #echo 'making the northern grd from the raster' #grdraster 3 -R70/100/40/45 -G/tmp/atfN.grd -I2.5m -V #echo 'making the northeast raster' #grdraster 6 -R100/130/40/45 -G/tmp/atfNE.grd -I2.5m -V #echo 'making the southeast raster' #grdraster 5 -R100/130/19/40 -G/tmp/atfSE.grd -I2.5m -V ##Then paste them together: #echo 'pasting the grds' #grdpaste /tmp/atfS.grd /tmp/atfN.grd \ #-G/tmp/atfW.grd -V #grdpaste /tmp/atfSE.grd /tmp/atfNE.grd \ #-G/tmp/atfE.grd -V #grdpaste /tmp/atfW.grd /tmp/atfE.grd -G/tmp/atf.grd -V ##Then make a -S10c -ent (intensity or illumination) file: #echo 'making the intensity map' #grdgradient /tmp/atf.grd -G/tmp/atf.int -A315 -V -Ne0.6 /bin/rm $filename.ps ###Basemap echo 'plotting basemap' psbasemap $region $projection -Ba5 -K > $filename.ps ###Topography echo 'plotting topo' grdimage \ /tmp/atf.grd \ -C/disk2/opt/igmt_1.0/colormaps/topo.cpt \ -I/tmp/atf.int \ -V -R -JM -O -K -U$filename>>$filename.ps ##Coastlines and whatnot echo 'plotting coastlines and whatnot' pscoast -R -JM -Di -N1 -I1/1/0/0/255 \ -S0/0/255 -O -K -V \ >> $filename.ps ##historic seismicity echo 'plotting historic seismicity with stars' nawk '{print $4, $3, 0.04*($5-4)}' \ histseis.txt | \ psxy -R -JM -Sa -G0/0/0 -W0.2/0/0/0 -L -O -K >>$filename.ps #our search parameters are: