small adjustment to the function that computes the renormalization sc…

…ale as a function of alpha_s
scarrazza · Feb 20, 2017 · 4b8630d · 4b8630d
1 parent 3702934
commit 4b8630d
Showing 1 changed file with 7 additions and 3 deletions.
diff --git a/src/Evolution/a_QCD.f b/src/Evolution/a_QCD.f
@@ -783,9 +783,13 @@ function muR2(as)
             snf = 0
          endif
 *
-         integral = integral + intBeta(nfi,ipt,as0,asthDown(nfi+snf))
-*
-         as0 = asthUp(nfi+snf)
+         if(nff.gt.nfi)then
+            integral = integral + intBeta(nfi,ipt,as0,asthDown(nfi+snf))
+            as0 = asthUp(nfi+snf)
+         else
+            integral = integral + intBeta(nfi,ipt,as0,asthUp(nfi))
+            as0 = asthDown(nfi)
+         endif
          nfi  = nfi + dnf
          goto 10
       endif