http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg11158.html
--- WIEN2k_14.orig/SRC_lapw5/main1.f    2014-09-19 13:30:20.757390463 +0200
+++ WIEN2k_14/SRC_lapw5/main1.f 2014-09-19 13:30:37.671774186 +0200
@@ -346,7 +346,7 @@
       if(addsub.eq.' ADD') fadd=1.d0
       if(addsub.eq.' SUB') fadd=-1.d0
                                   
-      IF(IUNIT.EQ.'ANG ') FAK=Bohr**3
+      IF(IUNIT.EQ.'ANG ') FAK=1.D0/Bohr**3      ! Change to avoid altering lapw5 output (fix to WIEN2k 14.1): Was 6.7551d0 in WIEN2k 13.1 (and likely in verisons < 13.1) but changed to increase accuracy
       R=VNORM(VX)                                                       
       RELX=R                                                            
       R=VNORM(VY)                                                       
@@ -373,7 +373,7 @@
 !
 !     Correct TCharge for units
       if(UseSTM)then
-        if(IUNIT.EQ.'ANG ')CTarget=CTarget/FAK
+        if(IUNIT.EQ.'ANG ')CTarget=CTarget*FAK  ! Change to avoid altering output
 !       Array for STM (LDM)
         allocate(STM(NPX,NPY))
       endif

