CarpetLib: Free ggf data structures in ggf, not in gf<T>

1 files changed, 10 insertions, 0 deletions

diff --git a/Carpet/CarpetLib/src/ggf.cc b/Carpet/CarpetLib/src/ggf.cc
index 84f867c46..7ec02ff30 100644
--- a/Carpet/CarpetLib/src/ggf.cc
+++ b/Carpet/CarpetLib/src/ggf.cc
@@ -52,10 +52,20 @@ ggf::ggf (const int varindex_, const operator_type transport_operator_,
   
   allggf.insert (this);
   d.insert (this);
+  
+  recompose_crop ();
+  for (int rl=0; rl<h.reflevels(); ++rl) {
+    recompose_allocate (rl);
+    recompose_free_old (rl);
+  } // for rl
 }
 
 // Destructors
 ggf::~ggf () {
+  for (int rl=0; rl<h.reflevels(); ++rl) {
+    recompose_free (rl);
+  } // for rl
+  
   d.erase (this);
   allggf.erase (this);
 }