% function [z,sigt,t,s,nn,nnv,mesh,header]=read_s3r3(file) function [z,sigt,t,s,nn,nnv,mesh,header]=read_s3r3(file) [pfid,message]=fopen([file]); datatype=file(length(file)-3:length(file)); mesh=fgets(pfid);mesh=blank(mesh(1:length(mesh)-1)); header=fgets(pfid); nnv=fscanf(pfid,'%f',1); data=fscanf(pfid,'%f',[1 inf])'; fclose(pfid); size(data); nn=ans(1)/nnv/5; string=[num2str(nn),'x',num2str(nnv),' ',datatype,' values read from ',file]; fprintf(1,string); fprintf(1,'\n'); data=reshape(data,[5 nn*nnv])'; z=data(:,2);sigt=data(:,3);t=data(:,4);s=data(:,5);