Changeset 9436 for issm/trunk/src/m/model/modelextract.m

Timestamp:

08/25/11 16:09:15 (14 years ago)

Author:

Mathieu Morlighem

Message:

Fixed some bad fields

File:

• issm/trunk/src/m/model/modelextract.m (modified) (3 diffs)

issm/trunk/src/m/model/modelextract.m

@@ -86 +86 @@
 
                 %get field
-                field=md1.(model_fields(i));
+                field=md1.(model_fields{i});
                 fieldsize=size(field);
 
                 %size = number of nodes * n
                 if fieldsize(1)==numberofnodes1
-                        md2.(model_fields(i))=field(pos_node,:);
+                        md2.(model_fields{i})=field(pos_node,:);
 
                 %size = number of elements * n
                 elseif fieldsize(1)==numberofelements1
-                        md2.(model_fields(i))=field(pos_elem,:);
+                        md2.(model_fields{i})=field(pos_elem,:);
                 end
         end
@@ -211 +211 @@
         if ~isnan(md1.spctemperature),
                 md2.spctemperature(nodestoflag2,1)=1;
-                if ~isnan(md1.observed_temperature)
-                        md2.spctemperature(nodestoflag2,2)=md2.observed_temperature(nodestoflag2); 
-                else
-                        md2.spctemperature(nodestoflag2,2)=zeros(length(nodestoflag2),2);
-                        disp(' ')
-                        disp('!! modelextract warning: spc values should be checked !!')
-                        disp(' ')
-                end
         end
 
@@ -253 +245 @@
         end
 
-        %reinitialize output parameters
-        md2.viscousheating=NaN;
-        md2.pressure_elem=NaN;
-        md2.stress=NaN;
-        md2.stress_surface=NaN;
-        md2.stress_bed=NaN;
-        md2.deviatoricstress=NaN;
-        md2.strainrate=NaN;
-
 %Keep track of pos_node and pos_elem
 md2.extractednodes=pos_node;