Fixed bug in crater_soften_accumulate

MintonGroup · Feb 22, 2017 · 396114b · 396114b
1 parent 0590ac5
commit 396114b
Showing 1 changed file with 2 additions and 4 deletions.
diff --git a/src/crater/crater_soften_accumulate.f90 b/src/crater/crater_soften_accumulate.f90
@@ -75,10 +75,8 @@ subroutine crater_soften_accumulate(user,surf,crater,domain,kdiff)
             call util_periodic(xpi,ypi,user%gridsize)
             areafrac = util_area_intersection(SOFTEN_SIZE * crater%frad,xbar,ybar,user%pix)
 
-            if (lradsq < fradsq) then 
-               kdiff(xpi,ypi) = 0.0_DP
-            else 
-               kdiff(xpi,ypi) = kappatmax * areafrac 
+            if (lradsq > fradsq) then 
+               kdiff(xpi,ypi) = kdiff(xpi,ypi) + kappatmax * areafrac 
             end if
 
          end if