From ad644b0d20d22a5b0decc3d7478dd92ac8cadf0d Mon Sep 17 00:00:00 2001
From: David Minton <daminton@purdue.edu>
Date: Fri, 23 Jul 2021 14:17:07 -0400
Subject: [PATCH] Removed obsolete diffusion code. This likely will need to be
 revisited.

---
 src/crater/crater_superdomain.f90 | 16 +---------------
 src/crater/module_crater.f90      |  3 +--
 2 files changed, 2 insertions(+), 17 deletions(-)

diff --git a/src/crater/crater_superdomain.f90 b/src/crater/crater_superdomain.f90
index 838046ab..9bb46f70 100644
--- a/src/crater/crater_superdomain.f90
+++ b/src/crater/crater_superdomain.f90
@@ -19,7 +19,7 @@
 !  Notes       :  
 !
 !**********************************************************************************************************************************
-subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,finterval,kdiff)
+subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,finterval)
    use module_globals
    use module_util
    use module_ejecta
@@ -35,7 +35,6 @@ subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,fin
    real(DP),dimension(:,:),intent(in)                  :: prod,nflux 
    type(domaintype),intent(in)                         :: domain
    real(DP),intent(in)                                 :: finterval
-   real(DP),dimension(:,:),intent(inout),optional      :: kdiff
 
    ! Internal variables
    integer(I4B),dimension(2,0:user%gridsize + 1,0:user%gridsize + 1) :: indarray
@@ -179,19 +178,6 @@ subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,fin
 
          end if
 
-         if (user%dosoftening .and. present(kdiff)) then
-            do j = 1, user%gridsize
-               do i = 1, user%gridsize
-                  xpi = i - crater%xlpx
-                  ypi = j - crater%ylpx
-                  if ((abs(xpi) > inc) .or. (abs(ypi) > inc)) cycle
-                  if (ejisray(xpi,ypi) == 0) cycle
-                  kdiff(i,j) = kdiff(i,j) + rayfrac * user%soften_factor / (PI * &
-                  user%soften_size**2) * crater%frad**(user%soften_slope - 2.0_DP)
-               end do
-            end do
-         end if
-
          deallocate(ejdistribution,ejisray)
 
        end do
diff --git a/src/crater/module_crater.f90 b/src/crater/module_crater.f90
index fa221ebe..e3038c0a 100644
--- a/src/crater/module_crater.f90
+++ b/src/crater/module_crater.f90
@@ -265,7 +265,7 @@ end function crater_critical_slope
    end interface
 
    interface
-      subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,finterval,kdiff)
+      subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,finterval)
       use module_globals
       type(usertype),intent(in)                           :: user
       type(surftype),dimension(:,:),intent(inout)         :: surf
@@ -274,7 +274,6 @@ subroutine crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,fin
       real(DP),dimension(:,:),intent(in)                  :: prod,nflux
       type(domaintype),intent(in)                         :: domain
       real(DP),intent(in)                                 :: finterval
-      real(DP),dimension(:,:),intent(inout),optional      :: kdiff
       end subroutine crater_superdomain
    end interface