Index: /issm/branches/trunk-larour-NatGeoScience2016/src/c/analyses/SealevelriseAnalysis.cpp
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/c/analyses/SealevelriseAnalysis.cpp	(revision 22162)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/c/analyses/SealevelriseAnalysis.cpp	(revision 22163)
@@ -62,4 +62,5 @@
 			element->InputUpdateFromConstant(0.0,SealevelNEsaRateEnum);
 			element->InputUpdateFromConstant(0.0,SealevelUEsaRateEnum);
+			element->InputUpdateFromConstant(0.0,SealevelRSLEustaticRateEnum);
 			counter++;
 		}
Index: /issm/branches/trunk-larour-NatGeoScience2016/src/c/classes/FemModel.cpp
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/c/classes/FemModel.cpp	(revision 22162)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/c/classes/FemModel.cpp	(revision 22163)
@@ -2458,5 +2458,5 @@
 		if(i<ns){
 		
-			if(VerboseConvergence())if(i%100==0)_printf0_("\r" << "      convolution progress: " << (double)i/(double)ns*100 << "%  ");
+			if(VerboseConvergence())if(i%100==0)_printf0_("\r" << "              convolution progress: " << (double)i/(double)ns*100 << "%  ");
 		
 			Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(i));
@@ -2513,5 +2513,5 @@
 	for(int i=0;i<nsmax;i++){
 		if(i<ns){
-			if(verboseconvolution)if(VerboseConvergence())if(i%100==0)_printf_("\r" << "      convolution progress: " << (double)i/(double)ns*100 << "%   ");
+			if(verboseconvolution)if(VerboseConvergence())if(i%100==0)_printf0_("\r" << "              convolution progress: " << (double)i/(double)ns*100 << "%   ");
 			Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(i));
 			element->SealevelriseNonEustatic(pRSLgo,RSLg_old,latitude,longitude,radius,eartharea);
@@ -2519,5 +2519,5 @@
 		if(i%loop==0)pRSLgo->Assemble();
 	}
-	if(verboseconvolution)if(VerboseConvergence())_printf_("\n");
+	if(verboseconvolution)if(VerboseConvergence())_printf0_("\n");
 	
 	/*Free ressources:*/
@@ -2637,5 +2637,5 @@
 	for(int i=0;i<nsmax;i++){
 		if(i<ns){
-			if(VerboseConvergence())if(i%100==0)_printf0_("\r" << "      convolution progress: " << (double)i/(double)ns*100 << "%  ");
+			if(VerboseConvergence())if(i%100==0)_printf0_("\r" << "              convolution progress: " << (double)i/(double)ns*100 << "%  ");
 			Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(i));
 			element->SealevelriseGeodetic(pUp,pNorth,pEast,RSLg,latitude,longitude,radius,xx,yy,zz,eartharea,horiz);
@@ -2656,5 +2656,5 @@
 		pEast->Assemble();
 	}
-	if(VerboseConvergence())_printf_("\n");
+	if(VerboseConvergence())_printf0_("\n");
 
 	/*Free ressources:*/
Index: /issm/branches/trunk-larour-NatGeoScience2016/src/c/cores/sealevelrise_core.cpp
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/c/cores/sealevelrise_core.cpp	(revision 22162)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/c/cores/sealevelrise_core.cpp	(revision 22163)
@@ -458,5 +458,5 @@
 		
 	}
-	if(VerboseConvergence()) _printf0_("\n   total number of iterations: " << count-1 << "\n");
+	if(VerboseConvergence()) _printf0_("\n              total number of iterations: " << count-1 << "\n");
 
 	xDelete<IssmDouble>(latitude);
@@ -803,5 +803,5 @@
 	int nv;
 	
-	if(VerboseSolution()) _printf0_("      computing earth mass transport\n");
+	if(VerboseSolution()) _printf0_("              computing earth mass transport\n");
 
 	/*This mass transport module for the Earth is because we might have thickness variations as spcs 
@@ -861,8 +861,8 @@
 	if(!xIsNan<IssmDouble>(eps_rel)){
 		if((ndS/nS)<eps_rel){
-			if(VerboseConvergence()) _printf0_(setw(50) << left << "      convergence criterion: norm(dS)/norm(S)" << ndS/nS*100 << " < " << eps_rel*100 << " %\n");
+			if(VerboseConvergence()) _printf0_(setw(50) << left << "              convergence criterion: norm(dS)/norm(S)" << ndS/nS*100 << " < " << eps_rel*100 << " %");
 		}
 		else{ 
-			if(VerboseConvergence()) _printf0_(setw(50) << left << "      convergence criterion: norm(dS)/norm(S)" << ndS/nS*100 << " > " << eps_rel*100 << " %\n");
+			if(VerboseConvergence()) _printf0_(setw(50) << left << "              convergence criterion: norm(dS)/norm(S)" << ndS/nS*100 << " > " << eps_rel*100 << " %");
 			converged=false;
 		}
@@ -870,8 +870,8 @@
 	if(!xIsNan<IssmDouble>(eps_abs)){
 		if(ndS<eps_abs){
-			if(VerboseConvergence()) _printf0_(setw(50) << left << "      convergence criterion: norm(dS)" << ndS << " < " << eps_abs << " \n");
+			if(VerboseConvergence()) _printf0_(setw(50) << left << "              convergence criterion: norm(dS)" << ndS << " < " << eps_abs << " ");
 		}
 		else{ 
-			if(VerboseConvergence()) _printf0_(setw(50) << left << "      convergence criterion: norm(dS)" << ndS << " > " << eps_abs << " \n");
+			if(VerboseConvergence()) _printf0_(setw(50) << left << "              convergence criterion: norm(dS)" << ndS << " > " << eps_abs << " ");
 			converged=false;
 		}