source: issm/oecreview/Archive/16133-16554/ISSM-16487-16488.diff@ 16556

Index: ../trunk-jpl/src/m/plot/processmesh.m
===================================================================
--- ../trunk-jpl/src/m/plot/processmesh.m       (revision 16487)
+++ ../trunk-jpl/src/m/plot/processmesh.m       (revision 16488)
@@ -20,55 +20,52 @@
        return;
 end
 
-if ~isempty(data)
-       %first load x,y, etc ... to speed up plot
+%first load x,y, etc ... to speed up plot
+if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
+       x=md.mesh.x;
+       if isprop(md.mesh,'x2d'), x2d=md.mesh.x2d; end
+       y=md.mesh.y;
+       if isprop(md.mesh,'y2d'), y2d=md.mesh.y2d; end
+else
+       x=md.mesh.long;
+       y=md.mesh.lat;
+end
 
-       if ~strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
-               x=md.mesh.x;
-               if isprop(md.mesh,'x2d'), x2d=md.mesh.x2d; end
-               y=md.mesh.y;
-               if isprop(md.mesh,'y2d'), y2d=md.mesh.y2d; end
-       else
-               x=md.mesh.long;
-               y=md.mesh.lat;
-       end
+z_field=getfieldvalue(options,'z',md.mesh.z);
+if ischar(z_field),
+       z=md.(z_field);
+elseif isnumeric(z_field),
+       z=z_field;
+else
+       z=md.mesh.z;
+end
 
-       z_field=getfieldvalue(options,'z',md.mesh.z);
-       if ischar(z_field),
-               z=md.(z_field);
-       elseif isnumeric(z_field),
-               z=z_field;
-       else
-               z=md.mesh.z;
-       end
+if isprop(md.mesh,'elements2d'), elements2d=md.mesh.elements2d; end
+elements=md.mesh.elements;
 
-       if isprop(md.mesh,'elements2d'), elements2d=md.mesh.elements2d; end
-       elements=md.mesh.elements;
-
-       %is it a 2d plot?
-       if (md.mesh.dimension==2);
+%is it a 2d plot?
+if (md.mesh.dimension==2);
+       is2d=1;
+elseif (md.mesh.dimension==3),
+       if getfieldvalue(options,'layer',0)>=1,
                is2d=1;
-       elseif (md.mesh.dimension==3),
-               if getfieldvalue(options,'layer',0)>=1,
-                       is2d=1;
-               else
-                       is2d=0;
-               end
        else
-               error(['dim = ' num2str(md.mesh.dimension) ' not supported yet']);
+               is2d=0;
        end
+else
+       error(['dim = ' num2str(md.mesh.dimension) ' not supported yet']);
+end
 
-       %layer projection? 
-       if getfieldvalue(options,'layer',0)>=1,
-               if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
-                       error('processmesh error message: cannot work with 3D meshes for now');
-               end
-               %we modify the mesh temporarily to a 2d mesh from which the 3d mesh was extruded. 
-               x=x2d;
-               y=y2d;
-               z=zeros(size(x2d));
-               elements=elements2d;
+%layer projection? 
+if getfieldvalue(options,'layer',0)>=1,
+       if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
+               error('processmesh error message: cannot work with 3D meshes for now');
        end
+       %we modify the mesh temporarily to a 2d mesh from which the 3d mesh was extruded. 
+       x=x2d;
+       y=y2d;
+       z=zeros(size(x2d));
+       elements=elements2d;
 end
 
 %units