def main(tmin, tmax, pmin, pmax): result = get_sounding('sounding.txt') T = result['T'] p = result['p'] RH = result['RH'] r = thermo.p_T_RH_to_r(p, T, RH) T_dew = thermo.T_d(r, p) make_skewT(tmin, tmax, pmax, pmin) tee, pee = skewIt(T, p, 30.) plot(tee, pee, 'k-') tee, pee = skewIt(T_dew, p, 30.) plot(tee, pee, 'k-') axis([tmin, tmax, pmax, pmin]) show()
def _main(tmin, tmax, rmin, rmax): result = _get_sounding('sounding.txt') T = result['T'] p = result['p'] RH = result['RH'] rt = thermo.p_T_RH_to_r(p, T, RH) r_star = thermo.r_star(p, T) rl = rt - r_star rl[rl < 0] = 0. r = rt - rl plot_rt_vs_theta_alpha(p, T, r, rl, 1.) show()
def main(tmin, tmax, rmin, rmax): result = get_sounding('sounding.txt') T = result['T'] p = result['p'] RH = result['RH'] rt = thermo.p_T_RH_to_r(p, T, RH) make_rt_vs_theta_l(270, 310, 0./1000., 6./1000., 84000.) r_star = thermo.r_star(p, T) rl = rt - r_star rl[rl < 0] = 0. r = rt - rl plot_rt_vs_theta_l(p, T, r, rl) axis([270,310,0,6]) show()
def main(tmin, tmax, rmin, rmax): result = get_sounding('sounding.txt') T = result['T'] p = result['p'] RH = result['RH'] rt = thermo.p_T_RH_to_r(p, T, RH) make_rt_vs_theta_l(270, 310, 0. / 1000., 6. / 1000., 84000.) r_star = thermo.r_star(p, T) rl = rt - r_star rl[rl < 0] = 0. r = rt - rl plot_rt_vs_theta_l(p, T, r, rl) axis([270, 310, 0, 6]) show()
def main(tmin, tmax, pmin, pmax): result = get_sounding('sounding.txt') T = result['T'] p = result['p'] RH = result['RH'] r = thermo.p_T_RH_to_r(p, T, RH) indexlim = 47 T_wsat = np.zeros(indexlim) for i in range(indexlim): print ConstantWsat.WsatShift(200, p[i], .01), ConstantWsat.WsatShift(330, p[i], .01), p[i], i T_wsat[i] = ConstantWsat.WsatTemp(p[i], .01) make_skewT(tmin, tmax, pmax, pmin) tee, pee = skewIt(T[0:indexlim,], p[0:indexlim,], 30.) plot(tee, pee, 'r') tee1, pee1 = skewIt(T_wsat[0:indexlim,], p[0:indexlim,], 30.) print "tee", tee, "tee1", tee1 plot(tee1, pee1, 'g') axis([tmin, tmax, pmax, pmin]) show()