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

correction of definition

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/trunk
SVN Revision: 2975
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 94fa4831
......@@ -6,8 +6,17 @@ c 03/03/99 by Thomas Forbriger (IfG Stuttgart)
c
c evaluate polynomial model description at any depth
c
c 11.2.2010: changed definition of a(3) return value
c I understand that a(2) should be the first derivative and a(3) should be the
c second derivative. This was not the case. I assume that this was a
c programming error. Currently the subroutine is only used in function
c mod_writeold and this uses only the return value a(1) and withx in the
c program moval.f which is not part of the inversion process. For this reason
c I modified the definition without hesitation.
c
c REVISIONS and CHANGES
c 03/03/99 V1.0 Thomas Forbriger
c 11/02/10 V1.1 modified definition of a(3)
c
c==============================================================================
c
......@@ -17,6 +26,9 @@ c mindex: index to model space
c pindex: index to model parameter
c depth: depth to calculate parameter value at
c a: value return-array
c a(1): value at depth
c a(2): first derivative at depth
c a(3): second derivative at depth
c
include 'glq_dim.inc'
include 'glq_model.inc'
......@@ -49,7 +61,14 @@ c go
enddo
a(1)=b(1)+b(2)*x+b(3)*x*x
a(2)=b(2)+2.*b(3)*x
a(3)=b(3)
c
c ATTENTION: should'nt this be
c a(3)=2.*b(3)
c
c a(3)=b(3)
c
c 11.2.2010 changed this to:
a(3)=2.*b(3)
c
return
end
......
......@@ -8,13 +8,14 @@ c
c REVISIONS and CHANGES
c 03/03/99 V1.0 Thomas Forbriger
c 24/03/99 V1.1 tell vp/vs
c 11/02/10 V1.2 added note regarding derivatives
c
c==============================================================================
c
program moval
c
character*79 version
parameter(version='MOVAL V1.1 MOdel parameter VALue at given depth')
parameter(version='MOVAL V1.2 MOdel parameter VALue at given depth')
c
include 'libs/glq_dim.inc'
c
......@@ -53,6 +54,10 @@ c
print *,'model name of file containing polynomial model'
print *,'-d depth output model coefficients at a given depth'
print *,'-o file reconvert to OLD format'
print *,' '
print *,'The program returns the value of the parameter as well'
print *,'as the first and second derivative of the polynomial'
print *,'defining the variation of the parameters with depth.'
stop
endif
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