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


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:
SVN Revision: 2369
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 3dc9595d
c this is <gresynoise.f>
c ----------------------------------------------------------------------------
c ($Id: gresynoise.f,v 1.22 2007-07-23 16:55:51 tforb Exp $)
c ($Id: gresynoise.f,v 1.23 2007-09-07 15:09:07 tforb Exp $)
c Copyright (c) 2007 by Thomas Forbriger (BFO Schiltach)
......@@ -40,6 +40,8 @@ c stacked; R component was Z componten plus last
c receivers R component
c R should be correct now
c 23/07/2007 V1.11 - report rms values for all traces
c 07/09/2007 V1.12 - correct label for scaling factors
c - apply scaling for horizontal components
c ============================================================================
......@@ -70,10 +72,10 @@ c
character*(*) version
&'GRESYNOISE V1.11 calculate noise seismograms')
&'GRESYNOISE V1.12 calculate noise seismograms')
character*(*) GRESYNOISE_CVS_ID
&'$Id: gresynoise.f,v 1.22 2007-07-23 16:55:51 tforb Exp $')
&'$Id: gresynoise.f,v 1.23 2007-09-07 15:09:07 tforb Exp $')
c dimensions
integer maxtr, maxsetsamp, maxom, maxu, maxset, maxsamp
......@@ -521,11 +523,12 @@ cc print *,i,io,Zfourier(io),Rfourier(io),phase(io)
c scale trace and stack
c R-component is scaled to be used as N- and E-component equally
cc print *,'rec # ', i, ' scal: ',scal
do io=1,Znom
cc print *,i,io,is,Zfourier(io),Rfourier(io),scal
& scal**2*real(Zfourier(io)*conjg(Zfourier(io)))
......@@ -708,7 +711,7 @@ c
94 stop 'ERROR closing response file'
50 format(//////)
51 format(3x,'trace: ',i3,1x,'offset: ',f10.2,' m',
& 2x,'rnd fac: ',f10.5,2x,'offs. fac: ',e10.3)
& 2x,'rnd fac: ',f10.5,2x,'scaling: ',e10.3)
52 format(5x,'trace: ',i3,' Z:',e12.5,' R:',e12.5)
53 format(5x,'trace: ',i3,' Z:',e12.5)
54 format(5x,'trace: ',i3,15x,' R:',e12.5)
......@@ -885,6 +888,32 @@ c
print *,' duration of the output is (N+1)*T/2, '
& ,'where N is the number of sets used.'
print *,' '
print *,'The horizontal component written by this program is the'
print *,'radial component of the PSV-case. The simulation problem'
print *,'to be solved for H/V studies requires N- and E-component'
print *,'signals, where we assume that signals from all azimuths'
print *,'contribute at the same level. For noise signals, the'
print *,'power of the resulting signal equals the integral over'
print *,'the power of the contributing signals.'
print *,' '
print *,'The power of the vertical component as an average over'
print *,'contributions from all azimuths is:'
print *,' PZ = int_0^(2*pi) pz(phi) d phi / (2*pi),'
print *,'where pz(phi) is the power of the wave contributing from'
print *,'azimuth phi. For the horizontal components this is'
print *,' PN = int_0^(2*pi) cos(phi)**2 pr(phi) d phi / (2*pi)'
print *,'for the north component and'
print *,' PE = int_0^(2*pi) sin(phi)**2 pr(phi) d phi / (2*pi)'
print *,'for the east component respectively if pr(phi) is the'
print *,'radial component in the PSV case.'
print *,' '
print *,'We take pz(phi)=pz, pr(phi)=pr constant for all phi.'
print *,'The power scaling factor for PZ then is 1 and for PN and'
print *,'PE is 1/2. Amplitudes scale with the square root of',
& ' power.'
print *,'To use the radial component written by this program as'
print *,'an effective horizontal component for both N and E,'
print *,'we must scale its amplitude by 1/sqrt(2)=0.707.'
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