% INDMAT.M %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Index matrices
global coldata
im_x1 = 1:2 * coldata.Nt(1);
im_x2 = im_x1(end) + [1:2 * coldata.Nt(2)];

im_x21 = im_x2(end) + [1:2 * coldata.Nt(3)];
im_x22 = im_x21(end) + [1:2 * coldata.Nt(4)];
im_x23 = im_x22(end) + [1:2 * coldata.Nt(5)];
im_x24 = im_x23(end) + [1:2 * coldata.Nt(6)];

% im_x31 = im_x24(end) + [1:2 * coldata.Nt(7)];
% im_x32 = im_x31(end) + [1:2 * coldata.Nt(8)];
% im_x33 = im_x32(end) + [1:2 * coldata.Nt(9)];
% im_x34 = im_x33(end) + [1:2 * coldata.Nt(10)];
% im_x35 = im_x34(end) + [1:2 * coldata.Nt(11)];
% im_x36 = im_x35(end) + [1:2 * coldata.Nt(12)];

im_xr = im_x24(end) + [1:2];
im_xt = im_xr(end) + [1:2];
im_xS1 = im_xt(end) + [1:2];

im_m1 = im_xS1(end) + [1:coldata.Nt(1)];
im_m2 = im_m1(end) + [1:coldata.Nt(2)];

im_m21 =  im_m2(end) + [1:coldata.Nt(3)];
im_m22 = im_m21(end) + [1:coldata.Nt(4)];
im_m23 = im_m22(end) + [1:coldata.Nt(5)];
im_m24 = im_m23(end) + [1:coldata.Nt(6)];
im_mr = im_m24(end) + 1;
im_mt = im_mr(end) + 1;
im_mS1 = im_mt(end) + 1;
% im_mS2 = im_mS1(end) + 1;
% 
% im_T1 = im_mt(end) + [1:coldata.Nt(1)];
% im_T2 = im_T1(end) + [1:coldata.Nt(2)];
% 
% im_T21  = im_T2(end) + [1:coldata.Nt(3)];
% im_T22 = im_T21(end) + [1:coldata.Nt(4)];
% im_T23 = im_T22(end) + [1:coldata.Nt(5)];
% im_T24 = im_T23(end) + [1:coldata.Nt(6)];
% 
% 
% im_T31 = im_T24(end) + [1:coldata.Nt(7)];
% im_T32 = im_T31(end) + [1:coldata.Nt(8)];
% im_T33 = im_T32(end) + [1:coldata.Nt(9)];
% im_T34 = im_T33(end) + [1:coldata.Nt(10)];
% im_T35 = im_T34(end) + [1:coldata.Nt(11)];
% im_T36 = im_T35(end) + [1:coldata.Nt(12)];
% 
% im_Tr = im_T36(end) + 1;





















