#!/bin/sh ##### # clean up ##### rm -f *PLT rm -f *egn rm -f *dsp rm -f sdisp96* rm -f *sac rm -f tmp* ##### # create the velocity model interactively ##### mkmod96 << EOF simple.mod Simple Crustal Model 0 40 6 3.5 2.7 0 0 0 0 1 1 0 8 4.7 3.3 0 0 0 0 1 1 EOF ##### # create the dfile use sprep96 -h for the format ##### cat > dfile << EOF 2500 0.5 2048 0 8 2550 0.5 2048 0 8 2600 0.5 2048 0 8 2650 0.5 2048 0 8 2700 0.5 2048 0 8 2750 0.5 2048 0 8 2800 0.5 2048 0 8 2850 0.5 2048 0 8 EOF ##### # setup multimode surface wave dispersion and synthetics ##### HS=10.0 HR=0.0 sprep96 -M simple.mod -d dfile -L -R -NMOD 100 -HS ${HS} -HR ${HR} ##### # get the dispersion and the eigenfunctions ##### sdisp96 sregn96 slegn96 ##### # make the synthetic for the station for a given focal mechanism and then convert to SAC ##### DIP=45 RAKE=45 STK=45 AZ=45 BAZ=225 MW=4.5 spulse96 -p -V -l 4 -d dfile | \ fmech96 -D ${DIP} -R ${RAKE} -S ${STK} -A ${AZ} -ROT -MW ${MW} | \ f96tosac -B ##### # do interactive multiple filter analysis to get group velocity dispersion ##### do_pom B*Z00.sac do_pom B*T00.sac ##### # The disperions curves are in the files GRN21T.dsp and GRN21Z.dsp # combine these to create a SURF96 dispersion file ##### cat rayl.dsp love.dsp > disp.d ##### # get the theoretical dispersion for this model ##### sdpegn96 -R -S -C -D disp.d -XLOG -PER -YMIN 2.5 -YMAX 5.0 sdpegn96 -L -S -C -D disp.d -XLOG -PER -YMIN 2.5 -YMAX 5.0 ##### # Look at the observed and predicted disperion ##### plotxvig < SREGNC.PLT plotxvig < SLEGNC.PLT