I added changes to allow calculation of reflected light phase curves. Most notable changes are confined to atmosphere_4d and clouds_4d in justdoit.py. I also assumed symmetry across phase = 0 deg (secondary eclipse) in fluxes.py > get_reflected_3d in order to calculate reflected phase intensity from 0-360 degrees.
I added changes to allow calculation of reflected light phase curves. Most notable changes are confined to atmosphere_4d and clouds_4d in justdoit.py. I also assumed symmetry across phase = 0 deg (secondary eclipse) in fluxes.py > get_reflected_3d in order to calculate reflected phase intensity from 0-360 degrees.