%% Construct Actual QDC for one year
%
%%

function year_actual_qdc(year)

% end date of period for calculation actual QD:
ds = datenum([year 1 1 0 0 0]); 

for dt=ds:10:ds+365
  [y, m, d] = datevec(dt);
  period = [y,m,d];
  disp(period);
  tic;
  qday_db(period);
  toc
end


