TIME DEPENDENT HADRONIC MODELING OF FLAT SPECTRUM RADIO QUASARS