From dca7b04dba177a5e89ce6181f605dbe9e10f60a6 Mon Sep 17 00:00:00 2001 From: daminton Date: Tue, 6 Dec 2016 15:21:44 +0000 Subject: [PATCH] Tweaked crater exterior funcdtion --- src/crater/crater_form_exterior.f90 | 6 +----- 1 file changed, 1 insertion(+), 5 deletions(-) diff --git a/src/crater/crater_form_exterior.f90 b/src/crater/crater_form_exterior.f90 index 9fb6af6d..ee4dc547 100644 --- a/src/crater/crater_form_exterior.f90 +++ b/src/crater/crater_form_exterior.f90 @@ -43,8 +43,7 @@ subroutine crater_form_exterior(user,surfi,crater,domain,lradsq,newelev,rd,delta lrad = sqrt(lradsq) ! change digital elevation map - !cform = (crater%rheight * ((crater%frad / lrad)**RIMDROP)) - cform = (crater%rheight * ((crater%frad / lrad)**rd)) + cform = (1.3_DP * crater%rheight * ((crater%rad / lrad)**rd)) ! elevation changes if (lrad <= crater%frim) then @@ -55,9 +54,6 @@ subroutine crater_form_exterior(user,surfi,crater,domain,lradsq,newelev,rd,delta end if surfi%dem = surfi%dem + elchange deltaMi = elchange - do l = 1, user%numlayers - surfi%baseline(l) = surfi%baseline(l) + elchange - end do !change ejecta coverage !surfi%ejcov = max(surfi%ejcov + elchange,0.0_DP)