Hi Rob,
I understand you are facing issues while redistributing the sum of rainfall from the least rainy days to the most rainy days. You can follow the following steps to do the same:
Step 1: Calculate the Annual Rainfall for Each Location
annual_rainfall = sum(P, 3);
Step 2: Estimate the True Number of Rainy Days
estimated_rainy_days = uint16(annual_rainfall ./10);
Step 3: Sort the rainfall data.
[P_sorted, I] = sort(P, 3, 'descend');
Step 4: Redistribute the Rainfall
total_days = 365;
for i = 1:size(P_sorted, 1)
for j = 1:size(P_sorted, 2)
% Calculate the sum of rainfall for the least rainy days
least_rainy_sum = sum(P_sorted(i, j, (estimated_rainy_days(i, j)+ 1):total_days));
% Calculate the average rainfall to be redistributed
least_rainy_average = least_rainy_sum / estimated_rainy_days(i, j);
% Redistribute the rainfall to the most rainy days
for k = 1:estimated_rainy_days(i, j)
P_sorted(i, j, k) = P_sorted(i, j, k) + least_rainy_average;
end
end
end
I hope this helps.
