From fe93d4f89721a3bf031273c7c705723dedd86e06 Mon Sep 17 00:00:00 2001 From: Austin Blevins Date: Wed, 21 Jun 2023 14:38:17 -0400 Subject: [PATCH] push to the cluster for debugging why ages are negative in some cases --- src/regolith/regolith_streamtube.f90 | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/regolith/regolith_streamtube.f90 b/src/regolith/regolith_streamtube.f90 index 77e2f74c..0f4a2863 100644 --- a/src/regolith/regolith_streamtube.f90 +++ b/src/regolith/regolith_streamtube.f90 @@ -339,16 +339,20 @@ subroutine regolith_streamtube(user,surf,crater,domain,ejb,ejtble,xp,yp,xpi,ypi, factor = newlayer%totvolume / newlayer%meltvolume newlayer%meltvolume = newlayer%totvolume distvol(:) = distvol(:) * factor + newlayer%age(:) = newlayer%age(:) * factor end if 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 newlayer%distvol(1+domain%rcnum) = 0.0_SP + newlayer%age(:) = 0.0_SP end if if (sum(newlayer%distvol) > newlayer%totvolume) then factor = newlayer%totvolume / sum(newlayer%distvol) newlayer%distvol(:) = newlayer%distvol(:) * factor + newlayer%age(:) = newlayer%age(:) * factor + end if newlayer%meltvolume = sum(newlayer%distvol)