inverse.txt

射线追踪程序

#!/bin/sh
(
PATH=$PATH:`pwd`/bin
test "."$1 = "." -o "."$2 = "." -o "."$3 = "." -o "."$4 = "." \
     && echo "$0 name OA1|OA2|1 RT|NRT [B|N]"
test "."$1 = "." -o "."$2 = "." -o "."$3 = "." -o "."$4 = "." && exit

(
test "."$2 = ".F" && echo OA=\"`cat OA.file`\"
test "."$2 = ".OA1" && echo "OA=\"0 1 2\""
test "."$2 = ".OA2" && echo "OA=\"10 11 12 13 14 15 16 17\""
test "."$2 = ".1" && echo OA=

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

rm -f final.geo final.xyz
rm -f $1.par $1.geo $1.var
test $RT = "TODO" &&\
     raytrace $NORM -m $MINR -i $INC \
           -n $1 -p $1.dat -Z $MAXDEPTH

layer2mod $1
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 $PARAM = "NODES" && \
       ray2Nmatrix -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
)
lsq -v -n $1 $damp #-r
block2dat $1
cat $1.geo >> final.geo
cat $1.xyz >> final.xyz
for layer in `awk '{ print $6 }' < final.geo | grep -v layer0 | sed 's/layer//' | sort | uniq`

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;
            print;
            la=$1;
            lo=$2;
            ld=$3; }' > final.geo.sort

for layer in `awk '{ print $6 }' < final.geo | sed 's/layer//' | sort | uniq`
    do ./mkgrdpl.tri $layer final.geo
    done
) 2>&1 | tee inverse.log