diff --git a/src/print.cpp b/src/print.cpp
index c9aceb3e8..aa43cb58b 100644
--- a/src/print.cpp
+++ b/src/print.cpp
@@ -1,4 +1,3 @@
-// -*- coding: windows-1252 -*-
 #include "Utils.h"
 #include "Phreeqc.h"
 #include "phqalloc.h"
@@ -1478,7 +1477,7 @@ print_species(void)
 			   "Molality", "Activity", "Molality", "Activity", "Gamma", "cm3/mol"));
 #else
 	output_msg(sformatf("   %-13s%12s%12s%10s%10s%10s%11s\n\n", "Species",
-			   "Molality", "Activity", "Molality", "Activity", "Gamma", "cm³/mol"));
+			   "Molality", "Activity", "Molality", "Activity", "Gamma", "cm³/mol"));
 #endif
 /*
  *   Print list of species
@@ -1646,7 +1645,7 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 				output_msg(sformatf("\t%11.3e  sigma, C/m2\n",
 #else
-				output_msg(sformatf("\t%11.3e  sigma, C/m²\n",
+				output_msg(sformatf("\t%11.3e  sigma, C/m²\n",
 #endif
 						   (double) (charge * F_C_MOL /
 									 (charge_ptr->Get_specific_area() *
@@ -1657,7 +1656,7 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 				output_msg(sformatf("\tundefined  sigma, C/m2\n"));
 #else
-				output_msg(sformatf("\tundefined  sigma, C/m²\n"));
+				output_msg(sformatf("\tundefined  sigma, C/m²\n"));
 #endif
 			}
 			if (use.Get_surface_ptr()->Get_type() == cxxSurface::CCM)
@@ -1679,7 +1678,7 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  specific area, m2/mol %s\n",
 #else
-						   "\t%11.3e  specific area, m²/mol %s\n",
+						   "\t%11.3e  specific area, m²/mol %s\n",
 #endif
 						   (double) charge_ptr->Get_specific_area(),
 						   comp_ptr->Get_phase_name().c_str()));
@@ -1687,7 +1686,7 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  m2 for %11.3e moles of %s\n\n",
 #else
-						   "\t%11.3e  m² for %11.3e moles of %s\n\n",
+						   "\t%11.3e  m² for %11.3e moles of %s\n\n",
 #endif
 						   (double) (charge_ptr->Get_grams() *
 									 charge_ptr->Get_specific_area()),
@@ -1700,7 +1699,7 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  specific area, m2/mol %s\n",
 #else
-						   "\t%11.3e  specific area, m²/mol %s\n",
+						   "\t%11.3e  specific area, m²/mol %s\n",
 #endif
 						   (double) charge_ptr->Get_specific_area(),
 						   comp_ptr->Get_rate_name().c_str()));
@@ -1708,7 +1707,7 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  m2 for %11.3e moles of %s\n\n",
 #else
-						   "\t%11.3e  m² for %11.3e moles of %s\n\n",
+						   "\t%11.3e  m² for %11.3e moles of %s\n\n",
 #endif
 						   (double) (charge_ptr->Get_grams() *
 									 charge_ptr->Get_specific_area()),
@@ -1721,13 +1720,13 @@ print_surface(void)
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  specific area, m2/g\n",
 #else
-						   "\t%11.3e  specific area, m²/g\n",
+						   "\t%11.3e  specific area, m²/g\n",
 #endif
 						   (double) charge_ptr->Get_specific_area()));
 #ifdef NO_UTF8_ENCODING
 				output_msg(sformatf("\t%11.3e  m2 for %11.3e g\n\n",
 #else
-				output_msg(sformatf("\t%11.3e  m² for %11.3e g\n\n",
+				output_msg(sformatf("\t%11.3e  m² for %11.3e g\n\n",
 #endif
 						   (double) (charge_ptr->Get_specific_area() *
 									 charge_ptr->Get_grams()),
@@ -1942,28 +1941,28 @@ print_surface_cd_music(void)
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  sigma, plane 0, C/m2\n",
 #else
-						   "\t%11.3e  sigma, plane 0, C/m²\n",
+						   "\t%11.3e  sigma, plane 0, C/m²\n",
 #endif
 						   (double) charge_ptr->Get_sigma0()));
 				output_msg(sformatf(
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  sigma, plane 1, C/m2\n",
 #else
-						   "\t%11.3e  sigma, plane 1, C/m²\n",
+						   "\t%11.3e  sigma, plane 1, C/m²\n",
 #endif
 						   (double) charge_ptr->Get_sigma1()));
 				output_msg(sformatf(
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  sigma, plane 2, C/m2\n",
 #else
-						   "\t%11.3e  sigma, plane 2, C/m²\n",
+						   "\t%11.3e  sigma, plane 2, C/m²\n",
 #endif
 						   (double) charge_ptr->Get_sigma2()));
 				output_msg(sformatf(
 #ifdef NO_UTF8_ENCODING
 						   "\t%11.3e  sigma, diffuse layer, C/m2\n\n",
 #else
-						   "\t%11.3e  sigma, diffuse layer, C/m²\n\n",
+						   "\t%11.3e  sigma, diffuse layer, C/m²\n\n",
 #endif
 						   (double) charge_ptr->Get_sigmaddl()));
 			}
@@ -1972,7 +1971,7 @@ print_surface_cd_music(void)
 #ifdef NO_UTF8_ENCODING
 				output_msg(sformatf("\tundefined  sigma, C/m2\n"));
 #else
-				output_msg(sformatf("\tundefined  sigma, C/m²\n"));
+				output_msg(sformatf("\tundefined  sigma, C/m²\n"));
 #endif
 			}
 			output_msg(sformatf("\t%11.3e  psi, plane 0, V\n",
@@ -2231,7 +2230,7 @@ print_totals(void)
 				   "Specific Conductance (uS/cm, ", tc_x, "oC)  = ", (int) SC));
 #else
 		output_msg(sformatf("%35s%3.0f%7s%i\n",
-				   "Specific Conductance (µS/cm, ", tc_x, "°C)  = ", (int) SC));
+				   "Specific Conductance (µS/cm, ", tc_x, "°C)  = ", (int) SC));
 #endif
 	}
 /* VP: Density Start */
@@ -2241,7 +2240,7 @@ print_totals(void)
 #ifdef NO_UTF8_ENCODING
 		output_msg(sformatf("%45s%9.5f", "Density (g/cm3)  = ",
 #else
-		output_msg(sformatf("%45s%9.5f", "Density (g/cm³)  = ",
+		output_msg(sformatf("%45s%9.5f", "Density (g/cm³)  = ",
 #endif
 			   (double) dens));
 		if (state == INITIAL_SOLUTION && use.Get_solution_ptr()->Get_initial_data()->Get_calc_density())
@@ -2266,7 +2265,7 @@ print_totals(void)
 				   " (solute contributions limited to 200 oC)"));
 #else
 			output_msg(sformatf("%18s\n",
-				   " (solute contributions limited to 200 °C)"));
+				   " (solute contributions limited to 200 °C)"));
 #endif
 		}
 		else output_msg(sformatf("\n"));
@@ -2296,7 +2295,7 @@ print_totals(void)
 #ifdef NO_UTF8_ENCODING
 	output_msg(sformatf("%45s%6.2f\n", "Temperature (oC)  = ",
 #else
-	output_msg(sformatf("%45s%6.2f\n", "Temperature (°C)  = ",
+	output_msg(sformatf("%45s%6.2f\n", "Temperature (°C)  = ",
 #endif
 			   (double) tc_x));