Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 27058)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 27059)
@@ -7482,8 +7482,13 @@
 	}
 
+		cout << "DEBUG: Sealevelchange convo:";
 	xDelete<IssmDouble>(G);
+		cout << " step1";
 	xDelete<IssmDouble>(Grot);
+		cout << ", step2";
 	for (int l=0;l<SLGEOM_NUMLOADS;l++) xDelete<IssmDouble>(Gsub[l]);
+		cout << ", step3";
 	xDelete<IssmDouble*>(Gsub);
+		cout << ", step4\n";
 
 	return;
Index: /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp	(revision 27058)
+++ /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp	(revision 27059)
@@ -340,4 +340,6 @@
 		}
 
+		cout << "DEBUG: Sealevelchange convolution cleared\n";
+
 		/*retrieve sea level average  and ocean area:*/
 		for(Object* & object : femmodel->elements->objects){
@@ -345,4 +347,6 @@
 			element->SealevelchangeOceanAverage(loads, oceanareas, subelementoceanareas, sealevelpercpu, slgeom);
 		}
+
+		cout << "DEBUG: Sealevelchange Oceanavaverage cleared\n";
 
 		loads->AssembleSealevelLoads();
@@ -360,4 +364,6 @@
 		ConserveOceanMass(femmodel,loads,barycontrib->Total()/totaloceanarea - oceanaverage,slgeom);
 
+		cout << "DEBUG: Sealevelchange Conserve ocean mass cleared\n";
+
 		//broadcast sea level loads 
 		loads->BroadcastSealevelLoads();
@@ -366,4 +372,5 @@
 		if(slcconvergence(loads->vsealevelloads,oldsealevelloads,eps_rel,eps_abs))break;
 
+		cout << "DEBUG: Sealevelchange convergence cleared\n";
 		//early return?
 		if(iterations>=max_nonlinear_iterations)break;