Dirac comb for Riemann zeta function zeros

Clear[n, k, t, A, nn];
nn = 200;
nnn = 60;
A = Table[
   Table[If[Mod[n, k] == 0, 1/(n/k)^(1/2 + I*t - 1), 0], {k, 1,
     nn}], {n, 1, nn}];
MatrixForm[A];
f[t_] = D[RiemannSiegelTheta[t], t]
g1 = ListLinePlot[
  Table[Total[
    1/Table[n*(f[t] + HarmonicNumber[nn]), {n, 1, nn}]*
     Total[Transpose[Re[Inverse[A]*Zeta[1/2 + I*t]]]]], {t, 1/1000,
    nnn, N[1/10]}], DataRange -> {0, nnn}, PlotRange -> {-0, 3},
  PlotStyle -> Thickness[0.003]]