Changeset 19121

Timestamp:

02/19/15 08:38:57 (10 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed SIA module with new md structure

File:

• issm/trunk-jpl/src/m/contrib/hack/sia.m (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/hack/sia.m

@@ -1 +1 @@
 function [velx,vely,vel]=sia(md)
-%BALVEL - computation of Shallow Ice velocities
+%SIA - computation of Shallow Ice velocities
 %
 %   This routine uses the model of SIA to compute the velocities
@@ -18 +18 @@
 summer=[1;1;1];
 hel=md.geometry.thickness(md.mesh.elements)*summer/3;
-Bel=md.B(md.mesh.elements)*summer/3;
+Bel=md.materials.rheology_B(md.mesh.elements)*summer/3;
 
 Ael=Bel.^(-3);