fixed bugs causing nan's

MintonGroup · Mar 2, 2023 · 6766f77 · 6766f77
1 parent 9d3c1db
commit 6766f77
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 0 deletions.
diff --git a/src/regolith/regolith_streamtube_lineseg.f90 b/src/regolith/regolith_streamtube_lineseg.f90
@@ -65,6 +65,7 @@ subroutine regolith_streamtube_lineseg(user,surfi,thetast,ri,rip1,zmin,zmax,erad
    end if
 
    vol = 0.0_DP
+   vsh = 0.0_DP
 
    do N=M,2,-1
 

diff --git a/src/regolith/regolith_subpixel_streamtube.f90 b/src/regolith/regolith_subpixel_streamtube.f90
@@ -104,6 +104,8 @@ subroutine regolith_subpixel_streamtube(user,surfi,deltar,ri,rip1,eradi,newlayer
    mvl = 0.0_DP
    mvr = 0.0_DP
    recyratio2 = 0.0_DP
+   vsh = 0.0_DP
+   vsh2 = 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
@@ -54,6 +54,8 @@ subroutine regolith_traverse_streamtube(user,surfi,deltar,ri,rip1,eradi,erado,ne
 
    !executable code
 
+   vsh = 0.0_DP
+
    erad = (eradi + erado)/2.0
    rzmax = erad * sqrt(3.0)/4.0
    if (ri .eq. 0) then
-Original file line number
+Diff line change
@@ Expand Up @@
        end if
        vol = 0.0_DP
+       vsh = 0.0_DP
        do N=M,2,-1
@@ Expand Down @@