# History file 'sal-inv3.h' to fit the SEG/EAGE 3-D Salt Model by a model # with interfaces. Input data for the salt interface are taken from the # GOCAD surfaces, data from file 'fault_a.ts' are also used. # Dimension of array RAM in include file 'ram.inc' should be # MRAM=20000000 or greater for this history file. # ===================================================================== # Input files required chk.pl: "data/sal/" "sal-m3.dat" chk.pl: "data/sal/" "sal-icb.cal" chk.pl: "model/" "sob22.dat" # Output of 'sal-ts.h' chk.pl: "" "top.pts" chk.pl: "" "hrz_1.pts" chk.pl: "" "gp2.pts" chk.pl: "" "bottom.pts" chk.pl: "" "fault_a.pts" # Output of 'sal-salt.h' chk.pl: "" "salt-t1.pts" chk.pl: "" "salt-t2.pts" chk.pl: "" "salt-b1.pts" # Output of 'sal-vel.h' chk.pl: "" "s-vc810r.out" chk.pl: "" "s-sc810r.out" # Initial model and its consistency check # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='sal-m3.dat' MODIN='sal-m3.dat' N1=20 N2=20 N3=10 LFREE=1 MODLOG='modchk.out' #modchk: # Form of the files with matrices # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ FORMM='unformatted' MATFORM='unformatted' # Inversion of interfaces # ~~~~~~~~~~~~~~~~~~~~~~~ ICLASS=1 # 646 model parameters for surfaces # Calculating matrices for inversion M1='m1.out' M2='m2.out' SOBOLEV='sob22.dat' MODSOB='modsob.out' SOBW00=1 invsoft: MODSOB= SOBW00= GM1='gm1.out' GM2='gm2.out' GM3=' ' DM1='dm1.out' MPAR=0 PTS='top.pts' INDFUN=1 ERRMUL=2.000 invpts: # 4 points M2IN='m2.out' # PTS='hrz_1.pts' INDFUN=2 ERRMUL=24.52 invpts: # 601 points PTS='gp2.pts' INDFUN=3 ERRMUL=42.19 invpts: # 1780 points PTS='bottom.pts' INDFUN=4 ERRMUL=2.000 invpts: # 4 points PTS='salt-t1.pts' INDFUN=5 ERRMUL=28.16 invpts: # 793 points PTS='salt-b1.pts' INDFUN=6 ERRMUL=24.02 invpts: # 577 points PTS='fault_a.pts' INDFUN=7 ERRMUL=31.91 invpts: # 1018 points PTS='salt-t2.pts' INDFUN=8 ERRMUL=14.04 invpts: # 197 points M2IN=' ' PTS=' ' INDFUN= ERRMUL= # total # 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='modsob.out' MATOUT='sm2.out' COEF1= COEF2=400000000 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='sal-m3.tmp' modmod: # Indices of complex blocks relevant to the velocities # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=68 N2=68 N3=20 D1=200 D2=200 D3=200 O1=0 O2=0 O3=40 N1=46 N2=46 N3=14 D1=300 D2=300 D3=300 O1=0 O2=0 O3=100 N1=34 N2=34 N3=20 D1=400 D2=400 D3=200 O1=200 O2=200 O3=40 N1=23 N2=23 N3=21 D1=600 D2=600 D3=200 O1=200 O2=200 O3=100 N1=17 N2=17 N3=20 D1=800 D2=800 D3=200 O1=400 O2=400 O3=40 N1=11 N2=11 N3=20 D1=1340 D2=1340 D3=200 O1=40 O2=40 O3=40 N1=17 N2=17 N3=39 D1=800 D2=800 D3=100 O1=400 O2=400 O3=40 ISRF=3 VEL='sal-s3.out' grid: ISRF=4 VEL='sal-s4.out' grid: ISRF= CAL='sal-icb.cal' GRD2='sal-s3.out' GRD3='sal-s4.out' GRD4='sal-icb.out' SALTV1=1000 SALTV2=1501 SALTV5=4449 SALTD3=152.4 SALTD4=1 GRD1='s-va200.out' GRD1='s-vc200.out' GRD1='s-va200r.out' GRD1='s-vc200r.out' GRD1='s-va300.out' GRD1='s-va420.out' GRD1='s-va620.out' GRD1='s-va820.out' GRD1='s-va134.out' GRD1='s-vc810r.out' grdcal: # Inversion of velocities # ~~~~~~~~~~~~~~~~~~~~~~~ MODEL='sal-m3.tmp' MODIN='sal-m3.tmp' ICLASS=2 NEGPAR=1 # Calculating matrices for inversion M1='m1.out' M2='m2.out' SOBOLEV='sob22.dat' MODSOB='modsob.out' SOBW01=1 invsoft: MODSOB= SOBW01= GM1='gm1.out' GM2='gm2.out' GM3=' ' DM1='dm1.out' GRDICB='sal-icb.out' MPAR=1 POWERM=-1 GRD='s-sa200.out' ERRMUL=304.11 #92480 points GRD='s-sc200.out' ERRMUL=304.11 #92480 points GRD='s-sa200r.out' ERRMUL=304.11 #92480 points GRD='s-sc200r.out' ERRMUL=304.11 #92480 points GRD='s-sa300.out' ERRMUL=172.12 #29624 points GRD='s-sa420.out' ERRMUL=152.05 #23120 points GRD='s-sa620.out' ERRMUL=105.40 #11109 points GRD='s-sa820.out' ERRMUL=76.026 #5780 points GRD='s-sa134.out' ERRMUL=50.408 #2541 points GRD='s-sc810r.out' ERRMUL=106.16 #11271 points invpts: GRD= GRDICB= 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='modsob.out' MATOUT='sm2.out' COEF1= COEF2=25000000 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='sal-m3.out' modmod: # Sobolev norm of the slowness in the updated model # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MATIN1='modsob.out' MATIN2='gm5.out' MATOUT='gm6.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='gm5.out' MATIN2='gm6.out' MATOUT='gm7.out' SYMMETRY='diag' MATT1=1 MATT2=0 matmul: MATIN1='gm7.out' MATOUT='s-sns.out' SPARSE=-1 MATFUN='sqrt' MATFORM='formatted' matfun: # List of output files # ~~~~~~~~~~~~~~~~~~~~ # 'sal-m3.tmp' ... Inverted model after inversion of interfaces. # 'sal-m3.out' ... Inverted model after inversions of interfaces and # of velocities. # 's-sns.ou@' ... Sobolev norm of the slowness in the model.