http://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=Lisstcheck 2026-06-01T22:28:40Z Revision history for this page on the wiki MediaWiki 1.34.0 http://wiki-pnb.eri.ucsb.edu/index.php?title=Lisstcheck&diff=854&oldid=prev 2015-10-27T21:34:58Z

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

Eriks