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

added slowness limit to boxcar

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: 1161
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 3cfabd3c
c this is <grenorm.f>
c------------------------------------------------------------------------------
c
c $Id: grenorm.f,v 1.7 2002-04-05 20:08:27 forbrig Exp $
c $Id: grenorm.f,v 1.8 2002-12-05 20:40:02 forbrig Exp $
c
c 28/06/99 by Thomas Forbriger (IfG Stuttgart)
c
......@@ -20,13 +20,14 @@ c 23/08/01 V1.6 allow amplitude normalize
c 27/08/01 V1.7 allow amplitude norm with taper file
c 27/03/02 V1.8 support calculation of modulus
c 05/04/02 V1.8a tell more about what we are doing
c 05/12/02 V1.9 added slowness limit to boxcar
c
c==============================================================================
c
program grenorm
c
character*79 version
parameter(version='GRENORM V1.8a GREens function will be NORMalized')
parameter(version='GRENORM V1.9 GREens function will be NORMalized')
c
c data
integer maxslo, maxfreq
......@@ -167,7 +168,7 @@ c
print *,'If -A is selected, the modulus is taken from all'
print *,'coefficients.'
print *,' '
print *,'$Id: grenorm.f,v 1.7 2002-04-05 20:08:27 forbrig Exp $'
print *,'$Id: grenorm.f,v 1.8 2002-12-05 20:40:02 forbrig Exp $'
stop
endif
c
......@@ -332,22 +333,28 @@ c----------------------------------------------------------------------
c
else
print *,' using default taper (boxcar)...'
print *,' slowness limit: ',slo_limit*1.e3
if (opt_amp) then
if (verbose) print *,' normalize to amplitude...'
if (verbose) then
print *,' normalize to amplitude...'
endif
do i=1,nom
intval=0.
do j=1,nslo
intval=max(intval,abs(green(i,j)))
if (slo(j).ge.slo_limit)
& intval=max(intval,abs(green(i,j)))
enddo
intval=max(intval,lowlim)
ofrefac(i)=intval
enddo
else
if (verbose) print *,' normalize to integral...'
if (verbose) then
print *,' normalize to integral...'
endif
do i=1,nom
intval=0.
do j=1,nslo
intval=intval+abs(green(i,j))
if (slo(j).ge.slo_limit) intval=intval+abs(green(i,j))
enddo
intval=max(intval,lowlim)
ofrefac(i)=intval
......@@ -388,7 +395,9 @@ c coslen=float(max(i-k,l-i)+1)
frefac(i)=ofrefac(i)
enddo
endif
if (verbose) print *,' rescale amplitudes...'
do i=1,nom
if (debug) print *,' i,fac',i,1./frefac(i)
do j=1,nslo
green(i,j)=green(i,j)/frefac(i)
enddo
......
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