WIP: Add eshift to restart file to preserve continuity

PHOTOX · Dec 1, 2023 · 1926f2e · 1926f2e
1 parent 51ccb46
commit 1926f2e
Showing 1 changed file with 2 additions and 0 deletions.
diff --git 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