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

refmetx produces SU data 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:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/branches/su1
SVN Revision: 3692
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 78a2dcce
......@@ -39,6 +39,7 @@
# discard f2c fallback option
# check environment
# discard package creation
# 14/01/2011 V1.4 implemented refmetx with libfapidxx interface
#
# ============================================================================
#
......@@ -101,6 +102,7 @@ subclean:
cd main; make clean
refmet: main/refmet.o $(addprefix sub/, $(REFMETOBS))
refmetx: main/refmet.o $(addprefix sub/, $(REFMETOBS))
resus: main/resus.o $(addprefix sub/, $(RESUSOBS))
resusnoise: main/resusnoise.o $(addprefix sub/, $(RESUSNOBS))
refmat: main/refmat.o $(addprefix sub/, $(REFMATOBS))
......@@ -120,6 +122,12 @@ refmatsh refmat resus refmet:
-L$(LOCLIBDIR) $(FFLAGS)
/bin/mv -fv $@ $(LOCBINDIR)
refmetx:
$(FC) $^ -o $@ -lemod -ltf \
-lfapidxx -ldatrwxx -lsffxx -lgsexx -ltime++ -laff \
-L$(LOCLIBDIR) $(FFLAGS)
/bin/mv -fv $@ $(LOCBINDIR)
resusnoise:
$(FC) $^ -o $@ -lsff -lemod -ltf -lgsl -lgslcblas \
-L$(LOCLIBDIR) $(FFLAGS)
......
......@@ -112,12 +112,13 @@ c 04/03/2010 V2.12 line source requires scaling with units of wavenumber to
c provide seismograms in units of 1m. Wavenumber
c is given in units of 1.e-3 rad / m
c V2.12a use factor 1.e3 rather than 1.e-3
c 14/01/2011 V2.13 implement libfapidxx interface
c
c======================================================================
PROGRAM refmet
character*70 version
parameter(version='REFMET V2.12a Reflectivity Method')
parameter(version='REFMET V2.13 Reflectivity Method')
character*79 cvsid
parameter(cvsid='$Id$')
......@@ -744,18 +745,21 @@ c Hankel functions
double precision cl_hankelswitch
c commandline
integer cl_maxopt, cl_lastarg, iargc
parameter(cl_maxopt=9)
character*2 cl_optid(cl_maxopt)
parameter(cl_maxopt=10)
character*3 cl_optid(cl_maxopt)
character*40 cl_optarg(cl_maxopt)
logical cl_optset(cl_maxopt), cl_opthasarg(cl_maxopt)
c command line options and parameter values
integer cl_vlevel
character*80 cl_comsel
character*80 cl_comsel, cl_fileformat
logical cl_debug, cl_rprogress, cl_comeach
c here are the keys to our commandline options
data cl_optid/2h-d,2h-v,2h-o,2h-c,2h-s,2h-1,2h-2,2h-p,2h-l/
data cl_opthasarg/.FALSE.,2*.TRUE.,.FALSE.,.TRUE.,2*.FALSE.,.TRUE.,.FALSE./
data cl_optarg/1h-,1h1,10hrefmet.out,1h-,2h ,2*1h-,2h0.,1h-/
data cl_optid/2h-d,2h-v,2h-o,2h-c,2h-s,2h-1,2h-2,2h-p,2h-l,
& '-ty'/
data cl_opthasarg/.FALSE.,2*.TRUE.,.FALSE.,.TRUE.,2*.FALSE.,
& .TRUE.,.FALSE.,.true./
data cl_optarg/1h-,1h1,10hrefmet.out,1h-,2h ,2*1h-,2h0.,1h-,
& 'sff'/
C**********************************************************************C
C 1. Herdfunktion (als Funktionsdefinition) C
......@@ -813,6 +817,7 @@ c set options
cl_hankel2=cl_optset(7)
read(cl_optarg(8), *) cl_hankelswitch
cl_linesrc=cl_optset(9)
cl_fileformat=cl_optarg(10)
cl_rprogress=.false.
if (cl_vlevel.lt.0) cl_rprogress=.true.
cl_vlevel=abs(cl_vlevel)
......@@ -823,6 +828,10 @@ c get filename
& print *,'WARNING: additional parameters are ignored'
call getarg((cl_lastarg+1), mainfile)
c
c select fapid file format
call sff_select_output_format(cl_fileformat, i)
if (i.ne.0) stop 'ERROR: selecting output file format'
c======================================================================
c
c read configuration from file
......
......@@ -38,11 +38,15 @@ c
c
print *, version
print *, 'Usage: refmet [-d] [-v level] [-o basename] [-c]'
print *, ' [-s select] [-1] [-2] [-p u] file'
print *, ' [-s select] [-1] [-2] [-p u] [-ty f] file'
print *, ' or: refmet -help'
print *, ' or: refmet -xhelp'
if (iargc().lt.1) stop 'ERROR: missing parameters'
call getarg(1, mainfile)
if (mainfile.eq.'-help') then
if (mainfile.eq.'-xhelp') then
call sff_help_details
stop
else if (mainfile.eq.'-help') then
call refmet_intro
print *,' '
print *,'commandline parameters are:'
......@@ -82,6 +86,7 @@ c
print *,' J0(w*u*r) -> 2.*cos(w*u*r)/max(1.d-50,w*u)'
print *,' J1(w*u*r) -> 0'
print *,' J2(w*u*r) -> 0'
print *,'-ty f select output file format f'
print *,' '
print *,'file Is the name of the main configuration file.'
print *,' It contains the names of the three file'
......@@ -112,6 +117,8 @@ c call other info routine
call refmet_modinf
call refmet_sourinf(hfktstr)
call refmet_rcvinf
print *,' '
call sff_help_formats
stop
endif
return
......
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