from rocketcea.cea_obj import CEA_Obj
from pylab import *
Pc = 200.0
eps = 20.0
xL = []
ispL = []
for x in range(10, 100, 5):
propName = 'HYD%g'%x
ispObj = CEA_Obj(propName=propName)
xL.append( x )
IspFroz, Cstar, Tcomb = ispObj.getFrozen_IvacCstrTc( Pc=Pc, eps=eps, frozenAtThroat=1 )
ispL.append( IspFroz )
plot(xL, ispL, linewidth=4)
grid(True)
title( 'Hydrazine Ideal Performance vs. Ammonia Dissociation\nPc=%g psia, Area Ratio=%g'%(Pc, eps) )
xlabel( '% Ammonia Dissociation' )
ylabel( 'IspVac (sec)' )
savefig('amm_dissociation.png', dpi=120)
show()