From 1eb6726d70e11aa4959d4356c101494085bd7010 Mon Sep 17 00:00:00 2001 From: Austin Blevins Date: Fri, 16 Jun 2023 14:38:01 -0400 Subject: [PATCH] fixed a bug that caused an exception when using debug flags; age still bugged for QMC runs --- src/regolith/regolith_streamtube_head.f90 | 1 + src/regolith/regolith_streamtube_lineseg.f90 | 1 + src/regolith/regolith_subpixel_streamtube.f90 | 2 ++ src/regolith/regolith_traverse_streamtube.f90 | 1 + 4 files changed, 5 insertions(+) diff --git a/src/regolith/regolith_streamtube_head.f90 b/src/regolith/regolith_streamtube_head.f90 index 3057cde7..42984663 100644 --- a/src/regolith/regolith_streamtube_head.f90 +++ b/src/regolith/regolith_streamtube_head.f90 @@ -61,6 +61,7 @@ subroutine regolith_streamtube_head(user,surfi,deltar,totmare,tots,age_collector zend = z zmin = zstart zmax = 2.0 * deltar + ratio = 0.0_DP if (zend >= zmax) then ! Stream tube's head is inside the 1st layer. tots = tots + vsgly diff --git a/src/regolith/regolith_streamtube_lineseg.f90 b/src/regolith/regolith_streamtube_lineseg.f90 index cf72a475..05d6dda1 100644 --- a/src/regolith/regolith_streamtube_lineseg.f90 +++ b/src/regolith/regolith_streamtube_lineseg.f90 @@ -66,6 +66,7 @@ subroutine regolith_streamtube_lineseg(user,surfi,thetast,ri,rip1,zmin,zmax,erad vol = 0.0_DP vsh = 0.0_DP + ratio = 0.0_DP do N=M,2,-1 diff --git a/src/regolith/regolith_subpixel_streamtube.f90 b/src/regolith/regolith_subpixel_streamtube.f90 index 7b015b7d..f7ebe46b 100644 --- a/src/regolith/regolith_subpixel_streamtube.f90 +++ b/src/regolith/regolith_subpixel_streamtube.f90 @@ -106,6 +106,8 @@ subroutine regolith_subpixel_streamtube(user,surfi,deltar,ri,rip1,eradi,newlayer recyratio2 = 0.0_DP vsh = 0.0_DP vsh2 = 0.0_DP + ratio = 0.0_DP + ratio2 = 0.0_DP ! Two cases: subpixel is inside the first layer, and its volume is simply the landing ejecta blanket. if (zend>=zmax) then diff --git a/src/regolith/regolith_traverse_streamtube.f90 b/src/regolith/regolith_traverse_streamtube.f90 index f777aabc..abfd7096 100644 --- a/src/regolith/regolith_traverse_streamtube.f90 +++ b/src/regolith/regolith_traverse_streamtube.f90 @@ -55,6 +55,7 @@ subroutine regolith_traverse_streamtube(user,surfi,deltar,ri,rip1,eradi,erado,ne !executable code vsh = 0.0_DP + ratio = 0.0_DP erad = (eradi + erado)/2.0 rzmax = erad * sqrt(3.0)/4.0