From 1790412071ff939e149426d7a20135a7b5db93ea Mon Sep 17 00:00:00 2001 From: David Minton Date: Thu, 26 Sep 2019 16:58:21 -0400 Subject: [PATCH] fixed bug where degradation state was in the wrong units --- src/crater/crater_get_degradation_state.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/crater/crater_get_degradation_state.f90 b/src/crater/crater_get_degradation_state.f90 index e94f293e..84fafae3 100644 --- a/src/crater/crater_get_degradation_state.f90 +++ b/src/crater/crater_get_degradation_state.f90 @@ -82,7 +82,7 @@ function crater_get_degradation_state(user,surf,crater,dd) result(Kval) bowl = bowl / nbowl outer = outer / nouter dd = max((rim - bowl) / crater%fcrat,SMALLFAC) - Kval = - 1._DP / Bconst * log(real(dd,kind=DP) / ddinit) + Kval = - (1._DP / Bconst * log(real(dd,kind=DP) / ddinit)) * crater%frad**2 return end function crater_get_degradation_state