written by Kris Kuhlman (klkuhlm at sandia dot gov)
The data (.dat) and Matlab (.m) files in this directory comprise a set
of MATLAB scripts that are the implementation of a universal kriging
program that computes estimates of potential via cokriging, taking
into account boundary conditions
Copyright (c) 2010 Kristopher L. Kuhlman (klkuhlm at sandia dot gov)
Permission is hereby granted, free of charge, to any person obtaining a copy
of this software and associated documentation files (the "Software"), to deal
in the Software without restriction, including without limitation the rights
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
copies of the Software, and to permit persons to whom the Software is
furnished to do so, subject to the following conditions:
The above copyright notice and this permission notice shall be included in
all copies or substantial portions of the Software.
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
THE SOFTWARE.
this function is part of the set of MATLAB scripts
comprising the implementation of a universal kriging
program that computes estimates of potential via
cokriging, taking into account boundary conditions.
Kuhlman, K.L., and E. Pardo-Igúzquiza, 2010. Universal
cokriging of hydraulic heads accounting for boundary
conditions, Journal of Hydrology, 384(1–2), 14–25.
http://dx.doi.org/10.1016/j.jhydrol.2010.01.002
The files in this directory are described below.
lhsu.m
------
MATLAB function for performing latin hypercube sampling
(modified from Jasper Vrugt's SCEM-UA code)
generate_synthetic_data.m
--------------------------
MATLAB script that creates synthetic data used in toy examples in JoH
paper.
test0{1,2}_{{ch,nf}vecs,head,{ns,ew}locs}.dat
---------------------------------------------
output from generate_synthetic_data.m
ukrig_fd_gradient.m, ukrig_true_gradient.m
------------------------------------------
driver MATLAB scripts that setup the universal kriging system, solve it
and plot the results. Calls covar_*.m and solve_*_gradient_system.m
solve_fd_gradient_system.m, solve_true_gradient_system.m
--------------------------------------------------------
MATLAB functions that build up the A matrix and b vectors to be solved
for the universal kriging systems in both the finite difference
gradient approximation and the true calculus-based gradient approach.
covar_derivs.m, covar_fcns.m
----------------------------
MATLAB functions defining the function and first directional
derivatives of the Gaussian, Cauchy, and Matern variogram models.
These functions are called by the solve_*_gradient_system.m Matlab
scripts.
compare_finite_diff_true_covariances.m
---------------------------------------
MATLAB script that compares the finite-difference and calculus
versions of the three covariance functions for first and second
derivatives. Used to create one of the figures in the JoH paper.
sample_steady_20.dat
--------------------
A random sample of 20 locations and flow model-generated heads from
eulogio_case/STEADY.DAT that were used to compare to kriged heads in
Figures 13, 14, and 15 of the JoH paper.
eulogio_case
------------
directory containing SS flow and Spain examples in JoH
paper.
eulogio_case/MAY1970.DAT
------------------------
data file containing three fixed-width columns of floating-point data.
The columns are X, Y and head (in meters). This data are plotted in
Figure 16 of the JoH paper.
eulogio_case/domain_boundary.dat
--------------------------------
data file containing two fixed-width columns of floating-point data
that are X,Y coordinates of the boundary of the domain shown in Figure
16 of the JoH paper.
eulogio_case/locs.dat
---------------------
data file containing two space-delimited columns of floating-point data
that are X,Y coordinates of the locations on the boundary of the
domain shown in Figure 16 of the JoH paper where boundary normals were
imposed (locations of blue arrows in Figure 16).
eulogio_case/vecs.dat
---------------------
data file containing two space-delimited columns of floating-point data
that are X,Y components of the boundary vectors shown in Figure 16 of
the JoH paper where boundary normals were imposed (direction of blue
arrows in Figure 16)
eulogio_case/STEADY.DAT
-----------------------
data file containing three fixed-width columns, two integer row/column
indices (for a 20x20 domain) and a floating-point head value. The
domain has a constant head in columns 0 and 21, and no-flow at rows 0
and 21.