Commit 0527c059 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

[WP] (i20_gremlin_example): make consistent use of basenames

A consitent way to use file name patterns is essential to support the addition
of further future examples.
parent 1fcc39dd
......@@ -4,6 +4,10 @@
# Copyright (c) 2016 by Thomas Forbriger (BFO Schiltach)
#
# rules and command to process toy example for gremlin
#
# basenames of examples:
# - layerHQ
# - layerLQ
#
# REVISIONS and CHANGES
# 18/06/2016 V1.0 Thomas Forbriger
......@@ -79,8 +83,8 @@ layer.pdf: layerHQsyg.ps layer.fig
-o $(patsubst %.$(FILETYPE),%,$@) $<
/bin/mv -fv $(patsubst %.$(FILETYPE),%,$@).TZ $@
.PRECIOUS: refmetLQ.su
.PRECIOUS: refmetHQ.su
.PRECIOUS: layerLQ.su
.PRECIOUS: layerHQ.su
# calculate table of rms amplitude
%.rms.xxx: %.$(FILETYPE)
......@@ -93,7 +97,7 @@ layer.pdf: layerHQsyg.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
%.sim.$(FILETYPE): %.$(FILETYPE) %.noiserms.xxx
printf "fac -2.\nnoi %s\n" \
$$(cat $(word 2,$^)) \
| tidofi -cs -v -o -type $(FILETYPE) -Type $(FILETYPE) $@.xxx $<
......@@ -103,9 +107,9 @@ simulation%.$(FILETYPE): refmet%.$(FILETYPE) refmet%.noiserms.xxx
# plot simulated data
# -------------------
expHQ.xxx: ; echo 0.4 > $@
expLQ.xxx: ; echo 0.8 > $@
simulation%.ps: simulation%.$(FILETYPE) exp%.xxx
layerHQ.exp.xxx: ; echo 0.4 > $@
layerLQ.exp.xxx: ; echo 0.8 > $@
%.sim.ps: %.sim.$(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" $<
......@@ -122,26 +126,26 @@ layer%.tt.xxx: layer.tt.cfg; /bin/cp -vpd $< $@
layer%.VP.mod.xxx: layer.VP.mod.cfg; /bin/cp -vpd $< $@
# display arrival time picks
arrival%.ps: simulation%.$(FILETYPE) exp%.xxx layer%.tt.xxx
%.arrival.ps: %.sim.$(FILETYPE) %.exp.xxx %.tt.xxx
refractx -ty $(FILETYPE) -d $@/cps -St -0.01,0.15 -Ev \
-Lt 3 -Cc -EP T -ET T -Et T -Fa $(word 3,$^) \
-Sm 2 -Sr 1.60 -Eg -Se $$(cat $(word 2,$^)) -Sa 30. -Sc 2. \
-Tt "simulated data and first arrivals" $<
#
# display travel time polygonial and P-wave velocity model
vpmod%.ps: simulation%.$(FILETYPE) exp%.xxx layer%.ttc.xxx layer%.VP.mod.xxx
%.vpmod.ps: %.sim.$(FILETYPE) %.exp.xxx %.ttc.xxx %.VP.mod.xxx
refractx -ty $(FILETYPE) -d $@/cps -St -0.01,0.15 -Ev -TM 0.06,20.\
-Lt 3 -Cc -EP T -ET T -Et T -Fp $(word 3,$^) -Fm $(word 4,$^)\
-Sm 2 -Sr 1.60 -Eg -Se $$(cat $(word 2,$^)) -Sa 30. -Sc 2. \
-Tt "simulated data and travel time curve and model for refracted P-waves" $<
#
# mute seismograms by application of taper
muted%.$(FILETYPE): simulation%.$(FILETYPE) layer%.tap.xxx
%.mute.$(FILETYPE): %.sim.$(FILETYPE) %.tap.xxx
tidofi -v -o -type $(FILETYPE) -Type $(FILETYPE) \
$@ $< poft:$(word 2,$^)
# display muted seismograms
muted%.ps: muted%.$(FILETYPE) exp%.xxx layer%.tap.xxx
%.mute.ps: %.mute.$(FILETYPE) %.exp.xxx %.tap.xxx
refractx -ty $(FILETYPE) -d $@/cps -St -0.15,0.5 \
-Lt 3 -Cc -EP T -ET T -Et T -Ft $(word 3,$^) \
-Sm 2 -Sr 0.37 -Eg -Se $$(cat $(word 2,$^)) -Sa 15. -Sc 5. \
......@@ -151,12 +155,12 @@ muted%.ps: muted%.$(FILETYPE) exp%.xxx layer%.tap.xxx
# calculate Fourier-Bessel expansion coefficients
# -----------------------------------------------
%.grn: muted%.$(FILETYPE)
%.grn: %.$(FILETYPE)
gredax $< $@ -ty $(FILETYPE) -n 100 -s 3.5 -f 80. -t 0. -T 10. \
-1 -v -o
# display Fourier-Bessel expansion coefficients
%grn.ps: %.grn
%.grn.ps: %.grn
export PGPLOT_PS_HEIGHT=7070 ; \
export PGPLOT_PS_WIDTH=10000 ; \
unset PGPLOT_IDENT ; \
......@@ -183,4 +187,10 @@ layerLQ.VP.ini.pmod.xxx layerHQ.VP.ini.pmod.xxx: layer.VP.ini.pmod.cfg
%.pmod.ps: %.pmod.xxx
mop -p -d $@/cps -T -l 4 -clw 3 $<
# ======================================================================
# start gremlin
# -------------
.PHONY:
# ----- END OF Makefile -----
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