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

[WP] (i20_gremlin_example): provide plots of simulated data

parent eee02da4
......@@ -88,11 +88,20 @@ layer.pdf: layersyg.ps layer.fig
# 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" \
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
# plot simulated data
# -------------------
expHQ.xxx: ; echo 0.4 > $@
expLQ.xxx: ; echo 0.8 > $@
simulation%.ps: simulation%.$(FILETYPE) exp%.xxx
refractx -ty $(FILETYPE) -d $@/ps -St -0.15,0.3 -Ev \
-Sm 2 -Sr 0.37 -Eg -Se $$(cat $(word 2,$^)) -Sa 15. -Sc 5. \
-Tt "simulated data with gaussian, uncorrelated noise added" $<
# ----- END OF Makefile -----
Layer on top of half-space (low Q version)
number of layers: 2
earth radius (km): 0.000 reference frequency (Hz): .000
Zb alpha beta rho Qalpha Qbeta Rb
-----------------------------------------------------------------------------
.000 .33180 .00000 .001300 1000.000 -1.000
0.0100000 0.6800 0.40000 1.700000 15.000 15.000
halfspace: 1.7000 1.00000 2.000000 15.000 15.000
par units description
----------------------------------------------------------------------------
Zb km bottom-of-layer depth
alpha km/s P-wave layer velocity
beta km/s S-wave layer velocity
rho g/cm^3 layer density
Qalpha quality factor related to alpha
Qbeta quality factor related to beta
Rb km radius corresponding to bottom of layer
Control file for refmet (calculate synthetic seismograms, low Q version)
Model taken from
layerLQ.mod.cfg
Source taken from
zforce.src.cfg
Receivers taken from
lin.rcv.cfg
umin, uwil, uwir, umax, Nu
0., 0., 4.5, 6.5, 2000
fmin, fwil, fwir, fmax, dt, T
0., 0., 80., 130., 0.0010, 3.0
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