source: issm/oecreview/Archive/12678-13393/ISSM-12959-12960.diff@ 21726

u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/utils/consistency/checkfield.m

@@ -135 +135 @@
         field2=reshape(field,prod(size(field)),1);
         if any(field2<lowerbound),
                 md = checkmessage(md,getfieldvalue(options,'message',...
-                        ['field ''' fieldname ''' should have values above ' num2str(lowerbound)]));
+                        ['field ''' fieldname ''' should have values above ' num2str(lowerbound(1,1))]));
         end
 end
 if exist(options,'>')
@@ -143 +143 @@
         field2=reshape(field,prod(size(field)),1);
         if any(field2<=lowerbound),
                 md = checkmessage(md,getfieldvalue(options,'message',...
-                        ['field ''' fieldname ''' should have values above ' num2str(lowerbound)]));
+                        ['field ''' fieldname ''' should have values above ' num2str(lowerbound(1,1))]));
         end
 end
 
@@ -153 +153 @@
         field2=reshape(field,prod(size(field)),1);
         if any(field2>upperbound),
                 md = checkmessage(md,getfieldvalue(options,'message',...
-                        ['field ''' fieldname ''' should have values below ' num2str(upperbound)]));
+                        ['field ''' fieldname ''' should have values below ' num2str(upperbound(1,1))]));
         end
 end
 if exist(options,'<')
@@ -161 +161 @@
         field2=reshape(field,prod(size(field)),1);
         if any(field2>=upperbound),
                 md = checkmessage(md,getfieldvalue(options,'message',...
-                        ['field ''' fieldname ''' should have values below ' num2str(upperbound)]));
+                        ['field ''' fieldname ''' should have values below ' num2str(upperbound(1,1))]));
         end
 end
 

u/astrid-r1b/morlighe/issmuci/trunk-jpl/../trunk-jpl/src/m/classes/mesh.m

@@ -86 +86 @@
                         md = checkfield(md,'mesh.elementonsurface','size',[md.mesh.numberofelements 1],'values',[0 1]);
                         md = checkfield(md,'mesh.vertexonbed','size',[md.mesh.numberofvertices 1],'values',[0 1]);
                         md = checkfield(md,'mesh.vertexonsurface','size',[md.mesh.numberofvertices 1],'values',[0 1]);
+                        if (md.mesh.dimension==3),
+                                md = checkfield(md,'mesh.z','>=',md.geometry.bed,'message','''mesh.z'' should be above the bed');
+                                md = checkfield(md,'mesh.z','<=',md.geometry.surface,'message','''mesh.z'' should be below the ice surface');
+                        end
                         if (md.mesh.dimension==2),
                                 md = checkfield(md,'mesh.average_vertex_connectivity','>=',9,'message','''mesh.average_vertex_connectivity'' should be at least 9 in 2d');
                         else