% Make meas
GlobalVars;
Ny=26;
Meas    =zeros(NVar,1);
MeasStd =zeros(NVar,1);
MeasU   =zeros(Ny,NVar);
i=1;

% Load matrix of 21 process measurements
ProcessData;

% Make measurement index matrix U and measurement uncertainy matrix MeasStd
Ix=IxStr( 1);   ix=Ix(3);  MeasU(i,ix)=1;   MeasStd(ix,1)=0.01;  i=i+1;  % Feed paraffine mol frac.    14-AI-015D
Ix=IxStr( 1);   ix=Ix(4);  MeasU(i,ix)=1;   MeasStd(ix,1)=0.01;  i=i+1;  % Feed naphthene mol frac.    14-AI-015E
Ix=IxStr( 1);   ix=Ix(5);  MeasU(i,ix)=1;   MeasStd(ix,1)=0.01;  i=i+1;  % Feed aromatics mol frac.    14-AI-015F
Ix=IxStr( 1);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % Feed temperature            14-TI-131
Ix=IxStr( 2);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % E-1401 csi temperature      14-TI-012
Ix=IxStr( 4);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % E-1401 cso temperature      14-TI-009
Ix=IxStr( 5);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % H-1401 outlet temperature   14-TIC-104
Ix=IxStr(15);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % R-1401 outlet temperature   14-TI-001
Ix=IxStr(16);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % H-1402 outlet temperature   14-TIC-098
Ix=IxStr(26);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % R-1402 outlet temperature   14-TI-002
Ix=IxStr(27);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % H-1403 outlet temperature   14-TIC-100
Ix=IxStr(37);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % R-1403 outlet temperature   14-TI-003
Ix=IxStr(38);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % H-1404 outlet temperature   14-TIC-102
Ix=IxStr(48);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % R-1404 outlet temperature   14-TI-004
Ix=IxStr(50);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % E-1401 hso temperature      14-TI-011
Ix=IxStr(51);   ix=Ix(8);  MeasU(i,ix)=1;   MeasStd(ix,1)=0.2;   i=i+1;  % D-1405 pressure             14-PIC-125
Ix=IxStr(52);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % D-1405 outlet temperature   14-TI-013
Ix=IxStr(54);   ix=Ix(1);  MeasU(i,ix)=1;   MeasStd(ix,1)=0.1;   i=i+1;  % Recic. gas H2 mole fraction 14-AI-025
Ix=IxStr(54);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % C-1401 inlet temperature    14-TI-132A
Ix=IxStr(55);   ix=Ix(7);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % C-1401 outlet temperature   14-TI-015
Ix=IxStr(55);   ix=Ix(8);  MeasU(i,ix)=1;   MeasStd(ix,1)=0.2;   i=i+1;  % C-1401 outlet pressure      14-PI-012
Ix=IxConstr(1); ix=Ix(1);  MeasU(i,ix)=1;   MeasStd(ix,1)=1;     i=i+1;  % Reformate RON               14-AI-005A
Ix=IxConstr(1); ix=Ix(3);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % Feed flow mass flow         14-FIC-004
Ix=IxConstr(1); ix=Ix(4);  MeasU(i,ix)=1;   MeasStd(ix,1)=1;     i=i+1;  % C-1401 outlet flow          14-FI-002
Ix=IxConstr(1); ix=Ix(5);  MeasU(i,ix)=1;   MeasStd(ix,1)=1;     i=i+1;  % D-1405 vapor product flow   14-FIC-086+14-FI-012
Ix=IxConstr(1); ix=Ix(6);  MeasU(i,ix)=1;   MeasStd(ix,1)=3;     i=i+1;  % Reformate product flow      14-FIC-013

Meas=PV(12,:)';
