Commit 294af201 authored by Daniel Wortmann's avatar Daniel Wortmann

Bugfixes for forcetheorem part

parent 1a19b942
......@@ -92,7 +92,7 @@ CONTAINS
!Now output the results
call closeXMLElement('Forcetheorem_Loop_DMI')
CALL openXMLElementPoly('Forcetheorem_DMI',(/'qPoints','Angles '/),(/SIZE(this%evsum,2),SIZE(this%evsum,2)/))
CALL openXMLElementPoly('Forcetheorem_DMI',(/'qPoints','Angles '/),(/SIZE(this%evsum,2),SIZE(this%evsum,1)/))
DO q=1,SIZE(this%evsum,2)
WRITE(attributes(1),'(i5)') q
WRITE(attributes(2),'(f12.7)') this%evsum(0,q)
......@@ -107,6 +107,7 @@ CONTAINS
END DO
ENDDO
CALL closeXMLElement('Forcetheorem_DMI')
CALL judft_end("Forcetheorem DMI")
END SUBROUTINE dmi_postprocess
SUBROUTINE dmi_dist(this,mpi)
......
......@@ -66,7 +66,7 @@ CONTAINS
TYPE IS (t_forcetheo_jij)
ALLOCATE(q_vectors(3,SIZE(forcetheo%qvec,2)))
q_vectors=forcetheo%qvec
TYPE IS (t_forcetheo) ! DEFAULT
CLASS IS (t_forcetheo) ! DEFAULT
ALLOCATE(q_vectors(3,1))
q_vectors(:,1)=noco%qss
END SELECT
......
......@@ -2154,7 +2154,7 @@ FUNCTION countStringTokens(line) RESULT(tokenCount)
INTEGER:: n,i
CHARACTER(len=256):: xpatha,valueString
n=evaluateFirstIntOnly(xmlGetAttributeValue(TRIM(ADJUSTL(path))//'/@count'))
n=xmlGetNumberOfNodes(TRIM(ADJUSTL(path))//'/q')
ALLOCATE(q(3,n))
DO i = 1, n
PRINT *, path,'/q[',i,']'
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment