Program: MKLSGRIDWDTDEMOTRI
mklsgridwdtdemotri.......demo of radius dependency of triangular grid mesh width
call: mklsgridwdtdemotri(gridbase,meshwdt,minr,maxr,dr);
gridbase: (string)
string denoting the base body to be used for grid
construction. This must be one of the values alloewd for
the lsopt.grid_base field.
meshwdt: (number) [km]
desired mesh width of the generated grid, in terms of
triangle circumcircle radius, measured in Kilometer.
The generated grid will have no triangles with
circumcircles greater than this.
minr: (number) [km]
Radius of the innermost layer.
maxr: (number) [km]
Radius of the outermost layer.
dr: (number) [km]
delta radius of the grid layers.
The generated grid will have a radius of the outermost
layer of 6371km (earth radius). The deeper layers will
have radii of 6371-n*dr, n=1...6371/dr.
result: none
This routine creates a multilayer triangular grid and iterates each layer
until its largest triangle's circumcircle is smaller than the predefined
mesh width.
Produces three plots:
Figure 1:
A 3D cutaway plot of the layers. Coloring is chosen according to the
number of iterations applied to each layer.
Figure 2:
Distributions of circumcircle radii in each layer, plotted against layer
radius. A Black line shows the maximum circumcircle radius in each layer.
A red line shos the desired mesh width. The red line should be above
the black line. Coloring of radius markers according to colors in Fig. 1
Figure 3:
A plane view of all points in the equator plane. Each layer is
represented by a circle with a dot at the positions of grid nodes.
If a layer has no nodes in the equator plane, only a circle (not dots) is
drawn. Coloring of radius according to colors in Figure 1.
Martin Knapmeyer, 26.01.2005
Read M-File Source Code
eof.