If I were to program it on python, how should I approach it. I want to def temperature(radius) and def scale_height(radius) which output both theta and dz/dR. rather than defining both temperature and scale height as variable theta =....., Tsc = (C_bw/4 ...........)-0.25. Use three of the variable to simulaneously sove for different radial compoent.
