Changeset 1759

Timestamp:

08/18/09 15:39:43 (16 years ago)

Author:

seroussi

Message:

fixed macayeal

Location:

issm/trunk/src/m/solutions/macayeal

Files:

• control.m (modified) (1 diff)
• diagnostic.m (modified) (1 diff)

issm/trunk/src/m/solutions/macayeal/control.m

@@ -40 +40 @@
 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;
+nodes_on_dirichlet=md.spcvelocity(:,1);
 nodes_on_icefront=zeros(nods,1); nodes_on_icefront(index_icefront)=1;
 nodes_on_iceshelf=md.gridoniceshelf;

issm/trunk/src/m/solutions/macayeal/diagnostic.m

@@ -249 +249 @@
 specified_velocity=zeros(2*nods,1);
 pos=find(node_on_dirichlet);
-specified_velocity(pos)=md.spcvelocity(pos,1)/md.yts;
-specified_velocity(nods+pos)=md.spcvelocity(pos,2)/md.yts;
+specified_velocity(pos)=md.spcvelocity(pos,4)/md.yts;
+specified_velocity(nods+pos)=md.spcvelocity(pos,5)/md.yts;