![]() | Sample 3 Array and Table |
PdependsOnH()

FUNCTION PdependsOnH()
SYSTEM(BackGround=995, FontSize=19)
DO step = 1, 81 ! from 1 to 81
height = (step - 1) * 100 ! from 0 to 8000
$HP(step, 1) = height
$HP(step, 2) = BarometricPressure(height)
ENDDO
F2()
F3()
F4()
END
FUNCTION F2() ! standard loop graph
DLG(R=2,C=4,BackG=99, AXis=3, X=1,MIN=0,MAX=8000,TIt='using standard loop', Y=1,MIN=300,MAX=1100)
DO i = 0, 10
height = 1000 * i
pressure = BarometricPressure(height)
LINE(AX=3, X=height, Y=pressure, Symbol='🞉', Draw=9)
ENDDO
END
FUNCTION F3() ! internal $ loop variable
DLG(R=2,C=4,BackG=990, AXis=4, X=1,TI='altitude/m',MAX=8000,TIt='using built-in loop variable $', Y=1,TI='pressure/hPa',MIN=300,MAX=1100)
LINE(AX=4, X=$HP($,1), Y=$HP($,2), Symbol='.', Draw=-9)
END
FUNCTION F4() ! show p(h) array and explain
DLG(Nr=5, L=1/2,T=1/2, Array=$HP, TItle=',height / m,barometric pressure / hPa', SCROLl=100)
DLG(Nr=6, L=0,W=1/2, B='Quit', XEQ='SYSTEM(QUIT=1)', E='#Explain')
END
FUNCTION BarometricPressure(h)
a = 0.0065 ! K/m
T0 = 15 + 273 ! K
M = 0.02896 ! kg/mol
g = 9.81 ! m/s2
R = 8.314 ! J/K/mol = kg*m/s2 /K/mol
P0 = 1013 ! hPa = h N/m2 = h m*kg/s2 /m2
P = P0 * ((1 - a*h/T0) ^ (M*g/R/a))
RETURN P
END