From a623e452ba0ab68f5ac4a31a9d9b2a860c959929 Mon Sep 17 00:00:00 2001 From: daminton Date: Tue, 6 Dec 2016 14:32:20 +0000 Subject: [PATCH] Updated obsolete crater_soften routine --- src/crater/crater_soften.f90 | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/src/crater/crater_soften.f90 b/src/crater/crater_soften.f90 index e192d96c..93bd4500 100644 --- a/src/crater/crater_soften.f90 +++ b/src/crater/crater_soften.f90 @@ -36,9 +36,9 @@ subroutine crater_soften(user,surf,crater,domain) integer(I4B),parameter :: MAXHITS = 1 integer(I4B) :: inc,incsq,N,xpi,ypi,iradsq,i,j - real(DP) :: xp,yp,fradsq,xbar,ybar,areafrac,kappatmax + real(DP) :: xp,yp,fradsq,xbar,ybar,areafrac,kappatmax,ss,sf - kappatmax = SOFTEN_FACTOR * crater%fcrat**SOFTEN_SLOPE + kappatmax = user%soften_factor * crater%fcrat**user%soften_slope inc = max(min(int(crater%fradpx + 2),PBCLIM * user%gridsize),2) crater%maxinc = max(crater%maxinc,inc) @@ -81,7 +81,6 @@ subroutine crater_soften(user,surf,crater,domain) end if end do end do !end area loopover - !read(*,*) call util_diffusion_solver(user,surf,2 * inc + 1,indarray,kappat,cumulative_elchange,maxhits)