Commit 2e3d419b 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: 2999
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 4b6c5000
c this is <gresy.f>
c------------------------------------------------------------------------------
c $Id: gresy.f,v 1.20 2010-03-01 15:23:17 tforb Exp $
c $Id: gresy.f,v 1.21 2010-03-02 13:35:57 tforb Exp $
c
c 26/06/97 by Thomas Forbriger (IfG Stuttgart)
c
......@@ -30,12 +30,13 @@ 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 02/03/10 V1.14 corrected (hopefully) line source integration
c
c==============================================================================
c
program gresy
character*79 version
parameter(version='GRESY V1.13 GREens function SYnthetics')
parameter(version='GRESY V1.14 GREens function SYnthetics')
c dimensions
integer maxtr, maxsamp, maxom, maxu
parameter(maxu=4500, maxtr=maxu, maxom=4100, maxsamp=maxom*2)
......@@ -145,7 +146,7 @@ c
print *,' amplitude unit: m**3/s if spectrum represents displacement'
print *,' waveform in m'
print *,' '
print *,'$Id: gresy.f,v 1.20 2010-03-01 15:23:17 tforb Exp $'
print *,'$Id: gresy.f,v 1.21 2010-03-02 13:35:57 tforb Exp $'
print *,' '
call refmet_rcvinf
stop
......@@ -412,7 +413,7 @@ c set slowness taper
& cos(arg)*du*tf_costap(iu,0,0,ltap,rtap)
endif
enddo
sdata(io)=0.5d0*sdata(io)
sdata(io)=2.d0*sdata(io)*om(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