Index: /issm/trunk/src/m/solutions/macayeal/control.m
===================================================================
--- /issm/trunk/src/m/solutions/macayeal/control.m	(revision 1752)
+++ /issm/trunk/src/m/solutions/macayeal/control.m	(revision 1753)
@@ -38,5 +38,5 @@
 rho_water=md.rho_water;
 g=md.g;
-index_icefront=md.segmentonneumann_diag; index_icefront=index_icefront(:,1:2); %we strip the last column, which holds the element number for the boundary segment
+index_icefront=md.pressureload; index_icefront=index_icefront(:,1:2); %we strip the last column, which holds the element number for the boundary segment
 nodes_on_boundary=md.gridonboundary;
 nodes_on_dirichlet=md.gridondirichlet_diag;
Index: /issm/trunk/src/m/solutions/macayeal/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/macayeal/diagnostic.m	(revision 1752)
+++ /issm/trunk/src/m/solutions/macayeal/diagnostic.m	(revision 1753)
@@ -36,8 +36,8 @@
 g=md.g;
 viscosity_overshoot=md.viscosity_overshoot;
-index_icefront=md.segmentonneumann_diag; index_icefront=index_icefront(:,1:2); %we strip the last column, which holds the element number for the boundary segment
+index_icefront=md.pressureload; index_icefront=index_icefront(:,1:2); %we strip the last column, which holds the element number for the boundary segment
 nodes_on_boundary=md.gridonboundary;
 nodes_on_icefront=zeros(nods,1); nodes_on_icefront(index_icefront)=1;
-node_on_dirichlet=md.gridondirichlet_diag;
+node_on_dirichlet=md.spcvelocity(:,1);
 nodes_on_icesheet=md.gridonicesheet;
 element_on_icesheet=md.elementonicesheet;
@@ -249,6 +249,6 @@
 specified_velocity=zeros(2*nods,1);
 pos=find(node_on_dirichlet);
-specified_velocity(pos)=md.dirichletvalues_diag(pos,1)/md.yts;
-specified_velocity(nods+pos)=md.dirichletvalues_diag(pos,2)/md.yts;
+specified_velocity(pos)=md.spcvelocity(pos,1)/md.yts;
+specified_velocity(nods+pos)=md.spcvelocity(pos,2)/md.yts;