diff --git a/src/crater/crater_soften_accumulate.f90 b/src/crater/crater_soften_accumulate.f90 index 4b116249..f50e2cf1 100644 --- 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