From cb23c74d9d6d438cd435ec560b8bc10e10b3055b Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 23 Feb 2017 12:53:35 +0000 Subject: [PATCH] Fixed unit problem in crater_soften_accumulate soften equation --- src/crater/crater_soften_accumulate.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/crater/crater_soften_accumulate.f90 b/src/crater/crater_soften_accumulate.f90 index 39b4d035..7314ef77 100644 --- a/src/crater/crater_soften_accumulate.f90 +++ b/src/crater/crater_soften_accumulate.f90 @@ -45,8 +45,8 @@ subroutine crater_soften_accumulate(user,surf,crater,domain,kdiff) !******** - kappatmax = user%soften_factor / (PI * (SOFTEN_SIZE * crater%frad)**(user%soften_slope - 2.0_DP)) - kappatmax = kappatmax - 0.84_DP / (PI * (SOFTEN_SIZE * crater%frad)**2.0) + kappatmax = user%soften_factor / (PI * SOFTEN_SIZE**2 * crater%frad**(user%soften_slope - 2.0_DP)) + kappatmax = kappatmax - 0.84_DP / (PI * (SOFTEN_SIZE * crater%frad)**2) inc = SOFTEN_SIZE * crater%fradpx + 2 crater%maxinc = max(crater%maxinc,inc) fradsq = crater%frad**2