Commit 893b379f authored by thomas.forbriger's avatar thomas.forbriger

ref/tools [TASK][API][!!!]: adopt to recent version of model file format

parent d32d52fb
......@@ -95,7 +95,7 @@ c new model
double precision qbeta(maxlayer)
real*8 nuref
c gemini
logical optani
logical optlegacy
integer maxsection, nsection, nl
parameter(maxsection=30)
double precision alpha(maxsection, 4), beta(maxsection, 4)
......@@ -126,7 +126,7 @@ c commandline
character*40 optarg(maxopt)
logical optset(maxopt), opthasarg(maxopt)
c here are the keys to our commandline options
data optid/'-d','-s','-r','-o','-S','-n','-f','-p','-N','-t'/
data optid/'-d','-s','-r','-o','-S','-n','-f','-p','-N','-l'/
data opthasarg/.FALSE.,2*.TRUE.,.FALSE.,.TRUE.,.FALSE.,3*.TRUE.,
& .false./
data optarg/'-','0.02','1000.','-','1.','-','0.','0.','-1','-'/
......@@ -147,6 +147,12 @@ c----------------------------------------------------------------------
print *,'transformed from elastic moduli to seismic layer'
print *,'velocities.'
print *,' '
print *,'In files providing anisotropic parameters with'
print *,'transverse isotropy, h-velocities and eta are ignored.'
print *,' '
print *,'gemini 2.2 (Jörg Dalkolmo, 1997) is available for download at'
print *,'http://www.spice-rtn.org/library/software/GEMINI.html'
print *,' '
print *,'-s stepsize Gives the maximum size of parameter function'
print *,' steps as a fraction of the absolute value.'
print *,' Default is 0.02 which means that every step'
......@@ -180,8 +186,10 @@ c----------------------------------------------------------------------
print *,' n=-2 for Rayleigh waves (Biswas)'
print *,' n=-5 for Love waves (Biswas and Knopoff)'
print *,' (default n=',optarg(9)(1:3),')'
print *,'-t earth model has transverse isotropy'
print *,' (h-velocities and eta are ignored)'
print *,'-l use reading function for legacy format'
print *,' for isotropic earth model data; this'
print *,' format is out-dated and not used with'
print *,' gemini 2.2'
print *,' '
print *,'infile GEMINI model file for spherical earth'
print *,'outfile Output containing model for flat earth'
......@@ -204,7 +212,7 @@ c set options
read(optarg(7), *) nuref
disperse=optset(8)
read(optarg(9), *) densexp
optani=optset(10)
optlegacy=optset(10)
if (disperse) then
if (nuref.le.0.d0) stop 'ERROR: You must set a reference frequency'
read(optarg(8), *) domper
......@@ -226,7 +234,7 @@ c set options
c----------------------------------------------------------------------
c read file
if (debug) print *,'DEBUG: filename ',filename(1:fin)
if (optani) then
if (.not.optlegacy) then
if (debug) print *,'DEBUG: call gemini_getani'
call gemini_getani(filename, lu, maxsection, eta,
& rb, qm, qk, rho, alpha, alphah, beta, betah,
......@@ -236,8 +244,6 @@ c read file
call gemini_getmod(filename, lu, maxsection,
& rb, qm, qk, rho, alpha, beta, nsection, iflso, nco, text)
endif
c call gemini_getmod(filename, lu, maxsection,
c & rb, qm, qk, rho, alpha, beta, nsection, iflso, nco, text)
R=rb(nsection)
c----------------------------------------------------------------------
c go through all sections of polynomal model and chop model
......
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