Hey @cgs5198, I assume that in the array dimensions of 144x73x841, 144 and 73 represent spatial dimensions (latitude and longitude), and 841 represents time in months.
Below I have provided implemenation to compute globally averaged yearly temperatures from your 3D array of temperature data
temperatureData = rand(144, 73, 841);
% Calculate the number of complete years
monthsPerYear = 12;
numCompleteYears = floor(size(temperatureData, 3) / monthsPerYear);
yearlyGlobalAverages = zeros(1, numCompleteYears);
for yearIdx = 1:numCompleteYears
startMonth = (yearIdx - 1) * monthsPerYear + 1;
endMonth = yearIdx * monthsPerYear;
yearlyData = temperatureData(:, :, startMonth:endMonth);
% Compute monthly global averages
monthlyGlobalAverages = squeeze(mean(mean(yearlyData, 1), 2)); % Average over lat and lon
% Compute the yearly average temperature
yearlyGlobalAverages(yearIdx) = mean(monthlyGlobalAverages);
end
disp('Yearly Global Averages:');
disp(yearlyGlobalAverages);
Hope it helps.
