Merge pull request #197 from PrincetonUniversity/fix_outputIsurf

fix bug in writing .sp.end, when Lconstraint.\neq.3 and Nvol=1 Isurf …
PrincetonUniversity · Dec 7, 2023 · 5a487ec · 5a487ec
2 parents c253ec3 + dfeb647
commit 5a487ec
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/global.f90 b/src/global.f90
@@ -1654,7 +1654,7 @@ subroutine wrtend
   write(iunit,'(" adiabatic   = ",257es23.15)') adiabatic(1:Mvol)
   write(iunit,'(" mu          = ",257es23.15)') mu(1:Mvol)
   write(iunit,'(" Ivolume     = ",257es23.15)') Ivolume(1:Mvol)
-  write(iunit,'(" Isurf       = ",257es23.15)') Isurf(1:Mvol-1)
+  write(iunit,'(" Isurf       = ",257es23.15)') Isurf(1:Mvol-1), 0.0
   write(iunit,'(" Lconstraint = ",i9        )') Lconstraint
   write(iunit,'(" pl          = ",257i23    )') pl(0:Mvol)
   write(iunit,'(" ql          = ",257i23    )') ql(0:Mvol)