Function temp_to_theta, t, p ; Function to convert T (K) and P (hPa) to theta ; using a reference pressure of 1000 hPa RoverCp=0.286 p0=1000. theta=t*(p0/p)^RoverCp return, theta END