xdemo1

su 的源代码库

#! /bin/sh
# This input file reproduces the output in Juhlin (1995) Figure 6a.
# 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 24.15 Hz to correpsond to alpha=10000 in
# Juhlin (1995)  
#
# original
#suea2df dt=0.0005 lt=0.4 nz=294 fx=0 nx=255  dx=6 verbose=1 \
#snfile=tmp hsz=00 vsx=200 snaptime=0.02,0.05,.1,.15,.2,.25.3,.35,.4 efile=mod \
#bc=2,1,1,0 qsw=0 asw=0 \
#sx=0 sz=12 favg=24.15 ts=0.05 wtype=dg \
#xl=1 \
# \
# zl=150,156 \
# vpl=2898,3320 \
# vsl=1398,1398 \
# rhol=2250,2250\
#>out2

modelfile=junk.unif             # input model file for unif2aniso
ninf=1                          # number of interfaces (surface counts)
x0=0,0,0                        # x-position(s) for  vp00,vs00,rho00, etc.
z0=0,150,156			# z-position(s) for  vp00,vs00,rho00, etc.
nz=294                          # size of z (depth) dimension  of model
nx=255                          # size of x (horizontal) dimension of model
dz=6                            # increment in z direction
dx=6                            # increment in x direction
vp00=2898,3320                  # P-wavespeed(s) at (z0,x0)
vs00=1398,1498                  # S-wavespeed(s) at (z0,x0)
rho00=2250,2250                 # density(s) at (z0,x0)
dt=0.0005                       # time sampling interval
lt=0.4	                        # latest time modeled
fx=0                            # first x value
verbose=1                       # =1 chatty, =0 silent
snfile="snaps.su"                # output file for snapshots
hsz=0.0		                # z-position of horizontal line of geophones
vsx=0.0		                # x-position of vertical line of geophones
snaptime=0.02,0.05,.1,.15,.2,.25.3,.35,.4	# times of snapshots
bc=2,1,1,0                      # boundary conditions
qsw=0                           # =1 put in attenuation
asw=0                           # =1 anisotropy
sx=0                            # x-position of sources
sz=0                            # z-position of sources
favg=24.15                      # average frequency
ts=.05                          # source duration
wtype=dg                        # waveform type


# build stiffness and density files
unif2aniso < $modelfile ninf=$ninf x0=$x0 z0=$z0 nz=$nz nx=$nx \
dx=$dx dz=$dz vp00=$vp00 vs00=$vs00 rho00=$rho00

# transpose stiffness and density
for i in c11 c13 c15 c33 c35 c55 rho
do
        mv ${i}_file tmp.file
        transp n1=294 < tmp.file > ${i}_file
done

rm tmp.file

# run suea2df
suea2df dt=$dt lt=$lt nz=$nz fx=$fx nx=$nx dx=$dx dz=$dz verbose=1 \
snfile=$snfile hsz=$hsz vsx=$vsx snaptime=$snaptime \
bc=$bc qsw=$qsw asw=$asw sx=$sx sz=$sz favg=$favg ts=$ts wtype=$wtype \
>out2


exit 0


suwind < ss.su key=duse min=1 max=1 | suwind j=3 | suwind key=offset max=420 | \
suxwigb f2=0 d2=18 perc=98
suwind < ss.su key=duse min=1 max=1 | suwind j=3 | suwind key=offset max=420 | \
sugain scale=-1 | supswigb f2=0 d2=18 nbpi=600 d1num=0.1 d2num=72 grid1=dash perc=98 \
 x2beg=414 x2end=0 wbox=1.75 hbox=2 labelsize=6 > F6a.ps
 
#awk < out2 '{if(NF==2) print $2,$1}' | a2b | xgraph n=101 pairs=1 
#awk < out2 '{if(NF==2) print $2}' | a2b n1=1 | suaddhead ns=101 | \
#sushw key=dt a=500 | suspecfx | suxwigb fx.su x1end=100