# Calculating the amplitudes of Gaussian packets # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Transform original Marmousi seismograms from binary to ASCII files N1=750 N2=96 N3=1 GRD='timess.out' NDIG=0 binasc: CAL='mul.cal' GRD1='timess.out' GRD2='mgpa-sm.dat' GRD3='timess.out' grdcal: echo.pl: "7500*0" ">time.out" copy.pl: "time.out" "times.out" append.pl: "times.out" "timess.out" append.pl: "times.out" "time.out" # Create files with constant initial parameters # R0, Y0 of Gaussian packets N1=12 N2=34 N3=36 echo.pl: '\$1=AR0' ">tmp.cal" CAL="tmp.cal" AR0=-0.25E-6 GRD1='fgbr22.out' grdcal: echo.pl: '\$1=AY0' ">tmp.cal" CAL="tmp.cal" AY0=0.25E-6 GRD1='fgby22.out' grdcal: # Calculating minimum and maximum values of R0, Y0 N1NEW=1 N2NEW=1 N3NEW=1 GNORM=-999 GRD='fgbr22.out' GRDNEW='fgbrmin.out' grdnorm: GRD='fgby22.out' GRDNEW='fgbymin.out' grdnorm: GNORM=999 GRD='fgby22.out' GRDNEW='fgbymax.out' grdnorm: GRD='fgbr22.out' GRDNEW='fgbrmax.out' grdnorm: # Calculating the amplitudes of Gaussian packets NX=116 NT=750 DX=25.0 DT=0.004 VSURF=1500.0 GWM=3.0 FMIN=5.0 FMAX=55.0 FREF=35.0 PARK=0.5 PARN=0.5 KAPPA=1.253314137 NXTP=14688 TIMESEC='times.out' SYNTSEC='synts.out' gpanal: