wing_param3

for wing=1:length(taper(1,:))
    figure(wing)
    clf
    set(gcf,'Position',[0 0 1250 750])
    
    subplot(3,2,1)
    plot3(latt_w(wing).XYZ(:,:,1)'/0.3048,latt_w(wing).XYZ(:,:,2)'/0.3048,latt_w(wing).XYZ(:,:,3)'/0.3048,'k')
    axis equal,hold on
    xlabel('Wing x-coordinate')
    ylabel('Wing y-coordinate')
    zlabel('Wing z-coordinate')
    title('3-D Wing configuration')
    axis tight
    grid on
    hold off
    
    cd output
    fname=strcat(Run_ID(wing).name,'-Cx');
    load(fname)
    cd ..
    
    cdi=iszero(results.CL)^2/(pi*(ref_w(wing).b_ref/0.3048)^2/(ref_w(wing).S_ref/0.3048^2));
    Oswald_e=iszero(results.CD)/cdi;
    subplot(3,2,3)
    hold on, grid on
    plot(results.ystation(:,1)/0.3048,results.ForcePerMeter(:,1)/4.44822162*0.3048);
    Title('Spanload on main wing');
    ylabel('lbf per foot');
    xlabel('Spanstation')
    hold off
    
    subplot(3,2,5)
    hold on, grid on
    plot(results.ystation(:,1)/0.3048,results.CL_local(:,1),'r');
    Title('Spanload on main wing');
    ylabel('CL Local');
    xlabel('Spanstation')
    hold off
       
    subplot(3,2,[2 4 6])
    axis off
text(0,1,'Tornado Computation Results ')
text(0,.95,'JID: '); text(0.25,0.95,Run_ID(wing).name)
text(0.4,.95,'Downwash matrix condition: '); text(0.8,0.95,num2str(results.dwcond))

text(0,.90,'Reference area: ');	text(0.25,0.90,num2str(ref_w(wing).S_ref/0.3048^2));
text(0,.85,'Root chord: ');text(0.25,.85,num2str(geometry(wing).c/0.3048));
text(0,.8,'Reference span: ');text(0.25,.8,num2str(ref_w(wing).b_ref/0.3048));  

text(0.4,.9,'Taper: ');text(0.8,.9,num2str(geometry(wing).T'));
text(0.4,.85,'Reference point pos: ');text(0.8,.85,num2str(geometry(wing).ref_point));
text(0.4,.8,'Oswalds efficiency: ');text(0.8,.8,num2str(Oswald_e));
   
text(0,.7,'Net Wind Forces: (lbf)');			
	text(0.0,.65,'Drag: ');		text(0.1,.65,num2str(iszero(results.D/4.44822162)));
   text(0.0,.6,'Side: ');			text(0.1,.6,num2str(iszero(results.C/4.44822162)));
   text(0.0,.55,'Lift: ');		    text(0.1,.55,num2str(iszero(results.L/4.44822162)));
   
text(0.35,.7,'Net Body Forces: (lbf)');			
   text(0.35,.65,'X: ');	    	text(0.4,.65,num2str(iszero(results.FORCE(1)/4.44822162)));
   text(0.35,.6,'Y: ');			text(0.4,.6,num2str(iszero(results.FORCE(2)/4.44822162)));
   text(0.35,.55,'Z: ');		    text(0.4,.55,num2str(iszero(results.FORCE(3)/4.44822162)));

text(0.7,.7,'Net Body Moments: (lbf-ft)');
   text(0.7,.65,'Roll: ');			text(0.8,.65,num2str(iszero(results.MOMENTS(1)/(4.44822162*0.3048))));
   text(0.7,.6,'Pitch: ');		    text(0.8,.6,num2str(iszero(results.MOMENTS(2)/(4.44822162*0.3048))));   
   text(0.7,.55,'Yaw: ');		    text(0.8,.55,num2str(iszero(results.MOMENTS(3)/(4.44822162*0.3048))));
   
text(0,.4,'CL ');	    text(0.1,.4,num2str(iszero(results.CL)))   
text(0,.35,'CD ');	    text(0.1,.35,num2str(iszero(results.CD))) 
text(0,.3,'CY ');	    text(0.1,.3,num2str(iszero(results.CY)))

text(0.35,.4,'CZ ');    text(0.45,.4,num2str(iszero(results.CZ)))
text(0.35,.35,'CX ');	text(0.45,.35,num2str(iszero(results.CX)))
text(0.35,.3,'CC ');	text(0.45,.3,num2str(iszero(results.CC)))

text(0.7,.4,'Cm ');		text(0.8,.4,num2str(iszero(results.Cm)))
text(0.7,.35,'Cn ');	text(0.8,.35,num2str(iszero(results.Cn)))
text(0.7,.3,'Cl ');		text(0.8,.3,num2str(iszero(results.Cl)))

text(0,.2,'STATE: ');
text(0,.15,'alpha: '); text(.15,.15,num2str(iszero(state.alpha*180/pi)));
text(0,.1,'beta: '); text(.15,.1,num2str(iszero(state.betha*180/pi)));
text(0,.05,'Airspeed: ');text(.15,.05,num2str(iszero(state.AS/0.3048)));
text(0,.0,'Density: ');text(.15,.0,num2str(iszero(state.rho)));

text(0.3,.15,'P: '); text(.45,.15,num2str(iszero(state.P)));
text(0.3,.10,'Q: '); text(.45,.10,num2str(iszero(state.Q)));
text(0.3,.05,'R: '); text(.45,.05,num2str(iszero(state.R)));


text(0.6,.1,'Rudder setting [deg]:'); text(.9,.1,num2str(geometry(wing).flap_vector*180/pi));
end
cd ..