Changeset 674 for issm/trunk/src/m/solutions/ice/Loadresults.m

Timestamp:

06/01/09 14:48:45 (16 years ago)

Author:

Mathieu Morlighem

Message:

results are now in md.results

File:

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

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

@@ -27 +27 @@
         end
 
-        md.vx=u_g(indx)*md.yts;
-        md.vy=u_g(indy)*md.yts;
-        md.vel=sqrt(md.vx.^2+md.vy.^2);
-        md.pressure=u_g(indp);
+        md.results.diagnostic.step=1;
+        md.results.diagnostic.time=0;
+        md.results.diagnostic.vx=u_g(indx)*md.yts;
+        md.results.diagnostic.vy=u_g(indy)*md.yts;
+        md.results.diagnostic.vel=sqrt(md.vx.^2+md.vy.^2);
+        md.results.diagnostic.pressure=u_g(indp);
 
 else
+
         m_dv=fem.m_dv;
         indx=1:6:m_dv.gridset.gsize;
@@ -43 +46 @@
         indp=indp(m_dv.tpart);
 
-        md.vx=u_g(indx)*md.yts;
-        md.vy=u_g(indy)*md.yts;
-        md.vz=u_g(indz)*md.yts;
-        md.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
-        md.pressure=u_g(indp);
+        md.results.diagnostic.step=1;
+        md.results.diagnostic.time=0;
+        md.results.diagnostic.vx=u_g(indx)*md.yts;
+        md.results.diagnostic.vy=u_g(indy)*md.yts;
+        md.results.diagnostic.vz=u_g(indz)*md.yts;
+        md.results.diagnostic.vel=sqrt(md.vx.^2+md.vy.^2+md.vz.^2);
+        md.results.diagnostic.pressure=u_g(indp);
 end