Index: ../trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp
===================================================================
--- ../trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18948)
+++ ../trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18949)
@@ -85,7 +85,7 @@
 
 		/*Get material parameters :*/
 		rho_ice=element->matpar->GetMaterialParameter(MaterialsRhoIceEnum);
-		rho_water=element->matpar->GetMaterialParameter(MaterialsRhoSeawaterEnum);
+		rho_water=element->matpar->GetMaterialParameter(MaterialsRhoFreshwaterEnum);
 
 		// loop over all vertices
 		for(v=0;v<numvertices;v++){
Index: ../trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 18948)
+++ ../trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 18949)
@@ -2252,7 +2252,7 @@
 
   /*Get material parameters :*/
   rho_ice=matpar->GetMaterialParameter(MaterialsRhoIceEnum);
-  rho_water=matpar->GetMaterialParameter(MaterialsRhoSeawaterEnum);
+  rho_water=matpar->GetMaterialParameter(MaterialsRhoFreshwaterEnum);
 
   /*Get other pdd parameters*/
   desfac=matpar->GetMaterialParameter(SurfaceforcingsDesfacEnum);
@@ -2284,8 +2284,8 @@
 	IssmDouble  rho_water,rho_ice,density;
 
 	/*material parameters: */
-	rho_water=matpar->GetMaterialParameter(MaterialsRhoIceEnum);
-	rho_ice=matpar->GetMaterialParameter(MaterialsRhoSeawaterEnum);
+	rho_water=matpar->GetMaterialParameter(MaterialsRhoSeawaterEnum);
+	rho_ice=matpar->GetMaterialParameter(MaterialsRhoIceEnum);
 	density=rho_ice/rho_water;
 	GetInputListOnVertices(&h[0],ThicknessEnum);
 	GetInputListOnVertices(&r[0],BedEnum);
