From 396114bc7dfc451eee9ab49e21de19036d6d91cc Mon Sep 17 00:00:00 2001 From: daminton Date: Wed, 22 Feb 2017 13:04:07 +0000 Subject: [PATCH] Fixed bug in crater_soften_accumulate --- src/crater/crater_soften_accumulate.f90 | 6 ++---- 1 file changed, 2 insertions(+), 4 deletions(-) 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