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


define rpprime 0.5
alias n2omoles r15
alias moles r14
alias volmoles r13
alias volume r12
alias n2otemp r11
alias voltemp r10
alias g r9
alias rprime r8
define MIN_PRESSURE 500

start:
yield

move n2omoles 0
move volmoles 0

l moles sensor TotalMoles
blez moles empty
l r0 sensor RatioNitrousOxide
mul n2omoles r0 moles

l r0 sensor RatioVolatiles
mul volmoles r0 moles

empty:

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

jal calcvolume
s db Setting volume

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

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