Lisstcheck
Jump to navigation
Jump to search
function lisstcheck(filename) %Written by Erik Stassinos 10/08/2011 %This function is designed to look at IOP fram lisst data. %file names are DH4 waped LISST names fp = fopen(filename) i=1; startbrak(i) = 1; count =1; deptharray =0; while ~feof(fp) %to get size of file databuff = fgets(fp); count = count +1; end frewind(fp) ind =1; clsind =1; for i = 1:count databuff = fgets(fp); %size(databuff) if strncmp(databuff,'{',1) % == 1 %index to find open bracket disp('found it first'); startbrack(ind) =i; ind = ind +1; end % if ind < 3 if i == (startbrack(ind-1)) + 37 %finding the depth disp('found depth') depth = databuff; depth = str2num(depth)*.01; deparr(ind) = depth; end % end if strncmp(databuff,'}',1) %== 1 %index to find close bracket endbrack(clsind) = i; disp('found close bracket') clsind = clsind +1; end end deparr; dlmwrite(strcat('lisstdepth',filename),deparr',''); %write the depth array %read just data loop %checking for indexes %endbrack %startbrack rp = fopen(filename); goodindex =1; np = fopen('lisstnobrack','w+') for rpi =1:count % loop to print without brackets num = fgets(rp) if strncmp(num,'{',1) == 1 disp('passing') continue end if strncmp(num,'}',1) == 1 continue end fprintf(np,'%s',num) lisstnumbers(rpi) = {num}; count rpi end fclose all ct =1; r=1; t=1; lisstdat = load('lisstnobrack'); %put data into parallel columns %40 columns in all while t <= length(lisstdat) holder(r,:) = lisstdat(t,:); r=r+1; if t == 40*ct %|| t ==1 cycle through 40 columns then restart to get data indexarr(ct,:) = holder'; ct = ct +1; clear holder r=1; end t = t +1; if t == length(lisstdat) dlmwrite('LISST_indexarr',indexarr,' ') return end end %save('indexarr', indexarr) return