From cacb4f6a5352a8647715b6dcd79bf18af480d2ba Mon Sep 17 00:00:00 2001 From: Austin Blevins Date: Thu, 23 Feb 2023 13:15:08 -0500 Subject: [PATCH] added volume to regolith mix subroutine --- src/regolith/regolith_mix.f90 | 2 ++ 1 file changed, 2 insertions(+) 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)