Fix error in strong-screening derivative w.r.t temperature

starkiller-astro · Jul 2, 2024 · 94020c1 · 94020c1
1 parent a3f6789
commit 94020c1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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