Index: ../trunk-jpl/src/m/contrib/gravity/vfsa.cpp
===================================================================
--- ../trunk-jpl/src/m/contrib/gravity/vfsa.cpp	(revision 18464)
+++ ../trunk-jpl/src/m/contrib/gravity/vfsa.cpp	(revision 18465)
@@ -129,12 +129,12 @@
 void* plouffT(void* vpthread_handle);
 void vec2gridsimple(Matrix *V,Matrix *V1,int nx, int ny);
 void reshape(Matrix* V,Matrix* V1,int nx,int ny);
-double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,int dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my);
+double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,double dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my);
 void GSLsquarefit(Matrix** pX,Matrix* A,Matrix* B);
 double signe(double a);
 void filtergrav(Matrix* A,Matrix* Ain,double ctr,double sd,int mx,int my);
 void newmodelgen(Matrix* m0,Matrix* m1,Matrix* bathy,Matrix* icethick,int mx,int my,double T,double ptval,double mmax,double mmax2,double ctr,double sd);
-double coolshed(double T0,int k,double c,double D);
+double coolshed(double T0,double k,double c,double D);
 void   LaunchThread(void* function(void*), void* usr,int num_threads);
 /*}}}*/
 
@@ -538,7 +538,7 @@
 		}
 	}
 }/*}}}*/
-double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,int dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my){/*{{{*/
+double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,double dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my){/*{{{*/
 	Matrix* m1=new Matrix(mx*my,4);
 	Matrix* m2=new Matrix(mx*my,3);
 	Matrix* g1=new Matrix(nx*ny,1);
@@ -690,7 +690,7 @@
 		}
 	}
 }/*}}}*/
-double coolshed(double T0,int k,double c,double D){/*{{{*/
+double coolshed(double T0,double k,double c,double D){/*{{{*/
 	double T1=T0*exp(-c*pow(k,1/D));
 	return T1;
 }/*}}}*/
Index: ../trunk-jpl/src/m/contrib/gravity/vfsa_mpi.cpp
===================================================================
--- ../trunk-jpl/src/m/contrib/gravity/vfsa_mpi.cpp	(revision 18464)
+++ ../trunk-jpl/src/m/contrib/gravity/vfsa_mpi.cpp	(revision 18465)
@@ -128,12 +128,12 @@
 void plouff(Matrix *g,Matrix *Pobs,Matrix *Pp,Matrix * mesh,Matrix *rho,int dx,int dy, int dn,int m,int n,int l,int my_rank,int num_procs);
 void vec2gridsimple(Matrix *V,Matrix *V1,int nx, int ny);
 void reshape(Matrix* V,Matrix* V1,int nx,int ny);
-double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,int dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my,int my_rank,int num_procs);
+double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,double dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my,int my_rank,int num_procs);
 void GSLsquarefit(Matrix** pX,Matrix* A,Matrix* B);
 double signe(double a);
 void filtergrav(Matrix* A,Matrix* Ain,double ctr,double sd,int mx,int my);
 void newmodelgen(Matrix* m0,Matrix* m1,Matrix* bathy,Matrix* icethick,int mx,int my,double T,double ptval,double mmax,double mmax2,double ctr,double sd);
-double coolshed(double T0,int k,double c,double D);
+double coolshed(double T0,double k,double c,double D);
 /*}}}*/
 
 int main(int argc,char *argv[]){/*{{{*/
@@ -620,7 +620,7 @@
 		}
 	}
 }/*}}}*/
-double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,int dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my,int my_rank,int num_procs){/*{{{*/
+double misfit(Matrix* m0,Matrix* evalid,Matrix* gobs,double dlevel,Matrix* Pobs,Matrix* xobs,Matrix* yobs,Matrix* Pp,Matrix* rho1, Matrix* rho2,int dx,int dy,int dn,int nx,int ny, int mx,int my,int my_rank,int num_procs){/*{{{*/
 	Matrix* m1=new Matrix(mx*my,4);
 	Matrix* m2=new Matrix(mx*my,3);
 	Matrix* g1=new Matrix(nx*ny,1);
@@ -741,7 +741,7 @@
 		}
 	}
 }/*}}}*/
-double coolshed(double T0,int k,double c,double D){/*{{{*/
+double coolshed(double T0,double k,double c,double D){/*{{{*/
 	double T1=T0*exp(-c*pow(k,1/D));
 	return T1;
 }/*}}}*/