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

resolved an error in the line source summation

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: 2996
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent ca2e2bd4
c this is <gresy.f>
c------------------------------------------------------------------------------
c $Id: gresy.f,v 1.19 2010-02-26 13:02:03 tforb Exp $
c $Id: gresy.f,v 1.20 2010-03-01 15:23:17 tforb Exp $
c
c 26/06/97 by Thomas Forbriger (IfG Stuttgart)
c
......@@ -29,12 +29,13 @@ c in m (rcv-file specifies 0.003, gresy receives 30km)
c I introduced verbose output and checked the program
c the problem disappeared for unkown reason...
c 26/02/10 V1.12 introduced line-source mode
c 01/03/10 V1.13 correction (linesource summation was simply wrong)
c
c==============================================================================
c
program gresy
character*79 version
parameter(version='GRESY V1.12 GREens function SYnthetics')
parameter(version='GRESY V1.13 GREens function SYnthetics')
c dimensions
integer maxtr, maxsamp, maxom, maxu
parameter(maxu=4500, maxtr=maxu, maxom=4100, maxsamp=maxom*2)
......@@ -144,7 +145,7 @@ c
print *,' amplitude unit: m**3/s if spectrum represents displacement'
print *,' waveform in m'
print *,' '
print *,'$Id: gresy.f,v 1.19 2010-02-26 13:02:03 tforb Exp $'
print *,'$Id: gresy.f,v 1.20 2010-03-01 15:23:17 tforb Exp $'
print *,' '
call refmet_rcvinf
stop
......@@ -407,10 +408,11 @@ c set slowness taper
if (radial) then
stop 'ERROR: radial component not defined for 2D'
else
sdata(io)=0.5d0*sdata(io)+green(io,iu)*
sdata(io)=sdata(io)+green(io,iu)*
& cos(arg)*du*tf_costap(iu,0,0,ltap,rtap)
endif
enddo
sdata(io)=0.5d0*sdata(io)
else
do iu=1,rtap
arg=slo(iu)*om(io)*r(i)
......
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