For a parameter vector, ak and a map-space domain-parameter vector, ma, return the lambda-parameterized quality diagnostic, Q, and the geographic domain-parameter vector ga.
Lambda is given by lam <1. Also, return the derivatives, qdak and qdma, of Q wrt ak and ma, and the derivatives gadak and gadma, of ga wrt ak and ma.
The domain averages of Q are accurately computed by bi-Gauss-Legendre quadrature over the positive quadrant of the domain (exploiting the symmetry) of the four constituent terms, v1, v2, v3, v4, from which the mean Q is computed using a quadratic formula of these constituents. The number of Gauss points in eaxh half-interval is ngh, and the nodes themselves are, in proportion to the half-interval, at xg. the normalized gauss weights are wg.
If a failure occurs, colmputations cease immediately and a failure flag, FF, is raised on return.
- Parameters
-
[in] | ngh | |
[in] | lam | Lambda |
[in] | xg | |
[in] | wg | |
[in] | ak | parameter vector |
[in] | ma | map-space domain-parameter vector |
[out] | q | lambda-parameterized quality diagnostic |
[out] | qdak | derivatives value |
[out] | qdma | derivatives value |
[out] | ga | geographic domain-parameter vector |
[out] | gadak | |
[out] | gadma | |
[out] | ff | error flag |
- Author
- R. J. Purser
Definition at line 602 of file pesg.f90.