Index: ../trunk-jpl/src/m/contrib/gravity/vfsa.cpp
===================================================================
--- ../trunk-jpl/src/m/contrib/gravity/vfsa.cpp	(revision 18453)
+++ ../trunk-jpl/src/m/contrib/gravity/vfsa.cpp	(revision 18454)
@@ -299,7 +299,6 @@
 	Matrix* m_min    = new Matrix(mx *my,3);
 	Matrix* m_new    = new Matrix(mx *my,3);
 	m_old->MatrixEqual(mesh_ini);
-	mesh_ini->Echo();
 
 	/* calculate initial misfit */
 	double E_old=misfit(m_old,evalid,gobs,dlevel,Pobs,xobs,yobs,Pp,rho1,rho2,dx,dy,dn,nx,ny,mx,my);
@@ -309,13 +308,14 @@
 	savefile1 << "nan   "<<  Ta<<"   "<< Tm<<"   "<< E_old<<"     "<< totaliter<<"         "<< Tred <<"  "<< endl;
 	savefile2 << totaliter<< endl;
 	for(int i=0;i<mx*my;i++){
-	savefile2 << m_old->GetValue(i,0)<<"   "<< m_old->GetValue(i,1)<<"   "<< m_old->GetValue(i,2)<<endl;
+		savefile2 << m_old->GetValue(i,0)<<"   "<< m_old->GetValue(i,1)<<"   "<< m_old->GetValue(i,2)<<endl;
 	}
-	
+	savefile2 << "111111111111111111111111111111111111111111111111111111111111111111111111111"<< endl;
+
 	/* beginning of the loop */
 
 	while(finished==0){
-		
+
 		iterT++;
 		totaliter++;
 
@@ -340,7 +340,7 @@
 		/* update model and calculate energy */
 
 		newmodelgen(m_old,m_new,bathy,icethick,mx,my,Tm,ptval,mmax,mmax2,ctr,sd);  /* new model */
-	   E_new=misfit(m_new,evalid,gobs,dlevel,Pobs,xobs,yobs,Pp,rho1,rho2,dx,dy,dn,nx,ny,mx,my); /* new energy */
+		E_new=misfit(m_new,evalid,gobs,dlevel,Pobs,xobs,yobs,Pp,rho1,rho2,dx,dy,dn,nx,ny,mx,my); /* new energy */
 		dE=E_new-E_old;                                        /* energy difference */
 
 		/* acceptance probability */
@@ -368,6 +368,7 @@
 				for(int i=0;i<mx*my;i++){
 					savefile2 << m_old->GetValue(i,0)<<"   "<< m_old->GetValue(i,1)<<"   "<< m_old->GetValue(i,2)<<endl;
 				}
+				savefile2 << "111111111111111111111111111111111111111111111111111111111111111111111111111"<< endl;
 			}
 		}
 		else{
@@ -382,6 +383,7 @@
 					for(int i=0;i<mx*my;i++){
 						savefile2 << m_old->GetValue(i,0)<<"   "<< m_old->GetValue(i,1)<<"   "<< m_old->GetValue(i,2)<<endl;
 					}
+					savefile2 << "111111111111111111111111111111111111111111111111111111111111111111111111111"<< endl;
 				}
 			}
 			else{
@@ -397,6 +399,8 @@
 	for(int i=0;i<mx*my;i++){
 		savefile2 << m_min->GetValue(i,0)<<"   "<< m_min->GetValue(i,1)<<"   "<< m_min->GetValue(i,2)<<endl;
 	}
+	savefile1.close;
+	savefile2.close;
 			/*}}}*/
 	return 0;
 }/*}}}*/
@@ -629,7 +633,6 @@
 	}
 	m1->ExtractColumn(m1col,1);
 	vec2gridsimple(m1col,m1gr,mx,my);
-	cout<<endl<<endl;
 	filtergrav(m1grsm,m1gr,ctr,sd,mx,my);
 	reshape(m1grsm,m1col,mx,my);
 	for (int i=0;i<mx*my;i++){
