trplall.txt

射线追踪程序

#!/bin/sh
INC=0.05
INC2=0.01
#REG=-3.5/-0.5/46.5/49.0
#REG=-4.0/6.0/44/50.0
REG=-5/10.0/41/50
#REG=-3/4.0/41.5/43.5
X=6
Y=4.5
T=5.5
#A1="-3 47"
#B1="-1.5 48.8"
#A2="-3.5 48.75"
#B2="-0.5 46.75"
LCR="6/255/255/255"

#./mkgrdpl

awk '{ print $4,$3 }' < all.sta > sta
\more all.lay |\
    awk -v layer=-1 -v Z=0 '{ LZ=Z; Z+=$5;
                             printf(" Layer %d %.0f-%.0fkm\n",layer,LZ,Z)
                             layer++;
                           }' > ti-layer

(
pstext -V -R0/1/0/1 -N -JM$T -X2 -Y12 -K \
       -U/0/-10.5/"Judenherc & al., 2001" << flag
flag

for i in `\ls dataf*grd | grep -v dataf0.grd |\
           sed -e 's/dataf//' -e 's/.grd//' | sort -n`
do SHIFT=-X$X
   test $i = 5 && SHIFT="-Y-6.25 -X-18"
   test $i = 9 && SHIFT="-Y-6.25 -X-18"
   test $i = 1 && SHIFT="-Y3"
   grdimage $SHIFT -O -JM$T -K -V -Cdeltav.cpt -R$REG dataf$i.grd
   #grdcontour -O -JM$T -K -V -C1 -R$REG dataf$i1.grd -G192
   #grdcontour -O -JM$T -K -V -C2 -R$REG dataf$i1.grd -G192
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/farmors.dat
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/fcentrals.dat
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/central.dat
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/rhenan.dat
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/fpyrs.dat
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/falpess.dat
   psxy -W2 -O -K -R -V -JM -M $HOME/Data/Cartes/centralV.dat
   psxy -Sc0.02 -G0 -O -K -R -V -JM -V sta
   BB=-B2g2::wens
   test $i = 1 && BB=-B2g2::WeNs
   pscoast -JM -O -Dh -R -N1 -V -K $BB -W1/0/0/0
   echo 0.1 0.05 9 0 1 1 `grep "Layer "$i" " ti-layer` |\
          pstext -N -O -K -R0/1/0/1 -V -JX$T
done
psscale -V -O -Cdeltav.cpt -D4/2/5.5/0.5 -X$X
) > layers.ps