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

../trunk-jpl/src/m/plot/processmesh.m

@@ -20 +20 @@
         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