Cd=: 0.65 d=: 0.03 gamma=: 1.4 R=: 53.343 rcr=: (2 % >:gamma) ^ gamma % <:gamma Tabs=: 459.67 + Tf =: 80 NB. Tf [deg F] chocked=: %: (gamma%<:gamma)*(rcr^2%gamma) - rcr^(>:gamma)%gamma P2=: 14.7 P1cr=: P2%rcr NB. Usage: mdot "0 P1 NB. where P1 is vector mdot=: 3 : 0 r=: P2%y md=: (6.30*Cd*d*d*y) % %: R*Tabs if. r > rcr do. md * %: (gamma%<:gamma)*(r^2%gamma) - r^(>:gamma)%gamma else. md * chocked end. ) mdot_sonic=: 3 : 0 r=: P2%y md=: (6.30*Cd*d*d*y) % %: R*Tabs md * chocked ) mdot_subsonic=: 3 : 0 r=: P2%y md=: (6.30*Cd*d*d*y) % %: R*Tabs md * %: (gamma%<:gamma)*(r^2%gamma) - r^(>:gamma)%gamma )