forced mix layers volume consistency

MintonGroup · Mar 1, 2023 · 776f477 · 776f477
1 parent feea6e6
commit 776f477
Showing 1 changed file with 6 additions and 0 deletions.
diff --git 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)