Changeset 18301 for issm/trunk/src/m/plot/plot_transient_movie.m

Timestamp:

07/30/14 14:22:43 (11 years ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 18299

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• src (modified) (1 prop)
• src/m/plot/plot_transient_movie.m (modified) (4 diffs)

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

@@ -13 +13 @@
                 field=getfieldvalue(options,'transient_movie_field');
         else
-                error('specify transient_movie_field in options list');
+                disp('List of available fields:');
+                F=fields(md.results.TransientSolution(1));
+                num = [];
+                for i=1:numel(F),
+                        if ~strcmp(F{i},'time') & ...
+                                ~strcmp(F{i},'step') & ...
+                                ~strcmp(F{i},'errlog') & ...
+                                ~strcmp(F{i},'outlog') & ...
+                                ~strcmp(F{i},'MaxIterationConvergenceFlag') & ...
+                                ~strcmp(F{i},'SolutionType'),
+                                disp(['   ' num2str(i) ': ' F{i} ]);
+                                num = [num i];
+                        end
+                end
+                choice=input(['please enter the field number? (between ' num2str(min(num)) ' and ' num2str(max(num)) ')  ']);
+                field =  F{choice};
         end
 
@@ -24 +39 @@
                 steps=1:length(results);
         end
-        results=md.results.TransientSolution;
+
+        %calculate caxis
+        if ~exist(options,'caxis'),
+                range = [Inf -Inf];
+                for i=steps
+                        [data datatype]=processdata(md,results(i).(field),options);
+                        range(1) = min(range(1),min(data));
+                        range(2) = max(range(2),max(data));
+                end
+                options=addfielddefault(options,'caxis',range);
+        end
+
+        %display movie
         nstep=1;
+        deltat = getfieldvalue(options,'pause',.5);
         for i=steps
 
@@ -33 +61 @@
                         [data datatype]=processdata(md,results(i).(field),options);
 
+                        clf;
                         titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
                         plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
@@ -49 +78 @@
                                 end
                         else
-                                pause(0.1)
+                                pause(deltat)
                         end
                         nstep=nstep+1;