Qmax 1 Vv theta sigma v_vm Qm D Vv tau_v 3600 Qmax 1 Vm theta sigma A v_mv Qv Vm tau_m 3600 mu Qm H chi v_vc C v_vh H c0 omega time_ alpha 2 24