Changeset 6815

Timestamp:

12/27/10 09:56:30 (14 years ago)

Author:

Mathieu Morlighem

Message:

Fixed plot_transient_movie.m

File:

• issm/trunk/src/m/model/plot/plot_transient_movie.m (modified) (2 diffs)

issm/trunk/src/m/model/plot/plot_transient_movie.m

@@ -16 +16 @@
         end
 
-        %process model
-        [x y z elements is2d]=processmesh(md,[],options);
+        if md.dim==2,
+                results=md.results.Transient2DSolution;
+        else
+                results=md.results.Transient3DSolution;
+        end
 
         %loop over the time steps
-        for i=1:length(md.results.transient)
-                eval(['data=md.results.transient(' num2str(i) ').' num2str(field) ';']);
+        for i=1:length(results)
 
                 %process data
-                [data datatype]=processdata(md,data,options);
-                titlestring=[field ' at time ' num2str(md.results.transient(i).time) ' year'];
+                [x y z elements is2d]=processmesh(md,results(i).(field),options);
+                [data datatype]=processdata(md,results(i).(field),options);
+
+                titlestring=[field ' at time ' num2str(results(i).time/md.yts) ' year'];
                 plot_unit(x,y,z,elements,data,is2d,datatype,options)
                 apply_options_movie(md,options,titlestring);
@@ -35 +39 @@
                                 frame=getframe(gcf);
                                 [images,map]=rgb2ind(frame.cdata,256,'nodither');
-                                images(1,1,1,length(md.results.transient))=0;
+                                images(1,1,1,length(results))=0;
                         else
                                 frame=getframe(gcf);