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.