Changed soften factor to be in terms of radius instead of diameter

MintonGroup · Dec 15, 2016 · 8c1fc33 · 8c1fc33
1 parent 3dedcc4
commit 8c1fc33
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/crater/crater_soften.f90 b/src/crater/crater_soften.f90
@@ -38,7 +38,7 @@ subroutine crater_soften(user,surf,crater,domain)
    integer(I4B) :: inc,incsq,N,xpi,ypi,iradsq,i,j
    real(DP) :: xp,yp,fradsq,xbar,ybar,areafrac,kappatmax,ss,sf
 
-   kappatmax = user%soften_factor * crater%fcrat**user%soften_slope
+   kappatmax = user%soften_factor * crater%frad**user%soften_slope
 
    inc = max(min(int(crater%fradpx + 2),PBCLIM * user%gridsize),2) 
    crater%maxinc = max(crater%maxinc,inc)