alias mixer d0
alias sensor d1
alias oxidizersensor d3
alias volsensor d4

# 0.33333 for O2, 0.5 for nitrous
define rpprime 0.33333333
alias oxidizermoles r15
alias moles r14
alias volmoles r13
alias volume r12
alias oxidizertemp r11
alias voltemp r10
alias g r9
alias rprime r8
alias beta r7
define MIN_PRESSURE 500

start:
yield

move oxidizermoles 0
move volmoles 0

l moles sensor TotalMoles
blez moles empty
l r0 sensor RatioOxygen		# change this to be the oxidizer
mul oxidizermoles r0 moles

l r0 sensor RatioVolatiles
mul volmoles r0 moles

empty:

mul r0 rpprime volmoles
sub r1 1 rpprime
mul r1 r1 oxidizermoles
sub g r1 r0
# beta = (1-r'')ma - r''mb

#sub g oxidizermoles volmoles # ratio of 0.5
#mul g rpprime g
# # formula gets more complex if r'' != 0.5

jal calcvolume
s db Setting volume

l oxidizertemp oxidizersensor Temperature
l voltemp volsensor Temperature

mul r3 volume 12.187
mul r2 g voltemp
mul r5 rpprime r3
sub r2 r2 r5 # gT_b - r''V12.2

sub r4 voltemp oxidizertemp
mul r4 g r4
sub r4 r4 r3 # (T_b-T_a)g - V12.2
div rprime r2 r4

mul r2 rprime oxidizertemp
sub r3 1 rprime
mul r3 r3 voltemp
mul r4 rprime oxidizertemp
add r3 r3 r4
div r2 r2 r3

mul r2 r2 100
#sub r2 100 r2  # if your volatiles are on Input 1
#s d5 Setting r2
s mixer Setting r2


l r0 sensor Pressure
slt r1 r0 MIN_PRESSURE
# need more pressure

l r0 oxidizersensor Pressure
sge r1 r0 10
l r0 volsensor Pressure
sge r2 r0 10
and r1 r1 r2

l r0 mixer On
and r1 r0 r1

s mixer On r1


j start

calcvolume:
l volume sensor Volume
j ra