function InitInvStokesM(Date)
		
global 	gB1InvStokesM...
	 	gB2InvStokesM...
	 	gB3InvStokesM...
	 	gB4InvStokesM...
	 	gB5InvStokesM...
	 	gB6InvStokesM...
	 	gC1InvStokesM...
	 	gC2InvStokesM...
	 	gC3InvStokesM;
		
global 	gAverageCimelFlag;

	% MMR108:5-7, Ref=184.65
Ref = 184.65;
B1Angles = [357.127; 43.061; 87.537] - [Ref; Ref; Ref];

	% MMR104:1-3, Ref=0
Ref = 0;
B2Angles = [91.893; 140.25; 4.202] - [Ref; Ref; Ref];

	% MMR108:1-3, Ref=184.65
Ref = 184.65;
B3Angles = [178.729; 222.001; 269.525] - [Ref; Ref; Ref];

	% MMR105:1-3, Ref=0
Ref = 0;
B4Angles = [93.695; 48.15; 2.235] - [Ref; Ref; Ref];

    % MMR105:5-7, Ref=0
Ref = 0;
B5Angles = [92.617; 46.856; 2.058] - [Ref; Ref; Ref];

	% MMR104:5-7, Ref=0
Ref = 0;
B6Angles = [273.035; 319.327; 2.216] - [Ref; Ref; Ref];

if Date <= 910716
		% MMR104:1-3, Ref=181.6
	Ref = 181.6;
	B2Angles = [115.6; 202.0; 250.9] - [Ref; Ref; Ref];
	
		% MMR104:5-7, Ref=95.35
	Ref = 95.35;
	B6Angles = [208.6; 299.5; 160.5] - [Ref; Ref; Ref];
	end;

if Date <= 910716 % (1 Degree FOV)

	gAverageCimelFlag = 1;

		% CimelQ:1-3, Ref=184.65
	Ref = 184.65;
	C1Angles = [123.308; 342.266; 80.697] - [Ref; Ref; Ref];
	
		% CimelE:1-3, Ref=90.2
	Ref = 90.2;
	C2Angles = [81.618; 121.136; 31.609] - [Ref; Ref; Ref];
	
		% CimelZ:1-3, Ref=95.35
	Ref = 95.35;
	C3Angles = [83.972; 118.719; 167.245] - [Ref; Ref; Ref];
	
else	% Date > 910716 (12 Degree FOV)

	gAverageCimelFlag = 0

		% CimelQ:1-3, Ref=2.5
	Ref = 2.5;
	C1Angles = [124.8; 343.7; 80.46] - [Ref; Ref; Ref];
	
		% CimelE:1-3, Ref=2.5
	Ref = 2.5;
	C2Angles = [3.; 40.8; 92.7] - [Ref; Ref; Ref];
	
		% CimelZ:1-3, Ref=2.5
	Ref = 2.5;
	C3Angles = [343.614; 21.377; 78.979] - [Ref; Ref; Ref];

end	% if Date <= 910716

gB1InvStokesM = InvStokesMatrix(B1Angles);
gB2InvStokesM = InvStokesMatrix(B2Angles);
gB3InvStokesM = InvStokesMatrix(B3Angles);
gB4InvStokesM = InvStokesMatrix(B4Angles);
gB5InvStokesM = InvStokesMatrix(B5Angles);
gB6InvStokesM = InvStokesMatrix(B6Angles);
gC1InvStokesM = InvStokesMatrix(C1Angles);
gC2InvStokesM = InvStokesMatrix(C2Angles);
gC3InvStokesM = InvStokesMatrix(C3Angles);

	% The test output values should be:
	%	 I =   6.6860
	%	 Q =   1.6914
	%	 U =  -2.2447
	% 	PP =  42.0373
	%	 X = -26.5004

TestAngles = [5.25; 95.00; 52.50];
TestInvStokesM = InvStokesMatrix(TestAngles);
TestStokesV = TestInvStokesM * [7.94; 5.41; 4.08]
TestP = DegOfPolarization(TestStokesV)
TestChi = AngOfPolarizedPlane(TestStokesV)
