diff --git a/src/crater/crater_form_exterior.f90 b/src/crater/crater_form_exterior.f90 index d36a7845..7c5de14d 100644 --- 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