fix bug in writing .sp.end, when Lconstraint.\neq.3 and Nvol=1 Isurf …

…was printing nothing and could not restart run
PrincetonUniversity · Nov 30, 2023 · dfeb647 · dfeb647
1 parent c253ec3
commit dfeb647
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)