Commit 71339bbc 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: 2275
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 7f04180b
c this is <gresynoise.f> c this is <gresynoise.f>
c ---------------------------------------------------------------------------- c ----------------------------------------------------------------------------
c ($Id: gresynoise.f,v 1.2 2007-05-10 08:04:43 tforb Exp $) c ($Id: gresynoise.f,v 1.3 2007-05-10 08:06:05 tforb Exp $)
c c
c Copyright (c) 2007 by Thomas Forbriger (BFO Schiltach) c Copyright (c) 2007 by Thomas Forbriger (BFO Schiltach)
c c
...@@ -21,7 +21,7 @@ c ...@@ -21,7 +21,7 @@ c
&'GRESYNOISE V1.0 calculate noise seismograms') &'GRESYNOISE V1.0 calculate noise seismograms')
character*(*) GRESYNOISE_CVS_ID character*(*) GRESYNOISE_CVS_ID
parameter(GRESYNOISE_CVS_ID= parameter(GRESYNOISE_CVS_ID=
&'$Id: gresynoise.f,v 1.2 2007-05-10 08:04:43 tforb Exp $') &'$Id: gresynoise.f,v 1.3 2007-05-10 08:06:05 tforb Exp $')
c c
c dimensions c dimensions
integer maxtr, maxsamp, maxom, maxu integer maxtr, maxsamp, maxom, maxu
...@@ -168,7 +168,7 @@ c ...@@ -168,7 +168,7 @@ c
hankel1=optset(7) hankel1=optset(7)
hankel2=optset(8) hankel2=optset(8)
read (optarg(9), *, err=99) scalexp read (optarg(9), *, err=99) scalexp
optscale=optset(10) optrndscale=optset(10)
c c
pi2=2.d0*pi pi2=2.d0*pi
c---------------------------------------------------------------------- c----------------------------------------------------------------------
...@@ -252,7 +252,7 @@ c scale down epicentral distances to meters ...@@ -252,7 +252,7 @@ c scale down epicentral distances to meters
r(i)=r(i)*1000. r(i)=r(i)*1000.
scalefac(i)=1. scalefac(i)=1.
enddo enddo
if (optscale) call tf_gsl_rnd_gaussian(scalefac, ntr) if (optrndscale) call tf_gsl_rng_ugaussian(scalefac, ntr)
c check stepwidth c check stepwidth
dom=Zom(2)-Zom(1) dom=Zom(2)-Zom(1)
do io=2,Znom do io=2,Znom
......
Supports Markdown
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