From f08287f008b6593ca3c6003a418b75639b47704f Mon Sep 17 00:00:00 2001 From: Georgia Stuart Date: Thu, 11 Feb 2021 10:34:08 -0600 Subject: [PATCH] Fix denominator in DiameterCritical for larger particles --- lpt_oil.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lpt_oil.F90 b/lpt_oil.F90 index ab7072e..5a23593 100755 --- a/lpt_oil.F90 +++ b/lpt_oil.F90 @@ -118,7 +118,7 @@ SUBROUTINE LPT_Oil_Buoyancy(DiameterEffective,VelocityTerminal) B2 = Y1 - A2*X1 ! Compute the critical diameter. - DiameterCritical = 10.D0**((B2-B1)/(A1-A1)) + DiameterCritical = 10.D0**((B2-B1)/(A1-A2)) ! Ellipsoidal shape, intermediate size range. IF(DiameterEffective.LE.DiameterCritical)THEN