Index: /issm/trunk-jpl/src/m/mech/robintemperature.m
===================================================================
--- /issm/trunk-jpl/src/m/mech/robintemperature.m	(revision 16341)
+++ /issm/trunk-jpl/src/m/mech/robintemperature.m	(revision 16342)
@@ -37,5 +37,5 @@
 zstar=sqrt(2*alphaT*thickness/accumrate);
 
-tprofile=surftemp+sqrt(2*thickness*pi/accumrate/alphaT)*(-heatflux)/2/rho/c*(erf(z/zstar)-erf(thickness/zstar));
+tprofile=surftemp+sqrt(2*thickness*pi/accumrate/alphaT)*(-heatflux)/2/rho/c.*(erf(z./zstar)-erf(thickness./zstar));
 
 % difference between surface and base temperature for check (Cuffey2010 p412):