From e3fc5e0924a8afb835d24f155dc71127cdc404d3 Mon Sep 17 00:00:00 2001 From: daminton Date: Mon, 5 Dec 2016 20:45:39 +0000 Subject: [PATCH] New model for crater form exterior --- src/crater/crater_form_exterior.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) 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