From dfeb64775b2de4b5ca371f4e77613cdfcdd26794 Mon Sep 17 00:00:00 2001 From: Joaquim Loizu Date: Thu, 30 Nov 2023 17:48:30 +0100 Subject: [PATCH] fix bug in writing .sp.end, when Lconstraint.\neq.3 and Nvol=1 Isurf was printing nothing and could not restart run --- src/global.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/global.f90 b/src/global.f90 index ef5056fe..c8529c78 100644 --- 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)