Changeset 24134

Timestamp:

09/06/19 04:06:01 (6 years ago)

Author:

rueckamp

Message:

CHG: extract function now keeps all TransientSolution time steps

File:

• issm/trunk-jpl/src/m/classes/model.m (modified) (1 diff)

issm/trunk-jpl/src/m/classes/model.m

@@ -653 +653 @@
                                         if isstruct(md1.results.(solutionfields{i}))
                                                 %get subfields
-                                                solutionsubfields=fields(md1.results.(solutionfields{i}));
-                                                for j=1:length(solutionsubfields),
-                                                        field=md1.results.(solutionfields{i}).(solutionsubfields{j});
+                                                % loop over time steps
+                                                for p=1:length(md1.results.(solutionfields{i}))
+                                                    current = md1.results.(solutionfields{i})(p);
+                                                    solutionsubfields=fields(current);
+                                                    for j=1:length(solutionsubfields),
+                                                        field=md1.results.(solutionfields{i})(p).(solutionsubfields{j});
                                                         if length(field)==numberofvertices1,
-                                                                md2.results.(solutionfields{i}).(solutionsubfields{j})=field(pos_node);
+                                                            md2.results.(solutionfields{i})(p).(solutionsubfields{j})=field(pos_node);
                                                         elseif length(field)==numberofelements1,
-                                                                md2.results.(solutionfields{i}).(solutionsubfields{j})=field(pos_elem);
+                                                            md2.results.(solutionfields{i})(p).(solutionsubfields{j})=field(pos_elem);
                                                         else
-                                                                md2.results.(solutionfields{i}).(solutionsubfields{j})=field;
+                                                            md2.results.(solutionfields{i})(p).(solutionsubfields{j})=field;
                                                         end
+                                                    end
                                                 end
                                         else