Index: /issm/trunk-jpl/src/m/mech/robintemperature.m
===================================================================
--- /issm/trunk-jpl/src/m/mech/robintemperature.m	(revision 17283)
+++ /issm/trunk-jpl/src/m/mech/robintemperature.m	(revision 17284)
@@ -35,7 +35,7 @@
 
 %create vertical coordinate variable
-zstar=sqrt(2*alphaT*thickness/accumrate);
+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):