I had the exponent with the wrong sign in crater_soften_accumulate

MintonGroup · Feb 23, 2017 · cbcbaed · cbcbaed
1 parent a5cfac8
commit cbcbaed
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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