Changeset 880 for trunk/src/proj_field.m

Timestamp:

Mar 4, 2015, 12:01:38 AM (9 years ago)

Author:

sommeria

Message:

various bug fixes

File:

• trunk/src/proj_field.m (modified) (2 diffs)

trunk/src/proj_field.m

@@ -509 +509 @@
         VarSize(nbvar)=mean((ProjData.([VarName 'Max'])-ProjData.([VarName 'Min']))/100);
     end
-    if  isempty(VarMesh) || isnan(VarMesh) % mesh not specified as input, estimate from the bounds
+    if  isempty(VarMesh)% || isnan(VarMesh) % mesh not specified as input, estimate from the bounds
         VarMesh=mean(VarSize);
         ord=10^(floor(log10(VarMesh)));%order of magnitude
@@ -524 +524 @@
         LowBound=VarMesh*ceil(ProjData.([VarName 'Min'])/VarMesh);
         UpperBound=VarMesh*floor(ProjData.([VarName 'Max'])/VarMesh);
+        if numel(indsel)<=1
+            errormsg='only one data point or less for histogram';
+            return
+        elseif isequal(LowBound,UpperBound)
+            errormsg='attempt histogram of uniform field: low bound = high bound';
+            return
+        end       
         ProjData.(VarName)=LowBound:VarMesh:UpperBound; % list of bin values
         ProjData.([VarName 'Histo'])=hist(double(FieldData.(VarName)(indsel,:)),ProjData.(VarName)); % histogram at predefined bin positions