Problems with smoke3d to_global() about fdsreader

[Sokdenumeros]

For the second problem, I just noticed that if I add back the loops as follows:
for axis in range(3):
dim = ('x', 'y', 'z')[axis]
# Add border points back again if needed
if np.isclose(subsmoke.mesh.coordinates[dim][-1], global_max[dim]):
temp_data_slices = [slice(s) for s in subsmoke_data.shape]
temp_data_slices[axis + 1] = slice(None)
subsmoke_data = np.concatenate((subsmoke_data, temp_data[dim][tuple(temp_data_slices)]), axis=axis + 1)

It can happen that if more than one dimension is close to the global_max, for the next axis it tries to concatenate arrays it will break. Because the shape of subsmoke_data will have been modified after the first concatenation.