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.