function DATA=EXTDIS(s,Node,DOF)
%******************************************************************************
% EXTRACT DISPLACEMENT OF A NODE FROM NLFEA OUTPUT FILE                       %
%******************************************************************************
%%                                                                            %
 IDATA=1; fid=fopen(s,'r');                                   % Open input file
 while 1                                       % Loop for all lines in the file
  tline=fgetl(fid);                                               % Read a line
  if strfind(tline,'TIME =')                                   % Find TIME line
   DATA(1,IDATA)=sscanf(tline,'TIME = %f');                    % Read time data
   while 1                                     % Loop for lines after TIME line
    tline=fgetl(fid);                                             % Read a line
    V=sscanf(tline,'%d %f %f %f');                     % Read displacement data
    if ~isempty(V) && V(1)==Node                           % Find matching node 
     DATA(2,IDATA)=V(DOF+1);                           % Store DOF displacement
     IDATA=IDATA+1;                                       % Increase data index 
     break;                          % Exit inner loop and search for TIME line
    end                                   % End of finding matching node branch
   end                                                % End of search node loop
  end                                         % End of finding TIME line branch
  if tline==-1, break; end          % Stop the program at the end of input file
 end                                           % End of reading input file loop
end                                                                           %