Changeset 17284

Timestamp:

02/14/14 17:03:27 (11 years ago)

Author:

cborstad

Message:

BUG: fix for elementwise division in case accumrate is a vector

File:

• issm/trunk-jpl/src/m/mech/robintemperature.m (modified) (1 diff)

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

@@ -35 +35 @@
 
 %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):