#/bin/sh -f
#
set modf = ma_45.fault
set modo = out
set modof = out.fault

#ns  nbf along strike
set ns = 9

#nd nbf along dip
set nd = 5


# fichier de donnees a inverser
set fic = ../DATA/dataF.xyz
set vec = /data/ult3/sar3/sar_chili/t_368/vecteur_moy

set fvecgps = /ult1/MODEL/INV2/GRAD/DATA/9592-vec.geo
#set fvecgps = /data/data4/gps_chili/cosismic/cos_anto_vec.geo


#INVERSION DONNEES 
#inv_gen_vec_plan << EOI
#2      parametrisation 1 ---> type invers3  ; 2 ---> type jbc
#2      UTM : 1 ; geographique :2
#19 -1      fuseau utm, et hemisphere (1:nord, -1:sud)
#m3.fault
#4       nb de fichiers de donnees
#$fic1
#1 1       facteurs multiplicatifs km (ou deg) et cm (1 1)
#1      incertitudes sur les donnees (si lu dans fichier:999)
#0       reference : 0 donnees absolues
#$vec1
#0 0 0
#0 0 0
#$fic2
#1 1       facteurs multiplicatifs km (ou deg) et cm (1 1)
#1      incertitudes sur les donnees (si lu dans fichier:999)
#0       reference : 0 donnees absolues
#$vec2
#0 0 0
#0 0 0
#$fic3
#1 1       facteurs multiplicatifs km (ou deg) et cm (1 1)
#1      incertitudes sur les donnees (si lu dans fichier:999)
#0       reference : 0 donnees absolues
#$vec3
#0 0 0
#0 0 0
#$fic4
#1 1       facteurs multiplicatifs km (ou deg) et cm (1 1)
#1      incertitudes sur les donnees (si lu dans fichier:999)
#0       reference : 0 donnees absolues
#$vec4
#0 0 0
#0 0 0
#EOI
#
#

invgrad <<EOI
$modf
2      UTM : 1 ; geographique :2
19      fuseau utm
$ns  nbf along strike
$nd nbf along dip
1   (normaliser)
100   (cm)
1       nb de fichiers de donnees
$fic
1      incertitudes sur les donnees (si lu dans fichier:999)
1       reference : 0 donnees absolues
$vec
1   poids
0  (nb de recherches aleatoires
0 1000 10  glissement min, max, increment (cm)
100  (nb iterations)
EOI
#
#/data/data4/data_chili/gps_deux/9592-x.geo
#/data/data4/data_chili/gps_deux/vecteur_x
#/data/data4/data_chili/gps_deux/9592-y.geo
#/data/data4/data_chili/gps_deux/vecteur_y
#
#
#
# ATTENENTION *********************
# si model direct
#\cp dir.fault out.fault

\cp out1.fault out.fault

defor_grid <<EOF
2
2
w84-19.prj
$modo
93 185
-70.7 -26.5
0.025
3  (dir t-s)
3
1
2  horiz vect
10
1
1 (1 dep verti)
3 (3 gmt)
0 (autre)
$vec
EOF

xyz2grd out_depts1.bin -Gout_depts1.grd -I0.025 -R-70.7/-68.4/-26.5/-21.9 -b -Z -V
set des = mod.ps
psbasemap -R-70.7/-68.4/-26.5/-21.9 -Jm1:2200000 -Ba1f.5 -K -P > ! $des
grdimage  out_depts1.grd -C../DATA/pal_data -R -Jm -O -K >> $des
grdcontour out_depts1.grd -C5 -A10 -O -K -Jm -R -V -G5>> $des
pscoast -R -Jm -O -K -W -Df >> $des
psxy <<! -Jm -R -O >> $des
!
xv $des &

COUP/coup.sh out_depts1.grd /data/ult5/ANTOFAGASTA/INV/COUP/m 

grdmath ../DATA/dataU.grd out_depts1.grd SUB = res.grd
COUP/coup.sh res.grd /data/ult5/ANTOFAGASTA/INV/COUP/r 


set des = res.ps
psbasemap -R-70.7/-68.4/-26.5/-21.9 -Jm1:2200000 -Ba1f.5 -K -P > ! $des
grdimage  res.grd -Cpal_r -R -Jm -O -K >> $des
grdcontour res.grd -C1 -A3 -O -K -Jm -R -V -G5>> $des
pscoast -R -Jm -O -K -W -Df >> $des
psxy <<! -Jm -R -O >> $des
!
xv $des &

COUP/plot_coup.sh  1
COUP/plot_coup.sh  2
COUP/plot_coup.sh  3
COUP/plot_coup.sh  4
COUP/plot_coup.sh  5
COUP/plot_coup.sh  6
COUP/plot_coup.sh  7
COUP/plot_coup.sh  8
COUP/plot_coup.sh  9


#mv out.fault $modof
#/ult1/MODEL/INV2/DATA/data1_c1.xy
#/ult1/MODEL/INV2/DATA/data2_c1.xy

# modelisation des  des donnees radars 
defor_fich <<EOF
2
2
w84-19.prj
$modo
data11.out
3
0 (autre grandeur
$vec
data11_calc.geo
0 (autre  fichier de poits
EOF


#
# calcul des points GPS
defor_fich <<EOF
2
2
w84-19.prj
out
$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
#
trcoor1 <<EOF
w84-19.prj
1   (km=1  m=2)
trfault.xy
trfault.utm
1
gps_z_calc.geo
gps_z_calc.utm
1
../DATA/9592-z.geo
gps_z_dat.utm
1
../DATA/9592-x.geo
gps_x_dat.utm
1
../DATA/9592-y.geo
gps_y_dat.utm
0
EOF
projxyz <<EOF
trfault.utm
193.6033      7397.9281
95
1000000
trf_95p.dz
nul
0
1
gps_z_calc.utm
193.6033      7397.9281
95
1000000
gps_z_calc.dz
nul
0
1
gps_z_dat.utm
193.6033      7397.9281
95
1000000
gps_z_calc.dz
nul
0
0
EOF
#
psxy /data/data4/sismo_chili/sismicutm/sismic-95p.dz -R0/300/-70/7 -Jx.06/.06 -Ba50f10/10f2/SnWe -Sc.1 -K -P > ! p.ps
psxy /data/data4/sismo_chili/sismicutm/main-95p.dz -R -Jx -Sa.25 -G255/0/255 -W10 -O -K>> p.ps
psxy trf_95p.dz -W15/255/0/0 -M -Jx -R -O -K >> p.ps
psxy trf_95p.dz -W15/255/0/0 -M -Jx -Sc.1 -R -O -K >> p.ps
psxy /data/data4/sismo_chili/sismicutm/topo/t2-95.dz -R0/300/-70000/7000 -Jx.06/0.00006 -O -K>> p.ps
psxy /data/data4/sismo_chili/sismicutm/topo/t1-95.dz -R0/300/-70000/7000 -Jx.06/0.00006 -W10 -O -K>> p.ps
psxy /data/data4/sismo_chili/sismicutm/topo/t1-95.dz -R0/100/-8000/6000 -Jx.06/0.00060  -W10 -O -K>> p.ps
psxy <<! -Jx5 -R-71.7/-69/-26/-22.5 -O -K  -Y5 -Ba1>>p.ps
!
psxy trfault.xy -W5 -M -R -Jx -O -K >>p.ps
psxy /data/data4/sismo_chili/sismic.dat -Jx -R -Sc.1 -O -K >> p.ps
psxy /data/data4/sismo_chili/mainshock.xy -R -Jx -Sa.25 -G255/0/255 -W10 -O -K>> p.ps
pscoast -Jx -R288.3/291/-26/-22.5 -O -K -W5 -Df >>p.ps
psxy <<! -Jx -R-71.7/-69/-26/-22.5 -O -K >> p.ps
!
psvelomeca $fvecgps -Jx -R -A0.1/.3/.2 -Se0.03/.9/1  -O -K>>p.ps
psvelomeca gps_h_calc.geo -Jx -R  -Se0.03/.9/1 -L -W10/255/0/0 -O -K>>p.ps
psxy <<! -Jx -R -O >>p.ps
!
xv p.ps &

xyz2grd sour1.out -Gresu.grd -R1/$nd/1/$ns -I1  -b -Z
\rm sour1.out

psbasemap -R1/$nd/1/$ns -Ba1 -Jx3 -K -P > ! resu.ps
grdimage resu.grd -Cpals -O -K -R -Jx >> resu.ps
grdcontour resu.grd -C50 -A200 -O -K -Jx >> resu.ps
pstext slip.txt -R -Jx -O -K >> resu.ps
psxy <<! -Jx -R -O >> resu.ps
!
psxy trfault.xy -R-72/-70/-25.5/-22.5 -Jx8 -K -Ba1 -P -M  > ! resu2.ps
pstext slip.txt -R -Jx -O -K >> resu2.ps
psxy /data/data4/sismo_chili/sismic.dat -Jx -R -Sc.1 -G255/0/0 -O -K >> resu2.ps
psxy /data/data4/sismo_chili/mainshock.xy -R -Jx -Sa.25 -G255/0/255 -W10 -O -K>> resu2.ps
pscoast -R288/290/-25.5/-22.5 -Jx -O -K -W -Df >> resu2.ps
psxy <<! -R  -Jx -O >> resu2.ps
!
xv resu.ps &
xv resu2.ps &