Index: /issm/trunk/src/c/objects/ElementResults/NodalValuesUnitConversion.cpp
===================================================================
--- /issm/trunk/src/c/objects/ElementResults/NodalValuesUnitConversion.cpp	(revision 4972)
+++ /issm/trunk/src/c/objects/ElementResults/NodalValuesUnitConversion.cpp	(revision 4973)
@@ -20,15 +20,22 @@
 	int i;
 	double yts;
+	double scale;
 
 	/*Get some values from parameters, that are often needed: */
 	parameters->FindParam(&yts,YtsEnum);
 
+	/*Get scaling factor: */
 	switch(enum_type){
-		case VxEnum:  for(i=0;i<num_nodal_values;i++)nodal_values[i]=nodal_values[i]*yts;break;
-		case VyEnum:  for(i=0;i<num_nodal_values;i++)nodal_values[i]=nodal_values[i]*yts;break;
-		case VzEnum:  for(i=0;i<num_nodal_values;i++)nodal_values[i]=nodal_values[i]*yts;break;
-		case VelEnum: for(i=0;i<num_nodal_values;i++)nodal_values[i]=nodal_values[i]*yts;break;
-		case MeltingRateEnum: for(i=0;i<num_nodal_values;i++)nodal_values[i]=nodal_values[i]*yts;break;
+		case VxEnum:  scale=yts;break; //m/yr
+		case VyEnum:  scale=yts;break; //m/yr
+		case VzEnum:  scale=yts;break; //m/yr
+		case VelEnum: scale=yts;break; //m/yr
+		case MeltingRateEnum: scale=yts;break; //m/yr
+		case MisfitEnum: scale=pow(yts,2);break; //(m/yr)^2
+		case MassFluxEnum: scale=pow(10,-12)*yts;break; // (GigaTon/year)
 		default: break;
 	}
+		
+	/*Scale : */
+	for(i=0;i<num_nodal_values;i++)nodal_values[i]=nodal_values[i]*scale; 
 }