Commit 7f8ac6b7 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: 3007
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 5bc3beb9
c this is <refmet.f> originally by J. Ungerer 1990 c this is <refmet.f> originally by J. Ungerer 1990
c====================================================================== c======================================================================
c $Id: refmet.f,v 1.18 2010-03-04 12:05:56 tforb Exp $ c $Id: refmet.f,v 1.19 2010-03-04 16:44:58 tforb Exp $
c c
c Reflectivity Method c Reflectivity Method
c c
...@@ -106,14 +106,15 @@ c 03/03/2010 V2.11 correction of line source kernel ...@@ -106,14 +106,15 @@ 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 04/03/2010 V2.12 line source requires scaling with units of wavenumber to
c provide seismograms in units of 1m. Wavenumber c provide seismograms in units of 1m. Wavenumber
c is given in units of 1.e-3 rad / m c is given in units of 1.e-3 rad / m
c V2.12a use factor 1.e3 rather than 1.e-3
c c
c====================================================================== c======================================================================
PROGRAM refmet PROGRAM refmet
character*70 version character*70 version
parameter(version='REFMET V2.12 Reflectivity Method') parameter(version='REFMET V2.12a Reflectivity Method')
character*79 cvsid character*79 cvsid
parameter(cvsid='$Id: refmet.f,v 1.18 2010-03-04 12:05:56 tforb Exp $') parameter(cvsid='$Id: refmet.f,v 1.19 2010-03-04 16:44:58 tforb Exp $')
c array dimension declaration c array dimension declaration
integer me, msl, mf, ms integer me, msl, mf, ms
...@@ -2104,7 +2105,7 @@ c no tangential component, no near field ...@@ -2104,7 +2105,7 @@ c no tangential component, no near field
else else
if (cl_linesrc) then if (cl_linesrc) then
c quick and dirty kernel substitution for line source c quick and dirty kernel substitution for line source
J0=2.d-3*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 J1=0.d0
J2=0.d0 J2=0.d0
else 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