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

correct scaling of line source

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: 3006
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent d27f4835
c this is <refmet.f> originally by J. Ungerer 1990
c======================================================================
c $Id: refmet.f,v 1.17 2010-03-03 10:55:32 tforb Exp $
c $Id: refmet.f,v 1.18 2010-03-04 12:05:56 tforb Exp $
c
c Reflectivity Method
c
......@@ -103,14 +103,17 @@ c only): tabluated velocity refers to real part of modulus
c 23/02/2009 V2.9 porting code to gfortran
c 02/03/2010 V2.10 implementing line source
c 03/03/2010 V2.11 correction of line source kernel
c 04/03/2010 V2.12 line source requires scaling with units of wavenumber to
c provide seismograms in units of 1m. Wavenumber
c is given in units of 1.e-3 rad / m
c
c======================================================================
PROGRAM refmet
character*70 version
parameter(version='REFMET V2.11 Reflectivity Method')
parameter(version='REFMET V2.12 Reflectivity Method')
character*79 cvsid
parameter(cvsid='$Id: refmet.f,v 1.17 2010-03-03 10:55:32 tforb Exp $')
parameter(cvsid='$Id: refmet.f,v 1.18 2010-03-04 12:05:56 tforb Exp $')
c array dimension declaration
integer me, msl, mf, ms
......@@ -2101,7 +2104,7 @@ c no tangential component, no near field
else
if (cl_linesrc) then
c quick and dirty kernel substitution for line source
J0=2.*cos(jarg)/max(1.d-50,w(f)*u)
J0=2.d-3*cos(jarg)/max(1.d-50,w(f)*u)
J1=0.d0
J2=0.d0
else
......
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