unified_Cfun.RdUnified Circadian Process (C) calculates C (circadian process)
unified_Cfun(tod, phi, tau = 24, A = 1)Time of day (in decimal hours)
Phase at beginning of the simulation (I think this should be 0 if t = tod)
Period of C process
Amplitute of process