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

improved dise functionality to compare a set of traces by a script

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: 246
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent f027ba1b
c this is <offli.f>
c------------------------------------------------------------------------------
c ($Source: /home/tforb/svnbuild/cvssource/CVS/thof/src/ts/hd/offli.f,v $)
c ($Id: offli.f,v 1.2 2000-05-30 08:46:42 thof Exp $)
c ($Id: offli.f,v 1.3 2000-09-19 16:05:13 thof Exp $)
c
c 30/05/2000 by Thomas Forbriger (IfG Stuttgart)
c
......@@ -9,15 +9,16 @@ c OFFset LIsting
c
c REVISIONS and CHANGES
c 30/05/2000 V1.0 Thomas Forbriger
c 20/09/2000 V1.1 introduced chaco option
c
c==============================================================================
c
program offli
c
character*(*) version
parameter(version='OFFLI V1.0 OFFset LIsting')
parameter(version='OFFLI V1.1 OFFset LIsting')
character*(*) OFFLI_CVS_ID
parameter(OFFLI_CVS_ID='$Id: offli.f,v 1.2 2000-05-30 08:46:42 thof Exp $')
parameter(OFFLI_CVS_ID='$Id: offli.f,v 1.3 2000-09-19 16:05:13 thof Exp $')
c
integer ierr,nsamp,msamp,nstack
parameter(msamp=100000)
......@@ -39,18 +40,18 @@ c
c commandline
integer maxopt, lastarg, iargc
character*80 argument, infile
parameter(maxopt=3)
parameter(maxopt=4)
character*2 optid(maxopt)
character*40 optarg(maxopt)
logical optset(maxopt), opthasarg(maxopt)
c
logical refmetform
logical refmetform,chacoform
c debugging
logical debug, verbose
c here are the keys to our commandline options
data optid/2h-d, 2h-v,2h-r/
data opthasarg/3*.FALSE./
data optarg/3*1h-/
data optid/2h-d, 2h-v,2h-r,2h-c/
data opthasarg/4*.FALSE./
data optarg/4*1h-/
c
c------------------------------------------------------------------------------
c basic information
......@@ -60,7 +61,7 @@ c
if (iargc().eq.1) call getarg(1, argument)
if ((argument(1:5).eq.'-help').or.(iargc().lt.1)) then
print *,version
print *,'Usage: offli filename [-r]'
print *,'Usage: offli filename [-r|-c]'
print *,' or: offli -help'
if (argument(1:5).ne.'-help') stop 'ERROR: wrong number of arguments'
print *,' '
......@@ -69,6 +70,7 @@ c
print *,'filename SFF data file to list offsets from'
print *,' '
print *,'-r write file appropriate for refmet'
print *,'-c write file appropriate for chaco'
print *,' '
print *,OFFLI_CVS_ID
stop
......@@ -83,6 +85,7 @@ c
debug=optset(1)
verbose=optset(2)
refmetform=optset(3)
chacoform=optset(4)
c
c------------------------------------------------------------------------------
c go
......@@ -108,6 +111,10 @@ c
do itrace=1,ntrace
print 53,offset(itrace)*1.e-3
enddo
elseif (chacoform) then
do itrace=1,ntrace
print 54,offset(itrace)
enddo
else
print 50,ntrace,infile(1:index(infile,' ')-1)
do itrace=1,ntrace
......@@ -121,6 +128,7 @@ c
52 format(a,/,a,/,'mode, v_red, Tl, Tr, N',/
& ,' 1, 0., 0., 0.,',i4,//,11x,'r, phi')
53 format(f12.6,', 0.')
54 format(f12.6,2(',0.'))
end
c
c ----- END OF offli.f -----
# ---------------------------------------
#
# $Header: /home/tforb/svnbuild/cvssource/CVS/thof/src/ts/wf/Makefile,v 1.4 2000-04-17 21:14:51 thof Exp $
# $Header: /home/tforb/svnbuild/cvssource/CVS/thof/src/ts/wf/Makefile,v 1.5 2000-09-19 16:05:13 thof Exp $
#
# Makefile fuer tools /src/ts/wf
#
......@@ -50,11 +50,13 @@ fesig: fesig.o
$(CC) $(CFLAGS) fesig.o -o fesig -ltf -lf2cstuff \
$(PGPLOTLIB) $(F2CLIB) -L$(LOCLIBDIR)
newprog fesig
tesiff teswf: %: %.o
evelo tesiff teswf: %: %.o
$(CC) $(CFLAGS) $< -o $@ -ltf -lf2cstuff $(F2CLIB) -L$(LOCLIBDIR)
newprog $@
evelo: evelo.o
$(CC) $(CFLAGS) $< -o $@ -ltf -lf2cstuff $(F2CLIB) -L$(LOCLIBDIR)
dise: %: %.o
$(CC) $(CFLAGS) $< -o $@ -ltf -lsffu -ltime -lf2cstuff $(F2CLIB) \
-L$(LOCLIBDIR)
newprog $@
......@@ -7,13 +7,14 @@ c DIfferential SEismograms
c
c REVISIONS and CHANGES
c 19/04/99 V1.0 Thomas Forbriger
c 20/09/00 V1.1 set matching source time
c
c==============================================================================
c
program dise
c
character*79 version
parameter(version='DISE V1.0 DIfferential SEismograms')
parameter(version='DISE V1.1 DIfferential SEismograms')
c
c commandline
integer maxopt, lastarg, iargc
......@@ -147,7 +148,8 @@ c
last=.false.
print *,'going to write A-B, A and B to ',
& outfile(1:index(outfile, ' ')),'...'
call sff_WOpen(lu, outfile, ierr)
call sff_WOpenS(lu, outfile, 'nil ',
& 'C', 0., 0., 0., '990420','000000.000', ierr)
if (ierr.ne.0) stop 'ERROR: opening output file'
c
call sff_PrepWID2(nout, 1./dt1, 'NSP ', 1999, 4, 20, 0, 0,
......
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