Put two more digits in output of eigenvalue sum in MAE calculations

forcetheorem/mae.F90

@@ -136,7 +136,7 @@ CONTAINS
 
     !Locals
     INTEGER:: n
-    CHARACTER(LEN=12):: attributes(3)
+    CHARACTER(LEN=16):: attributes(3)
     IF (this%directions_done==0) THEN
        RETURN
     ENDIF
@@ -151,9 +151,9 @@ CONTAINS
        DO n=1,SIZE(this%evsum)
           WRITE(attributes(1),'(f12.7)') this%theta(n)
           WRITE(attributes(2),'(f12.7)') this%phi(n)
-          WRITE(attributes(3),'(f12.7)') this%evsum(n)
+          WRITE(attributes(3),'(f16.9)') this%evsum(n)
           CALL writeXMLElementForm('Angle',(/'theta ','phi   ','ev-sum'/),attributes,&
-               RESHAPE((/5,3,6,12,12,12/),(/3,2/)))
+               RESHAPE((/5,3,6,12,12,16/),(/3,2/)))
        END DO
        CALL closeXMLElement('Forcetheorem_MAE')
     ENDIF