function md=loadfinalcontrolsol(md,solution);

md.cont_vx=solution.vx*md.yts;
md.cont_vy=solution.vy*md.yts;
md.cont_vel=solution.vel*md.yts;
md.cont_J=solution.J;
md.cont_parameters=solution.parameter;
