Fixed bug in crater_form_exterior function. Variable belchange was us…

…ed but not defined. Switched it to cform.
MintonGroup · Dec 14, 2016 · 3cd00fb · 3cd00fb
1 parent 86a0178
commit 3cd00fb
Showing 1 changed file with 3 additions and 3 deletions.
diff --git 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