A perturbative QCD based jet tomographic Monte Carlo model, CUJET2.0, is presented to predict jet quenching observables in relativistic heavy ion collisions at RHIC/BNL and LHC/CERN energies. This model generalizes the DGLV theory of flavor dependent radiative energy loss by including multi-scale running strong coupling effects. It generalizes CUJET1.0 by computing jet path integrations though more realistic 2+1D transverse and longitudinally expanding viscous hydrodynamical fields contrained by fits to low $p_T$ flow data. The CUJET2.0 output depends on three control parameters, $(alpha_{max},f_E,f_M)$, corresponding to an assumed upper bound on the vacuum running coupling in the infrared and two chromo-electric and magnetic QGP screening mass scales $(f_E mu(T), f_M mu(T))$ where $mu(T)$ is the 1-loop Debye mass. We compare numerical results as a function of $alpha_{max}$ for pure and deformed HTL dynamically enhanced scattering cases corresponding to $(f_E=1,2, f_M=0)$ to data of the nuclear modification factor, $R^f_{AA}(p_T,phi; sqrt{s}, b)$ for jet fragment flavors $f=pi,D, B, e$ at $sqrt{s}=0.2-2.76$ ATeV c.m. energies per nucleon pair and with impact parameter $b=2.4, 7.5$ fm. A $chi^2$ analysis is presented and shows that $R^pi_{AA}$ data from RHIC and LHC are consistent with CUJET2.0 at the $chi^2/d.o.f< 2$ level for $alpha_{max}=0.23-0.30$. The corresponding $hat{q}(E_{jet}, T)/T^3$ effective jet transport coefficient field of this model is computed to facilitate comparison to other jet tomographic models in the literature. The predicted elliptic asymmetry, $v_2(p_T;sqrt{s},b)$ is, however, found to significantly underestimated relative to RHIC and LHC data. We find the $chi^2_{v_2}$ analysis shows that $v_2$ is very sensitive to allowing even as little as 10% variations of the path averaged $alpha_{max}$ along in and out of reaction plane paths.