From 75e635b7a43d3dfcf39572c57de9c863ca8cdcf6 Mon Sep 17 00:00:00 2001 From: daminton Date: Wed, 7 Dec 2016 21:43:20 +0000 Subject: [PATCH] Changed subpixel diffusion to be in terms of radius like what is in the paper --- src/crater/crater_subpixel_diffusion.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/crater/crater_subpixel_diffusion.f90 b/src/crater/crater_subpixel_diffusion.f90 index 23986208..f4a20868 100644 --- a/src/crater/crater_subpixel_diffusion.f90 +++ b/src/crater/crater_subpixel_diffusion.f90 @@ -62,8 +62,8 @@ subroutine crater_subpixel_diffusion(user,surf,prod,nflux,domain,finterval) lambda_regolith(i) = dN(i) * 0.25_DP * PI * (nflux(2,i) + SQRT2 * user%pix)**2 if (user%dosoftening) then - kappat_bedrock(i) = user%soften_factor * nflux(1,i)**(user%soften_slope) - kappat_regolith(i) = user%soften_factor * nflux(2,i)**(user%soften_slope) + kappat_bedrock(i) = user%soften_factor * (0.5_DP * nflux(1,i))**(user%soften_slope) + kappat_regolith(i) = user%soften_factor * (0.5_DP * nflux(2,i))**(user%soften_slope) end if end do