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

do not expose misleading term "green" to the user

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.
updated copyright notice


SVN Path:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/trunk
SVN Revision: 3352
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent a095c9dc
......@@ -26,8 +26,6 @@
#
# grepg: display Fourier Bessel expansion coefficients
# like created by greda or syg
# grepg77: second make target linked againts g77 version of
# PGPLOT which may contain additional drivers (e.g.GIF)
#
# libraries/packages needed to compile the code
# f2c: Fortran to C converter from netlib
......
......@@ -2,11 +2,7 @@ c this is file <grepg.f>
c======================================================================
c $Id$
c
c GREPG.F
c
c Thomas Forbriger
c Institut fuer Geophysik
c Universitaet Stuttgart
c Copyright 1997,2010 by Thomas Forbriger (IfG Stuttgart)
c
c ----
c This program is free software; you can redistribute it and/or modify
......@@ -24,7 +20,7 @@ c along with this program; if not, write to the Free Software
c Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
c ----
c
c show amplitude plot of greens function matrix
c show amplitude plot of Fourier-Bessel expansion matrix
c
c 25/06/97 V1.0 first running version (from amppg.f)
c
......@@ -115,6 +111,7 @@ c affecting the CVS
c V2.38 provide adaptive scaling along time axis
c (option -ts, flag adaptivetimescaling)
c V2.39 04/12/09 use correct DIN notation for units
c V2.40 13/11/10 do not expose the term "green" to the user
c
c======================================================================
program grepg
......@@ -128,7 +125,7 @@ c
character titleformat*130
c program version:
character*79 version
parameter(version='GREPG V2.39 plot green amplitudes')
parameter(version='GREPG V2.40 plot Fourier-Bessel coefficients')
character*120 CVS_VERSION
parameter(CVS_VERSION=
......@@ -264,13 +261,13 @@ c give basic information
call getarg(1, filename)
if (filename(1:5).eq.'-help') then
print *,' '
print *,'Plot amplitudes of a greens function matrix.'
print *,'Plot amplitudes of Fourier-Bessel coefficients.'
print *,' '
print *,' filename any tf-type file'
print *,' you may read files with the following'
print *,' ''magic'' signatures:'
print *,' ',
& ' ''1234'' complex omega,p-spectra (aka green)'
& ' ''1234'' complex omega,p-spectra'
print *,' ',
& ' ''123S'' real weight values to be used as'
print *,' ',
......@@ -1526,7 +1523,7 @@ c complex green(maxfreq, maxslo)
real pi2
parameter(pi2=2.*3.141592653589)
c
print *,'read green file ',filename(1:index(filename, ' '))
print *,'read coefficient file ',filename(1:index(filename, ' '))
open(lu, file=filename, form='unformatted', status='old', err=99)
c check byte sex
......@@ -1552,7 +1549,7 @@ c check byte sex
print *,'bytesex not matching - we will have to swap!'
stop 'ERROR: do not know how to do that...'
else
print *,'bytesex read (',incmagic,') ist not green'
print *,'bytesex read (',incmagic,') is not coefficient file'
isgreen=.false.
call tf_bytesex(pcmagic, inmagic, cpu, match)
if (match.eq.1) then
......@@ -1568,7 +1565,7 @@ c check byte sex
isgreen=.false.
isgrereso=.false.
iskspec=.false.
print *,'bytesex read (',incmagic,') ist not grereso green'
print *,'bytesex read (',incmagic,') ist not grereso file'
call tf_bytesex(gcmagic, inmagic, cpu, match)
if (match.eq.1) then
print *,'matching bytesex (gabor output) - good...'
......@@ -1582,7 +1579,7 @@ c check byte sex
isgreen=.false.
isgrereso=.false.
iskspec=.false.
print *,'bytesex read (',incmagic,') ist not grereso green'
print *,'bytesex read (',incmagic,') ist not grereso file'
call tf_bytesex(kcmagic, inmagic, cpu, match)
if (match.eq.1) then
print *,'matching bytesex (wavenumber spectrum) - good...'
......@@ -1595,7 +1592,7 @@ c check byte sex
stop 'ERROR: do not know how to do that...'
else
print *,'bytesex read (',incmagic,
& ') ist not wavenumber green'
& ') ist not wavenumber coefficient file'
isgreen=.false.
isgabor=.false.
iskspec=.false.
......@@ -1618,7 +1615,7 @@ c
c read(lu, err=98, end=97) ((green(i,j), i=1,nom), j=1,nslo)
read(lu, err=98, end=97) ((green(j,i), i=1,nom), j=1,nslo)
close(lu, err=96)
print *,'green file read and closed'
print *,'coefficient file read and closed'
c convert gabor time values to 1.e3s unit (will be scaled back to 1s in the
c main program)
if (isgabor) then
......@@ -1660,10 +1657,10 @@ c read(lu, err=98, end=97) ((green(i,j), i=1,nom), j=1,nslo)
endif
return
99 stop 'ERROR: opening green file'
98 stop 'ERROR: reading green file'
97 stop 'ERROR: reading green file - unexpected end'
96 stop 'ERROR: closing green file'
99 stop 'ERROR: opening coefficient file'
98 stop 'ERROR: reading coefficient file'
97 stop 'ERROR: reading coefficient file - unexpected end'
96 stop 'ERROR: closing coefficient file'
88 stop 'ERROR: reading weight file'
87 stop 'ERROR: reading weight file - unexpected end'
86 stop 'ERROR: closing weight file'
......@@ -1713,7 +1710,7 @@ c
print *,'Wavenumber spectrum (magic number: ',kcmagic,')'
print *,' as used by ''flgevas''.'
print *,' '
print *,' In the case of wavenumber Green files verbose'
print *,' In the case of wavenumber coefficient files verbose'
print *,' messages regarding slowness concern wavenumber.'
call grepg_fourierinfo
return
......
......@@ -3,7 +3,7 @@ c======================================================================
c
c GREPG.F
c
c Thomas Forbriger
c Copyright 1997 by Thomas Forbriger
c Institut fuer Geophysik
c Universitaet Stuttgart
c
......
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