Here it is from Dokchitser in Magma:
> L:=LSeries(HilbertClassField(QuadraticField(145)) : Method:="Artin");
> L`prod;
[
<L-series of Riemann zeta function, 1>,
<L-series of Artin representation of Number Field with defining polynomial x^8 -
636*x^6 + 135214*x^4 - 11109740*x^2 + 290532025 over the Rational Field with
character ( 1, 1, 1, -1, -1 ) and conductor 5, 1>,
<L-series of Artin representation of Number Field with defining polynomial x^8 -
636*x^6 + 135214*x^4 - 11109740*x^2 + 290532025 over the Rational Field with
character ( 1, 1, -1, -1, 1 ) and conductor 145, 1>,
<L-series of Artin representation of Number Field with defining polynomial x^8 -
636*x^6 + 135214*x^4 - 11109740*x^2 + 290532025 over the Rational Field with
character ( 1, 1, -1, 1, -1 ) and conductor 29, 1>,
<L-series of Artin representation of Number Field with defining polynomial x^8 -
636*x^6 + 135214*x^4 - 11109740*x^2 + 290532025 over the Rational Field with
character ( 2, -2, 0, 0, 0 ) and conductor 145, 2>
]
> L5:=L`prod[5][1];
> CheckFunctionalEquation(L5); // LCfRequired(L5) demands 161 terms
1.57772181044202361082345713057E-30
> [<p,-Integers()!Coefficient(EulerFactor(L5,p),1)> : p in PrimesUpTo(100)];
[ <2, 0>, <3, 0>, <5, -1>, <7, 0>, <11, 0>, <13, 0>, <17, 0>, <19, 0>, <23, 0>,
<29, -1>, <31, 0>, <37, 0>, <41, 0>, <43, 0>, <47, 0>, <53, 0>, <59, -2>, <61, 0>,
<67, 0>, <71, -2>, <73, 0>, <79, 0>, <83, 0>, <89, 0>, <97, 0> ]
These are the same as your ap(p), essentially.