diff --git a/src/regolith/regolith_mix.f90 b/src/regolith/regolith_mix.f90
index 7ea0a444..06b7da90 100644
--- a/src/regolith/regolith_mix.f90
+++ b/src/regolith/regolith_mix.f90
@@ -65,6 +65,7 @@ subroutine regolith_mix(surfi,mixing_depth,domain)
       newlayer%distvol(:) = newlayer%distvol(:) + poppedarray(i)%thickness * poppedarray(i)%distvol(:)
       newlayer%ejm       = newlayer%ejm + poppedarray(i)%thickness * poppedarray(i)%ejm
       newlayer%ejmf      = newlayer%ejmf + poppedarray(i)%thickness * poppedarray(i)%ejmf
+      newlayer%meltvolume = newlayer%meltvolume + poppedarray(i)%thickness * poppedarray(i)%meltvolume
    end do
 
    ! Get average values of composition and melt fraction
@@ -74,6 +75,7 @@ subroutine regolith_mix(surfi,mixing_depth,domain)
    newlayer%distvol(:) = newlayer%distvol(:) / newlayer%thickness
    newlayer%ejm = newlayer%ejm / newlayer%thickness
    newlayer%ejmf = newlayer%ejmf / newlayer%thickness
+   newlayer%meltvolume = newlayer%meltvolume / newlayer%thickness
    
    call util_push_array(surfi%regolayer, newlayer)
    !call util_destroy_list(poppedlist_top)