Changeset 990

Timestamp:

06/13/09 20:02:40 (16 years ago)

Author:

Mathieu Morlighem

Message:

grid position was not updated in ice

File:

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

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

@@ -134 +134 @@
         %update grids
         disp(sprintf('%s','   updating grid positions...'));
-        grids_t=UpdateGridPosition(md,fem.m_t.grids,new_bed,new_thickness);
-        grids_dh=UpdateGridPosition(md,fem.m_dh.grids,new_bed,new_thickness);
-        grids_dv=UpdateGridPosition(md,fem.m_dv.grids,new_bed,new_thickness);
-        grids_p=UpdateGridPosition(md,fem.m_p.grids,new_bed,new_thickness);
+        if fem.ismacayealpattyn,
+                fem.m_dh.grids=UpdateGridPosition(md,fem.m_dh.grids,new_bed,new_thickness);
+        end
+        if fem.ishutter,
+                fem.m_dhu.grids=UpdateGridPosition(md,fem.m_dhu.grids,new_bed,new_thickness);
+                fem.m_ss.grids=UpdateGridPosition(md,fem.m_ss.grids,new_bed,new_thickness);
+        end
+        if fem.isstokes,
+                fem.m_bs.grids=UpdateGridPosition(md,fem.m_bs.grids,new_bed,new_thickness);
+                fem.m_ds.grids=UpdateGridPosition(md,fem.m_ds.grids,new_bed,new_thickness);
+        end
+        fem.m_dv.grids=UpdateGridPosition(md,fem.m_dv.grids,new_bed,new_thickness);
+        fem.m_p.grids=UpdateGridPosition(md,fem.m_p.grids,new_bed,new_thickness);
+        fem.m_t.grids=UpdateGridPosition(md,fem.m_t.grids,new_bed,new_thickness);
+        fem.m_m.grids=UpdateGridPosition(md,fem.m_m.grids,new_bed,new_thickness);
 
         %update time and counter