From 5ae11ab8cc1cee44a3a23de122b30c08c029ee46 Mon Sep 17 00:00:00 2001 From: daminton Date: Mon, 5 Dec 2016 20:47:08 +0000 Subject: [PATCH] New model for collapse to achieve initial morphology --- src/crater/crater_emplace.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/crater/crater_emplace.f90 b/src/crater/crater_emplace.f90 index 8293c8df..77f54b31 100644 --- a/src/crater/crater_emplace.f90 +++ b/src/crater/crater_emplace.f90 @@ -75,6 +75,7 @@ subroutine crater_emplace(user,surf,crater,domain,ejbmass) crater%maxinc = max(crater%maxinc,inc) fradsq = crater%frad**2 deltaMtot = ejbmass + incsq = inc**2 ! This loop may not be parallelizable because of the linked list operation inside crater_form_interior do j=-inc,inc ! Do the loop in pixel space do i=-inc,inc @@ -93,7 +94,7 @@ subroutine crater_emplace(user,surf,crater,domain,ejbmass) call util_periodic(xpi,ypi,user%gridsize) lradsq = (crater%xl - xp)**2 + (crater%yl - yp)**2 - if (lradsq <= fradsq) then + if (lradsq <= crater%rad**2) then call crater_form_interior(user,surf(xpi,ypi),crater,lradsq,newelev,deltaMi) deltaMtot = deltaMtot + deltaMi end if