Commit 050077d2 authored by thomas.forbriger's avatar thomas.forbriger Committed by thomas.forbriger
Browse files

proceeding

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: 2286
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent acd0d796
c this is <gresynoise.f>
c ----------------------------------------------------------------------------
c ($Id: gresynoise.f,v 1.7 2007-05-15 16:06:37 tforb Exp $)
c ($Id: gresynoise.f,v 1.8 2007-05-21 10:08:44 tforb Exp $)
c
c Copyright (c) 2007 by Thomas Forbriger (BFO Schiltach)
c
......@@ -13,17 +13,43 @@ c 09/05/2007 V1.0 Thomas Forbriger
c 10/05/2007 V1.1 introduced randomization
c 11/05/2007 V1.2 allow selection of number of output samples
c 15/05/2007 V1.3 reordered loops; phase is of correct type now
c 21/05/2007 V1.4i provide longer seismograms by keep multiple copies of
c noisymized Fourier coefficients
c not yet finished
c
c ============================================================================
c
c - calculate impulse response (subsurface's Green's function)
c for each receiver
c the sampling of Fourier coefficients is already defined by the parameters
c used together with syg; we can refine the sampling in the time domain by
c zero-padding the spectra; we can extend the time window of seismograms by
c stringing multiple compies of noisymize coefficients
c
c we need space for
c maxset sets of Fourier coefficients
c maxom samples of each set of Fourier coefficients
c 2*maxom samples in the Fourier transform array
c (maxset+1)*maxom samples for the output time series
c
c actual numbers:
c nset selected through command line
c must be: nset <= maxset
c nom number of Fourier coefficients to store in each set
c selected through Green's coefficient file
c nsamp number of samples to obtain from each set
c largest values with nsamp <= 2*maxom
c noutsamp number of output samples
c (nset+1)*nsamp/2
c
program gresynoise
c
character*(*) version
parameter(version=
&'GRESYNOISE V1.3 calculate noise seismograms')
&'GRESYNOISE V1.4i calculate noise seismograms')
character*(*) GRESYNOISE_CVS_ID
parameter(GRESYNOISE_CVS_ID=
&'$Id: gresynoise.f,v 1.7 2007-05-15 16:06:37 tforb Exp $')
&'$Id: gresynoise.f,v 1.8 2007-05-21 10:08:44 tforb Exp $')
c
c dimensions
integer maxtr, maxsamp, maxom, maxu
......@@ -510,7 +536,7 @@ c prepare wid2line
if (ierr.ne.0) stop 'ERROR: preparing WID2 line'
call sff_WTraceI(lu, wid2line, nsamp, Rfdata, Ridata, last,
& cs, c1, c2, c3, nstack, ierr)
if (debug) print *,'did writing'
if (debug) print *,'written'
if (ierr.ne.0) stop 'ERROR: writing trace'
c
stop
......
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