/** Calculate boiling point of water (and atmosphere density) by elevation. */

use StandardAtmosphere.frink
use water.frink

for elevation = concat[0 ft to 15000 ft step 500 ft, 16000 ft to 29000 ft step 1000 ft]
{
   [t, p, d] = StandardAtmosphere.getTPD[elevation]
   b = Water.boilingPoint[p]
   println[format[elevation, "feet", 0] + "\t" + format[elevation, "m", 0] + "\t" + format[F[b], 1, 2] + " F" + "\t"+ format[C[b],1,2] + " C" + "\t" + format[p, "millibar", 0]]
}