# Calculating velocity differences in models P1 and P1I # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Input files required # ~~~~~~~~~~~~~~~~~~~~ chk.pl: "forms/" "sub.cal" chk.pl: "data/p1/" "p1-vp.dat" # output of 'p1-inv.h' chk.pl: "" "p1-vs.out" chk.pl: "" "p1-mod.dat" chk.pl: "" "p1-mods.dat" # output of 'p1i-inv.h' chk.pl: "" "p1i-mod.dat" # Plot of the gridded velocity data # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ VREF=5.45 CREF=0.1666667 VCIRC=-1.5 HSIZE=24 YSIGN=-1 ROTATE=-90 VOFFSET=13 O1=0. D1=0.1 N1=474 O3=0. D3=0.1 N3=61 GRD='p1-vp.dat' PS='p1-vp.ps' grdps: VREF=3.14 CREF=0.096225 VCIRC=-0.866 GRD='p1-vs.out' PS='p1-vs.ps' grdps: # Plot of the gridded velocities in the smooth model # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='p1-mod.dat' VREF=5.45 CREF=0.1666667 VCIRC=-1.5 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 MPAR=1 VEL='p1-vep.out' grid: GRD='p1-vep.out' PS='p1-vep.ps' grdps: MODEL='p1-mods.dat' VREF=3.14 CREF=0.096225 VCIRC=-0.866 MPAR=2 VEL='p1-ves.out' grid: GRD='p1-ves.out' PS='p1-ves.ps' grdps: # Plot of the gridded velocities in the model with interfaces # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='p1i-mod.dat' VREF=5.45 CREF=0.1666667 VCIRC=-1.5 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 MPAR=1 VEL='p1i-vep.out' grid: GRD='p1i-vep.out' PS='p1i-vep.ps' grdps: VREF=3.14 CREF=0.096225 VCIRC=-0.866 MPAR=2 VEL='p1i-ves.out' grid: GRD='p1i-ves.out' PS='p1i-ves.ps' grdps: O3= D3= N3= # Subtracting the velocities # ~~~~~~~~~~~~~~~~~~~~~~~~~~ # Smooth model - model with intefaces CAL='sub.cal' N1=474 N2=61 N3= GRD1='p1-vep.out' GRD2='p1i-vep.out' GRD3='p1-difvp.out' grdcal: GRD1= GRD2= GRD3= VREF=0 CREF=0.1666667 VCIRC=0.4 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 GRD='p1-difvp.out' PS='p1-difvp.ps' grdps: N1=474 N2=61 N3= GRD1='p1-ves.out' GRD2='p1i-ves.out' GRD3='p1-difvs.out' grdcal: GRD1= GRD2= GRD3= VREF=0 CREF=0.1666667 VCIRC=0.23 #=0.4/sqrt(3) O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 GRD='p1-difvs.out' PS='p1-difvs.ps' grdps: # Smooth model - data GRD1='p1-vep.out' GRD2='p1-vp.dat' GRD3='p1-dvp.out' N1=474 N2=61 N3= grdcal: GRD1= GRD2= GRD3= VREF=0 CREF=0.1666667 VCIRC=0.4 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 GRD='p1-dvp.out' PS='p1-dvp.ps' grdps: GRD1='p1-ves.out' GRD2='p1-vs.out' GRD3='p1-dvs.out' N1=474 N2=61 N3= grdcal: GRD1= GRD2= GRD3= VREF=0 CREF=0.1666667 VCIRC=0.23 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 GRD='p1-dvs.out' PS='p1-dvs.ps' grdps: # Model with interfaces - data GRD1='p1i-vep.out' GRD2='p1-vp.dat' GRD3='p1i-dvp.out' N1=474 N2=61 N3= grdcal: GRD1= GRD2= GRD3= VREF=0 CREF=0.1666667 VCIRC=0.4 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 GRD='p1i-dvp.out' PS='p1i-dvp.ps' grdps: GRD1='p1i-ves.out' GRD2='p1-vs.out' GRD3='p1i-dvs.out' N1=474 N2=61 N3= grdcal: GRD1= GRD2= GRD3= VREF=0 CREF=0.1666667 VCIRC=0.23 O1=0. D1=0.1 N1=474 O2= D2= N2= O3=0. D3=0.1 N3=61 GRD='p1i-dvs.out' PS='p1i-dvs.ps' grdps: