From 776f4770517f096e7f1cc5a4d545f47c9f8f30a5 Mon Sep 17 00:00:00 2001 From: Austin Blevins Date: Wed, 1 Mar 2023 16:05:37 -0500 Subject: [PATCH] forced mix layers volume consistency --- src/regolith/regolith_mix.f90 | 6 ++++++ 1 file changed, 6 insertions(+) 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)