Index: /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp	(revision 25680)
@@ -69,5 +69,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,&extvar_enum,1);
 	}
Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 25680)
@@ -168,5 +168,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving levelset results\n");
 		int outputs[1] = {MaskIceLevelsetEnum};
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],1);
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 25680)
@@ -1050,6 +1050,4 @@
 
 	if(isFS && !(isSSA || isHO || isL1L2 || isMLHO)){
-		if(VerboseSolution()) _printf0_("   computing velocities\n");
-
 		femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum);
 
@@ -1075,6 +1073,4 @@
 	}
 	else if(!isFS && (isSSA || isHO || isL1L2 || isMLHO)){
-		if(VerboseSolution()) _printf0_("   computing velocities\n");
-
 		femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum);
 		if(newton>0)
Index: /issm/trunk-jpl/src/c/cores/adjointbalancethickness2_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/adjointbalancethickness2_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/adjointbalancethickness2_core.cpp	(revision 25680)
@@ -33,5 +33,4 @@
 	/*Save results*/
 	if(save_results || true){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		int outputs[1] = {AdjointEnum};
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],1);
Index: /issm/trunk-jpl/src/c/cores/adjointbalancethickness_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/adjointbalancethickness_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/adjointbalancethickness_core.cpp	(revision 25680)
@@ -33,5 +33,4 @@
 	/*Save results*/
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		int outputs[1] = {AdjointEnum};
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],1);
Index: /issm/trunk-jpl/src/c/cores/balancethickness2_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/balancethickness2_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/balancethickness2_core.cpp	(revision 25680)
@@ -33,5 +33,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		const int numoutputs = 1;
 		int outputs[numoutputs] = {ThicknessEnum};
Index: /issm/trunk-jpl/src/c/cores/balancethickness_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/balancethickness_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/balancethickness_core.cpp	(revision 25680)
@@ -36,5 +36,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		int outputs = ThicknessEnum;
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs,1);
Index: /issm/trunk-jpl/src/c/cores/balancevelocity_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/balancevelocity_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/balancevelocity_core.cpp	(revision 25680)
@@ -32,5 +32,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		int outputs[3] = {DrivingStressXEnum,DrivingStressYEnum,VelEnum};
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3);
Index: /issm/trunk-jpl/src/c/cores/bedslope_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/bedslope_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/bedslope_core.cpp	(revision 25680)
@@ -34,5 +34,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving bedslopes results\n");
 		if(domaintype!=Domain2DverticalEnum){
 			int outputs[2] = {BedSlopeXEnum,BedSlopeYEnum};
Index: /issm/trunk-jpl/src/c/cores/damage_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/damage_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/damage_core.cpp	(revision 25680)
@@ -40,5 +40,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving damage results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
 	}
Index: /issm/trunk-jpl/src/c/cores/depthaverage_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/depthaverage_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/depthaverage_core.cpp	(revision 25680)
@@ -16,6 +16,4 @@
 	int inputenum,input_average_enum;
 
-	if(VerboseSolution()) _printf0_("   depth averaging solution...\n");
-
 	/*Get parameters*/
 	femmodel->parameters->FindParam(&domaintype,DomainTypeEnum);
@@ -30,4 +28,6 @@
 	}
 
+	if(VerboseSolution()) _printf0_("   depth averaging "<<EnumToStringx(inputenum)<<"\n");
+
 	/*Special method for Penta, otherwise call solution sequence*/
 	if(elementtype==PentaEnum){
Index: /issm/trunk-jpl/src/c/cores/esa_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/esa_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/esa_core.cpp	(revision 25680)
@@ -99,5 +99,4 @@
 
 		if(save_results){
-			if(VerboseSolution()) _printf0_("   saving results\n");
 			femmodel->parameters->FindParam(&requested_outputs,&numoutputs,EsaRequestedOutputsEnum);
 			femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
Index: /issm/trunk-jpl/src/c/cores/extrudefrombase_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/extrudefrombase_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/extrudefrombase_core.cpp	(revision 25680)
@@ -15,6 +15,4 @@
 	int elementtype,domaintype;
 
-	if(VerboseSolution()) _printf0_("   extruding solution from base...\n");
-
 	/*Get parameters*/
 	femmodel->parameters->FindParam(&domaintype,DomainTypeEnum);
@@ -24,7 +22,10 @@
 	if(domaintype==Domain2DhorizontalEnum) return;
 
+	int inputenum;
+	femmodel->parameters->FindParam(&inputenum,InputToExtrudeEnum);
+	if(VerboseSolution()) _printf0_("   extruding "<<EnumToStringx(inputenum)<<" from base...\n");
+
 	/*Special method for Penta, otherwise call solution sequence*/
 	if(elementtype==PentaEnum){
-		int inputenum; femmodel->parameters->FindParam(&inputenum,InputToExtrudeEnum);
 		InputExtrudex(femmodel,inputenum,-1);
 	}
Index: /issm/trunk-jpl/src/c/cores/gia_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/gia_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/gia_core.cpp	(revision 25680)
@@ -58,5 +58,4 @@
 
 		if(save_results){
-			if(VerboseSolution()) _printf0_("   saving results\n");
 			int outputs[2] = {UGiaEnum,UGiaRateEnum};
 			femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],2);
Index: /issm/trunk-jpl/src/c/cores/groundingline_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/groundingline_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/groundingline_core.cpp	(revision 25680)
@@ -22,4 +22,5 @@
 
 	/*Move grounding line*/
+	if(VerboseSolution()) _printf0_("   computing new grounding line position\n");
 	GroundinglineMigrationx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters);
 
@@ -33,5 +34,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving groundingline results\n");
 		int outputs[3] = {SurfaceEnum,BaseEnum,MaskOceanLevelsetEnum};
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],3);
Index: /issm/trunk-jpl/src/c/cores/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/hydrology_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/hydrology_core.cpp	(revision 25680)
@@ -187,5 +187,4 @@
 	}
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving hydrology results \n");
 		if(hydrology_model==HydrologydcEnum && ThawedNodes==0){
 			if(VerboseSolution()) _printf0_("   No thawed node hydro is skiped \n");}
Index: /issm/trunk-jpl/src/c/cores/masstransport_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/masstransport_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/masstransport_core.cpp	(revision 25680)
@@ -85,5 +85,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
 	}
Index: /issm/trunk-jpl/src/c/cores/movingfront_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/movingfront_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/movingfront_core.cpp	(revision 25680)
@@ -48,8 +48,10 @@
 	}
 
+
 	/* start the work from here */
+	if(VerboseSolution()) _printf0_("   computing calving and undercutting\n");
 	Calvingx(femmodel);
 	FrontalForcingsx(femmodel);
-	if(VerboseSolution()) _printf0_("   computing level set transport\n");
+	if(VerboseSolution()) _printf0_("   computing new ice front position\n");
 
 	/* smoothen slope of lsf for computation of normal on ice domain*/
Index: /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp	(revision 25680)
@@ -65,6 +65,4 @@
 		int     numoutputs;
 		char **requested_outputs = NULL;
-
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		femmodel->parameters->FindParam(&requested_outputs,&numoutputs,SealevelriseRequestedOutputsEnum);
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
Index: /issm/trunk-jpl/src/c/cores/smb_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/smb_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/smb_core.cpp	(revision 25680)
@@ -103,5 +103,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving smb results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
 	}
Index: /issm/trunk-jpl/src/c/cores/steadystate_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/steadystate_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/steadystate_core.cpp	(revision 25680)
@@ -78,5 +78,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
 	}
Index: /issm/trunk-jpl/src/c/cores/stressbalance_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/stressbalance_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/stressbalance_core.cpp	(revision 25680)
@@ -94,7 +94,5 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving stressbalance results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
-		if(VerboseSolution()) _printf0_("   results saved\n");
 	}
 
Index: /issm/trunk-jpl/src/c/cores/surfaceslope_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/surfaceslope_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/surfaceslope_core.cpp	(revision 25680)
@@ -38,5 +38,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("saving surface slopes results:\n");
 		if(domaintype!=Domain2DverticalEnum){
 			int outputs[2] = {SurfaceSlopeXEnum,SurfaceSlopeYEnum};
Index: /issm/trunk-jpl/src/c/cores/thermal_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/thermal_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/thermal_core.cpp	(revision 25680)
@@ -50,5 +50,4 @@
 
 	if(save_results){
-		if(VerboseSolution()) _printf0_("   saving thermal results\n");
 		femmodel->RequestedOutputsx(&femmodel->results,requested_outputs,numoutputs);
 	}
Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 25679)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 25680)
@@ -190,10 +190,8 @@
 	/* Using Hydrology dc  coupled we need to compute smb in the hydrology inner time loop*/
 	if(issmb) {
-		if(VerboseSolution()) _printf0_("   computing smb\n");
 		smb_core(femmodel);
 	}
 
 	if(ishydrology){
-		if(VerboseSolution()) _printf0_("   computing hydrology\n");
 		int hydrology_model;
 		hydrology_core(femmodel);
@@ -203,15 +201,12 @@
 
 	if(isstressbalance && (step%sb_coupling_frequency==0 || step==1) ) {
-		if(VerboseSolution()) _printf0_("   computing stress balance\n");
 		stressbalance_core(femmodel);
 	}
 
 	if(isdamageevolution) {
-		if(VerboseSolution()) _printf0_("   computing damage\n");
 		damage_core(femmodel);
 	}
 
 	if(ismovingfront)	{
-		if(VerboseSolution()) _printf0_("   computing moving front\n");
 		movingfront_core(femmodel);
 	}
@@ -220,5 +215,4 @@
 
 	if(ismasstransport){
-		if(VerboseSolution()) _printf0_("   computing mass transport\n");
 		bmb_core(femmodel);
 		masstransport_core(femmodel);
@@ -227,10 +221,8 @@
 
 	if(isgroundingline){
-		if(VerboseSolution()) _printf0_("   computing new grounding line position\n");
 		groundingline_core(femmodel);
 	}
 
 	if(isgia){
-		if(VerboseSolution()) _printf0_("   computing glacial isostatic adjustment\n");
 		#ifdef _HAVE_GIA_
 		gia_core(femmodel);
@@ -461,4 +453,5 @@
 		delete hdl_setadjoint;
 
+		if(VerboseSolution()) _printf0_("   computing gradient\n");
 		tape_codi.evaluate();