... | ... | @@ -7,38 +7,23 @@ Location |
|
|
**gremlin** is located at
|
|
|
[src/green/gremlin1](/Seitosh/Seitosh/tree/master/src/green/gremlin1).
|
|
|
|
|
|
Purpose
|
|
|
-------
|
|
|
Joint inversion of Fourier-Bessel expansion coefficients and refracted P-wave
|
|
|
travel times.
|
|
|
|
|
|
### Array size adjustment
|
|
|
Array size adjustment
|
|
|
---------------------
|
|
|
The Fortran 77 program uses fixed size array.
|
|
|
Since it has to store the complete matrix of partial derivatives, memory requirements may become significant.
|
|
|
Since it has to store the complete matrix of partial derivatives, memory
|
|
|
requirements may become significant.
|
|
|
Array sizes therefore are to be carefully adjusted.
|
|
|
In particular the number of slowness values and frequencies in the input data (Fourier Bessel expansion coefficients) must not exceed the array size.
|
|
|
In case you experience problems, please use the program
|
|
|
[[gredim|trunk/src/green/gremlin1#gredim]]
|
|
|
|
|
|
## Problems
|
|
|
### Problems with broken library header files
|
|
|
If there are any problems with libraries(XXX macro: "br")
|
|
|
1) check your remcmmnt mechanism (see [wiki:docs/installation/remcmmnt])
|
|
|
|
|
|
printenv | grep REMCMMNT
|
|
|
|
|
|
if the output isn't
|
|
|
|
|
|
TF_REMCMMNT=cat
|
|
|
|
|
|
adjust settings in ~/.TFbashrc (see [wiki:docs/installation] point 1). (see also ticket:72)
|
|
|
2.) Call
|
|
|
|
|
|
make reinstall
|
|
|
|
|
|
in your directories (e.g. src/libs)(XXX macro: "br")
|
|
|
3.) Call
|
|
|
|
|
|
make install
|
|
|
|
|
|
in the same directories (e.g. src/libs)
|
|
|
In particular the number of slowness values and frequencies in the input data
|
|
|
(Fourier Bessel expansion coefficients) must not exceed the array size.
|
|
|
In case you experience problems, please use the program `gredim`.
|
|
|
|
|
|
Known Problems
|
|
|
--------------
|
|
|
### Problems with array size
|
|
|
If you encounter problems due to data size being too large ('data exceeds array bounds'), you may like to adjust
|
|
|
the array size used in gremlin.
|
... | ... | @@ -48,7 +33,8 @@ You can reduce the number of slowness values by an appropriate parameter to |
|
|
[[grepg|trunk/src/green/grepg]] or you might like to resample your frequency axis by applying
|
|
|
[[greto|trunk/src/green/tools#greto]] to the expansion coefficient data file.
|
|
|
|
|
|
## Documentation
|
|
|
Documentation
|
|
|
-------------
|
|
|
### Application
|
|
|
Documentation regarding the application of gremlin is collected at
|
|
|
wiki:docs/processing/Inversion[[OfExpansionCoefficients]].
|
... | ... | @@ -59,25 +45,33 @@ This is the main inversion program. |
|
|
Its application is documented at wiki:docs/processing/Inversion[[OfExpansionCoefficients]].
|
|
|
|
|
|
#### gredim
|
|
|
This is a program to check the fixed array dimensions in the Fortran 77 code of gremlin and the required memory.
|
|
|
This is a program to check the fixed array dimensions in the Fortran 77 code
|
|
|
of gremlin and the required memory.
|
|
|
The program just reports the current memory requirement.
|
|
|
It first reports the array dimensions set in
|
|
|
[[libs/glq_dim.inc|source:trunk/src/green/gremlin1/libs/glq_dim.inc]]
|
|
|
and then it calculates the number of bytes required for each major array and the total number of bytes in RAM to run [[gremlin|trunk/src/green/gremlin1#gremlin]].
|
|
|
This total must be significantly smaller than the total RAM available in your system.
|
|
|
The meaning of the dimensions set are defined by comments in [[libs/glq_dimstd.inc|source:trunk/src/green/gremlin1/libs/glq_dimstd.inc]].
|
|
|
and then it calculates the number of bytes required for each major array and
|
|
|
the total number of bytes in RAM to run
|
|
|
[[gremlin|trunk/src/green/gremlin1#gremlin]].
|
|
|
This total must be significantly smaller than the total RAM available in your
|
|
|
system.
|
|
|
The meaning of the dimensions set are defined by comments in
|
|
|
[[libs/glq_dimstd.inc|source:trunk/src/green/gremlin1/libs/glq_dimstd.inc]].
|
|
|
|
|
|
### Specific functions
|
|
|
#### gremlin
|
|
|
=#### med: edit model within gremlin=
|
|
|
The command `med` saves the current model to a file called `edit.p.mod` and then calls the text editor to allow model modification.
|
|
|
##### med: edit model within gremlin
|
|
|
The command `med` saves the current model to a file called `edit.p.mod` and
|
|
|
then calls the text editor to allow model modification.
|
|
|
Currently the `vim` editor is hardwired in the source code.
|
|
|
After closing the editor, the modified file is read for further use in gremlin.
|
|
|
This function is implemented at:
|
|
|
source:trunk/src/green/gremlin1/gremlin.f@3599:720-721#L705
|
|
|
|
|
|
=#### dgr: display expansion coefficients used as data=
|
|
|
The command `dgr` displays the modulus (amplitude) of the expansion coefficients which are read from file and used within gremlin as data in four different displays:
|
|
|
##### dgr: display expansion coefficients used as data
|
|
|
The command `dgr` displays the modulus (amplitude) of the expansion
|
|
|
coefficients which are read from file and used within gremlin as data in four
|
|
|
different displays:
|
|
|
1. expansion coefficients as read from file, full frequency and slowness range
|
|
|
1. expansion coefficients as read from file, selected frequency and slowness range is displayed only
|
|
|
1. expansion coefficients after application of scaling function selected by *green data mode*, full frequency and slowness range
|
... | ... | @@ -85,7 +79,7 @@ The command `dgr` displays the modulus (amplitude) of the expansion coefficients |
|
|
|
|
|
### Tips and tricks
|
|
|
#### gremlin
|
|
|
=#### What can I do if I do not have refracted P-wave traveltimes available?=
|
|
|
##### What can I do if I do not have refracted P-wave traveltimes available?
|
|
|
If you do not like to use the travel-time information at all, just set bala
|
|
|
to 1.
|
|
|
- first enter command "spa"
|
... | ... | |