Commit 33469f48 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: 2291
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 0a3f3b97
c this is <gresynoise.f>
c ----------------------------------------------------------------------------
c ($Id: gresynoise.f,v 1.12 2007-05-24 15:12:41 tforb Exp $)
c ($Id: gresynoise.f,v 1.13 2007-05-24 15:18:12 tforb Exp $)
c
c Copyright (c) 2007 by Thomas Forbriger (BFO Schiltach)
c
......@@ -50,7 +50,7 @@ c
&'GRESYNOISE V1.5 calculate noise seismograms')
character*(*) GRESYNOISE_CVS_ID
parameter(GRESYNOISE_CVS_ID=
&'$Id: gresynoise.f,v 1.12 2007-05-24 15:12:41 tforb Exp $')
&'$Id: gresynoise.f,v 1.13 2007-05-24 15:18:12 tforb Exp $')
c
c dimensions
integer maxtr, maxsetsamp, maxom, maxu, maxset, maxsamp
......@@ -385,13 +385,11 @@ c
else
print *,'I will use the Bessel function'
endif
last=.false.
c
c start receiver loop
c -------------------
do i=1,ntr
print *,'working on trace ',i,' at ',r(i),'m'
if (i.eq.ntr) last=.true.
c calculate bessel integral
c----------------------------------------------------------------------
c frequency loop
......@@ -598,12 +596,6 @@ c writing the same trace twice with equivalence buffers requires copying
do i=1,nsamp
Zfdata(i)=Rfdata(i)
enddo
if (optnew) call sff_New(lu, Rseisname, ierr)
if (ierr.ne.0) stop 'ERROR: deleting seismogram file'
call sff_WOpenFS(lu, Rseisname,
& free, nfree, srctype, cs, c1 ,c2 ,c3,
& date, time, ierr)
if (ierr.ne.0) stop 'ERROR: opening seismogram file'
c prepare wid2line
if (debug) print *,'go and write',nsamp
cchannel='N'
......@@ -618,12 +610,6 @@ c prepare wid2line
c east component
c --------------
last=.true.
if (optnew) call sff_New(lu, Rseisname, ierr)
if (ierr.ne.0) stop 'ERROR: deleting seismogram file'
call sff_WOpenFS(lu, Rseisname,
& free, nfree, srctype, cs, c1 ,c2 ,c3,
& date, time, ierr)
if (ierr.ne.0) stop 'ERROR: opening seismogram file'
c prepare wid2line
if (debug) print *,'go and write',nsamp
cchannel='E'
......
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