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.