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

implemented libfapidxx reading

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: 3653
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 6654ca1c
......@@ -69,7 +69,7 @@
# ============================================================================
#
all: greda gabor phadi hivexpanco hivgabor grereso
all: greda gabor phadi hivexpanco hivgabor grereso gaborx phadix gredax
#----------------------------------------------------------------------
# standard edit targets
......
......@@ -34,6 +34,7 @@ c 25/07/2010 V1.5 Zacharias found a bug: There was a parameter
c missing in a call to pgpt - I wonder how this
c went unnoticed over so many years...
c 27/07/2010 V1.6 correct plot labels
c 30.12.2010 V1.7 implemented libfapidxx
c
c==============================================================================
c
......@@ -41,15 +42,15 @@ c
c
character*(*) version
parameter(version=
& 'PHADI V1.6 calculate dispersion relation from phase differences')
& 'PHADI V1.7 calculate dispersion relation from phase differences')
character*(*) PHADI_CVS_ID
parameter(PHADI_CVS_ID=
& '$Id$')
c input dataset
character*80 filename
character*80 filename, informat
integer maxtraces, totmaxsamples
parameter(maxtraces=300, totmaxsamples=2000000)
integer lu
integer lu, ierr
parameter(lu=12)
real data(totmaxsamples)
integer idata(totmaxsamples)
......@@ -98,7 +99,7 @@ c constants
c commandline
integer maxopt, lastarg, iargc
character*80 argument
parameter(maxopt=22)
parameter(maxopt=23)
character*3 optid(maxopt)
character*80 optarg(maxopt)
logical optset(maxopt), opthasarg(maxopt)
......@@ -107,12 +108,12 @@ c debugging
c here are the keys to our commandline options
data optid/3h-DE, 2h-v, '-P', '-Pf', '-Po', '-d', '-p', '-pf', '-po',
& '-pt', '-D', '-Df', '-Do', '-Dv', '-f', '-Dt',
& '-DO', '-Dn', '-Dr', '-Dp','-De', '-Dl'/
& '-DO', '-Dn', '-Dr', '-Dp','-De', '-Dl', '-ty'/
data opthasarg/3*.FALSE.,3*.TRUE.,.FALSE.,2*.TRUE.,2*.FALSE.,
& 2*.TRUE.,.FALSE.,5*.TRUE.,2*.FALSE.,.true./
& 2*.TRUE.,.FALSE.,5*.TRUE.,2*.FALSE.,2*.true./
data optarg/3*1h-,'10.,20.','phasor.out','x11','-','10.,20.',
& 'phase.out',2*'-','10.,20.','disan.out','-','10.,20.',
& '0.', '1.e-4', '10', '1,1',2*'-','0.'/
& '0.', '1.e-4', '10', '1,1',2*'-','0.','sff'/
c
c------------------------------------------------------------------------------
c basic information
......@@ -128,7 +129,13 @@ c
print *,' [-D] [-Df f1,f2] [-Do file] [-Dv]'
print *,' [-Dt val] [-DO eps] [-Dn N]'
print *,' [-Dr Nl,Nr] [-Dp] [-De] [-Dl f]'
print *,' [-ty format]'
print *,' or: phadi -help'
print *,' or: phadi -xhelp'
if (argument(1:6).eq.'-xhelp') then
call sff_help_details
stop
endif
if (argument(1:5).ne.'-help') stop 'ERROR: wrong number of arguments'
print *,' '
print *,'calculate dispersion relation from phase differences'
......@@ -139,6 +146,7 @@ c
print *,'-DE print debuggin output'
print *,'-d device PGPLOT output device'
print *,'-f f1,f2 default frequency interval'
print *,'-ty format select file format (see list below)'
print *,' '
print *,'-P execute phasor test'
print *,'-Pf f1,f2 frequency interval for phasor test'
......@@ -172,6 +180,8 @@ c
print *,' '
call pgp_showdevices
print *,' '
call sff_help_formats
print *,' '
print *,PHADI_CVS_ID
stop
endif
......@@ -215,11 +225,14 @@ c
fitplot=optset(20)
errorplot=optset(21)
read(optarg(22), *) disanthreshlim
informat=optarg(23)
call getarg(1, filename)
c
c------------------------------------------------------------------------------
c go
call sff_select_input_format(informat, ierr)
if (ierr.ne.0) stop 'ERROR: selecting input file format'
call sffu_simpleread(lu, filename, maxtraces, totmaxsamples, data, idata,
& toffset, tracedt, roffset, innsamples, firstsample, ntraces, verbose)
c
......
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