xdemo3

su 的源代码库

#! /bin/sh
# This input file reproduces the output in Juhlin (1995) Figure 8b.
# Note that the source in Table 1 in Juhlin (1995) is given as g(t), but should be f(t)
# be f(t), i.e. the derivative of the Gaussian (dg)
# Since the dg source waveform is defined by its avg frequency content
# favg should be set to 17.077 Hz to correpsond to alpha=5000 in
# Juhlin (1995). Also the time step, dt, is not 0.8 ms, but 0.5 ms
#
# The unfiltered output contains high frequency numerical noise (I do not
# quite understand why) which is nicely removed usinf a highcut filter
#
# Note that the geometry of the model in Figure 8a is very schematic 
# and does not correspond to the model actually used in producing the 
# synthetics (the things we can get away with in our publications).

suea2df dt=0.0005 lt=1.2 zmax=2000 xmin=0 xmax=4000 dx=5 verbose=1 \
sx=-500 sz=20 favg=17.077 ts=0.08 wtype=dg \
snfile=tmp hsz=00 vsx=20 snaptime=0.19,0.59,0.99 efile=mod \
bc=2,10,10,10 qsw=0 asw=1 msw=1 \
xl=0,3500 \
0 \
 zl=0,5,500,505 \
 vpl=3320,2898,2898,3320 \
 vsl=1398,1398,1398,1398 \
 rhol=2250,2250,2250,2250 \
 epl=0,0.2566,0.2566,0\
 dsl=0,0.4025,0.4025,0\
 anl=0,-30,-30,0 \
3500 \
 zl=2000,2005,2500,2505 \
 vpl=3320,2898,2898,3320 \
 vsl=1398,1398,1398,1398 \
 rhol=2250,2250,2250,2250 \
 epl=0,0.2566,0.2566,0\
 dsl=0,0.4025,0.4025,0\
 anl=0,30,30,0 \
>out2

suwind < ss.su key=duse min=1 max=1 | suwind j=10 | suwind key=offset max=4000 | \
sufilter f=50,75 amps=1.,0. | suresamp dt=0.004 nt=301 | suxwigb interp=1 f2=0 d2=50 perc=50

suwind < ss.su key=duse min=1 max=1 | suwind j=10 | suwind key=offset max=4000 | \
sugain scale=1 | sufilter f=50,75 amps=1.,0. | suresamp dt=0.004 nt=301 | \
supswigb f2=0 d2=50 nbpi=600 d1num=0.2 d2num=1000 perc=90 \
 x2beg=0 x2end=4000 wbox=1.75 hbox=2.25 labelsize=6 > F8b.ps