Changeset 460

Timestamp:

05/15/09 17:13:22 (16 years ago)

Author:

Mathieu Morlighem

Message:

take care of 3d runs

File:

• issm/trunk/src/m/classes/public/loadresultsfromdisk.m (modified) (1 diff)

issm/trunk/src/m/classes/public/loadresultsfromdisk.m

@@ -30 +30 @@
         u_g=results{4};
 
-        %Used to recover velocities
-        indx=1:2:gsize;
-        indy=2:2:gsize;
-        indx=indx(part);
-        indy=indy(part);
+        if strcmpi(md.type,'2d'),
+                %Used to recover velocities
+                indx=1:2:gsize;
+                indy=2:2:gsize;
+                indx=indx(part);
+                indy=indy(part);
 
-        %Recover velocity
-        md.vx=u_g(indx)*md.yts;
-        md.vy=u_g(indy)*md.yts;
-        md.vel=sqrt(md.vx.^2+md.vy.^2);
+                %Recover velocity
+                md.vx=u_g(indx)*md.yts;
+                md.vy=u_g(indy)*md.yts;
+                md.vel=sqrt(md.vx.^2+md.vy.^2);
+        else
+                %Used to recover velocities
+                indx=1:3:gsize;
+                indy=2:3:gsize;
+                indz=3:3:gsize;
+                indx=indx(part);
+                indy=indy(part);
+                indz=indz(part);
+
+                %Recover velocity
+                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);
+        end
 
 elseif strcmpi(analysis_type,'control'),