Program: MKLSGRIDNEWTRI
mklsgridnewtri.........create new triangular grid
call: mygrid=mklsgridnewtri(lsopt);
mygrid=mklsgridnewtri(lsopt,radiuslist);
lsopt: (struct)
LOCSMITH options structure as returned by MKLSOPT
radiuslist: (numeric array) [km]
List of radius values: grid layers with the radii specified
in this list are produced instead of evaluatin the radius
limits given in LSOPT.
If no RADIUSLIST is given, the respective definitions in
LSOPT will be used.
It is assumed that RADIUSLIST is a unique list of radius
values.
result: mygrid: a structure variable representing a triangular grid
with properties as defined in LSOPT.
The structure has the following fields:
.type: (string)
defines the grid type.
this is the same as LSOPT.GRID_TYPE
.layers: (numeric array)
Array of layer structures. Each element of the
.layers field is a complete LAYER structure that
defines a triangular grid for a single radius.
.layers(1) is the innermost layer.
For details of LAYER structures see the
MKLSGRIDNEWTRILAYER routine.
in case of problems, the grid variable will remain empty.
Generates a grid according to specifications in LSOPT.
The initial layers will be iterated until the maximum circumcircle radius is
smaller than the node spacing defined by LSOPT.GRID_INIT_WDT.
(To save time, the necessary number of iterations is estimated from precomputed
values given by MKLSCONSTANTS, grid mesh width is not evaluated during construction.
This means that the initial mesh width may be realized only approximatively.)
NOTE that the value of the grid_init_wdt field overrides the value of the
iter_circmin field: the inital grid will be iterated until the
grid_init_wdt value of circumcircle radii is met, iter_circmin is
ignored (in fact, it is set to 0 for this purpose)!
Martin Knapmeyer, 08.09.2004, 26.01.2005, 05.07.2005
Read M-File Source Code
eof.