Changeset 1094 for trunk/src/plot_field.m

Timestamp:

Mar 25, 2021, 9:37:04 AM (3 years ago)

Author:

sommeria

Message:

phys_polar debugged

File:

• trunk/src/plot_field.m (modified) (4 diffs)

trunk/src/plot_field.m

@@ -768 +768 @@
             test_interp_X=0; %default, regularly meshed X coordinate
             test_interp_Y=0; %default, regularly meshed Y coordinate
-            if isfield(Data,'VarAttribute')
-                if numel(Data.VarAttribute)>=CellInfo{icell}.CoordIndex(end) && isfield(Data.VarAttribute{CellInfo{icell}.CoordIndex(end)},'units')
-                    x_units=Data.VarAttribute{CellInfo{icell}.CoordIndex(end)}.units;
-                end
-                if numel(Data.VarAttribute)>=CellInfo{icell}.CoordIndex(end-1) && isfield(Data.VarAttribute{CellInfo{icell}.CoordIndex(end-1)},'units')
-                    y_units=Data.VarAttribute{CellInfo{icell}.CoordIndex(end-1)}.units;
-                end
-            end
+%             if isfield(Data,'VarAttribute')
+%                 if numel(Data.VarAttribute)>=CellInfo{icell}.CoordIndex(end) && isfield(Data.VarAttribute{CellInfo{icell}.CoordIndex(end)},'units')
+%                     x_units=Data.VarAttribute{CellInfo{icell}.CoordIndex(end)}.units;
+%                 end
+%                 if numel(Data.VarAttribute)>=CellInfo{icell}.CoordIndex(end-1) && isfield(Data.VarAttribute{CellInfo{icell}.CoordIndex(end-1)},'units')
+%                     y_units=Data.VarAttribute{CellInfo{icell}.CoordIndex(end-1)}.units;
+%                 end
+%             end
             if numel(Coord_y)>2
                 DCoord_y=diff(Coord_y);
@@ -790 +790 @@
                 DCoord_x_min=min(DCoord_x);
                 DCoord_x_max=max(DCoord_x);
-                if sign(DCoord_x_min)~=sign(DCoord_x_max);% =1 for increasing values, 0 otherwise
+                if sign(DCoord_x_min)~=sign(DCoord_x_max)% =1 for increasing values, 0 otherwise
                     errormsg=['errror in plot_field.m: non monotonic dimension variable ' Data.ListVarName{VarRole.coord(2)} ];
                     return
@@ -821 +821 @@
     end
     %define coordinates as CoordUnits, if not defined as attribute for each variable
+%     if isfield(Data,'VarAttribute')&& numel(Data.VarAttribute)>=1 && isfield(Data.VarAttribute{1},'unit')
+%         y_units=Data.VarAttribute{1}.unit;
+%     end
     if isfield(Data,'CoordUnit')
         if isempty(x_units)
@@ -827 +830 @@
         if isempty(y_units)
             y_units=Data.CoordUnit;
+        end
+    elseif isfield(Data,'VarAttribute')
+        if numel(Data.VarAttribute)>=CellInfo{icell}.CoordIndex(end) && isfield(Data.VarAttribute{CellInfo{icell}.CoordIndex(end)},'units')
+            x_units=Data.VarAttribute{CellInfo{icell}.CoordIndex(end)}.units;
+        end
+        if numel(Data.VarAttribute)>=CellInfo{icell}.CoordIndex(end-1) && isfield(Data.VarAttribute{CellInfo{icell}.CoordIndex(end-1)},'units')
+            y_units=Data.VarAttribute{CellInfo{icell}.CoordIndex(end-1)}.units;
         end
     end