pg.meshtools.cellDataToNodeData does not take scalar fields

[florian-wagner]

Problem description

pg.meshtools.cellDataToNodeData does not take scalar numpy arrays

Steps to reproduce

Tell us how to reproduce this issue. Ideally, you could paste the code that produces the error:

import pygimli as pg
import numpy as np
grid = pg.createGrid(x=[-1.0, 0.0, 1.0, 4.0], y=(-1.0, 0.0, 1.0, 4.0))
celldata = np.ones(grid.cellCount())
nodedata = pg.meshtools.cellDataToNodeData(grid, celldata)

Expected behavior

Nodes take averages of neighboring cells.

Actual behavior

---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
<ipython-input-63-6c40652ca48f> in <module>()
      1 grid = pg.createGrid(x=[-1.0, 0.0, 1.0, 4.0], y=(-1.0, 0.0, 1.0, 4.0))
      2 celldata = np.ones(grid.cellCount())
----> 3 nodedata = pg.meshtools.cellDataToNodeData(grid, celldata)

~/miniconda3/envs/agu17/lib/python3.6/site-packages/pygimli/meshtools/mapping.py in cellDataToNodeData(mesh, data, style)
     70             return pg.cellDataToPointData(mesh, data)
     71         elif mesh.dim() == 2:
---> 72             return np.array([pg.cellDataToPointData(mesh, data[:, 0]),
     73                              pg.cellDataToPointData(mesh, data[:, 1])]).T
     74         elif mesh.dim() == 3:

IndexError: too many indices for array