How can i define sid_allyears = zeros(361,361,432); in my script

% open binary file for 
% Polar Pathfinder mean monthly 25 km EASE-Grid Sea Ice Motion Vectors, Version 
    folderpath = 'C:\Users\sk\Documents\MATLAB\nsidc0116_ice_motion_vector_v3\';

years = 1979:2014;
sid_allyears = zeros(361,361,432);
yearindex = 0;
for j = years;
    disp(j);
    yearindex = yearindex+1;

      stryear = num2str(j);
      filepath = fullfile(folderpath,stryear);

        files = dir(fullfile(filepath,'\icemotion.grid.month.*.n.v3.bin') );
        files = {files.name};

        r_years = []; % (361,361,12)
        d_years = []; % (361,361,12)

        mm  = 0 ;
        for i=1:numel(files);
            disp(files{i}); 
            mm = mm+1;
            filename = fullfile(filepath,files{i});

            fid = fopen(filename);
            rawdata = fread(fid,[361, inf],'int16');
            fclose(fid);

            % reshape it into a [3x361x361] matrix
            tdata = reshape(rawdata,[3, 361, 361]);

            % change 3d matrix to simple matrix, and divide by 10 to get cm/second 
            u = squeeze(tdata(1,:,:)./10.);
            v = squeeze(tdata(2,:,:)./10.);
            error = tdata(3,:,:)./10;

            [TH,R] = cart2pol(v,u);
            D = rad2deg(TH);

            r_years(:,:,mm) = R;
            d_years(:,:,mm) = D;  

        end
        **sid_allyears = ?????**

end

[Image Analyst]

Other than how do I format my code, what is your question? Is the zeros() function not preallocating space for your array?

[Sophia]

I have allocated the size but not sure how to define it in the code, I tried

*sid_allyears(:,:,:yearindex) = r_years;*

and thought once it would be 361*361*12*36 i can change it to 361*361*432, but its not working out

tr = zeros(361,361,432);
yrlst = 1979:2014;
mnthlst = 1:12;

for i=1:size( sid_r_allyears,1)
    for j=1:size( sid_r_allyears,2)
        ni=0;
        for k=1:length(mnthlst) 

              for l=1:length(yrlst) 
                  ni=ni+1;
                  tr(i,j,ni) =  sid_r_allyears(i,j,k,l);
              end
          end
      end
  end

ind=10:16; % Winter average over month Oct to following year's April
for k=1:(ni/12)-1
    meanr = nanmean(tr(i,j,ind));
    tmeanr(i,j,k) = meanr;
    ind = ind  + 12;
end

[Walter Roberson]

sid_allyears(:,:,:,yearindex) = r_years;

[Sophia]

Thanks Walter, it is working but the problem lies in the next part now. Since i have 432 files, ind values showing are 430-436, but since i don't have the complete winter data values for 2014(i just have oct to dec for 2014 not jan to apr 2015), so i want this ind to show 418-424 instead.

ind=10:16; % Winter average over month Oct to following year's April
for k=1:(ni/12)-1
    meanr = nanmean(tr(i,j,ind));
    tmeanr(i,j,k) = meanr;
    ind = ind  + 12;
end
The values in * *meanr* = nanmean(tr(i,j,ind));* are showing all zeros