From b94a8a62f72fdd5a526c30cfb1d350327cac020f Mon Sep 17 00:00:00 2001 From: David Minton Date: Fri, 23 Jul 2021 13:06:21 -0400 Subject: [PATCH] Checked out from main branch --- src/crater/crater_emplace.f90 | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/src/crater/crater_emplace.f90 b/src/crater/crater_emplace.f90 index 6e58bce2..79b91a44 100644 --- a/src/crater/crater_emplace.f90 +++ b/src/crater/crater_emplace.f90 @@ -64,7 +64,7 @@ subroutine crater_emplace(user,surf,crater,domain,deltaMtot) real(DP),intent(out) :: deltaMtot ! Internal variables - real(DP) :: lradsq,newelev + real(DP) :: lradsq,newelev, x_relative, y_relative integer(I4B) :: xpi,ypi,i,j,inc,incsq,iradsq real(DP) :: xp,yp,fradsq,deltaMi,rimheight logical :: lastloop @@ -94,10 +94,13 @@ subroutine crater_emplace(user,surf,crater,domain,deltaMtot) ! periodic boundary conditions call util_periodic(xpi,ypi,user%gridsize) - lradsq = (crater%xl - xp)**2 + (crater%yl - yp)**2 + x_relative = (crater%xl - xp) + y_relative = (crater%yl - yp) + + lradsq = x_relative**2 + y_relative**2 if (lradsq > crater%frad**2) cycle - call crater_form_interior(user,surf(xpi,ypi),crater,lradsq,newelev,deltaMi) + call crater_form_interior(user,surf(xpi,ypi),crater,x_relative, y_relative,newelev,deltaMi) deltaMtot = deltaMtot + deltaMi ! do porosity computation if (user%doporosity)