barometric_formula.nbt 502 B

1234567891011121314
  1. # This script calculates the air pressure at a specified
  2. # height above sea level using the barometric formula.
  3. let p0: Pressure = 1 atm
  4. let t0: Temperature = 288.15 K
  5. dimension TemperatureGradient = Temperature / Length
  6. let lapse_rate: TemperatureGradient = 0.65 K / 100 m
  7. fn air_pressure(height: Length) -> Pressure =
  8. p0 · (1 - lapse_rate · height / t0)^5.255
  9. print("Air pressure 1500 m above sea level: {air_pressure(1500 m) -> hPa}")
  10. assert_eq(air_pressure(1500 m), 845.586 hPa, 0.1 hPa)