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


define rpprime 0.33333333 # oxi / vol
#define rpprime 0.66666667
alias oximoles r15
alias moles r14
alias volmoles r13
alias volume r12
alias oxitemp r11
alias voltemp r10
alias g r9
alias rprime r8
define MIN_PRESSURE 500

start:
yield

move oximoles 0
move volmoles 0

l moles sensor TotalMoles
blez moles empty
l r0 sensor RatioOxygen
mul oximoles r0 moles

l r0 sensor RatioVolatiles
mul volmoles r0 moles

empty:

sub r1 1 rpprime
mul r1 r1 oximoles
mul r2 rpprime volmoles
sub g r1 r2 # g aka beta

jal calcvolume
#s db Setting volume

l oxitemp oxisensor 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 oxitemp
mul r4 g r4
sub r4 r4 r3 # (T_b-T_a)g - V12.2
div rprime r2 r4

s db Setting rprime

mul r2 rprime oxitemp
sub r3 1 rprime
mul r3 r3 voltemp
mul r4 rprime oxitemp
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

#s mixer On r1


j start

calcvolume:
l volume sensor Volume
j ra