From 3cd00fb5604059bcc360a7f9e596d4946e80e906 Mon Sep 17 00:00:00 2001 From: daminton Date: Wed, 14 Dec 2016 04:44:45 +0000 Subject: [PATCH] Fixed bug in crater_form_exterior function. Variable belchange was used but not defined. Switched it to cform. --- src/crater/crater_form_exterior.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/crater/crater_form_exterior.f90 b/src/crater/crater_form_exterior.f90 index ee4dc547..83d196a7 100644 --- a/src/crater/crater_form_exterior.f90 +++ b/src/crater/crater_form_exterior.f90 @@ -43,16 +43,16 @@ subroutine crater_form_exterior(user,surfi,crater,domain,lradsq,newelev,rd,delta lrad = sqrt(lradsq) ! change digital elevation map - cform = (1.3_DP * crater%rheight * ((crater%rad / lrad)**rd)) + cform = (1.3_DP * crater%rheight * ((crater%rad / lrad)**rd)) ! elevation changes if (lrad <= crater%frim) then - hcorr = belchng * ((crater%frim - lrad) / (crater%frim - crater%rad))**20 + hcorr = cform * ((crater%frim - lrad) / (crater%frim - crater%rad))**20 elchange = hcorr + cform else elchange = cform end if - surfi%dem = surfi%dem + elchange + surfi%dem = surfi%dem + elchange deltaMi = elchange !change ejecta coverage