clear all close all clc load PalletteLakeTemp.mat load PalletteLakeDO.mat wtemp=PalletteLakeTemp; wDO=PalletteLakeDO; wdepthstart=1.0564897943986; wdepthtotal=19; depths = length(wtemp(1,:)); %number of depths days = length(wtemp(:,1)); %number of days %wDO = flipud(wDO); %wtemp = flipud(wtemp); for i = 1:days; TDO3(i,1) = interp1(nonduplicate(nonduplicate(wDO(i,:))),wtemp(i,:),3); end % load TDO3_final %% TDO3 7 day TDO3_7day = movmean(TDO3,7); %% TDO3 30 day TDO3_30day = movmean(TDO3,30); %% Calculate the yearly TDO3 max values for i = 1:37 TDO3_max_7(i,1) = max(TDO3_7day(365*(i-1)+1:365*(i-1)+263,1)); TDO3_max_30(i,1) = max(TDO3_30day(365*(i-1)+1:365*(i-1)+263,1)); TDO3_noav(i,1) = max(TDO3(365*(i-1)+1:365*(i-1)+263,1)); end %% Calculate the daily dosage for i =1:length(TDO3); if TDO3(i)>=17, dosage_day(i,1) = TDO3(i,1)-17; else dosage_day(i,1) = 0; end end %% Calculate the yearly dosage for i = 1:37 dosage_year(i,1) = sum(dosage_day(365*(i-1)+1:365*(i-1)+263),1); end %% calculate the minimum depth between max and min for i = 1:length(TDO3) wdepths(i,:) = 0:wdepthstart:wdepthstart*wdepthtotal; end for i = 1:days; depthDO(i,1) = interp1(nonduplicate(nonduplicate(wDO(i,:))),wdepths(i,:),3); depthTEMP(i,1) = interp1(nonduplicate(nonduplicate(wtemp(i,:))),wdepths(i,:),17); distance(i,1) = depthTEMP(i,1)-depthDO(i,1); end for i = 1:37 distance_min(i,1) = nanmin(distance(365*(i-1)+1:365*(i-1)+263,1)); end %% make file alldata(:,1) = 1979:2015; alldata(:,2) = TDO3_noav; alldata(:,3) = TDO3_max_7; alldata(:,4) = TDO3_max_30; alldata(:,5) = dosage_year; alldata(:,6) = distance_min;