https://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=CERBbt_E 2026-04-20T04:19:47Z Revision history for this page on the wiki MediaWiki 1.34.0 https://wiki-pnb.eri.ucsb.edu/index.php?title=CERBbt_E&diff=748&oldid=prev 2014-04-10T23:22:58Z

<p>Created page with &quot;&lt;pre&gt; % AC9BT % FEB 98,FEB 99 % Stephen Wolfe % Input: none % Look for: list_prr % Output: z*.lcd.?, z*.lcd.strip % Assumptions: fixed for either single casts or yo-yo&#039;s. Ha...&quot;</p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> <br /> % AC9BT<br /> % FEB 98,FEB 99<br /> % Stephen Wolfe<br /> % Input: none<br /> % Look for: list_prr<br /> % Output: z*.lcd.?, z*.lcd.strip<br /> % Assumptions: fixed for either single casts or yo-yo's. Have to<br /> % pick something for at least 1 profile of each yo-yo<br /> % for the profile numbers to work correctly. Can't tell<br /> % if you skip a down-up pair.<br /> % <br /> % Purpose: pick tops and bottoms <br /> % <br /> % Calls: bt_ac9<br /> % Maintain: {filename, columnarraylist, disp_text}<br /> % Uses: perfect for prr processing, ref: prr.ins<br /> % Example: prrbt<br /> % Notes: This will run on NT if the call to bbopstrip is<br /> % commented out and the files are already stripped<br /> % Should label direction in castid table as well as<br /> % profile number. Work from left to right. If you<br /> % choose not to save points for a profile, you can<br /> % repick them and choose save and then go on.<br /> % TO DO: The archive scheme does not work due to changes in<br /> % naming convention. Doesn't create the log directory<br /> % if it's not there either.<br /> %Erik Stassinos Apr 2013<br /> clf;<br /> list_prr=[];<br /> prr_files=[];<br /> not_prr_files=[];<br /> tmpfile=fullfile(pwd,'goprrbt_tmp');<br /> % check for list_prr in pwd<br /> [s,m]=unix('ls list_prr');<br /> if isempty(m)<br /> msg=sprintf('Check your list_prr.');<br /> disp(msg)<br /> return<br /> end<br /> clear s m<br /> % load list_prr<br /> fid=fopen('list_prr.txt','r');<br /> while fid,<br /> check=fgetl(fid);<br /> if check == -1, break, end<br /> list_prr=[list_prr; check];<br /> end<br /> list_prr;<br /> clear fid check<br /> % check presence of prr files<br /> badfilectr=0;<br /> for i=1:size(list_prr,1)<br /> [s,m]=unix(['ls ',list_prr(i,:)]);<br /> if s<br /> badfilectr=badfilectr+1;<br /> end<br /> end<br /> clear i<br /> if badfilectr<br /> if badfilectr==1<br /> msg=sprintf('%i missing prr file.',badfilectr);<br /> disp(msg)<br /> error('Your prrlist includes a file not in pwd.')<br /> else<br /> msg=sprintf('%i missing prr files.',badfilectr);<br /> disp(msg)<br /> error('Your prrlist includes files not in pwd.')<br /> end<br /> end<br /> clear s m badfilectr msg<br /> % loop through <br /> for i=1:size(list_prr,1)<br /> prr_lcd=list_prr(i,:);<br /> % strip em<br /> prr_strip=[prr_lcd '.strip'];<br /> %eval(['!sbtprr_mob ',prr_lcd])<br /> %eval(['!bbopstrip -sz 1mer_time 1depth ' prr_lcd ' 'prr_strip]);<br /> %prrstrip(1,prr_lcd,'1mer_time','depth_ed') %changed 9-18-2012<br /> %prrstrip(1,prr_lcd,'1mer_time','1Depth') %need to change to accomodate screwed up LCD conversion 12-3-12<br /> % prrstrip(1,prr_lcd,'1mer_time','1Depth')<br /> prrstrip(1,prr_lcd,'Master_FrameNumber','LuZDepth')<br /> [fileA,headA]=read_strip(prr_strip);<br /> %eval(['!rm ''prr_strip']);<br /> <br /> [bt_pts]=bt_ac9(fileA);<br /> bt_pts<br /> <br /> if ~isempty(bt_pts),<br /> %create header, castid matrix<br /> %castid index 1mer_time 1depth<br /> %zub11008adt1 7.5000000e+01 2.1710000e+04 2.0443990e+02<br /> %zub11008adb1 1.2810000e+03 3.7655000e+05 1.3007170e+00 <br /> %castdir={ 'dt','db','ub','ut' };<br /> %castdir=[castdir castdir castdir castdir]; % good for 4 yoyo's<br /> %if size(bt_pts,1)&gt;size(castdir,2),<br /> % error('prrbt: increase size of castdir.');<br /> %end<br /> dd=[];<br /> idir=[];<br /> dd=diff(bt_pts(:,3));<br /> dd=[dd; 0];<br /> id=1;<br /> fid=fopen( 'bt.pts', 'w' );<br /> fprintf(fid,'castid index 1mer_time 1depth\n');<br /> for i=1:2:size(bt_pts,1),<br /> if dd(i)&gt;0,<br /> if idir==1,<br /> id=id+1;<br /> end<br /> idir=1;<br /> castdir1(i:i+1,:)=['dt'; 'db'];<br /> else<br /> idir=2;<br /> castdir1(i:i+1,:)=['ub'; 'ut'];<br /> end<br /> bt_castname=[ prr_lcd(1:9),castdir1(i,:),num2str(id) ];<br /> fprintf(fid,'%s %1.7e %1.7e %1.7e \n',bt_castname,bt_pts(i,:));<br /> bt_castname=[ prr_lcd(1:9),castdir1(i+1,:),num2str(id) ];<br /> fprintf(fid,'%s %1.7e %1.7e %1.7e \n',bt_castname,bt_pts(i+1,:));<br /> if idir==2,<br /> id=id+1;<br /> end<br /> end<br /> fclose(fid);<br /> <br /> %do this by calling a script that takes care of the awk and archive<br /> eval([ '!awk &lt; ' prr_lcd ' &gt; ' tmpfile...<br /> ' ''$1==&quot;&lt;sampled_parameters&gt;&quot;{system(&quot;cat bt.pts&quot;)}{print $0}'''])<br /> eval([ '!mv ' tmpfile ' ' prr_lcd ]);<br /> <br /> %archive<br /> %cruisename=eval(['echo ' prr_lcd ' | cut -b 3-6'])<br /> %cruisename=prr_lcd(3:6);<br /> %%cruisename=prr_lcd(11:17); commented 9-18-2012<br /> %archive=[ pwd '/ARCHIVE/' ];<br /> %=&gt;archive=[ '/home/data92/bbop/ARCHIVE/' ];<br /> %cp $lcd_file /home/data92/bbop/ARCHIVE/$cruisename/prr<br /> % eval(['!cp ' prr_lcd ' ' archive cruisename '/prr']);<br /> %echo $lcd_file has been archived to /home/data92/bbop/ARCHIVE/$cruisename/prr<br /> %=&gt;eval(['!echo ' prr_lcd ' has been archived to 'archive cruisename '/prr']);<br /> <br /> prr_files=[prr_files; prr_lcd];<br /> else<br /> not_prr_files=[not_prr_files; prr_lcd];<br /> end<br /> <br /> %cleanup<br /> %clear prr_lcd bt_pts fileA headA prr_strip;<br /> <br /> %eval('!rm bt.pts');<br /> end<br /> prr_files<br /> not_prr_files<br /> <br /> &lt;/pre&gt;</div>

Eriks