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

added definition of Fourier transformation to online help

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: 3670
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 279778bc
...@@ -103,6 +103,8 @@ c spectra is sample index ...@@ -103,6 +103,8 @@ c spectra is sample index
c 13/11/10 V3.20 - do not expose misleading term "green" to the c 13/11/10 V3.20 - do not expose misleading term "green" to the
c user c user
c - use GSL instead of numerical recipes c - use GSL instead of numerical recipes
c 12/01/2011 V3.20b - added definition of Fourier transformation
c online help text
c c
c============================================================================== c==============================================================================
c c
...@@ -112,7 +114,7 @@ c first we declare some general variables ...@@ -112,7 +114,7 @@ c first we declare some general variables
c c
character*79 version,CVSID character*79 version,CVSID
parameter(version= parameter(version=
&'GREDA V3.20 Calculate Fourier-Bessel expansion coefficients') &'GREDA V3.20b Calculate Fourier-Bessel expansion coefficients')
parameter(CVSID= parameter(CVSID=
& '$Id$') & '$Id$')
c c
...@@ -335,7 +337,9 @@ c ...@@ -335,7 +337,9 @@ c
print *,' ' print *,' '
print *,'Parameters defining the expansion:' print *,'Parameters defining the expansion:'
print *,'-1 Use the Hankel function H^(1)_0 instead of J_0' print *,'-1 Use the Hankel function H^(1)_0 instead of J_0'
print *,' (see remark on Fourier transformation below)'
print *,'-2 Use the Hankel function H^(2)_0 instead of J_0' print *,'-2 Use the Hankel function H^(2)_0 instead of J_0'
print *,' (see remark on Fourier transformation below)'
print *,'-q f,e Damping factor for inner product.' print *,'-q f,e Damping factor for inner product.'
print *,' The inner product damping will be calculated' print *,' The inner product damping will be calculated'
print *,' as rho=delta_r_min*f*omega**e.' print *,' as rho=delta_r_min*f*omega**e.'
...@@ -393,6 +397,22 @@ c ...@@ -393,6 +397,22 @@ c
print *,'rho should be chosen somewhere around delta_r_min*omega,' print *,'rho should be chosen somewhere around delta_r_min*omega,'
print *,'where delta_r_min is the minimal offset difference.' print *,'where delta_r_min is the minimal offset difference.'
print *,' ' print *,' '
print *,'Definition of the Fourier transformation:'
print *,'The Fourier transformation used in this program and in'
print *,'related programs (like gremlin, syg, and gresy) is'
print *,'defined as'
print *,' '
print *,' U(omega) = int_-infnity^+infnity u(t) exp(-i*omega*t) dt'
print *,' '
print *,'Theoretical descriptions of wave propagation often use'
print *,'exp(i*omega*t) as transform kernel instead of exp(-i*omega*t)'
print *,'in order to make positive wavenumbers equivalent to wave'
print *,'propagation in positive coordinate direction. The Fourier'
print *,'coefficients calculated by this program consequently are'
print *,'the complex conjugates of those used in theory. Where'
print *,'H^(2)_0 is used in theory, you have to use H^(1)_0'
print *,'in greda.'
print *,' '
print *,'$Id$' print *,'$Id$'
call pwo_cvsid call pwo_cvsid
print *,' ' print *,' '
......
...@@ -47,12 +47,14 @@ c the problem disappeared for unkown reason... ...@@ -47,12 +47,14 @@ c the problem disappeared for unkown reason...
c 26/02/10 V1.12 introduced line-source mode c 26/02/10 V1.12 introduced line-source mode
c 01/03/10 V1.13 correction (linesource summation was simply wrong) c 01/03/10 V1.13 correction (linesource summation was simply wrong)
c 02/03/10 V1.14 corrected (hopefully) line source integration c 02/03/10 V1.14 corrected (hopefully) line source integration
c 12/01/2011 V1.14b - added definition of Fourier transformation
c online help text
c c
c============================================================================== c==============================================================================
c c
program gresy program gresy
character*79 version character*79 version
parameter(version='GRESY V1.14 GREens function SYnthetics') parameter(version='GRESY V1.14b GREens function SYnthetics')
c dimensions c dimensions
integer maxtr, maxsamp, maxom, maxu integer maxtr, maxsamp, maxom, maxu
parameter(maxu=10000, maxtr=maxu, maxom=4100, maxsamp=maxom*2) parameter(maxu=10000, maxtr=maxu, maxom=4100, maxsamp=maxom*2)
...@@ -151,10 +153,28 @@ c ...@@ -151,10 +153,28 @@ c
print *,' produced by ''gremlin'' ' print *,' produced by ''gremlin'' '
print *,'-S suppress zero frequency and zero slowness' print *,'-S suppress zero frequency and zero slowness'
print *,'-1 use Hankel 1 instead of Bessel' print *,'-1 use Hankel 1 instead of Bessel'
print *,' (see remark on Fourier transformation below)'
print *,'-2 use Hankel 2 instead of Bessel' print *,'-2 use Hankel 2 instead of Bessel'
print *,' (see remark on Fourier transformation below)'
print *,'-L simulate seismograms from line source (2D)' print *,'-L simulate seismograms from line source (2D)'
print *,'-R calculate radial component' print *,'-R calculate radial component'
print *,' ' print *,' '
print *,'Definition of the Fourier transformation:'
print *,'The Fourier transformation used in this program and in'
print *,'related programs (like gremlin, syg, and greda) is'
print *,'defined as'
print *,' '
print *,' U(omega) = int_-infnity^+infnity u(t) exp(-i*omega*t) dt'
print *,' '
print *,'Theoretical descriptions of wave propagation often use'
print *,'exp(i*omega*t) as transform kernel instead of exp(-i*omega*t)'
print *,'in order to make positive wavenumbers equivalent to wave'
print *,'propagation in positive coordinate direction. The Fourier'
print *,'coefficients calculated by this program consequently are'
print *,'the complex conjugates of those used in theory. Where'
print *,'H^(2)_0 is used in theory, you have to use H^(1)_0'
print *,'in gresy.'
print *,' '
print *,'Input file units:' print *,'Input file units:'
print *,' ' print *,' '
print *,' Slowness unit: s/m' print *,' Slowness unit: s/m'
......
...@@ -30,6 +30,8 @@ c 14/03/02 V1.3 do not calculate static case with flgevas ...@@ -30,6 +30,8 @@ c 14/03/02 V1.3 do not calculate static case with flgevas
c 15/03/02 V1.4 iragc was declared twice c 15/03/02 V1.4 iragc was declared twice
c 26/03/02 V1.5 now optionally support complex angular frequencies c 26/03/02 V1.5 now optionally support complex angular frequencies
c 04/03/10 V1.6 comment regarding the source amplitude was wrong! c 04/03/10 V1.6 comment regarding the source amplitude was wrong!
c 12/01/2011 V1.6b - added definition of Fourier transformation
c online help text
c c
c============================================================================== c==============================================================================
c c
...@@ -38,7 +40,7 @@ c ...@@ -38,7 +40,7 @@ c
c variables c variables
character*79 version character*79 version
parameter (version= parameter (version=
& 'SYG V1.6 calculate synthetic greens matrix') & 'SYG V1.6b calculate synthetic greens matrix')
c c
integer mu, mom integer mu, mom
parameter (mu=10000, mom=4100) parameter (mu=10000, mom=4100)
...@@ -107,6 +109,20 @@ c go ...@@ -107,6 +109,20 @@ c go
print *,' -s sigma set imaginary part of angular frequency' print *,' -s sigma set imaginary part of angular frequency'
print *,' to -sigma.' print *,' to -sigma.'
print *,' ' print *,' '
print *,'Definition of the Fourier transformation:'
print *,'The Fourier transformation used in this program and in'
print *,'related programs (like gremlin, greda, and gresy) is'
print *,'defined as'
print *,' '
print *,' U(omega) = int_-infnity^+infnity u(t) exp(-i*omega*t) dt'
print *,' '
print *,'Theoretical descriptions of wave propagation often use'
print *,'exp(i*omega*t) as transform kernel instead of exp(-i*omega*t)'
print *,'in order to make positive wavenumbers equivalent to wave'
print *,'propagation in positive coordinate direction. The Fourier'
print *,'coefficients calculated by this program consequently are'
print *,'the complex conjugates of those used in theory.'
print *,' '
c 4.3.2010 notice: the following comment was wrong! c 4.3.2010 notice: the following comment was wrong!
c gr_refsub reads the source amplitude and uses it for scaling c gr_refsub reads the source amplitude and uses it for scaling
c print *,'Notice: the source amplitude is ignored.' c print *,'Notice: the source amplitude is ignored.'
......
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