Commit 767ecc08 authored by thomas.forbriger's avatar thomas.forbriger Committed by thomas.forbriger
Browse files

correction

This is a legacy commit from before 2015-03-01.
It may be incomplete as well as inconsistent.
See COPYING.legacy and README.history for details.


SVN Path:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/trunk
SVN Revision: 2414
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 8783afce
c this is <nyquist_check.f>
c ----------------------------------------------------------------------------
c ($Id: nyquist_check.f,v 1.1 2007-07-05 15:11:17 tforb Exp $)
c ($Id: nyquist_check.f,v 1.2 2007-10-22 13:55:35 tforb Exp $)
c
c Copyright (c) 2007 by Thomas Forbriger (BFO Schiltach)
c
......@@ -8,6 +8,7 @@ c Check Nyquist criterion for trapezoid rule
c
c REVISIONS and CHANGES
c 05/07/2007 V1.0 Thomas Forbriger
c 22/10/2007 V1.1 use correct units
c
c ============================================================================
c
......@@ -15,9 +16,9 @@ c
c
real r, p, o
c
c r: receiver offset in km
c p: slowness interval in s/km
c o: maximum frequency in 1/s
c r: receiver offset in m
c p: slowness interval in s/m
c o: maximum angular frequency in 1/s
c
c This function will check whether the slowness sampling is adequate for
c the oscillating Bessel functions in the slowness expansion integral when
......@@ -30,10 +31,11 @@ c
parameter(pi2=2.*3.14159265358979)
d=o*r*p/pi2
if (d.gt.0.25) then
print *,'WARNING: The trapezoid stepsize is ',d,'*2*pi.'
print *,' Receiver offset: ',r,' km'
print *,'WARNING: The trapezoid stepsize is ',d,'*2*pi = ',
& 2*d,'*pi.'
print *,' Receiver offset: ',1.e-3*r,' km'
print *,' Maximum frequency: ',o/pi2,' Hz'
print *,' Slowness stepsize: ',p,' s/km'
print *,' Slowness stepsize: ',p*1.e3,' s/km'
if (d.gt.0.9) then
print *,'!! The stepsize for the argument of the Bessel'
print *,'!! function is larger than 1.8*pi. The'
......
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