diff --git a/src/regolith/regolith_mix.f90 b/src/regolith/regolith_mix.f90
index 8b4c1b0f..bea3958b 100644
--- a/src/regolith/regolith_mix.f90
+++ b/src/regolith/regolith_mix.f90
@@ -81,6 +81,12 @@ subroutine regolith_mix(user,surfi,mixing_depth,domain)
    newlayer%ejmf = newlayer%ejmf / newlayer%thickness
    newlayer%meltvolume = newlayer%meltvolume / newlayer%thickness
    newlayer%totvolume = newlayer%totvolume / newlayer%thickness
+
+   !test forcing melt fraction to equal melt volume / total volume
+   newlayer%totvolume = newlayer%thickness * user%pix * user%pix
+   newlayer%meltfrac = newlayer%meltvolume / newlayer%totvolume
+   newlayer%meltdist(:) = newlayer%distvol(:) / newlayer%totvolume
+   newlayer%ejmf = newlayer%ejm / newlayer%totvolume
    
    call util_push_array(surfi%regolayer, newlayer)
    !call util_destroy_list(poppedlist_top)