Virtual Potential Temperature

[GGDLW]

Hello,

I am trying to get virtual potential temperature (Thetav). I have Tv and pressure values.
Would it be possible to get it by using the potential temperature equation substituting Tv instead of Ta?

[dcamron]

Yes, you can do that!

>>> from metpy.calc import potential_temperature, virtual_temperature, virtual_potential_temperature
>>> from metpy.units import units
>>> Ta = 283 * units.K
>>> r = 12 * units('g/kg')
>>> p = 800 * units.mbar
>>> virtual_potential_temperature(p, Ta, r)
<Quantity(303.804317, 'kelvin')>
>>> Tv = virtual_temperature(Ta, r)
>>> Tv
<Quantity(285.039709, 'kelvin')>
>>> potential_temperature(p, Tv)
<Quantity(303.804317, 'kelvin')>