g_i) + (\alpha_i  + \frac{5}{3}\tau \epsilon_n g_i)k_{\perp}^2 \right) & + & \nonumber \\
\omega\left( \tau \epsilon_n g_i k_y^2(\eta_i - \frac{7}{3} + \frac{5}{3}(1+\tau)\epsilon_n g_i + \frac{5}{3}k_{\perp}^2\alpha_i )- k_{\parallel}^2 (1 +  \frac{5}{3} \tau) \right) & - & \nonumber \\
k_{\parallel}^2  k_y \tau \left( (\eta_i -  \frac{2}{3}) +  \frac{5}{3} (1+\tau)\epsilon_n g_i \right) & = & 0. \nonumber \\\end{aligned}$$ Here the fields for the linear ITG mode are expanded according to $\tilde{\xi}_{DW}(x,y,t)  =  \xi_{DW} e^{i(k_x x + k_y y + k_{\parallel} z - \omega t)} + c.c.$ for ($\tilde{n}_i,\tilde{\phi},\tilde{T}_i, \tilde{v}_{i ||}$). The dispersion relation for the linear ITG mode including parallel ion motions effects is of 3rd order, however if $k_{\parallel}^2$ effects are neglected the dispersion relation is transformed into the previous result [@a23]. The effects of parallel ion motion on the real frequency and the growth rates are small. The important differences are that the real frequencies and growth rates are significantly changed for small toroidicity and in addition the $\eta_i$-threshold is decreased. In this particular case it is of interest to monitor the zonal flow generation from toroidal and slab like ITG mode turbulence. In the limit of small $\epsilon_n=2 L_n/R$ and small Finite Larmor Radius (FLR) effects an approximate solution may 