Updated pymangle so it works correctly with 4D arrays in anticipation of
updating the code for 4D (3-spatial + 1 dust type) arrays for the grid
physical properties.