

% ComputeVCMParametersATBC('UT2',2,2,OPL/2,0,-1) after redefining
% UTPupilBeforeDL for UT2 U= 26.3927 in the code  (flat VCM case, as shown above)


%!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
% CHECK for UTPupilBeforeDL for UT2 U= 26.3927 for flat VCM on
% AT-STS-G2
%!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


OPL=[60 50 40 30 20 10 ];

radius=[];
PIN=[];
POUT=[];

for i=1:length(OPL);
[RadiusVCMGIP, DistanceInputPupil, DistanceExitPupil, DistanceExitPupilFixedVCM]=ComputeVCMParametersATBC('UT2',2,2,OPL(i)/2,0,-1)
radius=[radius,RadiusVCMGIP];
PIN=[PIN,DistanceInputPupil];
POUT=[POUT,DistanceExitPupil];

end;


h=figure(2);
set(gca,'FontSize',13);
plot(OPL,radius*1000, 'k',OPL,radius*1000,'go');grid
xlabel('OPL in m');ylabel('DL2-VCM ROC in m');
title(' Theoretical DL2-ROC for AT3-G2-DL2-DDl2 with flat STS-VCM')
legend('Theoretical DL2-ROC ')