Changeset 21505 for issm/trunk-jpl/test/NightlyRun/test1501.m

Timestamp:

01/31/17 17:15:36 (8 years ago)

Author:

schlegel

Message:

CHG: update sawtooth test so that plotting works correctly

File:

• issm/trunk-jpl/test/NightlyRun/test1501.m (modified) (4 diffs)

issm/trunk-jpl/test/NightlyRun/test1501.m

@@ -108 +108 @@
         for t=starttime:endtime
                 thickness = [thickness (md.results.TransientSolution(t).Thickness)];
-                volume = [volume mean(md.results.TransientSolution(t).Thickness.value,2).*areas];
+                volume = [volume mean(md.results.TransientSolution(t).Thickness(md.mesh.elements),2).*areas];
                 massbal = [massbal (md.results.TransientSolution(t).SmbMassBalance)];
                 velocity = [velocity (md.results.TransientSolution(t).Vel)];
@@ -132 +132 @@
                 [data datatype]=processdata(md,results(i).(field),options);
 
-                titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
+                titlestring=[field ' at time ' num2str(results(i).time) ' year'];
                 plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
                 options=changefieldvalue(options,'title',titlestring);
@@ -147 +147 @@
                 [data datatype]=processdata(md,results(i).(field),options);
 
-                titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
+                titlestring=[field ' at time ' num2str(results(i).time) ' year'];
                 plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
                 options=changefieldvalue(options,'title',titlestring);
@@ -196 +196 @@
                 else
                         frame=getframe(gcf);
-                        images(:,:,1,count) = rgb2ind(frame.cdata,map,'nodither');
+                        images(:,:,1,count) = rgb2ind(frame.cdata,[map; map],'nodither');
                 end