Changeset 26901 for issm/trunk-jpl/src/m/contrib/defleurian/netCDF/restable.m

Timestamp:

02/25/22 00:11:37 (3 years ago)

Author:

bdef

Message:

CHG: update to nc export and checking procedure in runme

File:

• issm/trunk-jpl/src/m/contrib/defleurian/netCDF/restable.m (modified) (3 diffs)

issm/trunk-jpl/src/m/contrib/defleurian/netCDF/restable.m

@@ -15 +15 @@
                                 %we save the size of the current step for further treatment
                                 self.sizes=[self.sizes;fliplr(size(stepvar))];
-                                flatdata=reshape(stepvar, 1, []);
+                                %we need to transpose to follow the indexing
+                                flatdata=reshape(stepvar', 1, []);
                                 self.data=[self.data,flatdata];
                         end
@@ -32 +33 @@
                                         findim=[maxsize, rows];
                                         %first check if all steps are the same size
-                                        SameSize = sum(self.sizes - self.sizes(1))==0;
+                                        SameSize = sum(self.sizes - self.sizes(1, :))==0;
                                         if SameSize,
                                                 %same size for all steps, just reshape
@@ -40 +41 @@
                                                 startpoint=1;
                                                 datadim=ndims(self.data);
-                                                outdat=NaN*ones(findim);
+                                                outdat=nan(findim);
                                                 for r=1:rows
                                                         curlen = prod(self.sizes(r, :));
-                                                        outdat(1:self.sizes(r,1), 1:self.sizes(r,2), r) = reshape(self.data(startpoint:startpoint+curlen-1),self.sizes(r,:));
+                                                        outdat(1:self.sizes(r,1), 1:self.sizes(r,2), r) = reshape(self.data(startpoint:startpoint+ ...
+                                                                                                          curlen-1),self.sizes(r,:));
                                                         startpoint = startpoint+curlen;
                                                 end