diff --git a/src/lib_bbq.f90 b/src/lib_bbq.f90 index 9025494..428ad96 100644 --- 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