# History file 'sal-vel.h' to extract gridded velocities from file # 'Saltf@@' from the "SEG/EAGE 3-D Salt Model" CD-ROM, and to prepare # velocities and slownesses for the inversion and for the testing. # ====================================================================== # Input files required chk.pl: "forms/" "copy.pl" chk.pl: "forms/" "echo.pl" chk.pl: "forms/" "inv.cal" chk.pl: "data/sal/" "sal-vel.cal" chk.pl: "" "Saltf@@" # File 'Saltf@@' should be extracted from the GZIPped file 'saltf.gz' # located in directory 'bin' using the command # gzip -d -N saltf.gz # Swaping the input velocity file (if needed) # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ IN='Saltf@@' ESIZE=4 N1=676 N2=676 N3=210 # File 'Saltf@@' is big endian. swap: # Comment the above line on big endian computers (workstations)! # Refresh file 'Saltf@@' or comment the line before running this history # file second time even on a low endian computer (PC)! # Dense 'chosen' velocity and slowness cubes # to fit in 64 Mbytes of memory # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ IN='Saltf@@' ESIZE=4 NDIG=0 N1=676 N2=676 N3=210 D1=20 D2=20 D3=20 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4180 GRD='s-vc040.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=1 NO2=1 NO3=1 # 0-13480 x 0-13480 x 0-4160 binasc: N1=338 N2=338 N3=105 D1=40 D2=40 D3=40 O1=0 O2=0 O3=0 copy.pl: "s-vc040.out" "s-sc040.out" echo.pl: "@1=1/@1" ">cal.tmp" CAL='cal.tmp' GRD1='s-sc040.out' grdcal: N1=676 N2=676 N3=210 D1=20 D2=20 D3=20 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4180 GRD='s-vc041.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=2 NO2=1 NO3=1 # 20-13500 x 0-13480 x 0-4160 #binasc: GRD='s-vc042.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=1 NO2=2 NO3=1 # 0-13480 x 20-13500 x 0-4160 #binasc: GRD='s-vc043.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=2 NO2=2 NO3=1 # 20-13500 x 20-13500 x 0-4160 #binasc: GRD='s-vc044.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=1 NO2=1 NO3=2 # 0-13480 x 0-13480 x 20-4180 #binasc: GRD='s-vc045.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=2 NO2=1 NO3=2 # 20-13500 x 0-13480 x 20-4180 #binasc: GRD='s-vc046.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=1 NO2=2 NO3=2 # 0-13480 x 20-13500 x 20-4180 #binasc: GRD='s-vc047.out' N1NEW=338 N2NEW=338 N3NEW=105 ND1=2 ND2=2 ND3=2 NO1=2 NO2=2 NO3=2 # 20-13500 x 20-13500 x 20-4180 #binasc: GRD='s-vc060.out' N1NEW=226 N2NEW=226 N3NEW=210 ND1=3 ND2=3 ND3=1 NO1=1 NO2=1 NO3=1 # 0-13500 x 0-13500 x 0-4180 binasc: # Sparse 'averaged' slowness cubes # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=338 N2=338 N3=105 D1=40 D2=40 D3=40 O1=0 O2=0 O3=0 GNORM=1 GRD='s-sc040.out' N1NEW=68 N2NEW=68 N3NEW=20 D1NEW=200 D2NEW=200 D3NEW=200 O1NEW=0 O2NEW=0 O3NEW=40 # 0-13400 x 0-13400 x 40-3840 GRDNEW='s-sa200.out' #grdnorm: N1NEW=46 N2NEW=46 N3NEW=14 D1NEW=300 D2NEW=300 D3NEW=300 O1NEW=0 O2NEW=0 O3NEW=100 # 0-13500 x 0-13500 x 100-4000 GRDNEW='s-sa300.out' #grdnorm: N1NEW=34 N2NEW=34 N3NEW=20 D1NEW=400 D2NEW=400 D3NEW=200 O1NEW=200 O2NEW=200 O3NEW=40 # 200-13400 x 200-13400 x 40-3840 GRDNEW='s-sa420.out' #grdnorm: N1NEW=28 N2NEW=28 N3NEW=9 D1NEW=500 D2NEW=500 D3NEW=500 O1NEW=0 O2NEW=0 O3NEW=0 # 0-13500 x 0-13500 x 0-4000 GRDNEW='s-sa500.out' #grdnorm: N1NEW=23 N2NEW=23 N3NEW=20 D1NEW=600 D2NEW=600 D3NEW=200 O1NEW=300 O2NEW=300 O3NEW=40 # 300-13500 x 300-13500 x 40-3840 GRDNEW='s-sa620.out' #grdnorm: N1NEW=23 N2NEW=23 N3NEW=11 D1NEW=600 D2NEW=600 D3NEW=400 O1NEW=300 O2NEW=300 O3NEW=0 # 300-13500 x 300-13500 x 0-4000 GRDNEW='s-sa640.out' #grdnorm: N1NEW=23 N2NEW=23 N3NEW=9 D1NEW=600 D2NEW=600 D3NEW=500 O1NEW=300 O2NEW=300 O3NEW=0 # 300-13500 x 300-13500 x 0-4000 GRDNEW='s-sa650.out' #grdnorm: N1NEW=18 N2NEW=18 N3NEW=41 D1NEW=800 D2NEW=800 D3NEW=100 O1NEW=-50 O2NEW=-50 O3NEW=0 # -50-13550 x -50-13550 x 0-4000 GRDNEW='s-sa810.out' #grdnorm: N1NEW=18 N2NEW=18 N3NEW=21 D1NEW=800 D2NEW=800 D3NEW=200 O1NEW=-50 O2NEW=-50 O3NEW=0 # -50-13550 x -50-13550 x 0-4000 GRDNEW='s-sa820.out' #grdnorm: N1NEW=18 N2NEW=18 N3NEW=17 D1NEW=800 D2NEW=800 D3NEW=250 O1NEW=-50 O2NEW=-50 O3NEW=0 # -50-13550 x -50-13550 x 0-4000 GRDNEW='s-sa825.out' grdnorm: N1NEW=18 N2NEW=18 N3NEW=14 D1NEW=800 D2NEW=800 D3NEW=300 O1NEW=-50 O2NEW=-50 O3NEW=0 # -50-13550 x -50-13550 x 0-3900 GRDNEW='s-sa830.out' #grdnorm: N1NEW=18 N2NEW=18 N3NEW=11 D1NEW=800 D2NEW=800 D3NEW=400 O1NEW=-50 O2NEW=-50 O3NEW=0 # -50-13550 x -50-13550 x 0-4000 GRDNEW='s-sa840.out' #grdnorm: N1NEW=14 N2NEW=14 N3NEW=21 D1NEW=1000 D2NEW=1000 D3NEW=200 O1NEW=250 O2NEW=250 O3NEW=0 # 250-13250 x 250-13250 x 0-4000 GRDNEW='s-sa000.out' #grdnorm: N1NEW=11 N2NEW=11 N3NEW=21 D1NEW=1340 D2NEW=1340 D3NEW=200 O1NEW=40 O2NEW=40 O3NEW=0 # 40-13440 x 40-13440 x 0-4000 GRDNEW='s-sa134.out' #grdnorm: # Sparse 'averaged' velocity cubes # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CAL='inv.cal' N1=68 N2=68 N3=20 D1=200 D2=200 D3=200 O1=0 O2=0 O3=40 # 0-13400 x 0-13400 x 40-3840 GRD1='s-sa200.out' GRD2='s-va200.out' #grdcal: N1=46 N2=46 N3=14 D1=300 D2=300 D3=300 O1=0 O2=0 O3=100 # 0-13500 x 0-13500 x 100-4000 GRD1='s-sa300.out' GRD2='s-va300.out' #grdcal: N1=34 N2=34 N3=20 D1=400 D2=400 D3=200 O1=200 O2=200 O3=40 # 200-13400 x 200-13400 x 40-3840 GRD1='s-sa420.out' GRD2='s-va420.out' #grdcal: N1=28 N2=28 N3=9 D1=500 D2=500 D3=500 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4000 GRD1='s-sa500.out' GRD2='s-va500.out' #grdcal: N1=23 N2=23 N3=20 D1=600 D2=600 D3=200 O1=300 O2=300 O3=40 # 300-13500 x 300-13500 x 40-3840 GRD1='s-sa620.out' GRD2='s-va620.out' #grdcal: N1=23 N2=23 N3=11 D1=600 D2=600 D3=400 O1=300 O2=300 O3=0 # 300-13500 x 300-13500 x 0-4000 GRD1='s-sa640.out' GRD2='s-va640.out' #grdcal: N1=23 N2=23 N3=9 D1=600 D2=600 D3=500 O1=300 O2=300 O3=0 # 300-13500 x 300-13500 x 0-4000 GRD1='s-sa650.out' GRD2='s-va650.out' #grdcal: N1=18 N2=18 N3=41 D1=800 D2=800 D3=100 O1=-50 O2=-50 O3=0 # -50-13550 x -50-13550 x 0-4000 GRD1='s-sa810.out' GRD2='s-va810.out' #grdcal: N1=18 N2=18 N3=21 D1=800 D2=800 D3=200 O1=-50 O2=-50 O3=0 # -50-13550 x -50-13550 x 0-4000 GRD1='s-sa820.out' GRD2='s-va820.out' #grdcal: N1=18 N2=18 N3=17 D1=800 D2=800 D3=250 O1=-50 O2=-50 O3=0 # -50-13550 x -50-13550 x 0-4000 GRD1='s-sa825.out' GRD2='s-va825.out' grdcal: N1=18 N2=18 N3=14 D1=800 D2=800 D3=300 O1=-50 O2=-50 O3=0 # -50-13550 x -50-13550 x 0-3900 GRD1='s-sa830.out' GRD2='s-va830.out' #grdcal: N1=18 N2=18 N3=11 D1=800 D2=800 D3=400 O1=-50 O2=-50 O3=0 # -50-13550 x -50-13550 x 0-4000 GRD1='s-sa840.out' GRD2='s-va840.out' #grdcal: N1=14 N2=14 N3=21 D1=1000 D2=1000 D3=200 O1=250 O2=250 O3=0 # 250-13250 x 250-13250 x 0-4000 GRD1='s-sa000.out' GRD2='s-va000.out' #grdcal: N1=11 N2=11 N3=21 D1=1340 D2=1340 D3=200 O1=40 O2=40 O3=0 # 40-13440 x 40-13440 x 0-4000 GRD1='s-sa134.out' GRD2='s-va134.out' #grdcal: # Sparse 'chosen' velocity and slowness cubes # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=676 N2=676 N3=210 D1=20 D2=20 D3=20 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4180 GRD='s-vc100.out' N1NEW=136 N2NEW=136 N3NEW=42 ND1=5 ND2=5 ND3=5 NO1=1 NO2=1 NO3=1 # 0-13500 x 0-13500 x 0-4100 binasc: N1=136 N2=136 N3=42 D1=100 D2=100 D3=100 O1=0 O2=0 O3=0 CAL='inv.cal' GRD1='s-vc100.out' GRD2='s-sc100.out' grdcal: N1=676 N2=676 N3=210 D1=20 D2=20 D3=20 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4180 GRD='s-vc101.out' N1NEW=136 N2NEW=136 N3NEW=42 ND1=5 ND2=5 ND3=5 NO1=1 NO2=1 NO3=3 # 0-13500 x 0-13500 x 40-4140 binasc: N1=676 N2=676 N3=210 D1=20 D2=20 D3=20 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4180 GRD='s-vc200.out' N1NEW=68 N2NEW=68 N3NEW=20 ND1=10 ND2=10 ND3=10 NO1=1 NO2=1 NO3=3 # 0-13400 x 0-13400 x 40-3840 #binasc: N1=68 N2=68 N3=20 D1=200 D2=200 D3=200 O1=0 O2=0 O3=40 CAL='inv.cal' GRD1='s-vc200.out' GRD2='s-sc200.out' #grdcal: N1=676 N2=676 N3=210 D1=20 D2=20 D3=20 O1=0 O2=0 O3=0 # 0-13500 x 0-13500 x 0-4180 GRD='s-vc810.out' N1NEW=17 N2NEW=17 N3NEW=39 ND1=40 ND2=40 ND3=5 NO1=21 NO2=21 NO3=3 # 400-13200 x 400-13200 x 40-3840 binasc: N1=17 N2=17 N3=39 D1=800 D2=800 D3=100 O1=400 O2=400 O3=40 CAL='inv.cal' GRD1='s-vc810.out' GRD2='s-sc810.out' grdcal: # Homogeneous second velocity derivatives # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=136 N2=136 N3=42 D1=100 D2=100 D3=100 O1=0 O2=0 O3=40 GRD='s-vc101.out' GRD1= GRD2= GRD3= GRD11='s-vd11.tmp' GRD22='s-vd22.tmp' GRD33='s-vd33.tmp' grdfd: N1NEW=68 N2NEW=68 N3NEW=20 D1NEW=200 D2NEW=200 D3NEW=200 O1NEW=0 O2NEW=0 O3NEW=40 GRD='s-vd11.tmp' GRDNEW='s-vd211.out' #grdnew: GRD='s-vd22.tmp' GRDNEW='s-vd222.out' #grdnew: GRD='s-vd33.tmp' GRDNEW='s-vd233.out' #grdnew: N1NEW=17 N2NEW=17 N3NEW=39 D1NEW=800 D2NEW=800 D3NEW=100 O1NEW=400 O2NEW=400 O3NEW=40 GRD='s-vd11.tmp' GRDNEW='s-vd811.out' grdnew: GRD='s-vd22.tmp' GRDNEW='s-vd822.out' grdnew: GRD='s-vd33.tmp' GRDNEW='s-vd833.out' grdnew: # Removing artificial reflectors from the velocity cubes # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=68 N2=68 N3=20 D1=200 D2=200 D3=200 O1=0 O2=0 O3=20 CAL='sal-vel.cal' SALTREF=0.15 GRD1='s-vd211.out' GRD2='s-vd222.out' GRD3='s-vd233.out' GRD4='s-vc200.out' GRD5='s-vc200r.out' D1=100 D2=100 D3=100 # FD grid intervals for 'sal-vel.cal' #grdcal: D1=200 D2=200 D3=200 CAL='inv.cal' GRD1='s-vc200r.out' GRD2='s-sc200r.out' #grdcal: N1=17 N2=17 N3=39 D1=800 D2=800 D3=100 O1=400 O2=400 O3=40 CAL='sal-vel.cal' SALTREF=0.15 GRD1='s-vd811.out' GRD2='s-vd822.out' GRD3='s-vd833.out' GRD4='s-vc810.out' GRD5='s-vc810r.out' D1=100 D2=100 D3=100 # FD grid intervals for 'sal-vel.cal' grdcal: D1=800 D2=800 D3=100 CAL='inv.cal' GRD1='s-vc810r.out' GRD2='s-sc810r.out' grdcal: N1=68 N2=68 N3=20 D1=200 D2=200 D3=200 O1=0 O2=0 O3=40 CAL='sal-vel.cal' SALTREF=0.15 GRD1='s-vd211.out' GRD2='s-vd222.out' GRD3='s-vd233.out' GRD4='s-va200.out' GRD5='s-va200r.out' D1=100 D2=100 D3=100 # FD grid intervals for 'sal-vel.cal' #grdcal: D1=200 D2=200 D3=200 CAL='inv.cal' GRD1='s-va200r.out' GRD2='s-sa200r.out' #grdcal: # Computing average P wave slowness and average P wave velocity # (harmonic average) # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=338 N2=338 N3=105 D1=40 D2=40 D3=40 O1=0 O2=0 O3=0 N1NEW=1 N2NEW=1 N3NEW=1 GNORM=0 GRD='s-sc040.out' GRDNEW='s-avslo.out' grdnorm: N1=1 N2=1 N3=1 CAL='inv.cal' GRD1='s-avslo.out' GRD2='s-avvel.out' grdcal: # List of the output files # ~~~~~~~~~~~~~~~~~~~~~~~~ # original grid 20 x 20 x 20 m, 0-13500 x 0-13500 x 0-4180 m # # 'Chosen' velocities # 's-vc040.out' 40 x 40 x 40 m, 0-13480 x 0-13480 x 0-4160 m # 's-vc041.out' 40 x 40 x 40 m, 20-13500 x 0-13480 x 0-4160 m # 's-vc042.out' 40 x 40 x 40 m, 0-13480 x 20-13500 x 0-4160 m # 's-vc043.out' 40 x 40 x 40 m, 20-13500 x 20-13500 x 0-4160 m # 's-vc044.out' 40 x 40 x 40 m, 0-13480 x 0-13480 x 20-4180 m # 's-vc045.out' 40 x 40 x 40 m, 20-13500 x 0-13480 x 20-4180 m # 's-vc046.out' 40 x 40 x 40 m, 0-13480 x 20-13500 x 20-4180 m # 's-vc047.out' 40 x 40 x 40 m, 20-13500 x 20-13500 x 20-4180 m # 's-vc060.out' 60 x 60 x 20 m, 0-13500 x 0-13500 x 0-4180 m # 's-vc100.out' 100 x 100 x 100 m, 0-13500 x 0-13500 x 0-4100 m # 's-vc101.out' 100 x 100 x 100 m, 0-13500 x 0-13500 x 40-4140 m # 's-vc200.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # 's-vc810.out' 800 x 800 x 100 m, 400-13200 x 400-13200 x 40-3840 m # # 'Chosen' slownesses # 's-sc040.out' 40 x 40 x 40 m, 0-13480 x 0-13480 x 0-4160 m # 's-sc100.out' 100 x 100 x 100 m, 0-13500 x 0-13500 x 0-4100 m # 's-sc200.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # 's-sc810.out' 800 x 800 x 100 m, 400-13200 x 400-13200 x 40-3840 m # # 'Averaged' slownesses # 's-sa200.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # 's-sa300.out' 300 x 300 x 300 m, 0-13500 x 0-13500 x 100-4000 m # 's-sa420.out' 400 x 400 x 200 m, 200-13400 x 200-13400 x 40-3840 m # 's-sa500.out' 500 x 500 x 500 m, 0-13500 x 0-13500 x 0-4000 m # 's-sa620.out' 600 x 600 x 200 m, 300-13500 x 300-13500 x 40-3840 m # 's-sa640.out' 600 x 600 x 400 m, 300-13500 x 300-13500 x 0-4000 m # 's-sa650.out' 600 x 600 x 500 m, 300-13500 x 300-13500 x 0-4000 m # 's-sa810.out' 800 x 800 x 100 m, -50-13550 x -50-13550 x 0-4000 m # 's-sa820.out' 800 x 800 x 200 m, -50-13550 x -50-13550 x 0-4000 m # 's-sa825.out' 800 x 800 x 250 m, -50-13550 x -50-13550 x 0-4000 m # 's-sa830.out' 800 x 800 x 300 m, -50-13550 x -50-13550 x 0-3900 m # 's-sa840.out' 800 x 800 x 400 m, -50-13550 x -50-13550 x 0-4000 m # 's-sa000.out' 1000 x 1000 x 200 m, 250-13250 x 250-13250 x 0-4000 m # 's-sa134.out' 1340 x 1340 x 200 m, 40-13440 x 40-13440 x 0-4000 m # # 'Averaged' velocities # 's-va200.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # 's-va300.out' 300 x 300 x 300 m, 0-13500 x 0-13500 x 100-4000 m # 's-va420.out' 400 x 400 x 200 m, 200-13400 x 200-13400 x 40-3840 m # 's-va500.out' 500 x 500 x 500 m, 0-13500 x 0-13500 x 0-4000 m # 's-va620.out' 600 x 600 x 200 m, 300-13500 x 300-13500 x 40-3840 m # 's-va640.out' 600 x 600 x 400 m, 300-13500 x 300-13500 x 0-4000 m # 's-va650.out' 600 x 600 x 500 m, 300-13500 x 300-13500 x 0-4000 m # 's-va810.out' 800 x 800 x 100 m, -50-13550 x -50-13550 x 0-4000 m # 's-va820.out' 800 x 800 x 200 m, -50-13550 x -50-13550 x 0-4000 m # 's-va825.out' 800 x 800 x 250 m, -50-13550 x -50-13550 x 0-4000 m # 's-va830.out' 800 x 800 x 300 m, -50-13550 x -50-13550 x 0-3900 m # 's-va840.out' 800 x 800 x 400 m, -50-13550 x -50-13550 x 0-4000 m # 's-va000.out' 1000 x 1000 x 200 m, 250-13250 x 250-13250 x 0-4000 m # 's-va134.out' 1340 x 1340 x 200 m, 40-13440 x 40-13440 x 0-4000 m # # 'Chosen' velocities without artificial reflectors # 's-vc200r.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # 's-vc810r.out' 800 x 800 x 100 m, 400-13200 x 400-13200 x 40-3840 m # # 'Chosen' slownesses without artificial reflectors # 's-sc200r.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # 's-sc810r.out' 800 x 800 x 100 m, 400-13200 x 400-13200 x 40-3840 m # # 'Averaged' velocities without artificial reflectors # 's-va200r.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # # 'Averaged' slownesses without artificial reflectors # 's-sa200r.out' 200 x 200 x 200 m, 0-13400 x 0-13400 x 40-3840 m # # # Average slowness and velocity over the grid 40x40x40 m # 's-avslo.out' 's-avvel.out' #