From 20d2ac10f21759e49607f1e72715ba2e65244236 Mon Sep 17 00:00:00 2001
From: Josh Anibal <joanibal@umich.edu>
Date: Wed, 6 Jul 2022 18:29:45 -0400
Subject: [PATCH] Quick bug fix in getSolution

When you call getSolution with multiple families the values calculated by the earlier families are zeroed on the iteration.
I don't believe this was ever the intended behavior and frankly I'm surprised this was left on unfixed for so long.
---
 src/solver/surfaceIntegrations.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/solver/surfaceIntegrations.F90 b/src/solver/surfaceIntegrations.F90
index 35542b268..33e233dd7 100644
--- a/src/solver/surfaceIntegrations.F90
+++ b/src/solver/surfaceIntegrations.F90
@@ -959,7 +959,7 @@ subroutine getSolution(famLists, funcValues, globalValues)
        ! Extract the current family list
        nFam = famLists(iGroup, 1)
        famList => famLists(iGroup, 2:2+nFam-1)
-       funcValues = zero
+       funcValues(:, iGroup) = zero
        localVal = zero
 
        do sps=1, nTimeIntervalsSpectral