Index: /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 15738)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 15739)
@@ -241,5 +241,5 @@
 	else{
 		if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0)){
-			viscosity=0.5*pow((IssmDouble)10,(IssmDouble)14);
+			viscosity=0.5*pow(10.,14);
 		}
 		else{
@@ -253,5 +253,5 @@
 			if(A==0){
 				/*Maxiviscositym viscosity for 0 shear areas: */
-				viscosity=2.5*pow(10.,17.);
+				viscosity=2.5*pow(10.,17);
 			}
 			else{
@@ -308,5 +308,5 @@
 		if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0) && 
 				(epsilon[3]==0) && (epsilon[4]==0)){
-			viscosity3d=0.5*pow((IssmDouble)10,(IssmDouble)14);
+			viscosity3d=0.5*pow(10.,14);
 		}
 		else{
@@ -323,5 +323,5 @@
 			if(A==0){
 				/*Maxiviscosity3dm viscosity for 0 shear areas: */
-				viscosity3d=2.25*pow((IssmDouble)10,(IssmDouble)17);
+				viscosity3d=2.25*pow(10.,17);
 			}
 			else{
@@ -369,5 +369,5 @@
 
 	/*Get B and n*/
-	eps0=pow((IssmDouble)10,(IssmDouble)-27);
+	eps0=pow(10.,-27);
 	B=GetB();
 	n=GetN();
@@ -380,5 +380,5 @@
 		if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0) && 
 				(epsilon[3]==0) && (epsilon[4]==0) && (epsilon[5]==0)){
-			viscosity3d=0.5*pow((IssmDouble)10,(IssmDouble)14);
+			viscosity3d=0.5*pow(10.,14);
 		}
 		else{
@@ -396,5 +396,5 @@
 			if(A==0){
 				/*Maxiviscosity3dm viscosity for 0 shear areas: */
-				viscosity3d=2.25*pow((IssmDouble)10,(IssmDouble)17);
+				viscosity3d=2.25*pow(10.,17);
 			}
 			else{
@@ -449,5 +449,5 @@
 		if(A==0){
 			/*Maximum viscosity_complement for 0 shear areas: */
-			viscosity_complement=2.25*pow((IssmDouble)10,(IssmDouble)17);
+			viscosity_complement=2.25*pow(10.,17);
 		}
 		else{
@@ -458,5 +458,5 @@
 	}
 	else{
-		viscosity_complement=4.5*pow((IssmDouble)10,(IssmDouble)17);
+		viscosity_complement=4.5*pow(10.,17);
 	}
 
@@ -486,5 +486,5 @@
 	if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0) && 
 				(epsilon[3]==0) && (epsilon[4]==0)){
-		mu_prime=0.5*pow((IssmDouble)10,(IssmDouble)14);
+		mu_prime=0.5*pow(10.,14);
 	}
 	else{
@@ -520,5 +520,5 @@
 	if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0) && 
 				(epsilon[3]==0) && (epsilon[4]==0) && (epsilon[5]==0)){
-		mu_prime=0.5*pow((IssmDouble)10,(IssmDouble)14);
+		mu_prime=0.5*pow(10.,14);
 	}
 	else{
@@ -554,5 +554,5 @@
 
 	if((epsilon[0]==0) && (epsilon[1]==0) && (epsilon[2]==0)){
-		mu_prime=0.5*pow((IssmDouble)10,(IssmDouble)14);
+		mu_prime=0.5*pow(10.,14);
 	}
 	else{