Fix degradation calculate if energy is small

numaryu · Sep 27, 2024 · 1b21586 · 1b21586
1 parent 4e3833b
commit 1b21586
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/class_orbital.f90 b/src/class_orbital.f90
@@ -380,6 +380,8 @@ recursive subroutine calculate_degradation(self, ngen_in)
           nenergy1_max = egrid%grid_number(energy1_max)
           nenergy2_max = egrid%grid_number(energy2_max)
 
+          if (energy2_max < minval(self%energy_triplet)) cycle
+
           ! i: index of T2
           do i = nenergy2_max + 1, nenergy_max
              energy2 = egrid%val(i)