1) DFPT+U with ortho-atomic orbitals is disables,
because it is more complicated than it was expected.
First, Hubbard forces in pw.x must be extended to ortho-atomic orbitals
and then DFPT+U must be reconsidered.
2) Update the reference files for example18 for DFPT+U
3) Various minor changes