I'm trying to calculate the electron density distribution in a TiFe intermetallic. Using plot2d to plot this on a particular crystallographic plane works well, but plot3d does not. The output file RHO3D.xml contains the correct data point structure but the value at each point is exactly identical.
This seems to be a problem whenever plot3d is used: I've found the same thing on trying to plot the electron density gradient, electric field, etc. Can anyone shed any light on this or suggest a remedy?
I have noticed a couple of what I assume to be typos in the plot3d.f90 source file: the third axis is misnamed 'b' and the attribute 'larexunit' is set. Fixing those and recompiling does not solve the above problem. Perhaps there's something else that's beyond my fortran skills to debug?
Thanks in advance.