Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 18948)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 18949)
@@ -2253,5 +2253,5 @@
   /*Get material parameters :*/
   rho_ice=matpar->GetMaterialParameter(MaterialsRhoIceEnum);
-  rho_water=matpar->GetMaterialParameter(MaterialsRhoSeawaterEnum);
+  rho_water=matpar->GetMaterialParameter(MaterialsRhoFreshwaterEnum);
 
   /*Get other pdd parameters*/
@@ -2285,6 +2285,6 @@
 
 	/*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);
Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18948)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp	(revision 18949)
@@ -86,5 +86,5 @@
 		/*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