diff --git a/source/xnet_screening.f90 b/source/xnet_screening.f90 index 45631817..89678438 100644 --- a/source/xnet_screening.f90 +++ b/source/xnet_screening.f90 @@ -198,7 +198,7 @@ Subroutine screening(mask_in) dfhsdt9(0) = 0.0d0 dfhsdt9(1:izmax+2) = + cds(1) * gammaz(1:izmax+2) & & + cds(2) * gammaz(1:izmax+2)**cds(5) & - & + cds(3) * log(gammaz(1:izmax+2)) + & + cds(3) dfhsdt9(1:izmax+2) = -dfhsdt9(1:izmax+2) / t9t(izb) ! No screening term for 1-reactant reactions