mkgrdpl.tri.txt

射线追踪程序

#!/bin/sh
# ./mkgrdpl <couche> fichier-param
set -x
INC=0.025
#REG=-4.0/6.0/44/49.0
REG=-5/12.0/40/53
#REG=-2.5/4.0/42/46
MINRES=-10
X=12
Y=9
T=9.0
RES=-Dh

#DOFIL=1
MASK=1
FIL=`grep layer-$1" " FIL-NB.dat | awk '{ print 1.5*$2 }'`
RAD=`grep layer-$1" " FIL-NB.dat | awk '{ print $2"k" }'`

test "."$3 = ".G" || (
   ################# masque
   grep layer$1" " $2.sort | awk '{ print $1 " " $2 " " $7 }' > tmp.dat
   wc -l tmp.dat
   test "."$MASK = "." || nearneighbor -N1 -V -S$RAD -I$INC -R$REG -Gmask.grd tmp.dat
   #grdsample tmp.grd -Gmask.grd -I$INC
   #grdcontour -V -O -JM -K  -C1 -R tmp.grd
   

   ################# perturbations
   grep layer$1" " $2.sort |\
         awk -v minres=$MINRES '{ if ($7>minres) print $1 " " $2 " " $4 }' > tmp.dat
   wc -l tmp.dat
   blockmean  -I$INC -R$REG tmp.dat | triangulate -I$INC -R$REG -Gtmp.grd
   cp tmp.grd tmpdata.grd
   test "."$DOFIL = "." || grdfilter -V tmp.grd -Gtmpdata.grd -D4 -Fg$FIL -R$REG
   cp tmpdata.grd dataf$1.grd
   test "."$MASK = "." || grdmath -R -V tmpdata.grd mask.grd OR = dataf$1.grd

   test $1 = "0" && blockmean -I$INC -R$REG tmp.dat > dataf$1.grd

   ################# erreurs
   grep layer$1" " $2.sort | awk '{ print $1 " " $2 " " $5 }' > tmp.dat
   wc -l tmp.dat
   blockmean  -I$INC -R$REG tmp.dat | triangulate -I$INC -R$REG -Gtmp.grd
   cp tmp.grd tmpdata.grd
   #grdfilter tmp.grd -Gtmpdata.grd -D4 -Fg$FIL -R$REG 
   #grdsample tmpdata.grd -Gerror$1.grd -I$INC 
   cp tmpdata.grd error$1.grd
   test "."$MASK = "." || grdmath -R -V tmpdata.grd mask.grd OR = error$1.grd
   test $1 = "0" && blockmean -I$INC -R$REG tmp.dat > error$1.grd

   ################# diag res
   grep layer$1" " $2.sort | awk '{ print $1 " " $2 " " $8 }' > tmp.dat
   wc -l tmp.dat
   blockmean  -I$INC -R$REG tmp.dat | triangulate -I$INC -R$REG -Gtmp.grd
   cp tmp.grd tmpdata.grd
   #grdfilter tmp.grd -Gtmpdata.grd -D4 -Fg$FIL -R$REG 
   #grdsample tmpdata.grd -Gres$1.grd -I$INC 
   cp tmpdata.grd res$1.grd
   test "."$MASK = "." || grdmath -R -V tmpdata.grd mask.grd OR = res$1.grd
   test "."$1 = ".0" && blockmean -I$INC -R$REG tmp.dat > res$1.grd
   
   rm -f tmp* data
)


(
pstext  -R0/1/0/1 -N -JM8 -X2.5 -Y2.5 -K << flag
flag
#0 -0.1 8 0 1 1 Judenherc \& al., 1998
#flag

test "."$1 = ".0" && psxy -Sd0.3 -O -JM -K  -Cdeltav.cpt -R$REG dataf$1.grd 
test "."$1 = ".0" || grdimage -R$REG -O -JM -K -Cdeltav.cpt dataf$1.grd 
psxy -W3 -O -K -R  -JM -M -B1g1::WeNs $HOME/Data/Cartes/farmors.dat
psxy -Sc0.02 -G0 -O -K -R  -JM sta
pscoast -JM -O $RES -R -N1 -K -B1g1::WeNs -W1/0
pstext -O -K -R  -JM << flag
-3.5 49.0 7 0 1 1 Layer $1
flag
psscale -K  -O -Cdeltav.cpt -D8.5/4/8/0.5

test "."$1 = ".0" && psxy -Sd0.3 -Y9 -O -JM -K  -Cerror.cpt -R error$1.grd
test "."$1 = ".0" || grdimage -Y9 -O -JM -K  -Cerror.cpt -R error$1.grd
#grdcontour -O -JM -K  -C0.1 -R error$1.grd -G192
psxy -W3 -O -K -R  -JM -M -B1g1::WeNs $HOME/Data/Cartes/farmors.dat
psxy -Sc0.02 -G0 -O -K -R  -JM sta
pscoast -JM -O $RES -R -N1  -K -B1g1::WeNs -W1/0
pstext -O -K -R  -JM << flag
-3.5 49.0 7 0 1 1 Layer $1
flag
psscale -K  -O -Cerror.cpt -D8.5/4/8/0.5

test "."$1 = ".0" && psxy -Sd0.3 -X13 -O -JM -K  -Cres.cpt -R res$1.grd
test "."$1 = ".0" || grdimage -X13 -O -JM -K  -Cres.cpt -R res$1.grd
#grdcontour -O -JM -K  -C0.1 -R res$1.grd -G192
psxy -W3 -O -K -R  -JM -M -B1g1::WeNs $HOME/Data/Cartes/farmors.dat
psxy -Sc0.02 -G0 -O -K -R  -JM sta
pscoast -JM -O $RES -R -N1  -K -B1g1::WeNs -W1/0
pstext -O -K -R  -JM << flag
-3.5 49.0 7 0 1 1 Layer $1
flag
psscale -L -K  -O -Cres.cpt -D8.5/4/8/0.5

psxy -Y-9 -W3 -O -K -R  -JM -M -B1g1::WeNs $HOME/Data/Cartes/farmors.dat
psxy -Sc0.02 -G0 -O -K -R  -JM sta
pscoast -JM -O $RES -R -N1  -K -B1g1::WeNs -W1/0
grep layer$1" " $2.sort | awk '{ print $1 " " $2 " 10 0 1 6 " $7 }' |\
     pstext -O  -R -JM -K

pstext  -R -N -JM << flag
flag

) > layer-$1.ps