Program: MKLSGRIDNEWTRILAYER
mklsgridnewtrilayer.........construct new layer for triangular grid
call: layer=mklsgridnewtrilayer(lsopt,radius);
lsopt: LOCSMITH options structure as returned by MKLSOPT
radius: desired layer radius [km]
result: layer: a new layer structure, representing a single layer of
a triangular grid at the zero-th iteration.
The structure contains the following fields:
.id: (number)
layer ID.
this may be used to uniquely identify the layer.
.radius: layer radius (number) [km]
each layer consists of a number of points located on
a sphere with this radius
.times: (array structure)
See MKLSGRIDTIMESLIST for details.
.tri: (array structure)
triangle list structure as returned by MKLSGRIDNEWTRILIST.
See there for a description.
The triangle iteration counter will be initialized as
zero. All triangles will be active.
.nodes: (numeric array) [ID deg deg]
three column array containing the longitude and
latituide coordinates of all nodes used by this
layer.
nodes(i,1) is the ID of the i-th node in the list
nodes(i,2) is the longitude of the i-th node [deg]
nodes(i,3) is the latitude of the i-th node [deg]
.fit: (numeric array)
Array containing the maximum quality measure for each node
that is reached in the origin time intervals specified
by the .TIMES field. This is a 2D array consisting
of an ID that points to the corresponding node in
space and a quality value that corresponds
to the origin times given in the .TIMES field.
fit(i,1) is the ID of the node for which the quality
values in fit(i,2:end) are valid.
fit(i,2) is the quality value of the node with ID
fit(i,1) at origin times in the element of
.TIMES with the same node ID.
Upon initilaization, all quality values are set to 0.
.tested: (numeric array) [ID counter]
two column array showing with how many of the
available data values a node has been tested for
compatibility.
.tested(i,1) is the ID of th i-th node in the list
.tested(i,2) is the number of tests applied to the
node with ID tested(i,1).
Since this is a counter, elements for new nodes are
initialized with 0, and new noes may be identified by
searching for .tested(:,2)==0.
In case of problems, the result is empty (not an empty structure,
but simply empty: layer=[].
Martin Knapmeyer, 09.09.2004, 03.12.2004, 10.02.2005
Read M-File Source Code
eof.