# History file 'u2d-inv.h' to run 'invsoft.for' and 'invtt.for' # with demo data in 2-D model UNCONFORMITY # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Input files required: #chk.pl: "data/u2d/" "u2d-mod.dat" #chk.pl: "data/u2d/" "u2d-src.dat" #chk.pl: "data/u2d/" "u2d-rec.dat" #chk.pl: "data/u2d/" "u2d-rpa.dat" #chk.pl: "data/u2d/" "u2d-wri.dat" #chk.pl: "data/u2d/" "u2d-pts.dat" #chk.pl: "data/u2d/" "u2d-ftt.dat" #chk.pl: "data/u2d/" "u2d-crt1.dat" #chk.pl: "model/" "sob22.dat" # Initial model (model UNCONFORMITY by Cormier and Mellen 1984) MODEL='u2d-mod.dat' # Calculating L2 and Sobolev scalar products of basis functions # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Data describing Sobolev scalar products of basis functions SOBOLEV='sob22.dat' # Minimizing second derivatives M1='u2d-m1.out' # Output number of basis functions # L2 scalar products of surface, velocity and density basis functions MODL2='u2d-l2.out' MODSOB=' ' SOBW00=1 SOBW01=1 SOBW02=1 SOBW03=1 invsoft: 'u2d-inv.h' / # Sobolev scalar products of surface basis functions MODL2=' ' MODSOB='u2d-sobs.out' SOBW00=1 SOBW01=0 SOBW02=0 SOBW03=0 invsoft: 'u2d-inv.h' / # Sobolev scalar products of velocity basis functions MODL2=' ' MODSOB='u2d-sobv.out' SOBW00=0 SOBW01=1 SOBW02=1 SOBW03=0 invsoft: 'u2d-inv.h' / # Ray tracing # ~~~~~~~~~~~ # Source and receivers SRC='u2d-src.dat' REC='u2d-rec.dat' # Data filenames for CRT (for details see 'crt.for') DCRT='u2d-crt1.dat' CODE='u2d-crt1.dat' RPAR='u2d-rpa.dat' WRIT='u2d-wri.dat' CRTLOG='u2d-log.out' # Kind of ray parameters (for details see 'init.for') INIPAR=2 # geographic-like spherical coordinates # 2-D two-point ray tracing crt: 'u2d-inv.h' / # Travel-time derivatives with respect to basis functions # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ CRTOUT=' ' # default names of output files of program CRT PTS='u2d-pts.dat' # coordinates of shot and receiver points FTT='u2d-ftt.dat' # travel times to be inverted DIST=0.0001 # maximum distance between the point and rays VPOWER=1 # velocity power for the travel-time check sum M1=' ' # output file with number of model coefficients M2='u2d-m2.out' # output file with number of values to invert GM1='u2d-gm1.out' # matrix of the derivatives of travel times GM2='u2d-gm2.out' # vector of travel-time residuals GM3='u2d-gm3.out' # vector of travel times in the model DM1='u2d-dm1.out' # diagonal matrix of travel-time variances INVLOG='u2d-itt.out' # output log file # Generating system of equations for travel-time inversion invtt: 'u2d-inv.h' / # ====================================================================== # Output files of 'invsoft.for': # 'u2d-m1.out' Number of basis functions # 'u2d-l2.out' L2 scalar products of all basis functions # 'u2d-sobs.out' Sobolev scalar products of surface basis functions # 'u2d-sobv.out' Sobolev scalar products of velocity basis functions # Temporary files (may be deleted): # 'r0*.out' # 's0*.out' # Output files of 'invtt.for': # 'u2d-m2.out' # 'u2d-gm1.out' # 'u2d-gm2.out' # 'u2d-gm3.out' # 'u2d-dm1.out' # 'u2d-itt.out'