Index: /issm/trunk-jpl/src/m/mech/robintemperature.m
===================================================================
--- /issm/trunk-jpl/src/m/mech/robintemperature.m	(revision 16340)
+++ /issm/trunk-jpl/src/m/mech/robintemperature.m	(revision 16341)
@@ -21,5 +21,5 @@
 %
 %	Usage:
-%		tprofile=robintemperature(heatflux,accumrate,thickness,surftemp)
+%		tprofile=robintemperature(heatflux,accumrate,thickness,surftemp,z)
 
 %checks
@@ -37,6 +37,6 @@
 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 temperautre for check (Cuffey2010 p412):
+% difference between surface and base temperature for check (Cuffey2010 p412):
 % tprofile-surftemp