Changeset 18255

Timestamp:

07/17/14 08:26:59 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: can just call plotmodel(md,'data','transient_movie') now

File:

• issm/trunk-jpl/src/m/plot/plot_transient_movie.m (modified) (3 diffs)

issm/trunk-jpl/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
 
@@ -33 +48 @@
                         [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 +65 @@
                                 end
                         else
-                                pause(0.1)
+                                pause(0.5)
                         end
                         nstep=nstep+1;