Changeset 16342


Ignore:
Timestamp:
10/08/13 15:15:09 (11 years ago)
Author:
cborstad
Message:

ensure element-wise operations in temperature calculation

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r16341 r16342  
    3737zstar=sqrt(2*alphaT*thickness/accumrate);
    3838
    39 tprofile=surftemp+sqrt(2*thickness*pi/accumrate/alphaT)*(-heatflux)/2/rho/c*(erf(z/zstar)-erf(thickness/zstar));
     39tprofile=surftemp+sqrt(2*thickness*pi/accumrate/alphaT)*(-heatflux)/2/rho/c.*(erf(z./zstar)-erf(thickness./zstar));
    4040
    4141% difference between surface and base temperature for check (Cuffey2010 p412):
Note: See TracChangeset for help on using the changeset viewer.