That should work fine excepting you didn't allocate for the different means by year and overwrote the variable meanT on each iteration so all you have in the end is the last year average.
You could do something like:
yr1=min(yr); yr2=max(yr); % find first, last year in vector
meanT=zeros(yr2-yr1+1,1); % allocate vector for each year mean
k=0; % initialize mean storage counter
for j=yr1:yr2 % iterate over the year range found
k=k+1; % increment counter
meanT(k)=mean(temp(yr==j)); % compute mean for the year
end
NB: don't need find, the logical addressing vector works to return those locations w/o computing the indices explicitly.
This is fine for learning; the "Matlab way" is to use vector operations instead of loops -- investigate putting your data into a timetable and use the retime function to do this almost automagically.
