this is my sun data: (date, rise, set)

5/20/95 435.00 1910.00 5/21/95 434.00 1911.00 5/22/95 433.00 1912.00 5/23/95 432.00 1913.00 5/24/95 431.00 1914.00 5/25/95 431.00 1915.00 5/26/95 430.00 1916.00 5/27/95 430.00 1916.00 5/28/95 429.00 1917.00

and my other data: (date, measurements...) that goes every half hour for a week (then later a month to several years)

5/20/95 0:00 14.7 0.22 0.1 97.5 9.9 0.32 6.5 5/20/95 0:30 14.8 0.22 0.1 96.5 9.8 0.35 6.5 5/20/95 1:00 14.8 0.22 0.1 100.2 10.1 0.44 6.5 5/20/95 1:30 14.8 0.22 0.1 99.3 10.1 0.58 6.5 5/20/95 2:00 15 0.22 0.1 102.1 10.3 0.75 6.6 5/20/95 2:30 15 0.22 0.1 102.2 10.3 0.93 6.5 5/20/95 3:00 15 0.22 0.1 102.3 10.3 1.11 6.5 5/20/95 3:30 15.1 0.22 0.1 103 10.4 1.27 6.6 5/20/95 4:00 15.1 0.22 0.1 102.5 10.3 1.41 6.6 5/20/95 4:30 15 0.22 0.1 102 10.3 1.52 6.6 5/20/95 5:00 15 0.22 0.1 101.9 10.3 1.6 6.6 5/20/95 5:30 15 0.22 0.1 99.8 10.1 1.64 6.6 5/20/95 6:00 14.4 0.22 0.1 94.8 9.7 1.65 6.6 5/20/95 6:30 14.4 0.22 0.1 96.3 9.8 1.61 6.5 5/20/95 7:00 14.6 0.22 0.1 98.9 10.1 1.53 6.6

I need to run the following code following "day_2=..." for every day in the week and later the month and the year

%Hudson River Water Quality 1995-2000 data=xlsread('Hudson_test1.xlsx'); tm=data(:,1)+datenum(1904,1,1); [yr,mo,day,hr,mn]=datevec(tm); DO_p=data(:,5); DO=data(:,6); z=data(:,7);

dif_DO=DO(2:end)-DO(1:end-1); air_sea=(1-((DO_p(1:end-1)+DO_p(2:end))/200))0.50.5; do_flux=(dif_DO.*z(2:end))-air_sea;

sun=xlsread('rise_set_test.xls'); tm_sun=sun(:,1)+datenum(1904,1,1); [yr_s,mo_s,day_s]=datevec(tm_sun); %puts sunrise/set into time rise=(sun(:,2)); set=sun(:,3);

%put sunrise and sunset into proper matlab date format srh=floor(rise/100); srm=rise-100*srh; sr=(srh+srm/60)/24; sr_tm=sr+tm_sun; %puts in date format

ssh=floor(set/100); ssm=set-100*ssh; ss=(ssh+ssm/60)/24 ; ss_tm=ss+tm_sun;

day_2 = tm(2:96);

j_2 = day_2>=sr_tm(2,1) & day_2<=ss_tm(2,1); m_2 = day_2<=sr_tm(2,1) & day_2>=tm(49,1) | day_2>=ss_tm(2,1); do_dy_2 = sum(do_flux(j_2)); do_nt_2 = sum(do_flux(m_2)); resp_rt_2 = (do_nt_2.(-1))/(nnz(m_2)/2); ttl_resp_2 = resp_rt_2.24; day_resp_2 = resp_rt_2.*(nnz(j_2)/2); gp_2 = day_resp_2 + do_dy_2; NEM_2 = gp_2-ttl_resp_2;