Index: /issm/trunk-jpl/src/m/contrib/adhikari/sphericalharmonic_slr_solver/SHslr.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/adhikari/sphericalharmonic_slr_solver/SHslr.m	(revision 25447)
+++ /issm/trunk-jpl/src/m/contrib/adhikari/sphericalharmonic_slr_solver/SHslr.m	(revision 25448)
@@ -94,7 +94,7 @@
 	% Test whether mass is conserved (to ensure the solution is trustworthy!) {{{ 
 	% ocean-average sea level and eustatic value must be equal! 
-	avg_eustatic = -force_lm(1)/oce_lm(1); % see equation B16 
-	barystatic = sum(sealevel.*ocean)/nnz(ocean); % we can do this because of equal-area pixelization 
-	sol_diff = abs(abs(barystatic) - abs(avg_eustatic))/abs(barystatic)*100; % per cent 
+	barystatic = -force_lm(1)/oce_lm(1); % see equation B16 
+	gmsl = sum(sealevel.*ocean)/nnz(ocean); % we can do this because of equal-area pixelization 
+	sol_diff = abs(abs(gmsl) - abs(barystatic))/abs(gmsl)*100; % per cent 
 	if (sol_diff > 1) % per cent. 
 		error('Ocean-average (eustatic - sealevel) is NOT negligible. Returning.');