Index: /issm/trunk-jpl/src/m/contrib/hack/sia.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/hack/sia.m	(revision 19120)
+++ /issm/trunk-jpl/src/m/contrib/hack/sia.m	(revision 19121)
@@ -1,4 +1,4 @@
 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,5 +18,5 @@
 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);