Index: /issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.cpp	(revision 21855)
+++ /issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.cpp	(revision 21856)
@@ -12,14 +12,14 @@
 using namespace std;
 
-int InterpFromMeshToMesh2dx(IssmDouble** pdata_interp,int* index_data,IssmDouble* x_data,IssmDouble* y_data,int nods_data,int nels_data,
-			IssmDouble* data,int M_data,int N_data,IssmDouble* x_interp,IssmDouble* y_interp,int N_interp,Options* options){
+int InterpFromMeshToMesh2dx(double** pdata_interp,int* index_data,double* x_data,double* y_data,int nods_data,int nels_data,
+			double* data,int M_data,int N_data,double* x_interp,double* y_interp,int N_interp,Options* options){
 
 	/*Output*/
-	IssmDouble* data_interp=NULL;
+	double* data_interp=NULL;
 
 	/*Intermediary*/
-	IssmDouble xmin,xmax,ymin,ymax;
+	double xmin,xmax,ymin,ymax;
 	bool   isdefault;
-	IssmDouble defaultvalue;
+	double defaultvalue;
 	R2     r;
 	I2     I;
@@ -27,6 +27,6 @@
 	int    it;
 	int    i0,i1,i2;
-	IssmDouble areacoord[3];
-	IssmDouble aa,bb;
+	double areacoord[3];
+	double aa,bb;
 	long long dete[3];
 
@@ -46,5 +46,5 @@
 
 	/*Initialize output*/
-	data_interp=xNew<IssmDouble>(N_interp*N_data);
+	data_interp=xNew<double>(N_interp*N_data);
 
 	/*read background mesh*/
@@ -78,5 +78,5 @@
 	/*Loop over output nodes*/
 	for(i=0;i<N_interp;i++){
-		//if(i%100==0) _printf_("\r      interpolation progress: "<<setw(6)<<setprecision(2)<<IssmDouble(i)/IssmDouble(N_interp)*100.<<"%   ");
+		//if(i%100==0) _printf_("\r      interpolation progress: "<<setw(6)<<setprecision(2)<<double(i)/double(N_interp)*100.<<"%   ");
 
 		if(isdefault){
@@ -98,7 +98,7 @@
 
 			/*Area coordinates*/
-			areacoord[0]= reCast<IssmDouble>(dete[0])/reCast<IssmDouble>(tb.det);
-			areacoord[1]= reCast<IssmDouble>(dete[1])/reCast<IssmDouble>(tb.det);
-			areacoord[2]= reCast<IssmDouble>(dete[2])/reCast<IssmDouble>(tb.det);
+			areacoord[0]= reCast<double>(dete[0])/reCast<double>(tb.det);
+			areacoord[1]= reCast<double>(dete[1])/reCast<double>(tb.det);
+			areacoord[2]= reCast<double>(dete[2])/reCast<double>(tb.det);
 			/*3 vertices of the triangle*/
 			i0=Th->GetId(tb[0]);
Index: /issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.h
===================================================================
--- /issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.h	(revision 21855)
+++ /issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.h	(revision 21856)
@@ -8,6 +8,6 @@
 #include "../../classes/classes.h"
 
-int InterpFromMeshToMesh2dx(IssmDouble** pdata_interp,int* index_data,IssmDouble* x_data,IssmDouble* y_data,int nods_data,int nels_data,
-			IssmDouble* data,int M_data,int N_data,IssmDouble* x_interp,IssmDouble* y_interp,int N_interp,Options* options);
+int InterpFromMeshToMesh2dx(double** pdata_interp,int* index_data,double* x_data,double* y_data,int nods_data,int nels_data,
+			double* data,int M_data,int N_data,double* x_interp,double* y_interp,int N_interp,Options* options);
 
 #endif