jordansread / decadalntl Goto Github PK

%% open files, create time series

% Buffalo lake
[Bff_tme,Bff_lvl] = gFileOpenYYYY([rootDir 'Buffalo_lvl.txt']);
% Crystal lake
[Cry_tme,Cry_lvl] = gFileOpenYYYY([rootDir 'Crystal_lvl.txt']);
% Regional water table
[Reg_tme,Reg_lvl] = gFileOpenYYYY([rootDir 'Region_GW.txt']);
% Crystal water table
[CRw_tme,CRw_lvl] = gFileOpenYYYY([rootDir 'Crystal_GW.txt']);




%% interpolate to daily time series
warning off

Bff_lvl = interp1(Bff_tme,Bff_lvl,dates);

Cry_lvl = interp1(Cry_tme,Cry_lvl,dates);

Reg_lvl = interp1(Reg_tme,Reg_lvl,dates);

CRw_lvl = interp1(CRw_tme,CRw_lvl,dates);

warning on
%% create single daily time series
for j = 1:length(dates)
    vals = [Bff_lvl(j) Cry_lvl(j) Reg_lvl(j) CRw_lvl(j)];
    wtr_lvl(j) = mean(vals(~isnan(vals)));
end

for j = 1:length(years)

    stD = datenum(years(j),1,1);
    enD = datenum(years(j)+1,1,1);    
    useI = ge(dates,stD) & le(dates,enD);
    Stage(j) = mean(wtr_lvl(useI));
end