#/bin/sh -f
#
#set mod = m2n

#coupe
set xd = -72
set yd = -24.4
set xu =  195.7251115
set yu = 7298.1939733 

set azim = 77

set r1 = -R-70.7/-68.31/-26.7/-21.81
set r2 = -R-70.64/-69.2/-24.76/-22.7
set r3 = -R-70.67/-70/-23.91/-22.98

set modf = cyc.fault
set modo = out
set modof = out.fault

#ns  nbf along strike
set ns = 14

#nd nbf along dip
set nd = 8

# fichier de donnees a inverser


set vec1 = /data/ult3/sar3/sar_chili/t325_21355_20854/vecteur
set vec2 = /ult1/INTERFERO/96_4249-2455/MASSONNET/vecteur
set vec3 = /ult1/INTERFERO/367_20897_21398/vecteur

#nbd = 8677
set fic1 = ../DATA/4478_21856Fd_opeF2_sh-2.50.xyz
set rel1 = 3999

# nbd = 21571
set fic2 = ../DATA/4249_2455_opeF_cor.xyz
set rel2 =  9491
# centre : 9491
#nbd = 2446
#set fic2 = ../DATA/4249_2455_opeF2_cor.xyz
#set rel2 = 245
#set rel2 = 1068 
#(245 haut gauche 1068: centre)

# nbd = 15126
set fic3 = ../DATA/pha_20897_21398_topd_opePF1_cor.xyz
set rel3 = 10471
# nbd = 669
#set fic3 = ../DATA/pha_20897_21398_topd_opePF2_cor.xyz
#set rel3 = 461

#ATTENTION si decommente, detruit par  nearneighbor 
set grdd1 = data1.grd
set grdd2 = data2.grd
set grdd3 = data3.grd

set grdm1 = out_depts1.grd
set grdm2 = out_depts2.grd
set grdm3 = out_depts3.grd
set grdr1 = resi1.grd
set grdr2 = resi2.grd
set grdr3 = resi3.grd

set vec = /data/ult3/sar3/sar_chili/t_368/vecteur_moy

set fvecgps = /ult1/MODEL/INV2/GRAD/DATA/9592-vec.geo



\rm sour1.out
#INVERSION DONNEES 
invgrad <<EOI
$modf
2      UTM : 1 ; geographique :2
19      fuseau utm
$ns  nbf along strike
$nd nbf along dip
0   (normaliser)
3       nb de fichiers de donnees
$fic1
1      incertitudes sur les donnees (si lu dans fichier:999)
0      reference : 0 donnees absolues
$vec1
1   poids
$fic2
1      incertitudes sur les donnees (si lu dans fichier:999)
$rel2       reference : 0 donnees absolues
$vec2
1   poids
$fic3
1      incertitudes sur les donnees (si lu dans fichier:999)
$rel3      reference : 0 donnees absolues
$vec3
1   poids
0  (nb de recherches aleatoires
0 1000 10  glissement min, max, increment (cm)
100  (nb iterations)
EOI

\mv out1.fault out.fault
\mv data11.out data1.out 
\mv modl11.out modl1.out
\mv resi11.out resi1.out
\mv data12.out data2.out 
\mv modl12.out modl2.out
\mv resi12.out resi2.out
\mv data13.out data3.out 
\mv modl13.out modl3.out
\mv resi13.out resi3.out

projette.sh data1.out COUPES/data1
projette.sh data2.out COUPES/data2
projette.sh data3.out COUPES/data3
projette.sh modl1.out COUPES/modl1
projette.sh modl2.out COUPES/modl2
projette.sh modl3.out COUPES/modl3
projette.sh resi1.out COUPES/resi1
projette.sh resi2.out COUPES/resi2
projette.sh resi3.out COUPES/resi3

# modelisation des coupes continues
defor_fich <<EOF
2
2
w84-19.prj
$modo
COUPES/c1.xyd
3
0
$vec
COUPES/c1.xyz
1
COUPES/c2.xyd
3
0
$vec
COUPES/c2.xyz
1
COUPES/c3.xyd
3
0
$vec
COUPES/c3.xyz
1
COUPES/c4.xyd
3
0
$vec
COUPES/c4.xyz
1
COUPES/c5.xyd
3
0
$vec
COUPES/c5.xyz
1
COUPES/c6.xyd
3
0
$vec
COUPES/c6.xyz
1
COUPES/c7.xyd
3
0
$vec
COUPES/c7.xyz
1
COUPES/c8.xyd
3
0
$vec
COUPES/c8.xyz
1
COUPES/c9.xyd
3
0
$vec
COUPES/c9.xyz
0
EOF



plot_coup.sh


#TRACK 325
\rm out*bin
defor_grid <<EOF
2
2
w84-19.prj
$modo
97 197
-70.7 -26.7
0.025
1 (1 dep verti)
3 (3 gmt)
1  autre
2
10
1  autre
3   dirts
3 (3 gmt)
0 (autre)
$vec1
EOF
xyz2grd out_depts1.bin -Gout_depts1.grd -I0.025 -R-70.7/-68.31/-26.7/-21.81 -b -Z -V
nearneighbor data1.out -G$grdd1 -R-70.7/-68.31/-26.7/-21.81 -I0.025 -S0.03 -N1 -V
grdmath $grdd1 out_depts1.grd SUB = resi1.grd
grdtrack /data/ult5/ANTOFAGASTA/DATA/coup_cont.xy -Gresi1.grd > ! resi1_coup.xydz
grdtrack /data/ult5/ANTOFAGASTA/DATA/coup_cont.xy -Gout_depts1.grd > ! mod1_coup.xydz
coup_cont.sh $grdd1 COUPES/data1c
coup_cont.sh out_depts1.grd COUPES/modl1c
coup_cont.sh resi1.grd COUPES/resi1c


#track 96
\rm out*bin
defor_grid <<EOF
2
2
w84-19.prj
$modo
59 83
-70.64 -24.76
0.025
1 (1 dep verti)
3 (3 gmt)
1  autre
2
10
1  autre
3   dirts
3 (3 gmt)
0 (autre)
$vec2
EOF
xyz2grd out_depts1.bin -Gout_depts2.grd -R-70.64/-69.2/-24.76/-22.7 -I0.025 -b -Z -V
nearneighbor data2.out -G$grdd2 -R-70.64/-69.2/-24.76/-22.7 -I0.025 -S0.03 -N1 -V
grdmath $grdd2 out_depts2.grd SUB = resi2.grd
grdtrack /data/ult5/ANTOFAGASTA/DATA/coup_cont.xy -Gresi2.grd > ! resi2_coup.xydz
grdtrack /data/ult5/ANTOFAGASTA/DATA/coup_cont.xy -Gout_depts2.grd > ! mod2_coup.xydz
coup_cont.sh $grdd2 COUPES/data2c
coup_cont.sh out_depts2.grd COUPES/modl2c
coup_cont.sh resi2.grd COUPES/resi2c


#track 368
\rm out*bin
defor_grid <<EOF
2
2
w84-19.prj
$modo
28 38
-70.67 -23.91
0.025
1 (1 dep verti)
3 (3 gmt)
1  autre
2
10
1  autre
3   dirts
3 (3 gmt)
0 (autre)
$vec3
EOF
xyz2grd out_depts1.bin -Gout_depts3.grd -R-70.67/-70/-23.91/-22.98 -I0.025 -b -Z -V
nearneighbor data3.out -G$grdd3 -R-70.67/-70/-23.91/-22.98 -I0.025 -S0.03 -N1 -V
grdmath $grdd3 out_depts3.grd SUB = resi3.grd

grdtrack /data/ult5/ANTOFAGASTA/DATA/coup_cont.xy -Gresi3.grd > ! resi3_coup.xydz
grdtrack /data/ult5/ANTOFAGASTA/DATA/coup_cont.xy -Gout_depts3.grd > ! mod3_coup.xydz
coup_cont.sh $grdd3 COUPES/data3c
coup_cont.sh out_depts3.grd COUPES/modl3c
coup_cont.sh resi3.grd COUPES/resi3c


./plot_cont.sh

 

#determination de l'intersection avec le plan de faille
coup_faille << fin
$modof
2    geo
19   fuaseau
$xd $yd
1 $azim 1
projfault.geo
fin

projxyz << fin
projfault.geo
$xd $yd
$azim
10
projfault_p.geo
nul
0
0
fin

# calcul des points GPS
defor_fich <<EOF
2
2
w84-19.prj
$modo
$fvecgps
2
1
1
1
6
1
5
0
gps_h_calc.geo
gps_z_calc.geo
gps_x_calc.geo
gps_y_calc.geo
0
EOF
#


 
rms3 << fin > ! rms1.out
data1.out
modl1.out
0
fin
rms3 << fin > ! rms2.out
data2.out
modl2.out
0
fin
rms3 << fin > ! rms3.out
data3.out
modl3.out
0
fin

rms3 << fin > ! rms.out
data1.out
modl1.out
1
data2.out
modl2.out
1
data3.out
modl3.out
0
fin

#difference donnees inititial final --> si libre
dif3 << fin
data1.out
$fic1
difd_1
1
data2.out
$fic2
difd_2
1
data3.out
$fic3
difd_3
0
fin
head -1 difd_1 difd_2 difd_3




awk '{if ($1 == "TOTAL") print $7}' rms.out

set des = p.ps
set pal = pal_data
set grdd = data1.grd
set grdm = out_depts1.grd
set xy = data1.out
set dz0 = data1_p0.geo
set mz0 = data1_calcp0.geo
set dz1 = data1_p1.geo
set mz1 = data1_calcp1.geo
set dz2 = data1_p2.geo
set mz2 = data1_calcp2.geo

psbasemap -R-71.5/-68.4/-25.5/-22.5 -Y12 -Jm1:2400000 -Ba1f.5 -K -P > ! $des
#grdcontour $grdd -Ccontour -G5 -O -K -Jm -R -V -Wa20/0 >> $des
#grdcontour $grdm -Ccontour -O -K -Jm -R -V -Wa20/255/0/0  -F255/0/0>> $des
#psxy /data/data4/sismo_chili/sismic.dat -R -Jm -O -K -Sc.1 >> ! $des 
#awk '{print $3,$4 }' sismic_p.geo  |\
#psxy  -R -Jm -O -K -Sc.15 -G0 >> $des 
psvelomeca $fvecgps -Jm -R -A0.1/.3/.2 -Se0.03/.9/1  -O -K >> $des
psvelomeca gps_h_calc.geo -Jm -R  -Se0.03/.9/1 -L -W10/255/0/255 -O -K >> $des
psxy coup_cont_calc.xy -R -Jm -W5/255/0/0 -O -K >> $des
psxy $xy -R -Jm -C$pal -Sc.1 -O -K >> $des
psxy trfault.xy -Jm -R -M -W5  -O -K  >> $des
pscoast -R -Jm -O -K -W -Df >> $des
awk '{print $1,$2," 10 ",$4," 1 6 ",$7}' slip.txt |\
pstext -R -Jm -O -K -W255 >> $des

psxy <<!  -JX14.4/10 -R-71.5/-68.4/-70/5 -Ba1f.25/a10f2/SNW -Y-10 -O -K >> $des
!
awk '{print $3,$2,$2 }' sismic_p.geo |\
psxy   -Sc.05 -JX -R -O -K -G0>> $des

awk '{if ($1 != ">") print $3,$2 ; else print ">" }' projfault_p.geo |\
psxy  -JX -R  -W10/0/0/255 -O -K -M >> $des
awk '{print $1,$3,$3 }' coup_cont_calc.xy |\
psxy  -JX -R -W5/255/0/0 -O -K >> $des

awk '{print $3,$2,$2 }' $dz0 |\
psxy  -C$pal -Sc.1 -JX -R -O -K>> $des
awk '{print $3,$2,$2 }' $mz0 |\
psxy  -JX -R -Sc.05 -G150/0/0 -O -K >> $des
awk '{print $3,$2,$2 }' $dz1 |\
psxy  -C$pal -Sc.2 -JX -R -O -K>> $des
awk '{print $3,$2,$2 }' $mz1 |\
psxy  -JX -R -Sc.1 -G255/0/0 -O -K >> $des
awk '{print $3,$2,$2 }' $dz2 |\
psxy  -C$pal -Sc.1 -JX -R -O -K>> $des
awk '{print $3,$2,$2 }' $mz2 |\
psxy  -JX -R -Sc.05 -G255/200/200 -O -K >> $des



#awk '{if ($1 != "param") print $1,$3," 15 0.0 1 6 ",$8}' param.txt |\
#pstext -R -JX -O -K >> $des

awk '{if ($4 != "NaN") print $1,$4 }' resi_coup.xydz |\
psxy -R-71.5/-68.4/-10/10 -JX -O -K -W5/0/255/0 -Ba1f.25/a3f.5/E >> $des
awk '{if ($1 == "TOTAL") print "-69 -9 15 0.0  1 6 RMS ="$7}' rms.out |\
pstext -R -JX -O -K >> $des
psxy <<! -Jm -R -O >> $des
!
#xv $des &

set des = i1.ps
set grdd = $grdd1
set grdm = $grdm1
set grdr = $grdr1
set range = $r1
set pal = pal_data
set palr = pal_res
psbasemap $range  -Jm1:3500000 -Ba1f.5/WSNe -K  > ! $des
grdimage $grdd $range -Jm -O -K -C$pal >> $des
grdcontour $grdd -Canot  -G3 -O -K -Jm $range -V >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
psxy COUPES/c1.xyz $range -Jm -O -K >> $des
psxy COUPES/c2.xyz $range -Jm -O -K >> $des
psxy COUPES/c3.xyz $range -Jm -O -K >> $des
psxy COUPES/c4.xyz $range -Jm -O -K >> $des
psxy COUPES/c5.xyz $range -Jm -O -K >> $des
psxy COUPES/c6.xyz $range -Jm -O -K >> $des
psxy COUPES/c7.xyz $range -Jm -O -K >> $des
psxy COUPES/c8.xyz $range -Jm -O -K >> $des
psxy COUPES/c9.xyz $range -Jm -O -K >> $des
pstext <<! $range -Jm -O -K  >> $des
-69.8 -23.5 13 0.0 1 6 c1
-69.6 -23.5  13 0.0 1 6 c2
-70.2 -23.4  13 0.0 1 6 c3
-70.63 -23.27 13 0.0 1 6 c4
-70.65 -23.50  13 0.0 1 6 c5
-70.55 -23.10 13 0.0 1 6 c6
-70.6 -24.40  13 0.0 1 6 c7 
-70.6 -23.95  13 0.0 1 6 c8
-69.3 -23.7 13 0.0 1 6 c9
!
psscale -C$pal -D5.5/3/4/0.5 -O -K >> $des
pstext <<! -Jm $range -O -K >> $des
-69.5 -22 12 0.0 5 6 $fic1
!
head -1 difd_1 | awk '{print "-70 -22.2 12 0.0 5 6 shift=",$3 }'  |\
pstext -Jm $range -O -K >> $des

psxy <<! $range -Jm -Ba1f.5/NSew -K -O -X8.5 >> $des
!
grdimage $grdm $range -Jm -O -K -C$pal >> $des
grdcontour $grdm -Canot  -G3  -O -K -Jm $range -V >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
psscale -C$pal -D5.5/3/4/0.5 -O -K >> $des
psxy <<! $range -Jm -Ba1f.5/NSEw -K -O -X8.5 >> $des
!
grdimage $grdr $range -Jm -O -K -C$palr >> $des
grdcontour $grdr -A1 -G3 -O -K -Jm $range -V >> $des
psscale -C$pal -D10.3/3/4/0.5 -O -K >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
awk '{if ($1 == "TOTAL") print "-70 -22 15 0.0  1 6 RMS ="$7}' rms1.out |\
pstext $range -Jm -O -K >> $des
psscale -C$palr -D5.5/3/4/0.5 -O -K >> $des
psxy <<! -Jm $range -O >> $des
!
#xv $des &

set des = i2.ps
set grdd = $grdd2
set grdm = $grdm2
set grdr = $grdr2
set range = $r2
set pal = pal_data
set palr = pal_res
psbasemap $range  -Jm1:1500000 -Ba1f.5/WSNe -K  > ! $des
grdimage $grdd $range -Jm -O -K -C$pal >> $des
grdcontour $grdd -Canot  -G3  -O -K -Jm $range -V >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
psxy COUPES/c1.xyz $range -Jm -O -K >> $des
psxy COUPES/c2.xyz $range -Jm -O -K >> $des
psxy COUPES/c3.xyz $range -Jm -O -K >> $des
psxy COUPES/c4.xyz $range -Jm -O -K >> $des
psxy COUPES/c5.xyz $range -Jm -O -K >> $des
psxy COUPES/c6.xyz $range -Jm -O -K >> $des
psxy COUPES/c7.xyz $range -Jm -O -K >> $des
psxy COUPES/c8.xyz $range -Jm -O -K >> $des
psxy COUPES/c9.xyz $range -Jm -O -K >> $des
pstext <<! $range -Jm -O -K  >> $des
-69.8 -23.5 13 0.0 1 6 c1
-69.6 -23.5  13 0.0 1 6 c2
-70.2 -23.4  13 0.0 1 6 c3
-70.63 -23.27 13 0.0 1 6 c4
-70.65 -23.50  13 0.0 1 6 c5
-70.55 -23.10 13 0.0 1 6 c6
-70.6 -24.40  13 0.0 1 6 c7 
-70.6 -23.95  13 0.0 1 6 c8
-69.3 -23.7 13 0.0 1 6 c9
!
psscale -C$pal -D6.3/3/4/0.5 -O -K >> $des
pstext <<! -Jm $range -O -K >> $des
-70 -22.77 12 0.0 5 6 $fic2
!
head -1 difd_2 | awk '{print "-70.5 -22.81 12 0.0 5 6 shift=",$3 }'  |\
pstext -Jm $range -O -K >> $des

psxy <<! $range -Jm -Ba1f.5/NSew -K -O -X8.5 >> $des
!
grdimage $grdm $range -Jm -O -K -C$pal >> $des
grdcontour $grdm -Canot  -G3  -O -K -Jm $range -V >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
psscale -C$pal -D6.3/3/4/0.5 -O -K >> $des
psxy <<! $range -Jm -Ba1f.5/NSEw -K -O -X8.5 >> $des
!
grdimage $grdr $range -Jm -O -K -C$palr >> $des
grdcontour $grdr -A1 -G3 -O -K -Jm $range -V >> $des
psscale -C$pal -D10.3/3/4/0.5 -O -K >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
awk '{if ($1 == "TOTAL") print "-70.2 -22.78 15 0.0  1 6 RMS ="$7}' rms2.out |\
pstext $range -Jm -O -K >> $des
psscale -C$palr -D6.3/3/4/0.5 -O -K >> $des
psxy <<! -Jm $range -O >> $des
!
#xv $des &

set des = i3.ps
set grdd = $grdd3
set grdm = $grdm3
set grdr = $grdr3
set range = $r3
set pal = pal_data
set palr = pal_res
psbasemap $range  -Jm1:1000000 -Ba1f.5/WSNe -K  > ! $des
grdimage $grdd $range -Jm -O -K -C$pal >> $des
grdcontour $grdd -Canot  -G3  -O -K -Jm $range -V >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
psxy COUPES/c1.xyz $range -Jm -O -K >> $des
psxy COUPES/c2.xyz $range -Jm -O -K >> $des
psxy COUPES/c3.xyz $range -Jm -O -K >> $des
psxy COUPES/c4.xyz $range -Jm -O -K >> $des
psxy COUPES/c5.xyz $range -Jm -O -K >> $des
psxy COUPES/c6.xyz $range -Jm -O -K >> $des
psxy COUPES/c7.xyz $range -Jm -O -K >> $des
psxy COUPES/c8.xyz $range -Jm -O -K >> $des
psxy COUPES/c9.xyz $range -Jm -O -K >> $des
pstext <<! $range -Jm -O -K  >> $des
-69.8 -23.5 13 0.0 1 6 c1
-69.6 -23.5  13 0.0 1 6 c2
-70.2 -23.4  13 0.0 1 6 c3
-70.63 -23.27 13 0.0 1 6 c4
-70.65 -23.50  13 0.0 1 6 c5
-70.55 -23.10 13 0.0 1 6 c6
-70.6 -24.40  13 0.0 1 6 c7 
-70.6 -23.95  13 0.0 1 6 c8
-69.3 -23.7 13 0.0 1 6 c9
!
psscale -C$pal -D6.3/3/4/0.5 -O -K >> $des
pstext <<! -Jm $range -O -K >> $des
-70.4 -23 12 0.0 5 6 $fic2
!
head -1 difd_3 | awk '{print "-70.5 -23.05 12 0.0 5 6 shift=",$3 }'  |\
pstext -Jm $range -O -K >> $des

psxy <<! $range -Jm -Ba1f.5/NSew -K -O -X8.5 >> $des
!
grdimage $grdm $range -Jm -O -K -C$pal >> $des
grdcontour $grdm -Canot  -G3  -O -K -Jm $range -V >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
psscale -C$pal -D6.3/3/4/0.5 -O -K >> $des
psxy <<! $range -Jm -Ba.5f.1/NSEw -K -O -X8.5 >> $des
!
grdimage $grdr $range -Jm -O -K -C$palr >> $des
grdcontour $grdr -A1 -G3 -O -K -Jm $range -V >> $des
psscale -C$pal -D10.3/3/4/0.5 -O -K >> $des
psxy trfault.xy -Jm $range -M -W5  -O -K  >> $des
pscoast $range -Jm -O -K -W -Df >> $des
awk '{if ($1 == "TOTAL") print "-70.5 -23 15 0.0  1 6 RMS ="$7}' rms3.out |\
pstext $range -Jm -O -K >> $des
psscale -C$palr -D6.3/3/4/0.5 -O -K >> $des
psxy <<! -Jm $range -O >> $des
!
#xv $des &

xyz2grd sour1.out -Gsour1.grd -R1/$nd/1/$ns -I1 -Z -b  
grdsample sour1.grd -Gsour1S.grd -I.1

set des = plan.ps
psbasemap -R-71.6/-69.1/-25.5/-22.7  -Jm1:2200000 -Ba1f.5/WsNe -K -Y10 -P > !  $des
psxy trfault.xy -Jm -R -M -W5  -O -K  >> $des
pscoast -R -Jm -O -K -W -Df >> $des
psxy /data/topo1/data4/sismo_chili/mainshock.xy -Sc.2 -G255/0/0 -R -Jm -O -K  >> $des

psvelomeca $fvecgps -Jm -R -A0.1/.3/.2 -Se0.03/.9/1  -O -K >> $des
psvelomeca gps_h_calc.geo -Jm -R  -Se0.03/.9/1 -L -W10/255/0/255 -O -K >> $des
awk '{if ($1 != ">") print $3,$4 ; else print ">" }' projfault_p.geo |\
psxy -Jm -R  -Sc.12 -G0/0/255 -O -K -M >> $des
awk '{print $1,$2," 13 ",$4," 4 6 ",$7}' slip.txt |\
pstext -R -Jm -O -K  >> $des
awk '{print "-70 -25.45 15 0.0 4 6", $7,$8,$9,$10 }' moment.txt |\
pstext   -R -Jm -O -K  >> $des
awk '{if ($1 == "TOTAL") print "-69.5 -25.3 15 0.0  1 6 RMS ="$7}' rms.out |\
pstext  -R -Jm -O -K >> $des

psxy <<! -R-71.6/-69.1/-70/0 -Jx5.1/.07 -O -K -Ba1f.2/a10f2/NSWe -Y-6 >> $des
!
awk '{if ($1 != ">") print $3,$2 ; else print ">" }' projfault_p.geo |\
psxy  -Jx -R  -W10/0/0/255 -O -K >> $des
awk '{if ($1 != ">") print $3,$2 ; else print ">" }' projfault_p.geo |\
psxy  -Jx -R  -Sc.2 -G0/0/255 -O -K -M >> $des
awk '{print $3,$2,$2 }' sismic_p.geo |\
psxy   -Sc.05 -Jx -R -O -K -G0>> $des

grdimage sour1.grd -Cpal_p -R1/$nd/1/$ns -Jx0.65 -B1/nsew -K -O -X13.0 -Y12 >>  $des
grdcontour sour1.grd  -Ccont_plan  -O -R -Jx -K >> $des
awk '{print '$nd'-3,'$ns'-3,$1 }' null |\
psxy  -Jx -O -R -K -Sc.2 -G255/0/0 >> $des

grdimage sour1S.grd -Cpal_p -R1/$nd/1/$ns -Jx0.65 -B1/nsew -K -O  -Y-12 >>  $des
grdcontour sour1S.grd  -Ccont_plan  -O -R -Jx -K  >> $des
awk '{print '$nd'-3,'$ns'-3,$1 }' null |\
psxy  -Jx -O -R -K -Sc.2 -G255/0/0 >> $des

psscale -X-14  -Cpal_p -D10.5/17/6/.6 -O -K >> $des
psxy <<! -Jx -R -O  >> $des
!
xv $des &