Changeset 1067

Timestamp:

06/23/09 14:38:48 (15 years ago)

Author:

seroussi

Message:

no more m_dh.gridset because of hutter

File:

• issm/trunk/src/m/solutions/ice/transient3d.m (modified) (2 diffs)

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

@@ -35 +35 @@
 
 %initialize (velocity,pressure,...)
-solution.u_g=zeros(fem.m_dh.gridset.gsize,1);
-if ~isnan(md.vx), solution.u_g(1:6:fem.m_dh.gridset.gsize)=md.vx/md.yts; end
-if ~isnan(md.vy), solution.u_g(2:6:fem.m_dh.gridset.gsize)=md.vy/md.yts; end
-if ~isnan(md.vz), solution.u_g(3:6:fem.m_dh.gridset.gsize)=md.vz/md.yts; end
-solution.pressure=zeros(fem.m_dh.gridset.gsize,1); 
+solution.u_g=zeros(fem.m_dv.gridset.gsize,1);
+if ~isnan(md.vx), solution.u_g(1:6:fem.m_dv.gridset.gsize)=md.vx/md.yts; end
+if ~isnan(md.vy), solution.u_g(2:6:fem.m_dv.gridset.gsize)=md.vy/md.yts; end
+if ~isnan(md.vz), solution.u_g(3:6:fem.m_dv.gridset.gsize)=md.vz/md.yts; end
+solution.pressure=zeros(fem.m_dv.gridset.gsize,1); 
 if ~isnan(md.pressure), solution.pressure(1:6:end)=md.pressure; end
 solution.t_g=zeros(fem.m_t.gridset.gsize,1);
@@ -53 +53 @@
 %initialize geometry and time
 solution.mesh=struct('elements',md.elements,'grids',[md.x md.y md.z]);
-solution.thickness=zeros(fem.m_dh.gridset.gsize,1);solution.thickness(1:6:fem.m_dh.gridset.gsize)=md.thickness;
-solution.surface=zeros(fem.m_dh.gridset.gsize,1);solution.surface(1:6:fem.m_dh.gridset.gsize)=md.surface;
-solution.bed=zeros(fem.m_dh.gridset.gsize,1);solution.bed(1:6:fem.m_dh.gridset.gsize)=md.bed;
+solution.thickness=zeros(fem.m_p.gridset.gsize,1);solution.thickness(1:6:fem.m_p.gridset.gsize)=md.thickness;
+solution.surface=zeros(fem.m_p.gridset.gsize,1);solution.surface(1:6:fem.m_p.gridset.gsize)=md.surface;
+solution.bed=zeros(fem.m_p.gridset.gsize,1);solution.bed(1:6:fem.m_p.gridset.gsize)=md.bed;
 solution.time=0;