Commit eee02da4 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

[WP] (i20_gremlin_example): provide synthetic seismogram data

- calculate seismograms with reflectivity method
- add gaussian noise
parent 32ab1793
......@@ -10,6 +10,8 @@
#
# ============================================================================
#
# seismogram file type to be used
FILETYPE=su
all:
......@@ -23,8 +25,28 @@ edit: flist; vim $<
clean: ;
-find . -name \*.bak | xargs --no-run-if-empty /bin/rm -v
-/bin/rm -vf flist
-/bin/rm *.inctex *.bak *.ps *.syg *.pdf
-/bin/rm -vf *.inctex *.bak *.ps *.syg *.pdf *.d *.xxx
.PHONY: cleanall
cleanall: clean
-/bin/rm -fv *.$(FILETYPE)
#======================================================================
# dependency files
# ----------------
%.d: %.main.cfg
echo $< $@: \
$(shell tail $< -n +4 | head -n 1) \
$(shell tail $< -n +7 | head -n 1) \
$(shell tail $< -n +10 | head -n 1) \
> $@
include $(patsubst %.main.cfg,%.d,$(wildcard *.main.cfg))
# ======================================================================
#
# calculate synthetic Fourier-Bessel expansion coefficients
# ---------------------------------------------------------
%.syg: %.mod.cfg zforce.src.cfg
syg $^ $@ -P 200,200,100.,3.
......@@ -42,4 +64,35 @@ clean: ;
layer.pdf: layersyg.ps layer.fig
fig2dev -L pdf layer.fig $@
# ======================================================================
#
# calculate synthetic seismograms
# -------------------------------
%.$(FILETYPE): %.main.cfg
refmetx -ty $(FILETYPE) -v -4 -c -s TR \
-o $(patsubst %.$(FILETYPE),%,$@) $<
/bin/mv -fv $(patsubst %.$(FILETYPE),%,$@).TZ $@
.PRECIOUS: refmetLQ.su
.PRECIOUS: refmetHQ.su
# calculate table of rms amplitude
%.rms.xxx: %.$(FILETYPE)
sigval -type $(FILETYPE) -format "%OFF %RMS" $< \
| cut -c 3- | sort -n > $@
%.noiserms.xxx: %.rms.xxx
tail -n 1 $< | cut -f 2 -d ' ' > $@
# simulated data
# --------------
# add gaussian white noise at appropriate rms level
# discard samples at times larger than 1s
simulation%.$(FILETYPE): refmet%.$(FILETYPE) refmet%.noiserms.xxx
printf "fac 2.\nnoi %s\n" \
$$(cat $(word 2,$^)) \
| tidofi -cs -v -o -type $(FILETYPE) -Type $(FILETYPE) $@.xxx $<
tijerasxx --verbose --overwrite \
--iformat $(FILETYPE) --oformat $(FILETYPE) \
--secondsd 1. $@ $@.xxx
# ----- END OF Makefile -----
Schicht ueber Halbraum Referenz Modell
Layer on top of half-space
number of layers: 2
earth radius (km): 0.000 reference frequency (Hz): .000
......
lineare Auslage
linear profile of receivers
mode, Vred, Tl, Tr, Nr
2, 0., 0., 0., 20
2, 0., 0., 0., 72
rmin, rmax, phi
0.005, 0.100, 0.
0.001, 0.072, 0.
Synthetics main
Control file for refmet (calculate synthetic seismograms)
Model taken from
<MODEL>
layer.mod.cfg
Source taken from
zforce.src
zforce.src.cfg
Receivers taken from
lin.rcv
lin.rcv.cfg
umin, uwil, uwir, umax, Nu
0., 0., 6.5, 8.5, 1000
0., 0., 4.5, 6.5, 8000
fmin, fwil, fwir, fmax, dt, T
0., 0., 130., 130., 0.0010, 3.0
0., 0., 80., 130., 0.0010, 3.0
Just a single vertical force at the surface
typ, sig, src, Ts, Td, Zs
2, 2, 2, 0., 0., 0.000
2, 2, 1, 0., 0.02, 0.000
F0
1.
units
[m/s]
m/s
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment