From fb5552bcb1b4616a2ec2f8901263782002695255 Mon Sep 17 00:00:00 2001
From: acse-ej321 <89605848+acse-ej321@users.noreply.github.com>
Date: Thu, 24 Oct 2024 16:30:55 +0100
Subject: [PATCH] #220 more clarity in comments for clearning subinterval
 fields in indicate errors

---
 goalie/go_mesh_seq.py | 19 ++++++++++---------
 1 file changed, 10 insertions(+), 9 deletions(-)

diff --git a/goalie/go_mesh_seq.py b/goalie/go_mesh_seq.py
index 0d24c799..74b3d288 100644
--- a/goalie/go_mesh_seq.py
+++ b/goalie/go_mesh_seq.py
@@ -232,16 +232,17 @@ def indicate_errors(
                     indi.interpolate(abs(indi))
                     self.indicators[f][i][j].interpolate(ufl.max_value(indi, 1.0e-16))
 
-            # discard extra subinterval solution field
-            self.solutions[f][FWD_OLD].pop()
-            self.solutions[f][ADJ_NEXT].pop()
-            enriched_mesh_seq.solutions[f][FWD_OLD].pop()
-            enriched_mesh_seq.solutions[f][ADJ_NEXT].pop()
-
-            # delete extra mesh to reduce the memory footprint
-            enriched_mesh_seq.meshes.pop()
+            # discard current subinterval solution field
+            for f in self.fields:
+                self.solutions[f][FWD_OLD].pop(-1)
+                self.solutions[f][ADJ_NEXT].pop(-1)
+                enriched_mesh_seq.solutions[f][FWD_OLD].pop(-1)
+                enriched_mesh_seq.solutions[f][ADJ_NEXT].pop(-1)
+
+            # delete current subinterval enriched mesh to reduce the memory footprint
+            enriched_mesh_seq.meshes.pop(-1)
             enriched_mesh_seq.time_partition.num_subintervals -= 1
-            enriched_mesh_seq.time_partition.timesteps.pop()
+            enriched_mesh_seq.time_partition.timesteps.pop(-1)
 
         return self.solutions, self.indicators