Commit b6cd3dc3 authored by thomas.forbriger's avatar thomas.forbriger

refmet [WP]: regional exmple is operational

parent 92b5c2ae
......@@ -38,6 +38,7 @@ clean: ;
$(shell tail $< -n +7 | head -n 1) \
$(shell tail $< -n +10 | head -n 1) \
> $@
echo $@: $< >> $@
include $(patsubst %.main.cfg,%.d,$(wildcard *.main.cfg))
......@@ -58,15 +59,24 @@ iasp91.syg.xxx: refmet.iasp91flat.xxx refmet.source.cfg
# -------------------
#
# shallow seismic examples creating seismograms of a receiver profile
$(EXAMPLE): refmet.main.cfg
refmetx -v -4 -o example -ty $(FILETYPE) -c refmet.main.cfg
MAINFILE?=refmet.main.cfg
$(EXAMPLE): $(MAINFILE)
refmetx -v -4 -o example -ty $(FILETYPE) -c $<
.PHONY: example
example: $(EXAMPLE)
# ======================================================================
# display seismograms
$(addsuffix .ps,$(EXAMPLE)): %.ps: %
refractx -ty $(FILETYPE) -R 6371. \
-S3 20000.,50000. -Sa 8000. -Sc 2.e4 -Se 0.5 \
-d $@/ps $<
# ======================================================================
# display results
%.psp: %.ps; gv $<; /bin/rm -fv $<
%.pdf: %.ps; ps2pdf $<
%.pdp: %.pdf; $(PREVIEWER) $<
......
# Isotropic Earth model IASP91 with quality factors from PREM at 1 Hz
number of layers: 14
earth radius (km):-1* -6371.000 reference frequency (Hz): 0.000
Zb alpha beta rho Qalpha Qbeta Rb
-----------------------------------------------------------------------------
0.000 0.33180 0.00000 0.001300 1000.000 -1.000 6371.000
10.000 5.80000 3.36000 2.600000 1323.917 600.000 6361.000
20.000 5.80000 3.36000 2.600000 1323.917 600.000 6351.000
27.200 6.50000 3.75000 2.900000 1334.643 600.000 6343.800
34.400 6.50000 3.75000 2.900000 1334.643 600.000 6336.600
77.200 8.04244 4.47733 3.377336 1430.862 600.000 6293.800
120.000 8.04748 4.49244 3.372684 1423.228 600.000 6251.000
148.500 8.08959 4.50285 3.368810 193.275 80.000 6222.500
193.500 8.19167 4.51020 3.364816 197.524 80.000 6177.500
210.000 8.27709 4.51635 3.361474 201.103 80.000 6161.000
232.500 8.34106 4.54158 3.436521 360.400 143.000 6138.500
278.000 8.46516 4.60074 3.456825 361.711 143.000 6093.000
324.500 8.63306 4.68078 3.484294 363.435 143.000 6046.500
372.000 8.80461 4.76256 3.512361 365.141 143.000 5999.000
halfspace: 8.93510 4.82477 3.533709 366.402 143.000 5975.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
amplitude correction factors depending on source depth are
source depth transform. exp. explosion corr. double-couple corr.
100.km -1 1.01582 1.03189
100.km -1 1.01582 1.03189
600.km -1 1.09875 1.20726
600.km -1 1.09875 1.20726
This file was generated by:
CHOPMOD V1.10 spherical --> flat earth
Minimum allowed relative parameter stepsize: 0.0200
Minimum processed radius: 6000.000km
Scanning stepsize: 1.000km
A reference frequency of 0Hz means values are frequency independent.
No earth flattening approximation was applied!
Example control file for refmet
Filename for subsurface model:
refmet.iasp91_flat.cfg
Filename for source description:
refmet.source.cfg
Filename for receiver definition:
refmet.rcv.cfg
slowness sampling (s/km) pmin, pwil, pwir, pmax, Np:
0.,0.,0.5,0.7,500
temporal sampling (Hz und s) fmin, fwil, fwir, fmax, dt, length of trace:
0., 0., 10., 15., 0.025, 90.
horizontal double couple (as defined by Ungerer)
typ, outsig, srcsig, The, Thd, ZQ:
1, 2, 1, 0., 3., 50.
1, 2, 1, 0., 3., 34.6
M0, Mxx, Myy, Mzz, Mxy, Mxz, Myz
100., 0., 0., 0., -1., 0., 0.
......
Markdown is supported
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