Note: The output unit for <T> is hbar/2. ...this has to be changed in the future.
Note2: Switches have to be implemented to toggle the computation of these quantities.

...I hope I found all of them. :)

...otherwise the stepfunction is not automatically recalculated if
the MT radii change.

These files lead to too many problems for new users. If you need
this data to be stored, e.g., for performace, use the cdn.hdf file.

Minor fixes at I/O. HDF is now working with FullyFullyNoco but not functioning yet since the offdiagonal elements of the density are not stroed in hdf yet.

Started implementing the rotation of the Spin directions without rotating the spin quantization axis.

score is now only in plotdop.f90 and plot.f90=> TO BE RENAMED SOON. Next step: Make Testfiles free of score and plplot.

Removed plplot and score from inp.xml. They still need to be removed from the tests and in case of score from the code.

the previous  z(:,i) = CMPLX( z1(1,:,i) ,z1(2,:,i) ) avoids getting error message "3-rd index is 2 but upper bound
is 1" when doing array-bound-checking    

...Christoph asked for this because a too low precision breaks the respective requirements in Spex.