From 363f84f216e14fe14562a73106a2dd6bba64fdc8 Mon Sep 17 00:00:00 2001 From: Austin Blevins Date: Fri, 30 Jun 2023 15:21:47 -0400 Subject: [PATCH] modify regolith_streamtube --- src/regolith/regolith_streamtube.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/regolith/regolith_streamtube.f90 b/src/regolith/regolith_streamtube.f90 index a2220a6f..4088d894 100644 --- a/src/regolith/regolith_streamtube.f90 +++ b/src/regolith/regolith_streamtube.f90 @@ -343,7 +343,7 @@ subroutine regolith_streamtube(user,surf,crater,domain,ejb,ejtble,xp,yp,xpi,ypi, newlayer%distvol(:) = newlayer%distvol(:) + distvol(:) newlayer%distvol(1+domain%rcnum) = newlayer%meltvolume - sum(newlayer%distvol(1:domain%rcnum)) - if (newlayer%distvol(1+domain%rcnum) < 0.0_DP) then !pixel consists entirely of QMC melt + if (newlayer%distvol(1+domain%rcnum) < 0.0) then !pixel consists entirely of QMC melt newlayer%distvol(1+domain%rcnum) = 0.0_SP newlayer%age(:) = 0.0_SP end if @@ -358,11 +358,11 @@ subroutine regolith_streamtube(user,surf,crater,domain,ejb,ejtble,xp,yp,xpi,ypi, end if !conserve volume in the age array - if (sum(newlayer%age(:)) > 0.0_DP) then + if (sum(newlayer%age(:)) > 0.0) then agefactor = newlayer%distvol(1+domain%rcnum) / sum(newlayer%age(:)) newlayer%age(:) = newlayer%age(:) * agefactor else - newlayer%age(:) = 0.0_DP + newlayer%age(:) = 0.0_SP end if call util_push_array(surf(xpi,ypi)%regolayer,newlayer)