GitLab

Not writing out the explicit version of the file could
lead to problems, e.g., dependence on a non-existing
sym.out file.

The allocation was required for he call to r_inpXML

The problem only occured because stars%gamx is changed after it
is initialized. I therefore introduced an additional variable
gmaxInit in type t_stars.

Note: There are still bugs in this part of the code.

The atom positions were wrong because the array zat in
force/geo.f90 was not initialized.

The here fixed bugs are due to a missing division by 2*Pi wherever
the old cotra1 subroutine was replaced by a call to MATMUL.

The old code is not yet commented out, the result is just overwritten.

1. The total forces are now written out to out.xml iff l_f=T.
2. There were bugs in some array dimensionings in force_a21. Corrected.
3. Some cleanup in cdnval, force_a21, and force_a21_lo.

...as suggested by Daniel.

The tlmplm arrays in that routine only have a spin dimension of 1.
Depending on the variable isp these arrays either contain the values
for spin 1 or 2.

This commit also provides helper routines to write out data into out.xml.
Nothing is tested yet. The routines might be horribly wrong.

Based on v0.26 hybrid branch