Fix #110: correct coefficients in MS subroutine

nrc-cnrc · Aug 17, 2016 · a93c726 · a93c726
1 parent d70e11d
commit a93c726
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/HEN_HOUSE/src/egsnrc.mortran b/HEN_HOUSE/src/egsnrc.mortran
@@ -4089,7 +4089,7 @@ tau2    = tau*tau;
 epsilon = eloss/e0;
 epsilonp= eloss/e;
 "e       = e * (1 - epsilonp*epsilonp*((6+tau*(10+5*tau))/(tau+1)/(tau+2))/24);
-e       = e * (1 - epsilonp*epsilonp*(6+10*tau+5*tau2)/(24*tau2+48*tau+72));
+e       = e * (1 - epsilonp*epsilonp*(6+10*tau+5*tau2)/(24*tau2+72*tau+48));
 p2      = e*(e + rmt2);
 "p2i     = 1/p2;
 beta2   = p2/(p2 + rmsq);