Make sure we dont get very tiny numbers

rjfarmer · Jun 29, 2023 · 0780e02 · 0780e02
1 parent 0d7bbd7
commit 0780e02
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/lib_bbq.f90 b/src/lib_bbq.f90
@@ -552,6 +552,10 @@ subroutine do_clean1(species, xa, ierr)
          return
       end if
       xa(1: species) = xa(1: species)/xsum
+      do j = 1, species
+         if (xa(j) < tiny_x) xa(j) = tiny_x
+      end do    
+
    end subroutine do_clean1