Removed obsolete diffusion code. This likely will need to be revisited.

MintonGroup · Jul 23, 2021 · ad644b0 · ad644b0
1 parent de2c6ec
commit ad644b0
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 17 deletions.
diff --git 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
@@ -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