Bond Albedo #358
Description
The bond albedo we provide to users ignores the phase integral q (as pointed out by Vighnesh Nagpal at University of Chicago vnagpal@uchicago.edu -- advised by Diana Powell). Through some back and forth we decided that there should maybe be a more direct way of getting the full phase integral + associated bond albedo. The code snippet attached provides an initial schematic to get the phase integral from compute a series of albedos at varying phases -- following what I think Cahoy did in 2010 paper. Work is being done by Vighnesh to ensure that PICASO is properly reproducing e.g., Rayleigh phase integrals from the code snippet attached. Then we can move toward having a bond_albedo function in justdoit.py. Should also think about the 4d phase curve code because that should be fairly straightforward to add in the bond albedo plug in.
- provide a function to compute bond albedos for 1D
- provide a function to compute bond albedos for 3D
- add documentation about the different albedo returns