inv.txt

射线追踪程序

#!/bin/sh
PATH=$PATH:bin:$HOME/src/tomo3d
export PATH
PROJ_LIB=.
export PROJ_LIB

test "."$1 = "." -o "."$2 = "." \
    && exec echo "$0 name clean|F|OA1|OA2|1 RT|NRT [B|N] [dir layer]"

test "."$2 = ".clean" \
    && exec rm -f final.geo final.xyz   $1.sys $1.sta $1.blk $1.out \
                $1.hit $1.num $1.inv $1.ray $1.res $1.xyz $1-blk.gmt \
                sta FIL-NB.dat final.geo.sort inverse.log ti-layer

test "."$3 = "." -o "."$4 = "." \
    && exec echo "$0 name clean|F|OA1|OA2|1 RT|NRT [B|N] [dir layer]"

(
echo OA=-1

test "."$2 = ".F" && echo OA=\"`grep -v ^\# LAYER.file`\"
test "."$2 = ".F" && echo LAYERS=\"`grep -v ^\# OA.file`\"
test "."$2 = ".A" && echo OA=$5
test "."$2 = ".A" && echo LAYERS=$6
test "."$2 = ".OA1" && echo "OA=\"-1 0 1 2\""
test "."$2 = ".OA2" && echo "OA=\"-1 10 11 12 13 14 15 16 17\""

test "."$3 = ".RT" && echo RT=TODO
test "."$3 = ".NRT" && echo RT=NO

test "."$4 = ".N" && echo PARAM=NODES
test "."$4 = ".B" && echo PARAM=BLOCKS
) > $1.inv

. $1.cnf
. $1.inv

clear
echo Running inv for OA=$OA LAYERS=$LAYERS

test $RT = "TODO" &&\
     raytrace $NORM -m $MINR -i $INC \
           -n $1 -p $1.dat -Z $MAXDEPTH

test "."$2 = ".F" || LAYERS=`cat $1.lay | wc -l | awk '{ for (i=1; i<$1-1; i++) print i }'`

test "."$2 = ".1" && LAYERS=-1

test "."$2 = ".A" && LAYERS=$6
test "."$2 = ".A" && OA=$5

for layer in $LAYERS
do for dir in $OA
       do layer2mod $1 $layer $dir
          ray2Bmatrix $ACH -n $1 -m $MINH -s $SIG2
          test $PARAM = "NODES" && \
                ray2Nmatrix -n $1 -m $MINH -s $SIG2
          test "."$REBLK1 = ".1" && (
             reblk $1 $MAXHIT1 $MAXLAY &&\
                ray2Bmatrix $ACH -n $1 -m $MINH -s $SIG2
          )
          test "."$REBLK2 = ".1" && (
             reblk $1 $MAXHIT2 $MAXLAY &&\
                ray2Bmatrix $ACH -n $1 -m $MINH -s $SIG2
                test $PARAM = "NODES" &&\
                      ray2Nmatrix -n $1 -m $MINH -s $SIG2
          )
          rm -f $1.par $1.geo
          lsq -v -n $1 $damp
          block2dat $1
          cat $1.geo >> final.geo
          cat $1.xyz >> final.xyz
       done
done

awk '{ print $3 " " $2 }' < $1.dat | sort | uniq > sta

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

cat $1.lay | awk -v layer=-1 '{ printf("layer-%d %.0f\n",
                                layer,0.50*($2+$4));
                                layer++;
                              }' > FIL-NB.dat

rm -f dataf*.grd

sort -n final.geo |\
     awk -v la=0 -v lo=0 -v ld=-1 '{
            if (($1!=la)||($2!=lo)||($3!=ld)) print;
            la=$1;
            lo=$2;
            ld=$3; }' > final.geo.sort

test -f mkgrdpl.tri && (
    for layer in `awk '{ print $6 }' <final.geo|sed 's/layer//' | sort | uniq`
        do ./mkgrdpl.tri $layer final.geo
    done
)