# Constructing model P1I with interfaces # ====================================== # Input files required # ~~~~~~~~~~~~~~~~~~~~ chk.pl: "data/p1/" "p1-vp.dat" chk.pl: "data/p1/" "p1-vpvs.dat" chk.pl: "data/p1/" "p1i-mod0.dat" chk.pl: "data/p1/" "p1-if2.dat" chk.pl: "data/p1/" "p1-if9.dat" chk.pl: "data/p1/" "p1-if3.dat" chk.pl: "data/p1/" "p1-il4.dat" chk.pl: "data/p1/" "p1-il5.dat" chk.pl: "data/p1/" "p1-il6.dat" chk.pl: "data/p1/" "p1-il7.dat" chk.pl: "data/p1/" "p1-if6.dat" chk.pl: "data/p1/" "p1-if7.dat" chk.pl: "data/p1/" "p1-il8.dat" chk.pl: "data/p1/" "p1-if8.dat" chk.pl: "data/p1/" "p1-il9.dat" chk.pl: "data/p1/" "p1-il10.dat" chk.pl: "data/p1/" "p1-bnd.dat" chk.pl: "model/" "sob22.dat" chk.pl: "forms/" "inv.cal" chk.pl: "forms/" "mul.cal" chk.pl: "forms/" "echo.pl" chk.pl: "forms/" "append.pl" # Converting velocity to slowness # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CAL='inv.cal' N1=474 N2=61 N3= GRD1='p1-vp.dat' GRD2='p1-sp.out' grdcal: N1= N2= N3= GRD1= GRD2= CAL='mul.cal' N1=474 N2=61 N3= GRD1='p1-vpvs.dat' GRD2='p1-sp.out' GRD3='p1-ss.out' grdcal: GRD1= GRD2= GRD3= N1= N2= N3= # Constructing velocity model with interfaces # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Inversion of the interfaces # ~~~~~~~~~~~~~~~~~~~~~~~~~~~ ICLASS=1 # Initial model, name of the resulting model MODEL='p1i-mod0.dat' NEGPAR=0 # Calculating matrices for inversion M1='m1.out' M2='m2.out' SOBOLEV='sob22.dat' MODSOB='modsobi.out' SOBW00=1 invsoft: MODSOB= SOBW01= GM1='gm1.out' GM2='gm2.out' GM3=' ' DM1='dm1.out' PTS='p1-if2.dat' INDFUN=1 ERRMUL=3.464 invpts: # 12 points M2IN='m2.out' # PTS='p1-if9.dat' INDFUN=2 ERRMUL=1.732 invpts: # 3 points PTS='p1-if3.dat' INDFUN=3 ERRMUL=2.646 invpts: # 7 points PTS='p1-il4.dat' INDFUN=4 ERRMUL=2.646 invpts: # 7 points PTS='p1-il5.dat' INDFUN=5 ERRMUL=4.472 invpts: # 20 points PTS='p1-il6.dat' INDFUN=6 ERRMUL=2.646 invpts: # 7 points PTS='p1-il7.dat' INDFUN=7 ERRMUL=3.162 invpts: # 10 points PTS='p1-if6.dat' INDFUN=8 ERRMUL=2.236 invpts: # 5 points PTS='p1-if7.dat' INDFUN=9 ERRMUL=2. invpts: # 4 points PTS='p1-il8.dat' INDFUN=10 ERRMUL=1.732 invpts: # 3 points PTS='p1-if8.dat' INDFUN=11 ERRMUL=2.236 invpts: # 5 points PTS='p1-il9.dat' INDFUN=12 ERRMUL=2. invpts: # 4 points PTS='p1-il10.dat' INDFUN=13 ERRMUL=1.732 invpts: # 3 points M2IN= PTS= INDFUN= ERRMUL= # Matrix operations MATIN1='dm1.out' MATOUT='dm2.out' MATFUN='inv' matfun: MATIN1='dm2.out' MATIN2='gm1.out' MATOUT='dm2gm1.out' SYMMETRY= MATT1= MATT2=1 matmul: MATIN1='gm1.out' MATIN2='dm2gm1.out' MATOUT='sm1.out' SYMMETRY='sym' MATT1= MATT2= matmul: MATIN1='sm1.out' MATIN2='modsobi.out' MATOUT='sm2.out' # Smoothing interfaces by adding COEF2*modsobi.out COEF1= COEF2=0.010 matlin: MATIN1='sm2.out' MATOUT='sm3.out' matinv: MATIN1='dm2.out' MATIN2='gm2.out' MATOUT='gm3.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='gm1.out' MATIN2='gm3.out' MATOUT='gm4.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='sm3.out' MATIN2='gm4.out' MATOUT='gm5.out' SYMMETRY=' ' MATT1= MATT2= matmul: # Updating the model M1='m1.out' MODNEW='gm5.out' MODOUT='p1i-mod.tmp' modmod: M1= M2= MODNEW= MODOUT= MODIN= # Plot of the data for interfaces # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ PICDAT='picdat.tmp' PICTURE='p1i-pts.ps' echo.pl: "'PROJECTION' " "> picdat.tmp" echo.pl: "16.35 27.2 0. -0.5 0. 0. -0.5 0. / " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'ATTRIBUTES' " ">>picdat.tmp" echo.pl: "'MSZSF' 0.50 / Marker size scale factor. " ">>picdat.tmp" echo.pl: "'ITC' 0 / No texts drawn. " ">>picdat.tmp" echo.pl: "'LWIDTH' 0.50 / " ">>picdat.tmp" echo.pl: "/ " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'LINES' " ">>picdat.tmp" echo.pl: "'p1-bnd.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-if2.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-if3.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-if9.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il4.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il5.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il6.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il7.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-if6.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-if7.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il8.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-if8.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il9.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" echo.pl: "'POINTS' " ">>picdat.tmp" echo.pl: "'p1-il10.dat' " ">>picdat.tmp" echo.pl: " " ">>picdat.tmp" pictures: # Plot of the complex blocks indices together with the data for interfaces # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='p1i-mod.tmp' O2= D2= N2= O1=0. D1=0.01 N1=4731 O3=0. D3=0.01 N3=601 O1=0. D1=0.05 N1=947 O3=0. D3=0.05 N3=121 MPAR=1 VEL=' ' ICB='p1i-icb.out' grid: VREF=7.5 CREF=0.166667 VCIRC=-5.5 HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35 HOFFSET= GRD='p1i-icb.out' PS='p1i-icb.ps' SHOWPAGE=0 grdps: append.pl: "p1i-icb.ps" "p1i-pts.ps" # Indices of complex blocks relevant to the velocities # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='p1i-mod.tmp' O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 MPAR=1 VEL=' ' ICB='p1i-icb.out' grid: # Inversion of velocities # ~~~~~~~~~~~~~~~~~~~~~~~ MODEL='p1i-mod.tmp' MODIN='p1i-mod.tmp' ICLASS=2 # Calculating matrices for inversion M1='m1.out' M2='m2.out' SOBOLEV='sob22.dat' MODSOB='modsobp.out' SOBW01=1 invsoft: MODSOB= SOBW01= SOBOLEV='sob22.dat' MODSOB='modsobs.out' SOBW02=1 invsoft: MODSOB= SOBW02= GM1='gm1.out' GM2='gm2.out' GM3=' ' DM1='dm1.out' O1=0. D1=0.1 N1=474 O3=0. D3=0.1 N3=61 PTS= LIN= GRD='p1-sp.out' GRDICB='p1i-icb.out' INDFUN= MPAR=1 POWERM=-1 ERRMUL=170. # SQRT(28914) invpts: M2IN='m2.out' PTS= LIN= GRD='p1-ss.out' GRDICB='p1i-icb.out' INDFUN= MPAR=2 POWERM=-1 ERRMUL=170. # SQRT(28914) invpts: GRD= GRDICB= INDFUN= MPAR= POWERM= ERRMUL= # Matrix operations MATIN1='dm1.out' MATOUT='dm2.out' MATFUN='inv' matfun: MATIN1='dm2.out' MATIN2='gm1.out' MATOUT='dm2gm1.out' SYMMETRY= MATT1= MATT2=1 matmul: MATIN1='gm1.out' MATIN2='dm2gm1.out' MATOUT='sm1.out' SYMMETRY='sym' MATT1= MATT2= matmul: MATIN1='sm1.out' MATIN2='modsobp.out' MATOUT='sm2a.out' # Smoothing P-wave velocities by adding COEF2*modsobp.out COEF1= COEF2=0.001 matlin: MATIN1='sm2a.out' MATIN2='modsobs.out' MATOUT='sm2.out' # Smoothing S-wave velocities by adding COEF2*modsobs.out COEF1= COEF2=0.001 matlin: MATIN1='sm2.out' MATOUT='sm3.out' matinv: MATIN1='dm2.out' MATIN2='gm2.out' MATOUT='gm3.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='gm1.out' MATIN2='gm3.out' MATOUT='gm4.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='sm3.out' MATIN2='gm4.out' MATOUT='gm5.out' SYMMETRY=' ' MATT1= MATT2= matmul: # Updating the model M1='m1.out' MODNEW='gm5.out' MODOUT='p1i-mod.dat' modmod: M1= M2= MODNEW= MODOUT= MODIN= # Plot of the gridded velocities in the model # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='p1i-mod.dat' O1=0. D1=0.1 N1=474 O3=0. D3=0.1 N3=61 O2= D2= N2= O1=0. D1=0.01 N1=4731 O3=0. D3=0.01 N3=601 O1=0. D1=0.05 N1=947 O3=0. D3=0.05 N3=121 MPAR=1 VEL='p1i-vep.out' ICB= grid: VREF=5.45 CREF=0.1666667 VCIRC=-1.5 HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35 HOFFSET= GRD='p1i-vep.out' PS='p1i-vep.ps' SHOWPAGE= grdps: VREF=3.14 CREF=0.096225 VCIRC=-0.866 MPAR=2 VEL='p1i-ves.out' ICB= grid: GRD='p1i-ves.out' PS='p1i-ves.ps' SHOWPAGE= grdps: