Automated adjoints of coupled PDE-ODE systems