#!/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 is the big map: #Here are some variables to save some trouble set region = '-R74/130/19/45' set filename = 'atfall' set projection = '-JM6.8' set offsets = '-X1 -Y6' ###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 -Ba10 -K $offsets -P> $filename.ps ###Topography echo 'plotting topo' grdimage \ /tmp/atf.grd \ -C/disk2/opt/igmt_1.0/colormaps/topo.cpt \ -I/tmp/atf.int \ -V $region $projection -O -K >>$filename.ps ##Coastlines and whatnot echo 'plotting coastlines and whatnot' pscoast $region $projection -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 $region $projection -Sa -G0/0/0 -W0.2/0/0/0 -L -O -K \ >>$filename.ps #our search parameters are: