function [Jstar,CTstar,CPstar,Etastar,ifig]=PlotProp(J,CT,CP,string1,ifig)
ifig=ifig+1;
figure(ifig)
Eta=CT.*J./CP;
subplot(311); plot(J,CT); grid
ylabel('Thrust Coef, CT')
title(string1)
subplot(312); plot(J,CP); grid
ylabel('Power Coef, CP')
subplot(313); plot(J,Eta); grid
ylabel('Efficiency, eta')
xlabel('Advace ratio, J=V/(nD)')
j=find(Eta==max(Eta));
Jstar=J(j);
CTstar=CT(j);
CPstar=CP(j);
Etastar=CTstar.*Jstar./CPstar;

subplot(311); hold on; plot([Jstar,Jstar],[0,CTstar],'r:',[0,Jstar],[CTstar,CTstar],':'); 
string2=['CT*= ',num2str(CTstar),' for J*= ',num2str(Jstar)];
text(.1,CTstar+.006,string2)
hold off

subplot(312); hold on; plot([Jstar,Jstar],[0,CPstar],'r:',[0,Jstar],[CPstar,CPstar],':'); 
string3=['CP*= ',num2str(CPstar),' for J*= ',num2str(Jstar)];
text(.1,CPstar+.005,string3)
hold off

subplot(313); hold on; plot([Jstar,Jstar],[0,Etastar],'r:',[0,Jstar],[Etastar,Etastar],':'); 
string3=['Eta*= ',num2str(Etastar),' for J*= ',num2str(Jstar)];
text(.1,Etastar-.1,string3)
hold off
