Note

Download full notebook here

Computing Thermal Flux

Computing thermal flux is pretty much the same deal as the reflected light. There are just a couple tweaks we will have to make in order to make sure we have the full set of info.

import warnings
warnings.filterwarnings('ignore')
import numpy as np
import pandas as pd
import astropy.units as u
#picaso
from picaso import justdoit as jdi
from picaso import justplotit as jpi

#plotting
from bokeh.io import output_notebook
from bokeh.layouts import column, row
output_notebook()
from bokeh.plotting import show

Loading BokehJS ...

We will use new Hot Jupiter template to guide us through the exercise.

opa = jdi.opannection()

case1 = jdi.inputs()

case1.phase_angle(0,num_gangle=8, num_tangle=2)


#here we are going to have to specify gravity through R and M since we need it in the Flux calc
case1.gravity(mass=1, mass_unit=u.Unit('M_jup'), radius=1.2, radius_unit=u.Unit('R_jup'))

#here we are going to have to specify R as well
case1.star(opa, 4000,0.0122,4.437,radius=0.7, radius_unit = u.Unit('R_sun') )

#atmo
case1.atmosphere(filename = jdi.HJ_pt(), delim_whitespace=True)
case1.clouds(filename = jdi.HJ_cld(), delim_whitespace=True)

Return PICASO Full Ouput

df= case1.spectrum(opa, full_output=True,calculation='thermal') #note the new last key

wno, fpfs , fp = df['wavenumber'] , df['fpfs_thermal'], df['thermal']
full_output = df['full_output']

Analyzing Thermal Emission Output

All the functionality that we used for reflected light we will also be able to use for thermal emission

Mixing Ratios and Pressure Temperature plots

show(row(jpi.mixing_ratio(full_output), jpi.pt(full_output)))
#can also input any key word argument acceptable for bokeh.figure:
#show(jpi.mixing_ratio(full_output, plot_width=500, y_axis_type='linear',y_range=[10,1e-3]))

Standard Relative Flux Fp/Fs

The same function as reflected light will work with any of the thermal emission output

show(jpi.spectrum(wno,fpfs*1e6,plot_width=500,y_axis_type='log'))

Exploring Thermal Emission Spectrum

This is a useful plot to see the interplay between the computed flux and blackbodies of temperatures at various pressures along the PT profile.

When you specify a pressure here, it will find the corresponding temperature along the PT profile and plot that blackbody.

show(jpi.flux_at_top(full_output, pressures=[1e-3, 1e-1], plot_width=500, plot_height=400))