From cbcbaed8899ba3c92721394b0f9ecb5a609ecaca Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 23 Feb 2017 20:47:24 +0000 Subject: [PATCH] I had the exponent with the wrong sign in crater_soften_accumulate --- src/crater/crater_soften_accumulate.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/crater/crater_soften_accumulate.f90 b/src/crater/crater_soften_accumulate.f90 index 7b499bd0..da703e73 100644 --- a/src/crater/crater_soften_accumulate.f90 +++ b/src/crater/crater_soften_accumulate.f90 @@ -37,7 +37,7 @@ subroutine crater_soften_accumulate(user,surf,crater,domain,kdiff) real(DP) :: kappatmax,lrad,lradsq,xp,yp,fradsq,areafrac,xbar,ybar kappatmax = user%soften_factor / (PI * user%soften_size**2) * crater%frad**(user%soften_slope - 2.0_DP) - kappatmax = kappatmax - 0.84_DP / (PI * (user%soften_size * crater%frad)**2) + kappatmax = kappatmax - 0.84_DP / (PI * user%soften_size**2) * crater%frad**2 inc = int(user%soften_size * crater%frad / user%pix) + 2 crater%maxinc = max(crater%maxinc,inc) fradsq = crater%frad**2