# Controlled initial-value ray tracing in the model with # a LENticular inclusion followed by Interpolation of travel times # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Input files required #chk.pl: "data/len/" "len-mod.dat" #chk.pl: "data/len/" "leni-crt.dat" #chk.pl: "data/len/" "leni-rpa.dat" #chk.pl: "data/len/" "len-src.dat" #chk.pl: "crt/" "codep.dat" #chk.pl: "crt/" "writall.dat" # Model and source MODEL='len-mod.dat' SRC='len-src.dat' # Data filenames for CRT (for details see 'crt.for') DCRT='leni-crt.dat' RPAR='leni-rpa.dat' CODE='codep.dat' # just refracted P wave WRIT='writall.dat' # all rays stored for interpolation CRTLOG='leni-log.out' # Kind of ray parameters (for details see 'init.for') INIPAR=2 # geographic-like spherical coordinates # Grid dimensions N1=601 N2=1 N3=301 O1=0 O2=10 O3=-25 D1=0.1 D2=1 D3=0.1 # Output files of program 'mtt.for' NUM='len-num.out' MTT='len-mtt.out' MHIST='len-his.out' # Ray tracing of refracted wave crt: 'leni-mtt.h' / # Interpolation of travel times within ray cells mtt: 'leni-mtt.h' / # Plotting numbers of arrivals at the gridpoints VPLUS=2 CREF=0.166667 GRD='len-num.out' PS='len-num.ps' grdps: 'leni-mtt.h' / # Plotting travel times at the gridpoints, sorted according to their value VCIRC=1 VREF=0 CREF=0.166667 NUM='len-num.out' MGRD='len-mtt.out' GRD='len-t00.out' mgrd: 'leni-mtt.h' / GRD='len-t00.out' PS='len-t01.ps' grdps: 'leni-mtt.h' / # ====================================================================== # Temporary files (may be deleted) # '*.tmp' # 'r0*.out' # 's0*.out' # 't0*.out' # List of output files # 'len-num.out' Numbers of travel times # 'len-his.out' Multivalued ray histories # 'len-mtt.out' Multivalued travel times # 'len-t00.out' Travel times separated into N4 univalued slices # List of output colour PostScript figures # 'len-num.ps' Numbers of travel times at the gridpoints. # 'len-t**.ps' Travel times, sorted according to their size # (e.g., file 'leni-t01.ps' contains the travel times # of the fastest arrival).