fixed a bug where volumes didn't display for melt sheets

MintonGroup · Feb 23, 2023 · e8fbdbd · e8fbdbd
1 parent 0a64f69
commit e8fbdbd
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/regolith/regolith_interior.f90 b/src/regolith/regolith_interior.f90
@@ -73,8 +73,10 @@ subroutine regolith_interior(user,surf,crater,domain,incval,nmeltsheet,vmeltshee
             newlayer%meltvolume = vmeltsheet / nmeltsheet
             newlayer%totvolume = newlayer%meltvolume
             newlayer%meltdist(:) = 0.0_SP
+            newlayer%distvol(:) = 0.0_SP
             if(domain%currentqmc) then
                 newlayer%meltdist(domain%nqmc) = newlayer%meltfrac
+                newlayer%distvol(domain%nqmc) = newlayer%meltvolume
             end if
             call util_push_array(surf(xpi,ypi)%regolayer,newlayer)
         end do