Changeset 9436 for issm/trunk/src/m/model/modelextract.m
- Timestamp:
- 08/25/11 16:09:15 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/modelextract.m
r8833 r9436 86 86 87 87 %get field 88 field=md1.(model_fields (i));88 field=md1.(model_fields{i}); 89 89 fieldsize=size(field); 90 90 91 91 %size = number of nodes * n 92 92 if fieldsize(1)==numberofnodes1 93 md2.(model_fields (i))=field(pos_node,:);93 md2.(model_fields{i})=field(pos_node,:); 94 94 95 95 %size = number of elements * n 96 96 elseif fieldsize(1)==numberofelements1 97 md2.(model_fields (i))=field(pos_elem,:);97 md2.(model_fields{i})=field(pos_elem,:); 98 98 end 99 99 end … … 211 211 if ~isnan(md1.spctemperature), 212 212 md2.spctemperature(nodestoflag2,1)=1; 213 if ~isnan(md1.observed_temperature)214 md2.spctemperature(nodestoflag2,2)=md2.observed_temperature(nodestoflag2);215 else216 md2.spctemperature(nodestoflag2,2)=zeros(length(nodestoflag2),2);217 disp(' ')218 disp('!! modelextract warning: spc values should be checked !!')219 disp(' ')220 end221 213 end 222 214 … … 253 245 end 254 246 255 %reinitialize output parameters256 md2.viscousheating=NaN;257 md2.pressure_elem=NaN;258 md2.stress=NaN;259 md2.stress_surface=NaN;260 md2.stress_bed=NaN;261 md2.deviatoricstress=NaN;262 md2.strainrate=NaN;263 264 247 %Keep track of pos_node and pos_elem 265 248 md2.extractednodes=pos_node;
Note:
See TracChangeset
for help on using the changeset viewer.