Properly handle 2D fields
If any one dimension has only length 1, then ignore defining and declaring that dimension.
This helps in viewing the 2D case because ParaView then assumes that the fields are 2D and not 3D which they currently appear to be.
This isn’t too difficult (I’ve tested it when Ny=1) and just requires many if blocks around each definition of a dimension.