From 0780e0257b27d242c05174f4541b424dbe4376b7 Mon Sep 17 00:00:00 2001 From: Rob Farmer Date: Thu, 29 Jun 2023 17:07:10 +0200 Subject: [PATCH] Make sure we dont get very tiny numbers --- src/lib_bbq.f90 | 4 ++++ 1 file changed, 4 insertions(+) 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