From 679e6732ae94e59525785fda53d7f3d5faab1e85 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Wed, 25 Jul 2012 18:09:13 -0400 Subject: [PATCH] Moved line in ACE law 44 calculation to save a few flops. --- src/physics.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/physics.F90 b/src/physics.F90 index dc2f9b815f..70a70fcab1 100644 --- a/src/physics.F90 +++ b/src/physics.F90 @@ -1821,8 +1821,8 @@ contains ! Sampled correlated angle from Kalbach-Mann parameters r3 = prn() r4 = prn() - T = (TWO*r4 - ONE) * sinh(KM_A) if (r3 > KM_R) then + T = (TWO*r4 - ONE) * sinh(KM_A) mu_out = log(T + sqrt(T*T + ONE))/KM_A else mu_out = log(r4*exp(KM_A) + (ONE - r4)*exp(-KM_A))/KM_A