Program: MKLSGRIDDEMO

  mklsgriddemo.......demo program for grid handling
 
  call: mklsgriddemo('iteration',miniter,maxiter);
        mklsgriddemo('iteration',gridbase,maxiter);
        mklsgriddemo('resolution',gridbase,gridwdt,minr,maxr,layerdr);
        mklsgriddemo('continents',level,maxiter);
        mklsgriddemo('coasts',level,maxiter);
 
     The first parameter defines which demo to show:
 
     'iteration': shows the behaviour of tetrahedron, octahedron and
                   icsoahedron upon iteration. The user can define which
                   iteration levels to show, these are then plotted in a
                   matrix-like display with one columnt per body and lines
                   corresponding to iteration level.
 
                   additional parameters:
 
                   miniter: (number)
                            minimum iteration level
                   maxiter: (number)
                            maximum iteration level
 
     'iteration2': shows the n-th iteration for one of the base bodies.
 
                   additional parameters:
 
                   gridbase: (string)
                         name of the base body to be used. This may be one
                         of the names allowed for the grid_base field in
                         the LSOPT structure.
 
                   maxiter: (number)
                              iteration level. The base body will be
                              iterated to this level.
 
     'resolution': shows the spatial and vertical resolution of the initial
                   grid, depending on layer spacing and mesh width. layer
                   spacing can be set directly. Mesh width is controlled
                   via tha max size of riangle circumcircles, the grid is
                   iterated so fasr that all circumcircles fall below a
                   given threshold.
                   Draws 3D cutoff views of multilayer grids, statistics of
                   circumcircle distibutions and plane views of the
                   equatorial nodes of all layers.
 
                   additional parameters:
 
                   gridbase: (string)
                             base body to use for grid construction
                   gridwdt:  (number) [km]
                             desired lateral grid resolution
                      minr:  (number) [km]
                             radius of innermost layer
                      maxr:  (number) [km]
                             Radius of outermost layer
                   layerdr:  (number) [km]
                             layer spacing
 
     'continents': a single grid layer is refined locally in order to
                   produce a grid that is fine on land and coarse on
                   oceans.
                   Draws a node map and a nice 3D view.
 
                   additional parameters:
 
                   level: (number) [m]
                          elevation level in meters that is used as boundary
                   maxiter: (number)
                            max number of iteraions to be applied
 
 
     'coasts':     similar to 'continents', but this time the grid is
                   refined at coasts only.
 
                   additional parameters:
 
                   As in 'continents' demo
 
 
  result: none
 
  Martin Knapmeyer, 29.09.2004, 06.05.2008

Read M-File Source Code


eof.