--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/+multiblock/StaggeredSurface.m	Thu Mar 10 16:54:26 2022 +0100
@@ -0,0 +1,97 @@
+classdef StaggeredSurface < handle
+    properties
+        grid
+        surfs
+        gridGroup
+        subGrid
+
+        ZData
+        CData
+
+    end
+
+    methods
+        function obj = StaggeredSurface(g, gf, gridGroup, subGrid)
+
+            default_arg('gridGroup', 1);
+            default_arg('subGrid', 1);
+
+            obj.grid = g;
+            obj.gridGroup = gridGroup;
+            obj.subGrid = subGrid;
+
+            % coords = obj.grid.points();
+            % X = obj.grid.funcToPlotMatrices(coords(:,1));
+            % Y = obj.grid.funcToPlotMatrices(coords(:,2));
+            % V = obj.grid.funcToPlotMatrices(gf);
+            X = {};
+            Y = {};
+            V = {};
+
+            holdState = ishold();
+            hold on
+
+            surfs = cell(1, obj.grid.nBlocks);
+            gfIndex = 1;
+            for i = 1:g.nBlocks()
+
+                gi = g.grids{i}.gridGroups{gridGroup}{subGrid};
+
+                X{i} = grid.funcToPlotMatrix(gi, gi.coords(:,1));
+                Y{i} = grid.funcToPlotMatrix(gi, gi.coords(:,2));
+
+                Ni = gi.N();
+                gf_i = gf(gfIndex:gfIndex+Ni-1);
+                V{i} = grid.funcToPlotMatrix(gi, gf_i);
+
+                surfs{i} = surf(X{i}, Y{i}, V{i});
+                gfIndex = gfIndex + Ni;
+            end
+
+            if holdState == false
+                hold off
+            end
+
+            obj.surfs = [surfs{:}];
+
+            obj.ZData = gf;
+            obj.CData = gf;
+        end
+
+        function set(obj, propertyName, propertyValue)
+            set(obj.surfs, propertyName, propertyValue);
+        end
+
+        function obj = set.ZData(obj, gf)
+            obj.ZData = gf;
+
+            % V = obj.grid.funcToPlotMatrices(gf);
+            gfIndex = 1;
+            for i = 1:obj.grid.nBlocks()
+                gi = obj.grid.grids{i}.gridGroups{obj.gridGroup}{obj.subGrid};
+                Ni = gi.N();
+                gf_i = gf(gfIndex:gfIndex+Ni-1);
+                Vi = grid.funcToPlotMatrix(gi, gf_i);
+                obj.surfs(i).ZData = Vi;
+
+                gfIndex = gfIndex + Ni;
+            end
+        end
+
+        function obj = set.CData(obj, gf)
+            obj.CData = gf;
+
+            % V = obj.grid.funcToPlotMatrices(gf);
+            gfIndex = 1;
+            for i = 1:obj.grid.nBlocks()
+                gi = obj.grid.grids{i}.gridGroups{obj.gridGroup}{obj.subGrid};
+                Ni = gi.N();
+                gf_i = gf(gfIndex:gfIndex+Ni-1);
+                Vi = grid.funcToPlotMatrix(gi, gf_i);
+                obj.surfs(i).CData = Vi;
+
+                gfIndex = gfIndex + Ni;
+            end
+        end
+    end
+end