tauCO = 70*24*3600 # s (wright et al. 2018 -- calibrated for mid M dwarfs, likeley not correct for Brown Dwarfs)
offset = 300 # to allow for thermal transfer time # need to connect this to the convective overturn timescale or better yet dont be stupid about it
coldSideTemp = np.zeros_like(orbits) + BDTEMP
for oID, orbit in tqdm(enumerate(orbits[offset:]), total=orbits[offset:].shape[0]):
    delayedHotSideTemp = hst(oID*period)
    coldSideTemp[oID + offset] = delayedHotSideTemp
    coldSideTemp[oID + offset] = cool_blackbody(BDRAD, coldSideTemp[oID + offset], bdmass, orbitJump * period, internalEnergy)