Ignore:
Timestamp:
08/25/11 16:09:15 (14 years ago)
Author:
Mathieu Morlighem
Message:

Fixed some bad fields

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/model/modelextract.m

    r8833 r9436  
    8686
    8787                %get field
    88                 field=md1.(model_fields(i));
     88                field=md1.(model_fields{i});
    8989                fieldsize=size(field);
    9090
    9191                %size = number of nodes * n
    9292                if fieldsize(1)==numberofnodes1
    93                         md2.(model_fields(i))=field(pos_node,:);
     93                        md2.(model_fields{i})=field(pos_node,:);
    9494
    9595                %size = number of elements * n
    9696                elseif fieldsize(1)==numberofelements1
    97                         md2.(model_fields(i))=field(pos_elem,:);
     97                        md2.(model_fields{i})=field(pos_elem,:);
    9898                end
    9999        end
     
    211211        if ~isnan(md1.spctemperature),
    212212                md2.spctemperature(nodestoflag2,1)=1;
    213                 if ~isnan(md1.observed_temperature)
    214                         md2.spctemperature(nodestoflag2,2)=md2.observed_temperature(nodestoflag2);
    215                 else
    216                         md2.spctemperature(nodestoflag2,2)=zeros(length(nodestoflag2),2);
    217                         disp(' ')
    218                         disp('!! modelextract warning: spc values should be checked !!')
    219                         disp(' ')
    220                 end
    221213        end
    222214
     
    253245        end
    254246
    255         %reinitialize output parameters
    256         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 
    264247%Keep track of pos_node and pos_elem
    265248md2.extractednodes=pos_node;
Note: See TracChangeset for help on using the changeset viewer.