From d603c0a57249e94c85bb5073a9ab74e965168f19 Mon Sep 17 00:00:00 2001
From: David Bowler <david.bowler@ucl.ac.uk>
Date: Wed, 9 Sep 2026 16:00:56 +0900
Subject: [PATCH] Restore limits on band structure

---
 tools/PostProcessing/process_module.f90 | 4 ++--
 tools/PostProcessing/read_module.f90    | 2 +-
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/tools/PostProcessing/process_module.f90 b/tools/PostProcessing/process_module.f90
index 253566fc..77d26f4c 100644
--- a/tools/PostProcessing/process_module.f90
+++ b/tools/PostProcessing/process_module.f90
@@ -2001,13 +2001,13 @@ subroutine process_band_structure
        E_DOS_min = minval(eigenvalues(1,:,:))
        write(*,fmt='(2x,"Band structure lower limit set automatically: ",f12.5," Ha")') E_DOS_min
     else
-       write(*,fmt='(2x,"Band structure lower limit set by user: ",f12.5," Ha")') E_DOS_min
+       write(*,fmt='(2x,"Band structure lower limit set by user: ",e12.5," Ha")') E_DOS_min
     end if
     if(abs(E_DOS_max)<RD_ERR) then
        E_DOS_max = maxval(eigenvalues(n_bands_total,:,:))
        write(*,fmt='(2x,"Band structure upper limit set automatically: ",f12.5," Ha")') E_DOS_max
     else
-       write(*,fmt='(2x,"Band structure upper limit set by user: ",f12.5," Ha")') E_DOS_max
+       write(*,fmt='(2x,"Band structure upper limit set by user: ",e12.5," Ha")') E_DOS_max
     end if
     write(*,fmt='(2x,"Writing band structure files")')
     if(flag_proc_band_str==4) then
diff --git a/tools/PostProcessing/read_module.f90 b/tools/PostProcessing/read_module.f90
index ce5d33d7..f8038dae 100644
--- a/tools/PostProcessing/read_module.f90
+++ b/tools/PostProcessing/read_module.f90
@@ -273,7 +273,7 @@ subroutine read_input
     ! DOS
     flag_wf_range_Ef = fdf_boolean('IO.WFRangeRelative',.true.)
     flag_procwf_range_Ef = fdf_boolean('Process.WFRangeRelative',flag_wf_range_Ef)
-    if(i_job==6.or.i_job==7) then
+    if(i_job==6.or.i_job==7.or.i_job==8) then
        ! Add flag for window relative to Fermi level
        E_DOS_min = fdf_double('Process.min_DOS_E',E_wf_min)
        E_DOS_max = fdf_double('Process.max_DOS_E',E_wf_max)