From 431fcaaebcc080440935904c35e48f29ada97c62 Mon Sep 17 00:00:00 2001
From: David Minton <daminton@purdue.edu>
Date: Fri, 26 Jul 2019 23:10:38 -0400
Subject: [PATCH] Updated the complex crater d/D cutoff based on simple
 diffusion calculation calibrated for Kv

---
 src/crater/crater_tally_observed.f90 | 13 ++++++++-----
 1 file changed, 8 insertions(+), 5 deletions(-)

diff --git a/src/crater/crater_tally_observed.f90 b/src/crater/crater_tally_observed.f90
index 97acd465..f6b5a90d 100644
--- a/src/crater/crater_tally_observed.f90
+++ b/src/crater/crater_tally_observed.f90
@@ -200,11 +200,14 @@ subroutine crater_tally_observed(user,surf,domain,nkilled,onum,obsdist,obslist,o
       bowl = bowl / nbowl 
       outer = outer / nouter
       tmp_depthdiam(craternum) = (rim - bowl) / crater%fcrat
-      if (crater%fcrat < 20e3_DP) then
-         dd = DDCUTOFF
-      else
-         dd = DDCUTOFF - (crater%fcrat - 20e3_DP) * 2e-7
-      end if
+
+      !if (crater%fcrat < 20e3_DP) then
+      !   dd = DDCUTOFF
+      !else
+      !   dd = DDCUTOFF - (crater%fcrat - 20e3_DP) * 2e-7
+      !end if
+      dd = min(DDCUTOFF, 24.897 * crater%fcrat ** (-0.632545))
+
       if (((tmp_depthdiam(craternum) > dd).and.((outer - rim) / crater%fcrat < OCUTOFF)).and.&
           (nrim /= 0).and.(nbowl /= 0).and.(nouter /= 0)) then
          countable(craternum) = .true.