Index: /issm/trunk-jpl/src/m/contrib/gravity/vfsa_mpi.cpp
===================================================================
--- /issm/trunk-jpl/src/m/contrib/gravity/vfsa_mpi.cpp	(revision 18478)
+++ /issm/trunk-jpl/src/m/contrib/gravity/vfsa_mpi.cpp	(revision 18479)
@@ -147,5 +147,5 @@
 
 	/* Seed the random number generator {{{*/
-	srand (time(NULL));               /*}}}*/
+		srand (time(NULL));            /*}}}*/
 	/* Define the variables {{{*/
 
@@ -444,5 +444,4 @@
 	int    M,N;
 	double chisq;
-
 	/*Get system size*/
 	A->GetSize(&M,&N);
@@ -533,5 +532,4 @@
 
 	double *glocal=new double[n]();
-	double *gcalc=new double[n]();
 
 	for(int c=my_rank;c<n;c+=num_procs){
@@ -567,5 +565,5 @@
 					si*=-1;
 					for(int j=0;j<2;j++){
-						si*=-1;
+						sj*=-1;
 						s=si*sj;
 						for(int k=0;k<l;k++){
@@ -588,9 +586,5 @@
 	}
 
-	MPI_Allreduce(glocal,gcalc,n,MPI_DOUBLE,MPI_SUM,MPI_COMM_WORLD);
-	if(!my_rank){for(int c=0;c<n;c++){
-		g->SetValue(c,0,gcalc[c]);
-	}
-	}
+	MPI_Allreduce(glocal,g->GetPointer(),n,MPI_DOUBLE,MPI_SUM,MPI_COMM_WORLD);
 
 	delete xp;
@@ -603,5 +597,4 @@
 	delete U2;
 	delete []glocal;
-	delete []gcalc;
 }/*}}}*/
 void vec2gridsimple(Matrix *V,Matrix *V1,int nx, int ny){/*{{{*/