Commit 4779b1f1 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: 2947
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 10c94c15
c this is <damplo.f>
c------------------------------------------------------------------------------
c ($Id: damplo.f,v 1.2 2001-08-07 17:19:09 forbrig Exp $)
c ($Id: damplo.f,v 1.3 2009-12-04 12:22:11 tforb Exp $)
c
c 07/08/2001 by Thomas Forbriger (IMGF Frankfurt)
c
......@@ -8,15 +8,17 @@ c plot energy damping
c
c REVISIONS and CHANGES
c 07/08/2001 V1.0 Thomas Forbriger
c 04/12/2009 V1.1 use correct DIN notation for units
c
c==============================================================================
c
program damplo
c
character*(*) version
parameter(version='DAMPLO V1.0 plot energy damping')
parameter(version='DAMPLO V1.1 plot energy damping')
character*(*) DAMPLO_CVS_ID
parameter(DAMPLO_CVS_ID='$Id: damplo.f,v 1.2 2001-08-07 17:19:09 forbrig Exp $')
parameter(DAMPLO_CVS_ID=
& '$Id: damplo.f,v 1.3 2009-12-04 12:22:11 tforb Exp $')
c
c input dataset
character*80 filename
......@@ -184,7 +186,7 @@ c go
call pgbox(XOPT,1.,1,YOPT,1.,1)
call pgunsa
call pgslw(2)
call pglab('offset (m)','scaled relative rms',facstring)
call pglab('offset / m','scaled relative rms',facstring)
call pgpt(ntraces,roffset,energy, -4)
call pgmtxt('R',2.5,0.5,0.5,filename)
call pgclos
......
c this is <stuplo.f> by Thomas Forbriger 1996
c
c $Id: stuplo.f,v 1.29 2008-01-18 11:55:40 tforb Exp $
c $Id: stuplo.f,v 1.30 2009-12-04 12:22:11 tforb Exp $
c
c This is a simple plotting tool for seismic time series in
c SFF format
......@@ -60,6 +60,7 @@ c V1.25 29/06/07 provide station location in caption
c V1.26 02/07/07 provide option to set tick intervals
c V1.27 18/01/08 do not use write statement for annotations
c just assign to character variable line
c V1.28 04/12/09 use correct DIN notation for units
c
c======================================================================
program stuplo
......@@ -70,7 +71,7 @@ c
c version
character*77 version, creator
parameter(version=
& 'STUPLO V1.27 plot seismic time series (SFF format)')
& 'STUPLO V1.28 plot seismic time series (SFF format)')
parameter(creator='1996 by Thomas Forbriger (IfG Stuttgart)')
c parameter definitions
integer maxsamples, maxselect, lu, maxtraces, maxchain, maxstyle
......@@ -800,10 +801,10 @@ c
botlabel=' '
if ((trace.eq.ntraces).or.(.not.(optscalex))) then
if (optabstime) then
botlabel='time since midnight (sec)'
botlabel='time since midnight / sec'
if (opttbox) botlabel='time since midnight'
else
botlabel='time since first sample (sec)'
botlabel='time since first sample / sec'
if (opttbox) botlabel='time since first sample'
endif
if (optxlabel) botlabel=parxlabel
......
......@@ -7,11 +7,12 @@ c plot sum of spectra
c
c REVISIONS and CHANGES
c 08/04/99 V1.0 Thomas Forbriger
c 04/12/09 V1.1 use correct DIN notation for units
c
c==============================================================================
program splot
program susplo
character*70 version
parameter(version='SUSPLOT V1.0 plot stacked spectra')
parameter(version='SUSPLOT V1.1 plot stacked spectra')
c declare data variables
integer maxtraces, maxsamples
parameter(maxtraces=50)
......@@ -50,7 +51,7 @@ c
c give basic information
print *,version
if (iargc().ne.4) then
print *,'Usage: splot datafile tapfrac device frac'
print *,'Usage: susplo datafile tapfrac device frac'
print *,' '
print *,' datafile name of data file (sff format)'
print *,' tapfrac tapering fraction for cosine taper'
......@@ -226,7 +227,7 @@ c enddo
call pgupdt
c
c label data
call pglab('frequency [Hz]',' ',text)
call pglab('frequency / Hz',' ',text)
return
end
......
c this is <refract_cmdopt.f>
c------------------------------------------------------------------------------
cS
c ($Id: refract_cmdopt.f,v 2.6 2007-06-21 11:06:47 tforb Exp $)
c ($Id: refract_cmdopt.f,v 2.7 2009-12-04 12:24:12 tforb Exp $)
c
c 24/05/2000 by Thomas Forbriger (IfG Stuttgart)
c
......@@ -30,6 +30,7 @@ c 11/01/2001 V1.2 introduced plflag_hypoffset
c 18/01/2001 V1.3 corrected braces around units
c 17/06/2003 V1.4 introduced new option plflag_tracenum
c 09/09/2004 V1.5 introduced new option plflag_tracename
c 04/12/2009 V1.6 use correct DIN notation for units
c
c==============================================================================
c
......@@ -55,7 +56,8 @@ c declare local variables
integer i
c CVS Id
character*(*) refract_cmdopt_id
parameter (refract_cmdopt_id='$Id: refract_cmdopt.f,v 2.6 2007-06-21 11:06:47 tforb Exp $')
parameter (refract_cmdopt_id=
& '$Id: refract_cmdopt.f,v 2.7 2009-12-04 12:24:12 tforb Exp $')
c commandline
integer maxopt
parameter(maxopt=59)
......@@ -184,11 +186,11 @@ c seismogram scaling
c
c titles, label, legends
opt_Ttitle=version
opt_Txlabel='time (s)'
opt_Txlabel='time / s'
if (opt_Sradius.gt.0.) then
opt_Tylabel='offset (^)'
opt_Tylabel='offset / ^'
else
opt_Tylabel='offset (m)'
opt_Tylabel='offset / m'
endif
if (optset(16)) opt_Ttitle=optarg(16)
if (optset(17)) opt_Txlabel=optarg(17)
......
......@@ -9,6 +9,7 @@ c REVISIONS and CHANGES
c 18/11/98 V1.0 Thomas Forbriger
c 23/11/98 V1.1 allow stacking to arrival at mid of sread
c 23/02/99 V1.2 a factor of two was missing
c 04/12/09 V1.3 use correct DIN notation for units
c
c==============================================================================
c
......@@ -119,7 +120,8 @@ c
enddo
call pgenv(opt_minvel,opt_maxvel,0.,1.,0,1)
write(title, 50) 1./slowestim
call pglab('deconvolution velocity of air-coupled sound wave [m/s]',
call pglab(
& 'deconvolution velocity of air-coupled sound wave / m s\u-1',
& 'normalized rms amplitude of correlated signal', title)
call pgsave
call pgslw(3)
......@@ -356,14 +358,14 @@ c & minvalue, 'shade')
endif
if (opt_pickmode) then
call pglab(xlabel,
& 'deconvolution slowness [s/km]','pick your optimum')
& 'deconvolution slowness / s km\u-1','pick your optimum')
else
if (opt_usetitle) then
call pglab(xlabel,
& 'deconvolution slowness [s/km]',opt_title)
& 'deconvolution slowness / s km\u-1',opt_title)
else
call pglab(xlabel,
& 'deconvolution slowness [s/km]',
& 'deconvolution slowness / s km\u-1',
& 'correlation test of air-coupled sound wave')
endif
endif
......@@ -412,7 +414,7 @@ c
enddo
c
return
50 format('apparent arrival time [sec] at offset ',f6.1,'m')
50 format('apparent arrival time (in seconds) at offset ',f6.1,'m')
end
c
c----------------------------------------------------------------------
......
......@@ -12,13 +12,15 @@ c 28/02/99 V1.2 allow calculation of hilbert transform
c 01/03/99 V1.3 calculate spectral coefficients
c 25/06/01 V1.4 depart from GSE sampling rate standard formst if
c necessary
c 04/12/09 V1.5 corrected comment
c
c==============================================================================
c
program evelo
c
character*79 version
parameter(version='EVELO V1.4 calculate EnVELOpe of time series')
parameter(version=
& 'EVELO V1.5 calculate EnVELOpe of time series')
c
integer maxsamples,nsamp,nspa,i, npow, limit
parameter(maxsamples=310000)
......@@ -96,8 +98,8 @@ c
print *,' '
print *,'The output Fourier spectrum is normalized to be'
print *,'equivalent to the corresponding integral transform.'
print *,'Thus values are spectral amplitude densities'
print *,'with unit [amplitude]/Hz, where [amplitude] is'
print *,'Thus Fourier coeffients and spectral amplitudes are'
print *,'given in [amplitude]/Hz, where [amplitude] is'
print *,'the unit of the input samples.'
stop
endif
......
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