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

fidasex operates on SU data

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/branches/su1
SVN Revision: 3689
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 31923021
......@@ -25,6 +25,7 @@
# REVISIONS and CHANGES
# 08/04/2009 V1.0 Thomas Forbriger
# 15/11/2010 V1.1 reworked makefile; avoid makedoc.tcl and incdep
# 14/01/2011 V1.2 implemented fidasex with libfapidxx interface
#
# ============================================================================
#
......@@ -84,10 +85,18 @@ fidase.doc: $(patsubst %.o,%.f,$(FIDAOBS))
$(FC) -o $@ -c $< $(FFLAGS)
fidase: $(FIDAOBS)
$(FC) -o fidase $(FIDAOBS) \
$(FC) -o $@ $^ \
-ltf -lts -lsffu -ltime_trad -lsff \
$(PGPLOTLIB) $(LINLIB) \
-L$(LOCLIBDIR)
/bin/mv -fv $@ $(LOCBINDIR)
fidasex: $(FIDAOBS)
$(FC) -o $@ $^ \
-ltf -lts -lsffu -ltime_trad \
-lfapidxx -ldatrwxx -lsffxx -lgsexx -ltime++ -laff \
$(PGPLOTLIB) $(LINLIB) \
-L$(LOCLIBDIR)
/bin/mv -fv $@ $(LOCBINDIR)
# ----- END OF Makefile -----
......@@ -32,6 +32,7 @@ c 23/10/01 V1.4 destack option
c 06/12/02 V1.5 support flgevask inversion inv1d
c 09/09/04 V1.6 time window
c 05/02/10 V1.7 explain trace selection
c 14/01/11 V1.8 use libfapidxx interface to select file format
c
c==============================================================================
c
......@@ -39,7 +40,7 @@ c
c
character*79 version
parameter(version=
& 'FIDASE V1.7 FIt DAtaSets - make datasets homogeneous')
& 'FIDASE V1.8 FIt DAtaSets - make datasets homogeneous')
c
c common blocks
include 'fidase_para.inc'
......@@ -49,6 +50,7 @@ c common blocks
c options
logical overwrite, opttaper, optenfit, optdestack
logical optrenamestat, optwritesingle
character*80 optfileformat
c
real enfminoff, enfmindelta, enfexpo
logical optresamp
......@@ -61,15 +63,15 @@ c filenames
c commandline
integer maxopt, lastarg, iargc
character*80 argument
parameter(maxopt=12)
character*2 optid(maxopt)
parameter(maxopt=13)
character*3 optid(maxopt)
character*40 optarg(maxopt)
logical optset(maxopt), opthasarg(maxopt)
c here are the keys to our commandline options
data optid/2h-D,2h-v,2h-o,2h-t,2h-e,2h-d,2h-r,2h-E,2h-s,2h-R,2h-S,
& 2h-w/
data opthasarg/3*.FALSE.,.TRUE.,4*.TRUE.,3*.FALSE.,.TRUE./
data optarg/3*1h-,1h-,1h-,3hx11,1h-,2h1.,4*1h-/
& 2h-w,'-ty'/
data opthasarg/3*.FALSE.,.TRUE.,4*.TRUE.,3*.FALSE.,2*.TRUE./
data optarg/3*1h-,1h-,1h-,3hx11,1h-,2h1.,4*1h-,'sff'/
c
c------------------------------------------------------------------------------
c basic information
......@@ -78,13 +80,17 @@ c
c print *,'iargc: ',iargc()
if (iargc().eq.1) call getarg(1, argument)
c
if ((argument(1:5).eq.'-help').or.(iargc().lt.1)) then
if ((argument(1:6).eq.'-xhelp').or.(iargc().lt.1)) then
call sff_help_details
stop
elseif ((argument(1:5).eq.'-help').or.(iargc().lt.1)) then
print *,version
print *,'Usage: fidase [-D] [-v] [-o] [-t file] [-e min,mind]'
print *,' [-d device] [-r n,dt] [-E exp] [-s]'
print *,' [-R] [-S] [-w l,r]'
print *,' [-R] [-S] [-w l,r] [-ty f]'
print *,' file1 [t:list] file2 [t:list] ... target'
print *,' or: fidase -help'
print *,' or: fidase -xhelp'
if (iargc().lt.1) stop 'ERROR: missing arguments'
print *,' '
print *,'FIt DAtaSets - make datasets homogeneous'
......@@ -92,6 +98,7 @@ c
print *,'-D Print debug information.'
print *,'-v Be verbose.'
print *,'-o Overwrite existing output file.'
print *,'-ty f select data file format f'
print *,'-t file Read quad taper from file.'
print *,'-e min,mind Fit amplitudes to smooth energy decay.'
print *,' Seismograms with receiver offset larger than'
......@@ -141,6 +148,8 @@ c
print *,' t:9,8,10,14 will select traces 8, 9, 10 and 14'
print *,' '
print *,'$Id$'
print *,' '
call sff_help_formats
stop
endif
c
......@@ -167,6 +176,7 @@ c
optrenamestat=optset(10)
optwritesingle=optset(11)
optwinset=optset(12)
optfileformat=optarg(13)
if (optwinset) read (optarg(12), *, err=99) optwinl,optwinr
c
if ((iargc()-lastarg).lt.2) stop 'ERROR: too few filenames'
......@@ -175,6 +185,8 @@ c
c
c------------------------------------------------------------------------------
c go
call sff_select_format(optfileformat, i)
if (i.ne.0) stop 'ERROR: selecting file format'
call readdata(lastarg)
c
if (optdestack) call destack
......
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