https://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=BB9spectra 2026-05-23T23:52:56Z Revision history for this page on the wiki MediaWiki 1.34.0 https://wiki-pnb.eri.ucsb.edu/index.php?title=BB9spectra&diff=855&oldid=prev 2015-10-27T21:36:34Z

<p>Created page with &quot;&lt;pre&gt; function BB9spectra() %plot BB9 spectra 5m binned with depth pathname = pwd; dhnumber = pathname(end-2:end); %=====LISST DATA SET MODULE======== data_lst = load(&#039;LI...&quot;</p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> function BB9spectra()<br /> %plot BB9 spectra 5m binned with depth<br /> <br /> <br /> <br /> <br /> <br /> pathname = pwd;<br /> dhnumber = pathname(end-2:end);<br /> %=====LISST DATA SET MODULE========<br /> data_lst = load('LISST_indexarr');<br /> lst_depth = .01*data_lst(:,37);<br /> lst_date1 = num2str(data_lst(:,39));<br /> lst_date2 = num2str(data_lst(:,40));<br /> lst_day = str2num(lst_date1(:,1:3));<br /> lst_hrs = str2num(lst_date1(:,4:5));<br /> <br /> %when lisst minutes reset every 60 sec<br /> for ml =1:length(lst_date2)<br /> hld = lst_date2(ml,1:2);<br /> char(hld);<br /> hld = str2num(hld);<br /> if isempty(hld)<br /> lst_min(ml) =0;<br /> continue<br /> end<br /> lst_min(ml) = hld;%str2num(lst_date2(ml,1:2))<br /> <br /> end<br /> lst_min = (lst_min)';<br /> lst_sec = str2num(lst_date2(:,3:4));<br /> <br /> JD = lst_hrs.*3600 + lst_min.*60 + lst_sec; %this is in seconds<br /> %start from zero<br /> JD = JD-JD(1); %lisst time in seconds starting from zero<br /> %======END LISST MODULE========<br /> <br /> <br /> %=====BB9 MODULE========<br /> % since BB9 time is more accurate, use it for merging<br /> ECOname = strcat('archive_23_ECO.',dhnumber);<br /> %[BB9_data,BB9_HDR] = hdrload('archive_23_ECO.040');<br /> [BB9_data,BB9_HDR] = hdrload(ECOname);<br /> %BB9_data = load(ECOname);<br /> BB9_time = BB9_data(:,1); %(in ms) si<br /> BB9_time = BB9_time - BB9_time(1);<br /> BB9_time = BB9_time./1000;<br /> <br /> <br /> %====CTD MODULE=======<br /> CTDname = strcat('archive_22_SBE25CTD.',dhnumber)<br /> <br /> [CTD_dat,CTD_hdr] = hdrload(CTDname);<br /> <br /> %[CTD_dat,CTD_hdr] = hdrload('archive_22_SBE25CTD.040');<br /> ctd_time = CTD_dat(:,1); %ms<br /> ctd_time = ctd_time./1000; %seconds<br /> %ctd_depth = CTD_dat(:,2);<br /> %ctd_temp = CTD_dat(:,3);<br /> %ctd_cond = CTD_dat(:,4);<br /> %ctd_sal = CTD_dat(:,5);<br /> <br /> %====ACS MODULE=======<br /> ACSname = strcat('archive_21_ACS.',dhnumber)<br /> [ACS_dat,ACS_hdr] = hdrload(ACSname);<br /> %ACS_D = interp1(nmini30s_T,nmini30s_D,ctd_time);<br /> <br /> <br /> %=====MERGING MODULE=====<br /> <br /> BB9_new = interp1(BB9_time,BB9_data,ctd_time);<br /> sample =[];<br /> if issorted(data_lst(:,1)) == 0<br /> issorted(JD)<br /> issorted(ctd_time)<br /> issorted(data_lst(:,1))<br /> issorted(data_lst(:,37))<br /> <br /> %at this point CTD and BB9 data is parsed. BB data is matched to CTD time<br /> <br /> %======= FILTERING NAN IF TIME RECORDS GO OVER =====<br /> nn= isnan(BB9_new);<br /> non= find(nn(:,1) ==1);<br /> CTD_dat(non,:) = [];<br /> BB9_new(non,:) = [];<br /> %=====<br /> ctd_time = CTD_dat(:,1); %ms<br /> ctd_time = ctd_time./1000; %seconds<br /> ctd_depth = CTD_dat(:,2);<br /> ctd_temp = CTD_dat(:,3);<br /> ctd_cond = CTD_dat(:,4);<br /> ctd_sal = CTD_dat(:,5);<br /> <br /> CTD30s_D = ctd_depth;%interp1(nmini30s_T,nmini30s_D,ctd_time); %gets a depth dimention relative to 8Hz ctd time (same dimension as ctd data)<br /> ACS_D = interp1(ctd_time,ctd_depth,.001*ACS_dat(:,1));<br /> BB9_D = interp1(ctd_time,ctd_depth,BB9_time);<br /> <br /> <br /> CTD_dat(:,2) = CTD30s_D; %repair CTD depth column in data<br /> end<br /> <br /> <br /> %===============================<br /> %====BB9 spectra plotting=======<br /> <br /> <br /> <br /> sepparation_index = find(BB9_D == max(BB9_D)); %where to split cast<br /> sep_ind = sepparation_index(1);<br /> for i = 1:size(BB9_new,2)<br /> BB9_downcast(:,i) = smooth(BB9_new(1:sep_ind,i));<br /> end<br /> BB9_D = smooth(BB9_D(1:sep_ind,:));<br /> <br /> %Tref = 22.6 %(from factory cal)<br /> %header_b = ['Time(ms) c400.9 c404.1 c407.3 c410.3 c413.7 c417.1 c420.5 c424.2 c427.9 c431.7 c435.2 c438.6 c442.5 c446.3 c450.0 c454.3 c458.1 c462.2 c466.2 c470.5 c474.5 c478.5 c482.7 c486.7 c490.6 c494.6 c498.6 c502.4 c506.6 c510.7 c514.7 c518.9 c522.7 c526.8 c530.8 c534.6 c538.4 c542.3 c546.1 c549.9 c553.6 c557.4 c561.0 c564.7 c568.1 c571.3 c574.7 c577.8 c581.2 c584.6 c588.3 c591.8 c595.2 c598.9 c602.6 c606.3 c610.2 c613.9 c617.7 c621.5 c625.2 c629.2 c632.9 c636.8 c640.7 c644.3 c648.4 c652.4 c656.0 c660.0 c663.9 c667.8 c671.3 c674.9 c679.1 c682.7 c686.3 c689.7 c693.4 c697.0 c700.5 c703.9 c707.1 c710.7 c714.2 c717.4 c720.8 c724.1 c727.5 c730.5 c734.2 a402.6 a405.9 a408.9 a411.9 a415.3 a418.7 a422.4 a426.0 a429.5 a433.3 a436.8 a440.3 a444.0 a447.9 a452.0 a455.9 a459.9 a463.8 a468.2 a472.2 a476.4 a480.4 a484.5 a488.5 a492.3 a496.1 a500.1 a504.1 a508.3 a512.5 a516.5 a520.3 a524.4 a528.2 a532.3 a536.1 a539.7 a543.9 a547.6 a551.4 a555.2 a558.7 a562.8 a566.2 a569.0 a572.2 a575.6 a579.0 a582.4 a585.8 a589.3 a592.7 a596.4 a599.9 a603.6 a607.3 a611.3 a615.0 a618.9 a622.7 a626.4 a630.2 a634.0 a638.0 a641.8 a645.6 a649.6 a653.5 a657.5 a661.1 a665.0 a668.8 a672.6 a676.3 a680.1 a683.8 a687.1 a690.8 a694.7 a698.1 a701.5 a705.0 a708.4 a711.8 a715.2 a718.5 a721.7 a725.1 a728.4 a731.8 a735.0 iTemp(C) Conduct Depth XTemp Salinity'];<br /> header_b = ['Time(ms) Beta(412) BetaP(412) bbP(412) bb(412) Beta(440) BetaP(440) bbP(440) bb(440) Beta(488) BetaP(488) bbP(488) bb(488) Beta(510) BetaP(510) bbP(510) bb(510) Beta(532) BetaP(532) bbP(532) bb(532) Beta(595) BetaP(595) bbP(595) bb(595) Beta(660) BetaP(660) bbP(660) bb(660) Beta(676) BetaP(676) bbP(676) bb(676) Beta(715) BetaP(715) bbP(715) bb(715)'];<br /> BBwvln = [412,440,488,510,532,595,660,676,715]<br /> legendvals = regexp(header_b,'\s','split')<br /> channels = legendvals(4:4:end-1); %generating a list of ACS channels for TScal<br /> <br /> <br /> <br /> bwindow =5;<br /> disp('doing BB9 binning')<br /> [binned,stats,stdev,medi]=databinner(BB9_D(1:sep_ind),BB9_downcast,bwindow)<br /> %[binned,stats,stdev,medi]=databinner(ACS_D(1:sep_ind),ACS_downcast,bwindow)<br /> for i=1:length(legendvals)<br /> header_std{i} = char([' std',legendvals{i}])';<br /> end<br /> %header_std = char(header_std)';<br /> <br /> <br /> <br /> m35_index = find(abs(BB9_D(1:sep_ind) -35) == min(abs(BB9_D(1:sep_ind) - 35)));<br /> m50_index = find(abs(BB9_D(1:sep_ind) -50) == min(abs(BB9_D(1:sep_ind) - 50)));<br /> m80_index = find(abs(BB9_D(1:sep_ind) -80) == min(abs(BB9_D(1:sep_ind) - 80)));<br /> m95_index = find(abs(BB9_D(1:sep_ind) -95) == min(abs(BB9_D(1:sep_ind) - 95)));<br /> <br /> figure<br /> plot(BBwvln,BB9_downcast(m35_index(1),4:4:end-1),BBwvln,BB9_downcast(m50_index(1),4:4:end-1),BBwvln,BB9_downcast(m80_index(1),4:4:end-1),BBwvln,BB9_downcast(m95_index(1),4:4:end-1))%,400:3.8202:740,avg_spec_001May2012unfilb(:,2:91), 400:3.8202:740,avg_spec_001May2012filc(:,2:91),400:3.8202:740,avg_spec_001May2012unfilc(:,2:91))<br /> title(['BB9 279 spectra ',ECOname(end-7:end)])<br /> xlabel('400-740nm');<br /> ylabel('bbp data 1/m')<br /> legend('35m','50m','80m','95m')%,'May2012-unfilter','May2012-newfilter','May2012-newunfilter','May2012-newfilter runII','May2012-newunfilter runII')<br /> saveas(gcf, ['BB9_bbp',ECOname(end-7:end)], 'png')<br /> <br /> <br /> figure<br /> plot(BBwvln,BB9_downcast(m35_index(1),5:4:end),BBwvln,BB9_downcast(m50_index(1),5:4:end),BBwvln,BB9_downcast(m80_index(1),5:4:end),BBwvln,BB9_downcast(m95_index(1),5:4:end))%,400:3.8202:740,avg_spec_001May2012unfilb(:,2:91), 400:3.8202:740,avg_spec_001May2012filc(:,2:91),400:3.8202:740,avg_spec_001May2012unfilc(:,2:91))<br /> title(['BB9 279 spectra ',ECOname(end-7:end)])<br /> xlabel('400-740nm');<br /> ylabel('bb data 1/m')<br /> legend('35m','50m','80m','95m')%,'May2012-unfilter','May2012-newfilter','May2012-newunfilter','May2012-newfilter runII','May2012-newunfilter runII')<br /> saveas(gcf, ['BB9_bb',ECOname(end-7:end)], 'png')<br /> <br /> <br /> &lt;/pre&gt;</div>

Eriks