Esempi in Python per p_T_RH_to_r

Esempio n. 1

File: skewT.py Progetto: cpatrizio88/thermolib

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()

Esempio n. 2

File: qt_vs_theta_alpha.py Progetto: cpatrizio88/thermolib

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()

Esempio n. 3

File: qt_vs_theta_alpha.py Progetto: CCheng1231/A405

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()

Esempio n. 4

File: skewT.py Progetto: CCheng1231/A405

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()

Esempio n. 5

File: qt_vs_theta_l.py Progetto: CCheng1231/A405

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()

Esempio n. 6

File: qt_vs_theta_l.py Progetto: cpatrizio88/thermolib

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()

Esempio n. 7

File: nchaparr_skewT.py Progetto: nchaparr/A405

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()