% average Eta from daily averages


clear
nx=4320;


folderEta = ['~dmenemen/llc_4320/regions/global/Eta_daily'];
  
fnamesEta = dir([folderEta '/*.2*']);  


numFiles = length(fnamesEta);
numAvg = 27;

etat = 0;
for ii =42:42+numAvg
    
    fnamEta = [folderEta '/' fnamesEta(ii).name];
   
    [eta]= quikreadpcolor_llc(fnamEta,nx);
    etat = eta + etat;
   
end

Etaout = etat./numAvg;


save(['Eta_' num2str(numAvg) 'day_avg'],'Etaout');

