Commit 94255d50 authored by Gregor Michalicek's avatar Gregor Michalicek

Write out vacuum energy parameters to out.xml

parent e3e500a1
......@@ -502,6 +502,14 @@ CONTAINS
ENDIF
ENDIF
evac(ivac,jsp) = enpara%evac0(ivac,jsp) + vz0
attributes = ''
WRITE(attributes(1),'(i0)') ivac
WRITE(attributes(2),'(i0)') jsp
WRITE(attributes(3),'(f16.10)') vz(1,ivac,jsp)
WRITE(attributes(4),'(f16.10)') vz(vacuum%nmz,ivac,jsp)
WRITE(attributes(5),'(f16.10)') evac(ivac,jsp)
CALL writeXMLElementForm('vacuumEP',(/'vacuum','spin ','vzIR ','vzInf ','value '/),&
attributes(1:5),reshape((/6+4,4,4,5,5+30,8,1,16,16,16/),(/5,2/)))
ENDDO
IF (vacuum%nvac.EQ.1) THEN
evac(2,jsp) = evac(1,jsp)
......
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