Commit 084a311e authored by Gregor Michalicek's avatar Gregor Michalicek

Bugfix: The last Fermi energy wasn't read in properly since we

don't store the output density anymore. I now read it from the
last input density. ...this may be slightly inconsistent but
in the case of self convergence it should be irrelevant.
parent ff4216a5
......@@ -704,9 +704,9 @@ MODULE m_cdn_io
previousDensityIndex, jspins, date, time, distance, fermiEnergy, l_qfix)
archiveName = ''
WRITE(archiveName,'(a,i0)') '/cdn-', previousDensityIndex
l_exist = isDensityEntryPresentHDF(fileID,archiveName,DENSITY_TYPE_NOCO_OUT_const)
l_exist = isDensityEntryPresentHDF(fileID,archiveName,DENSITY_TYPE_NOCO_IN_const)
IF(l_exist) THEN
CALL peekDensityEntryHDF(fileID, archiveName, DENSITY_TYPE_NOCO_OUT_const,&
CALL peekDensityEntryHDF(fileID, archiveName, DENSITY_TYPE_NOCO_IN_const,&
iter, starsIndex, latharmsIndex, structureIndex, stepfunctionIndex,&
previousDensityIndex, jspins, date, time, distance, fermiEnergy, l_qfix)
eFermiPrev = fermiEnergy
......
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