code now really uses saturated, 2-phase ullage

introduced new get_P function, completely re-did Vdot solution
function. Added a 6th variable (U_tg).