New model for crater form exterior

MintonGroup · Dec 5, 2016 · e3fc5e0 · e3fc5e0
1 parent d9c94f7
commit e3fc5e0
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/crater/crater_form_exterior.f90 b/src/crater/crater_form_exterior.f90
@@ -44,11 +44,11 @@ subroutine crater_form_exterior(user,surfi,crater,domain,lradsq,newelev,rd,delta
 
    ! change digital elevation map
    !cform = (crater%rheight * ((crater%frad / lrad)**RIMDROP))
-   cform = (crater%rheight * ((crater%frad / lrad)**rd))
+   cform = (1.30*crater%rheight * ((crater%rad / lrad)**rd))
 
    ! elevation changes
    if (lrad <= crater%frim) then
-      hcorr = belchng * ((crater%frim - lrad) / (crater%frim - crater%frad))**20
+      hcorr = belchng * ((crater%frim - lrad) / (crater%frim - crater%rad))**20
       elchange = hcorr + cform 
    else 
       elchange = cform