Index: /issm/trunk-jpl/src/m/mech/robintemperature.py
===================================================================
--- /issm/trunk-jpl/src/m/mech/robintemperature.py	(revision 17835)
+++ /issm/trunk-jpl/src/m/mech/robintemperature.py	(revision 17836)
@@ -37,5 +37,6 @@
 	
 	tprofile=surftemp+npy.sqrt(2.*thickness*npy.pi/accumrate/alphaT)*(-heatflux)/2./rho/c*(erf(z/zstar)-erf(thickness/zstar))
-	
+
+	return tprofile	
 	# difference between surface and base temperature for check (Cuffey2010 p412):
 	# print tprofile-surftemp