psplot

su 的源代码库

#! /bin/sh
# generate synthetic common midpoint gathers

# convenient scratch directory location
S=.
CMP=150

# set parameters
outcmg=$S/cmgsynhti
outeps=$S/outsynti.eps
nxo=30 dxo=0.1 fxo=0.1
nxm=360 dxm=0.025 fxm=1.5
nt=401 dt=0.01 ft=0
v00=2.5 dvdz=0.0

#susynlvfti >$outcmg trans=1 nt=$nt dt=$dt ft=$ft ob=0 \
#	nxo=$nxo dxo=$dxo fxo=$fxo \
#	nxm=$nxm dxm=$dxm fxm=$fxm \
#	v00=$v00 dvdx=0 dvdz=$dvdz fpeak=15 \
#	ref="3.449,0.259;4.663,0.960;7.337,0.960" \
#	ref="3.366,0.366;4.357,1.357;7.643,1.357" \
#	ref="3.259,0.449;3.960,1.663;8.040,1.663" \
#	ref="3.135,0.500;3.550,1.855;8.500,1.855" \
#	ref="3.000,0.518;3.000,1.920;9.000,1.920" \
#	a=1.3333334 f=.4 l=.3

# Plot a CMP gather
suwind <$outcmg key=cdp min=$CMP max=$CMP count=$nxo |
supswigp -geometry 300x800+100+100 \
	n1=$nt d1=$dt f1=$ft label1="Time (s)" \
	n2=$nxo d2=$dxo f2=$fxo label2="Offset (km)" \
	x1beg=1.0 x1end=2.5 labelsize=24 \
	d1num=0.2 f1num=0.6 d2num=1 f2num=0 \
	title="synthetic CMP $CMP" \
	perc=99.7 grid1=solid>$outeps

echo "You may print or view the file $outeps with a PostScript previewer."

exit 0