raydens_plane.gmt

这是一个TOMO TO GMT的源码

#!/bin/csh
#
# script to plot ray-density tensor in planes
#
# using output from tomo2gmt
#
# for classified plotting
#
# classes are a combination of 
# ray distribution & weighted ray length
#
#						(c) [e-mail]

set plane=015.00
set depth=15.0

# !! be sure to set paths to your correct directory !!
set wpath=/ek1/tomoled/tomo2gmt/PACKAGE
set inpath=$wpath/example_raydens

set psfile=$wpath/rayden$plane.ps
set infile=$inpath/VP/rayden$plane.xyz

# set colors for each class
set c1=255/10/10
set c2=235/185/30
set c3=10/10/255
set c4=80/200/0

# create qual. files foreach class and total ray length class
nawk -f $wpath/GMT/splitrayden.awk $infile 

psbasemap -Jm6 -R-71.2730/-69.2166/-25.1260/-22.0547 -Ba30m -K -V -P -X5 -Y5 >$psfile

# plot inv. grid
psxy $inpath/VP/invgrd.xy -Jm -R -W0.5top -M -O -K -V -P >> $psfile

# Positionen der Landstationen einzeichnen
nawk '{print $1,$2,"0"}' ~/GMT/pos_land.xy |\
psxyz  -R -Jm -St0.3 -L -G0 -O -P  -K  >>$psfile

# Positionen der OBH einzeichnen
nawk '{print $1,$2,"0"}' ~/GMT/pos_obh.xy |\
psxyz  -R -Jm -St0.3 -G0 -O -P -K >>$psfile

# coastlines
pscoast -Jm -R  -A0/1 -Di -W3 -O -K -V >> $psfile

# work on class1 (circles)
if (-f qual1_1.out) psxy qual1_1.out -Jm -R -G$c1 -Sc -O -K -V -P >>$psfile

# on class1
if (-f qual1_1.out) psxy qual1.out -Jm -R -W20/$c1 -M -O -K -V -P >>$psfile

# on class2 (circles)
if (-f qual2_1.out) psxy qual2-1.out -Jm -R -G$c2 -Sc -O -K -V -P >>$psfile

# on class2 
psxy qual2.out -Jm -R -W20/$c2 -M -O -K -V -P >> $psfile

# on class3 
psxy qual3.out -Jm -R -W20/$c3 -M -O -K -V -P >> $psfile

# on class3 
psxy qual4.out -Jm -R -W20/$c4 -M -O -K -V -P >> $psfile

pstext -Jm -R -N -O -P -V <<END >>$psfile
-70.80 -21.70 12 0 0 1 ray density tensor PLANE$plane
END

gs $psfile

# cleanup
rm qual*