def p2h(p, T=293., P0=1000., m=28.966, unit_p='mbar'): """ Returns an elevation from barometric pressure Parameters ---------- p: {float, array} barometric pressure in mbar or torr specified with unit_p T: float, optional Temperature in K P0: float, optional Pressure at reference altitude in hPa (default = 1000.) m: float, optional average mass of gas molecules in u (default = 28.966) unit_p: {[mbar], torr}, optional Source ------ http://en.wikipedia.org/wiki/Barometric_formula """ if unit_p == 'torr': p = unit_conversion.torr2mbar(p) k = const.physical_constants['Boltzmann constant'][0] g = const.physical_constants['standard acceleration of gravity'][0] m *= 1 / const.physical_constants['Avogadro constant'][0] / 1000. h = (np.log(P0) - np.log(p)) * ((k * T) / (m * g)) return h
def p2h(p,T = 293.,P0 = 1000., m = 28.966, unit_p = 'mbar'): """ parameters: h: altitude from reference higth in meter (if reference hight is not sealevel, you probably want to adjust P0) optional parameters: T: Temperature in K P0: Pressure at reference altitude in hPa (default = 1000.) m: average mass of gas molecules in u (default = 28.966) unit_p: either 'mbar' or 'torr' source: http://en.wikipedia.org/wiki/Barometric_formula """ if unit_p == 'torr': p = unit_conversion.torr2mbar(p) k = const.physical_constants['Boltzmann constant'][0] g = const.physical_constants['standard acceleration of gravity'][0] m *= 1/const.physical_constants['Avogadro constant'][0]/1000. h = (np.log(P0) - np.log(p)) * ((k * T)/(m * g)) return h
def p2h(p, T=293., P0=1000., m=28.966, unit_p='mbar'): """ parameters: h: altitude from reference higth in meter (if reference hight is not sealevel, you probably want to adjust P0) optional parameters: T: Temperature in K P0: Pressure at reference altitude in hPa (default = 1000.) m: average mass of gas molecules in u (default = 28.966) unit_p: either 'mbar' or 'torr' source: http://en.wikipedia.org/wiki/Barometric_formula """ if unit_p == 'torr': p = unit_conversion.torr2mbar(p) k = const.physical_constants['Boltzmann constant'][0] g = const.physical_constants['standard acceleration of gravity'][0] m *= 1 / const.physical_constants['Avogadro constant'][0] / 1000. h = (np.log(P0) - np.log(p)) * ((k * T) / (m * g)) return h