 Solution

Use the electromagnetic field tensor convention stated above, and lower indices with the Minkowski metric $g=\mathrm{diag}(-1,1,1,1)$. Lowering a single time index changes its sign. Summing both orders of each antisymmetric index pair gives the Lorentz scalarFor $T^{00}$, $F_j{}^0=F^{j0}=-E_j/c$ and $g^{00}=-1$. Direct substitution in the printed electromagnetic stress-energy tensor yieldsFor a spatial index $i$, $g^{0i}=0$ and $F_j{}^i=F^{ji}=-{ϵ}_{ijk}{B}_k$. HenceThus the full top row isHere $u_{\mathrm{EM}}$ is the electromagnetic energy density and $\mathbf{S}$ the Poynting vector, using ${ϵ}_0{\mu }_0{c}^2=1$. By the symmetry $T^{0i}=T^{i0}$, these components describe both energy flux divided by $c$ and $c$ times the electromagnetic momentum density $\mathbf{g}=\mathbf{S}/c^2$. Keeping the factor $c$ is essential because the time coordinate is $ct$.