Changeset 16341

Timestamp:

10/08/13 14:44:08 (11 years ago)

Author:

cborstad

Message:

missed negative sign for heat flux

File:

• issm/trunk-jpl/src/m/mech/robintemperature.m (modified) (2 diffs)

issm/trunk-jpl/src/m/mech/robintemperature.m

@@ -21 +21 @@
 %
 %       Usage:
-%               tprofile=robintemperature(heatflux,accumrate,thickness,surftemp)
+%               tprofile=robintemperature(heatflux,accumrate,thickness,surftemp,z)
 
 %checks
@@ -37 +37 @@
 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