From b9f8b0a2abfae04d0f43858f15cd8e9e9a5d7839 Mon Sep 17 00:00:00 2001 From: mhirabay Date: Thu, 29 Dec 2016 14:09:07 +0000 Subject: [PATCH] Moved the porosity function to crater_emplace. --- src/crater/crater_emplace.f90 | 12 ++++++++++++ 1 file changed, 12 insertions(+) diff --git a/src/crater/crater_emplace.f90 b/src/crater/crater_emplace.f90 index 77f54b31..4b7e7660 100644 --- a/src/crater/crater_emplace.f90 +++ b/src/crater/crater_emplace.f90 @@ -52,6 +52,7 @@ subroutine crater_emplace(user,surf,crater,domain,ejbmass) use module_globals use module_util + use module_porosity use module_crater, EXCEPT_THIS_ONE => crater_emplace implicit none @@ -98,6 +99,17 @@ subroutine crater_emplace(user,surf,crater,domain,ejbmass) call crater_form_interior(user,surf(xpi,ypi),crater,lradsq,newelev,deltaMi) deltaMtot = deltaMtot + deltaMi end if + + ! do porosity computation if (user%doporosity) + ! It is still important to consider the physical meaning of frad and rad. + ! frad is the final crater, while rad is the transient crater. + ! which one should be reasonable here. + if (lradsq <= crater%frad**2) then + if (user%doporosity) then + call porosity_form_interior(user, surf(xpi,ypi), crater, lradsq) + end if + end if + end do end do