Question

Step I: Put the degree 24 polynomial into Magma, make it a number field, and call LSeries on it. This divides the $L$-function into a product of 7 distinct ones (Dokchitsers code, under an attribute called "prod" on the L-series object), given by Artin representations. So my plan was to compute zeros for each of these $L$-functions, with them being a subset of those of $\zeta_M$ naturally, and recombine. As representations for $SL(2,3)$ this is as $$1\oplus \omega\oplus\bar\omega\oplus 2\tau_2\oplus 2\tau_2\omega\oplus 2\tau_2\bar\omega\oplus 3\kappa_3.$$ Or as an $L$-function product $$\zeta_M=\zeta\cdot L(\omega)\cdot L(\bar\omega)\cdot L(\tau_2)^2\cdot L(\tau_2\omega)^2\cdot L(\tau_2\bar\omega)^2\cdot L(\kappa_3)^3.$$

EDIT: Oh I see now, you wanted it for $M$ w/o assuming analyticity of Artin constituent parts, but Magma automatically dissembles it, and assumes. But now at this point, I see you were trying to avoid this decomposition perhaps, but then I really don't see how you could work with $\zeta_M$ directly, as the conductor is too big, being $163^{16}$. From the standpoint of scientific evidence, it is probably enough that CheckFunctionalEquation in Magma for each of the above constituents gives an answer indistinguishable from zero, though the feel of zeros is also nice.

Another way to test individual analyticity is to use subfields. There is a unique (up to isomorphism) quartic subfield $K_4\subset M$, with $$\zeta_K=\zeta\cdot L(\kappa_3).$$ The conductor of $\zeta_K$ is small enough to compute with it directly (no decomposition), so the analyticity of $L(\kappa_3)=\zeta_K/\zeta$ can be checked (numerically) by seeing if Riemann zeros are a subset of those of $\zeta_K$. EDIT: Of course, $\kappa_3$ is induced from the a nontrivial linear character of the $Q_8$ subgroup of $SL(2,3)$, so analyticity already follows by induction.

There is also a unique (up to isomorphism) sextic subfield $L_6\subset M$ with $$\zeta_L=\zeta\cdot L(\omega)\cdot L(\bar\omega)\cdot L(\kappa_3),$$ so the linear parts can be isolated here by $\zeta_L/\zeta_K$, though maybe superfluous ("easy" theorem)? In that matter, already the cubic subfield $C_3\subset M$ has $\zeta_C=\zeta\cdot L(\omega)\cdot L(\bar\omega)$ if desired.

There is a unique (up to isomorphism) octic subfield $N_8\subset M$ with $$\zeta_N=\zeta\cdot L(\tau_2\omega)\cdot L(\tau_2\bar\omega)\cdot L(\kappa_3).$$ The integer ring has discriminant $163^4$ so $\zeta_N$ is likely still computable directly. By quotient $\zeta_N/\zeta_K$, this can tell about the product of conjugate 2-dimensional $L$-functions. EDIT: see below for another idea, using twists of $\zeta_N$.

Finally the unique duodecic subfield $R_{12}\subset M$ has $$\zeta_R=\zeta\cdot L(\omega)\cdot L(\bar\omega)\cdot L(\kappa_3)^3.$$ So this gives nothing new, and the discriminant is too large anyway.

Note that none of the parts has $L(\tau_2)$ directly, only $\zeta_M$ itself that is too hard to compute. EDIT: However, by Rankin-Selberg I think(?) it follows that the analyticity of $L(\tau_2)$ is equivalent to that after twisting by $\omega$ to get $L(\tau_2\omega)$.

Answer?: For that matter, recurring to the octic subfield $N_8$, instead of using just the Dedekind $\zeta$-function of $N_8$, twisting it by $\omega$ could be profitable, achieving $$L(\sigma_8\omega)=L(\omega)\cdot L(\tau_2\bar\omega)\cdot L(\tau_2)\cdot L(\kappa_3),$$ $$L(\sigma_8\bar\omega)=L(\bar\omega)\cdot L(\tau_2)\cdot L(\tau_2\omega)\cdot L(\kappa_3),$$ where $\sigma_8=1\oplus\tau_2\omega\oplus\tau_2\bar\omega\oplus\kappa_3$ is for the Dedekind of $N_8$. Also $\kappa_3\omega=\kappa_3$ nicely. Noting here the twisting gives $$L(\sigma_8)L(\sigma_8\omega)L(\sigma_8\bar\omega)=\zeta_M,$$ this could provide a way to compute $\zeta_M$, as each part of the left is known as analytic I presume. The conductor of $L(\sigma_8)$ is $163^4$, and that for the twists is $163^6$. Or one can avoid $\zeta_M$ alternatively, as my computation is $$L(\tau_2)^2={L(\sigma_8\omega)L(\sigma_8\bar\omega)\zeta^2\over\zeta_{N_8}\zeta_{L_6}},$$ where each factor on the right should be known (easily?) to be holomorphic away from the $\zeta$-pole. Note this exploits the linear characters of $SL(2,3)$, and you have none for your next case of $SL(2,5)$.

In all cases, zeros need to be computed, and the right tool is L-calc. But I don't know if it is really feasible to go too far for $L(\sigma_8\omega)$ of conductor $163^6$, without intensive effort.

Part II: Zeros of $L(\tau_2)$ computation (example): I compute its first few zeros for the 2-dimensional representation of real character. For this representation, with $10^5$ coeffs (taking 6sec in Magma), I exported these to Lcalc (a somewhat hackish tool of M. O. Rubinstein, included in Sage I suppose but I did it direct), which returns after 12 seconds, listing the imaginary parts of the first 10 zeros on the half-line (also proving RH up to that point by A Turing method):

 0
 0.99014365233
 1.38830360231
 2.35103235859
 3.45296629741
 4.32708276131
 4.73989005257
 5.42392092883
 5.99574967707
 6.70167394143

The first zero is at $s=1/2$ as the sign is $-1$. The second is about $s\approx 1/2+0.990i$

For reference, here are the L-calc settings I used, hackish as I say:

1
0
100000
0
2
.5
0.5 0.0
.5
0.5 0.0
51.884511447957879460656106859439682023
-1.0 0.0
0

As explained in their help, the first "1" says the coeffs are integers, the second "0" says they are not special, the 3rd "100000" says they are that many, the 4th "0" says not periodic, the 5th "2" says degree 2, the 6th ".5" and "0.5 0.0" say the form of the gamma factor, the "51.884" is $\sqrt{163^2/\pi^2}$ as the analytic conductor, the "-1.0 0.0" is the sign, and the final "0" says no poles. Then the 100000 coeffs are given as integers.