fixed a bug with melt distribution syntax

MintonGroup · Feb 23, 2023 · b6ec18e · b6ec18e
1 parent 8d9122c
commit b6ec18e
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/src/regolith/regolith_streamtube_head.f90 b/src/regolith/regolith_streamtube_head.f90
@@ -68,7 +68,7 @@ subroutine regolith_streamtube_head(user,surfi,deltar,totmare,tots,age_collector
       age_collector(:) = age_collector(:) + current(M)%age(:) * recyratio
       headmeltvol = current(M)%meltfrac * vsgly * recyratio
       meltinejecta = meltinejecta + headmeltvol
-      distvol(:) = distvol + (current(M)%meltdist(:)*vsgly*recyratio)
+      distvol(:) = distvol(:) + (current(M)%meltdist(:)*vsgly*recyratio)
       totvol = totvol + vsgly
    else ! head is not intersected with layers. 
 
@@ -83,7 +83,7 @@ subroutine regolith_streamtube_head(user,surfi,deltar,totmare,tots,age_collector
             age_collector(:) = age_collector(:) + current(N)%age(:) * recyratio
             headmeltvol = current(N)%meltfrac * tothead * recyratio
             meltinejecta = meltinejecta + headmeltvol
-            distvol(:) = distvol + (current(N)%meltdist(:)*tothead*recyratio)
+            distvol(:) = distvol(:) + (current(N)%meltdist(:)*tothead*recyratio)
             totvol = totvol + tothead
             !current => current%next
             !N = N - 1
@@ -97,7 +97,7 @@ subroutine regolith_streamtube_head(user,surfi,deltar,totmare,tots,age_collector
             age_collector(:) = age_collector(:) + current(N)%age(:) * recyratio
             headmeltvol = current(N)%meltfrac * tothead*recyratio
             meltinejecta = meltinejecta + headmeltvol
-            distvol(:) = distvol + (current(N)%meltdist(:)*tothead*recyratio)
+            distvol(:) = distvol(:) + (current(N)%meltdist(:)*tothead*recyratio)
             totvol = totvol + tothead
             exit
          end if

diff --git a/src/regolith/regolith_subpixel_streamtube.f90 b/src/regolith/regolith_subpixel_streamtube.f90
@@ -114,7 +114,7 @@ subroutine regolith_subpixel_streamtube(user,surfi,deltar,ri,rip1,eradi,newlayer
          vsh              = regolith_shock_damage(eradi,deltar,xmints,xsfints,0.0_DP,eradi)
          recyratio        = max(vseg-vsh,0.0 )/ (user%pix**2) / (surfi%regolayer(M)%thickness)
          meltinejecta     = surfi%regolayer(M)%meltfrac * vseg * recyratio
-         distvol          = distvol + (surfi%regolayer(M)%meltdist(:)*vseg*recyratio)
+         distvol(:)       = distvol(:) + (surfi%regolayer(M)%meltdist(:)*vseg*recyratio)
          totvol           = vseg
          age_collector(:) = age_collector(:) + surfi%regolayer(M)%age(:) * recyratio
          vol              = vol + sum(age_collector(:))
@@ -205,7 +205,7 @@ subroutine regolith_subpixel_streamtube(user,surfi,deltar,ri,rip1,eradi,newlayer
 
          !current => current%next
          meltinejecta = meltinejecta + mvl + mvr
-         distvol = distvol + (current(N)%meltdist(:)*((vsgly1*recyratio)+(vsgly2*recyratio2)))
+         distvol(:) = distvol(:) + (current(N)%meltdist(:)*((vsgly1*recyratio)+(vsgly2*recyratio2)))
          !distvol = distvol + (current(N)%meltdist(:)*(vsgly1+vsgly2))
          totvol = totvol + vsgly1 + vsgly2
          !N = N - 1