Changeset 6134

Timestamp:

10/04/10 08:47:19 (14 years ago)

Author:

Eric.Larour

Message:

some improvements to transient plot, minor

Location:

issm/trunk/src/m/model/plot

Files:

• plot_manager.m (modified) (1 diff)
• plot_transient_field.m (modified) (3 diffs)
• plotdoc.m (modified) (1 diff)

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

@@ -118 +118 @@
                 case 'transient_results',
                         plot_transient_results(md,options,subplotwidth,i);
-                case {'transient_thickness','transient_bed','transient_surface','transient_temperature','transient_melting','transient_vel','transient_vx','transient_vy','transient_vz','transient_pressure'}
-                        plot_transient_field(md,options,subplotwidth,i,data);
+
+                case 'transient_field',
+                        plot_transient_field(md,options,subplotwidth,i);
                         return;
 

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

@@ -8 +8 @@
 
 %Check that they are transient results
-if ~isfield(md.results,'transient'),
+if (~isfield(md.results,'Transient2DSolution') && ~isfield(md.results,'Transient3DSolution')),
         error('plot_transient_field error message: no transient results in the model');
 end
 
+if md.dim==2,
+        transient=md.results.Transient2DSolution;
+else
+        transient=md.results.Transient3DSolution;
+end
+
 %Figure out the iterations to plot and check if it is possible
-steps=getfieldvalue(options,'steps');
-maxiteration=size(md.results.transient,2);
+maxiteration=size(transient,2);
+steps=getfieldvalue(options,'steps',1:1:maxiteration);
+
 if max(steps)>maxiteration | min(steps)<1,
         error(['plot_transient_field error message: problem with the steps requested, must be an interger between 0 and ' num2str(maxiteration)]);
@@ -21 +28 @@
 
 %Figure out the field to plot
-if strncmpi(fliplr(data),fliplr('vx'),2),
-        field='vx';
-elseif strncmpi(fliplr(data),fliplr('vy'),2),
-        field='vy';
-elseif strncmpi(fliplr(data),fliplr('vz'),2),
-        field='vz';
-elseif strncmpi(fliplr(data),fliplr('vel'),3),
-        field='vel';
-elseif strncmpi(fliplr(data),fliplr('pressure'),8),
-        field='pressure';
-elseif strncmpi(fliplr(data),fliplr('melting'),6),
-        field='melting';
-elseif strncmpi(fliplr(data),fliplr('temperature'),11),
-        field='temperature';
-elseif strncmpi(fliplr(data),fliplr('bed'),3),
-        field='bed';
-elseif strncmpi(fliplr(data),fliplr('thickness'),9),
-        field='thickness';
-elseif strncmpi(fliplr(data),fliplr('surface'),7),
-        field='surface';
-end
+
+field=getfieldvalue(options,'field','Vel');
 
 %process mes only once
@@ -50 +38 @@
 
         %process data and change title if needed
-        [data datatype]=processdata(md,md.results.transient(steps(i)).(field),options);
-        options=changefieldvalue(options,'title',[field ' at time ' num2str(md.results.transient(steps(i)).time) ' a']);
+        [data datatype]=processdata(md,transient(steps(i)).(field),options);
+        options=changefieldvalue(options,'title',[field ' at time ' num2str(transient(steps(i)).time/md.yts) ' a']);
 
         %create plot of step i

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

@@ -40 +40 @@
 disp('                  - ''transient_vel'': this will display the velocity for the time steps requested in ''steps'' of a transient run');
 disp('                  - ''transient_vel'': vel can be by any field of the transient results (vx, vy, vz, vel, temperature, melting, pressure, bed, thickness, surface)');
+disp('                  - ''transient_field'': dynamic plot of results. specify ''steps'' option, as fell as ''field'' (defaults are all steps, for ''Vel'' field)');
 disp('                  - ''transient_movie'': this will display the time steps of a given field of a transient run');
 disp('                  - ''transient_movie_field'': field to be displayed when doing  transient_movie data display');