A C1-virtual element method of high order for the Brinkman equations in stream function formulation with pressure recovery. (22nd December 2021)