diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90
index 641cfb42..8eab7a58 100644
--- a/src/crater/crater_populate.f90
+++ b/src/crater/crater_populate.f90
@@ -316,7 +316,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
       if (user%doregotrack) then
          call crater_superdomain(user,surf,age,age_resolution,prod,nflux,domain,finterval)
          call regolith_depth_model(user,domain,finterval,nflux,p)
-         !call regolith_subcrater_mix(user,surf,domain,nflux,finterval,p)
+         call regolith_subcrater_mix(user,surf,domain,nflux,finterval,p)
          age = age - finterval * user%interval
       end if 
 
diff --git a/src/util/util_traverse_pop_array.f90 b/src/util/util_traverse_pop_array.f90
index 15bebc5e..01dc1c32 100644
--- a/src/util/util_traverse_pop_array.f90
+++ b/src/util/util_traverse_pop_array.f90
@@ -45,7 +45,7 @@ subroutine util_traverse_pop_array(regolayer,traverse_depth,poppedarray)
  
   
     ! Initialize popped array
-    !call util_init_array(poppedarray,initstat)
+    call util_init_array(poppedarray,initstat)
  
     
     !if (initstat) then
@@ -61,7 +61,7 @@ subroutine util_traverse_pop_array(regolayer,traverse_depth,poppedarray)
     end do
 
 
-    allocate(poppedarray,source=regolayer(maxi:N))
+    !allocate(poppedarray,source=regolayer(maxi:N))
 
 
     depth = regolayer(maxi)%depth