- when used for example with wallShearStress, the stress field is
  initially created as incompressible but later updated with the
  correct compressible/incompressible dimensions.

  If this field is sampled as a surface and stored on the registry
  the dimensions should be reset() and not '=' assigned, since that
  causes a dimension check which will obviously fail.

- add "point data" or "face data" note to the boundaryData output
  for easier diagnositics in the future.

- the generic constructor could be called twice in pointPatchField::New,
  which results in an attempt to transfer tokens twice.

- setup writer outside the data loop to ensure that the number of
  output fields is correct (VTK format).

- ignore 'interpolate' on sampled surfaces to ensure proper
  face sampling, never allow point sampling

BUG: incorrect debug-switch for sampledIsoSurface

- number of fields was not set.
- interpolated surfaces incorrectly written for all formats

Should not call map() since it triggers building the
(one-to-one) mapping which conflicts with the self-contained
mapping of cyclicAMI

BUG: setAlphaField: fix incompatibilities with BCs (#1962)

See merge request !418

Johan Roenby authored Dec 21, 2020 and Andrew Heather committed Dec 21, 2020

Co-authored-by: Johan Roenby <johan.roenby@gmail.com>
Co-authored-by: Henning Scheufler <Henning.Scheufler@dlr.de>