Changeset 32 for issm/trunk/src/m/solutions/ice/ModelProcessorPrognostic.m

Timestamp:

04/24/09 10:38:58 (16 years ago)

Author:

Mathieu Morlighem

Message:

initial input from ice1

File:

• issm/trunk/src/m/solutions/ice/ModelProcessorPrognostic.m (modified) (1 diff)

issm/trunk/src/m/solutions/ice/ModelProcessorPrognostic.m

@@ -116 +116 @@
 loads=struct('load',cell([0,1]));
 
-%Single point constraints: none;
-constraints=struct('constraint',cell(0,0));
+
+%Initialize constraints structure
+constraints=struct('constraint',cell(length(find(md.gridondirichlet_prog)),1));
+
+pos=find(md.gridondirichlet_prog);
+count=1;
+for i=1:length(find(md.gridondirichlet_prog)),
+
+                %constrain the thickness
+                constraints(count).constraint=spc;
+                constraints(count).constraint.grid=pos(i);
+                constraints(count).constraint.dof=1;
+                constraints(count).constraint.value=md.dirichletvalues_prog(pos(i)); %in m
+                count=count+1;
+end
 
 %Last thing, return a partitioning vector, and its transpose.