From 1926f2e84fe9aba8c6bcc70a74644b191565504f Mon Sep 17 00:00:00 2001 From: Daniel Hollas Date: Wed, 9 Mar 2022 14:50:28 +0000 Subject: [PATCH] WIP: Add eshift to restart file to preserve continuity --- src/sh_integ.F90 | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/sh_integ.F90 b/src/sh_integ.F90 index 504d1448..b3c2e8f7 100644 --- a/src/sh_integ.F90 +++ b/src/sh_integ.F90 @@ -453,6 +453,7 @@ end subroutine sh_set_energy_shift subroutine sh_write_wf(outunit) integer, intent(in) :: outunit integer :: ist1 + write (outunit, *) eshift do ist1 = 1, nstate write (outunit, *) cel_re(ist1), cel_im(ist1) end do @@ -466,6 +467,7 @@ subroutine sh_read_wf(inunit) allocate (cel_re(nstate), cel_im(nstate)) end if + read (inunit, *) eshift do ist1 = 1, nstate read (inunit, *) cel_re(ist1), cel_im(ist1) end do