https://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=Hs6_split.m 2026-05-24T12:20:02Z Revision history for this page on the wiki MediaWiki 1.34.0 https://wiki-pnb.eri.ucsb.edu/index.php?title=Hs6_split.m&diff=185&oldid=prev 2010-08-03T00:01:39Z

<p>Created page with &#039;&lt;pre&gt; function hs6_split() %Created by Erik Stassinos 5/26/2010 %This function takes the input single lcd and then slits it into up and %downcast matricies line =0; h1=0; id = …&#039;</p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> function hs6_split()<br /> %Created by Erik Stassinos 5/26/2010<br /> %This function takes the input single lcd and then slits it into up and<br /> %downcast matricies<br /> <br /> <br /> line =0;<br /> h1=0;<br /> id = fopen('list_hs6')<br /> if id == -1<br /> disp('cannot find hs6 list file')<br /> return<br /> end<br /> <br /> <br /> while feof(id) ~=1<br /> tu = fgets(id);<br /> h1 = h1+1;<br /> end<br /> <br /> frewind(id)<br /> for i = 1:h1%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br /> %clear h1;<br /> line =0;<br /> r = fgets(id);<br /> length(r);<br /> t = textscan(r,'%s');<br /> p = char(t{1,:})<br /> [in,mer,dep] = index1(p); <br /> in{1};<br /> in{2};<br /> in{3};<br /> in{4};<br /> <br /> m = strrep(p,'.lcd','.datmat');<br /> mat = load(m,'.ASC');<br /> width = length(mat(1,:));<br /> <br /> %format long<br /> clear ij<br /> clear downcastmatrix<br /> clear upcastmatrix<br /> <br /> for ij = in{1}:1:in{2};<br /> downcastmatrix(ij-(in{1}-1),:) = mat(ij,1:width); %starts from 1 and goes to in2<br /> end<br /> for r = in{3}:in{4};<br /> upcastmatrix(r-(in{3}-1),:) = mat(r,1:width);<br /> end<br /> dlmwrite(strcat(m,'.1'),downcastmatrix,'delimiter',' ','precision',6);<br /> dlmwrite(strcat(m,'.2'),upcastmatrix,'delimiter',' ','precision',6);<br /> <br /> <br /> name1 = strrep(p,'.lcd','.lcd.1')<br /> name2 = strrep(p,'.lcd','.lcd.2')<br /> nfile1 = fopen(strrep(p,'.lcd','.lcd.1'),'w+');<br /> nfile2 = fopen(strrep(p,'.lcd','.lcd.2'),'w+');<br /> oglcd = fopen(p)<br /> <br /> while strncmp(line,'&lt;data&gt;',6) ~= 1<br /> if strncmp(line,'&lt;data&gt;',6) == 1<br /> break<br /> end<br /> <br /> line = fgets(oglcd)<br /> fprintf(nfile1,'%s',line);<br /> fprintf(nfile2,'%s',line); <br /> end<br /> dn = fopen(strcat(m,'.1'),'r');<br /> up = fopen(strcat(m,'.2'),'r');<br /> <br /> for dc = 1:length(downcastmatrix(:,1))<br /> down = fgets(dn);<br /> fprintf(nfile1,'%s',down);<br /> end<br /> for uc = 1:length(upcastmatrix(:,1))<br /> upc = fgets(up);<br /> fprintf(nfile2,'%s',upc);<br /> end<br /> fprintf(nfile1,'%s\n','&lt;filters_used&gt;')<br /> fprintf(nfile2,'%s\n','&lt;filters_used&gt;')<br /> <br /> %=====update list_hs6======<br /> <br /> h = i*2<br /> r = h-1<br /> newl(r,:) = name1;<br /> newl(h,:) = name2<br /> <br /> end<br /> <br /> dlmwrite('list_hs6_s',newl,'');<br /> <br /> fclose all<br /> <br /> &lt;/pre&gt;<br /> <br /> <br /> &lt;big&gt;[[hs6_split]]&lt;big&gt;</div>

128.111.101.185