raydens_latsec.gmt

这是一个TOMO TO GMT的源码

#!/bin/csh
#
# script to plot ray-density tensor distribution
# in depth sections along constant latitude
#
#  3 sections on 1 page
#
# using output files from tomo2gmt
# for classified plotting
#
#					(c) [e-mail]

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

set psfile = $wpath/rayden_latsec.ps

# 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

# first section
set latsec=04
set lat=23.50
set infile=$inpath/VP/latsec.rayden.$latsec.xyz

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


psbasemap -Jx5/-0.05 -R-71.2730/-69.2166/-2.5/100 -Baf0.5:"longitude [degrees]":/a10f5:"depth [km]":WeNs -K -V -P -X5 -Y4 >$psfile

# mark block boundaries
psxy $inpath/VP/latsec_invgrd.xz -Jx -R -W0.5top -M -O -K -V -P >> $psfile

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

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

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

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

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

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

pstext -Jx -R -W255 -O -K -P -V <<END >>$psfile
-71.10 95 10 0 1 1 RAY DENSITY TENSOR along Latitude $lat S
END


# second section
set latsec=05
set lat=23.25
set infile=$inpath/VP/latsec.rayden.$latsec.xyz

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

# mark block boundaries
psxy $inpath/VP/latsec_invgrd.xz -Jx -R -Baf0.5:"longitude [degrees]":/a10f5:"depth [km]":WeNs -W0.5top -M -O -K -V -P -Y8 >> $psfile

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

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

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

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

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

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

pstext -Jx -R -W255 -O -K -P -V <<END >>$psfile
-71.10 95 10 0 1 1 RAY DENSITY TENSOR along Latitude $lat S
END

# third section
set latsec=06
set lat=23.00
set infile=$inpath/VP/latsec.rayden.$latsec.xyz

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

# mark block boundaries
psxy $inpath/VP/latsec_invgrd.xz -Jx -R -Baf0.5:"longitude [degrees]":/a10f5:"depth [km]":WeNs -W0.5top -M -O -K -V -P -Y8 >> $psfile

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

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

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

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

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

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

pstext -Jx -R -W255 -O -P -V <<END >>$psfile
-71.10 95 10 0 1 1 RAY DENSITY TENSOR along Latitude $lat S
END

gs $psfile

# cleanup
rm qual*