EXAMPLE OF RUNNING RANDALL REFLECTIVITY CODE Executable versions of the programs are contained in ~shearer/CLASS/227C/REFL mijkennett computes response for all 6 moment tensor excitations mijfkkennett computes response for 6 mom. tensor + 3 point source mgenkennett computes response for specific moment tensor xtsynth creates seismograms from output of above programs Here is an example of an input file to mgenkennett (in.kennett) ---------------------------------------------------------------------- 42 1.0 4.845 2.795 2.290 200. 50. 10000. 0.0001 10000. 0.0001 1.0 5.136 2.964 2.339 200. 50. 10000. 0.0001 10000. 0.0001 1.0 5.427 3.133 2.388 200. 50. 10000. 0.0001 10000. 0.0001 1.0 5.718 3.302 2.437 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.009 3.471 2.486 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.300 3.640 2.535 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.315 3.649 2.547 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.330 3.657 2.558 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.344 3.666 2.570 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.359 3.674 2.582 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.374 3.683 2.593 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.389 3.691 2.605 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.404 3.700 2.617 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.419 3.708 2.628 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.433 3.717 2.640 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.448 3.725 2.652 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.463 3.734 2.663 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.478 3.742 2.675 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.493 3.751 2.687 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.507 3.759 2.698 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.522 3.768 2.710 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.537 3.776 2.722 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.552 3.785 2.733 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.567 3.793 2.745 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.581 3.802 2.757 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.596 3.810 2.768 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.611 3.819 2.780 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.626 3.827 2.792 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.641 3.836 2.803 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.656 3.844 2.815 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.670 3.853 2.827 200. 50. 10000. 0.0001 10000. 0.0001 1.0 6.685 3.861 2.838 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.800 4.500 3.000 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.803 4.502 3.001 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.806 4.503 3.002 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.810 4.505 3.003 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.813 4.506 3.004 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.816 4.508 3.005 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.819 4.510 3.006 200. 50. 10000. 0.0001 10000. 0.0001 1.0 7.822 4.511 3.007 200. 50. 10000. 0.0001 10000. 0.0001 50.0 7.822 4.511 3.007 200. 50. 10000. 0.0001 10000. 0.0001 0.0 7.822 4.511 3.007 200. 50. 10000. 0.0001 10000. 0.0001 0.01 0.125 0.211766 0.220159 1.0 0. 10.0 10.0 1024 0. 30.0000 20.0000 20.0000 15 1 0.01 0.0 1.00000 0. 0. 1.00000 0. 1.00000 ------------------------------------------------------------------------ Explanation: # Number of layers in model (max 500) 42 # # The model: # thickness v_p v_s den Q_p Q_s absorption band parameters 1.0 4.845 2.795 2.290 200. 50. 10000. 0.0001 10000. 0.0001 1.0 5.136 2.964 2.339 200. 50. 10000. 0.0001 10000. 0.0001 # Note that the bottom layer has thickness 0. # Also note that the formatting does not matter. # Velocities and Q cannot equal 0. # # Slowness integration parameters: # pmin p2 p3 pmax pfac 0.01 0.125 0.211766 0.220159 1.0 # How does one choose these parameters? # pmin is the minimum slowness synthesized # p2 should be a little less than the slowness corresponding to the # highest model velocity # p3 should be a little higher than the slowness corresponding to he # lowest model velocity # pmax is the maximum slowness synthesized # pfac determines the slowness sampling; pfac=1 gives a nominal slowness # sampling, pfac>1 increases sampling, pfac<1 decreases sampling; # slowness sampling is calculated as dp=1./(2*fmax*rmax*pfac) # where fmax=max frequency and rmax=max distance synthesized # # Frequency parameters: # fmin fmax fnyq npts 0.0 10.0 10.0 1024 # Note that the sample interval dt = 1/(2*fnyq) # # Distance (in km) and azimuth (1st line) # distance>0 means compute synthetic at this one distance # distance=0 means look at next line for min distance, distance spacing, # and total number of distances (max 10?) (this option shown below) # distance<0 means the next line says how many distances there are and # the third line is a list of distances # azimuth doesn't matter since there is no assumed fault geometry 0. 30.0000 20.0000 20.0000 15 # # Number of sources (max 5 in same layer) 1 # Source depth (.01 insures that it isn't located on a layer interface) 0.01 # Receiver depth (must be above source, not on an interface) 0.0 # components of moment tensor (Mxx, Mxy, Mxz, Myy, Myz, Mzz) # This example is for an isotropic source 1.00000 0. 0. 1.00000 0. 1.00000 (last line is for mgenkennett only) Here is an example command file (do.kennett) to run mgenkennett and xtsynth: ---------------------------------------------------------------------- /net/rock/shearer/CLASS/227C/REFL/mgenkennett << ! $1 ! mkdir $2 mv mspec $2 cd $2 /net/rock/shearer/CLASS/227C/REFL/xtsynth -d << ! 8 0 ! rm t_* rm r_* cd .. ---------------------------------------------------------------------- This command file requires two arguments. The first is the name of the input file to mgenkennett (in.kennett in this case). The second is the SAC file where the seismograms are to be stored. Here is an example of what to say: do.kennett in.kennett SACDIR The output of mgenkennett (mspec) is moved to the SAC directory. Then the xtsynth program is run. Options for this program include: -d = generate displacement synthetics -v = generate velocity synthetics The program will automatically generate all three components of output (z=vertical, r=radial, t=transverse) in the SAC files. Note that the command file in this case deletes the radial and transverse files. xtsynth also asks for a reducing velocity and a time delay: 8 0 (reduction velocity, time delay) The resulting seismograms can be plotted using the SAC program.