#!/bin/csh gmtset BASEMAP_AXES = WeSn set z=1.0 /bin/rm tmp.ps ### Left panel. Both circles pswiggle out.dat -R-200/200/00/50 -JX1.8/-5.0 \ -Ba90f30/a5f1 -G200/50/50 \ -C-0.00 -M -W3 -Z$z -X1.5 -Y3 -K -P >tmp.ps #pswiggle out.dat -R -JX \ #-G200/250/250 -N -C-0.00 -M -W3 -Z$z -K -O >>tmp.ps pswiggle out.dat -R -JX -C-0.00 -M -W3 -Z$z -K -O >>tmp.ps pstext -N -R -J -O -K <>tmp.ps -330 25 15 90 5 6 Depth (km) END /bin/rm nrfs.dat grep ' 0.000 ' out.dat | sort -n >nrfs.dat awk '{print $1, $5}' nrfs.dat | \ psxy -R-200/200/0/200 -JX1.8/1.0 -Y5.6 -Ba90f30/a100f20 \ -K -O -W3/255/0/0 >>tmp.ps pstext -N -R -J -O -K <>tmp.ps -330 100 15 90 5 6 N. RFs/bin END ##--------------------- ### Middle panel. inner circles gmtset BASEMAP_AXES = weSn pswiggle out_inner.dat -R-200/200/00/50 -JX1.8/-5.0 \ -Ba90f30/a5f1 -G200/50/50 \ -C-0.00 -M -W3 -Z$z -X2.2 -Y-5.6 -O -K -P >>tmp.ps pswiggle out_inner.dat -R -JX -C-0.00 -M -W3 -Z$z -K -O >>tmp.ps pstext -N -R -J -O -K <>tmp.ps 0 55 15 0 5 6 Azimuth from North (deg.) END /bin/rm nrfs.dat grep ' 0.000 ' out_inner.dat | sort -n >nrfs.dat awk '{print $1, $5}' nrfs.dat | \ psxy -R-200/200/0/200 -JX1.8/1.0 -Y5.6 -Ba90f30/a100f20 -O \ -K -W3/255/0/0 >>tmp.ps ##--------------------- ### Right panel. outer circles pswiggle out_outer.dat -R-200/200/00/50 -JX1.8/-5.0 \ -Ba90f30/a5f1 -G200/50/50 \ -C-0.00 -M -W3 -Z$z -Y-5.6 -X2.2 -O -K -P >>tmp.ps #pswiggle out.dat -R -JX \ #-G200/250/250 -N -C-0.00 -M -W3 -Z$z -K -O >>tmp.ps pswiggle out_outer.dat -R -JX -C-0.00 -M -W3 -Z$z -K -O >>tmp.ps /bin/rm nrfs.dat grep ' 0.000 ' out_outer.dat | sort -n >nrfs.dat awk '{print $1, $5}' nrfs.dat | \ psxy -R-200/200/0/200 -JX1.8/1.0 -Y5.6 -Ba90f30/a100f20 -O -W3/255/0/0 >>tmp.ps ##--------------------- ps2jpeg150 tmp.ps