< y < 8.75$, and $-10 < z < 10$ all in kpc. The simulated $500~{\rm pc}\times 500~{\rm pc}$ region corresponds to $\sim 10''\times 10''$ if we put the galaxy at a distance of 10 Mpc from the Earth (e.g., at Virgo cluster).

We divide the computational region by 32 grids each in $x$ and $y$, and by 1280 grids in $z$. For each grid, we calculate the Faraday depth using Equation (\[eq:FD\]), and calculate the polarized intensity described in Section \[subsection:2.2\]. We integrate the intensity along the LOS and then we obtain the FDF by sorting the intensity as a function of Faraday depth. We also obtain Faraday depth as a function of $z$. The Faraday depth for a grid cell is calculated at the middle of the cell, and the polarized intensity of a cell is an average within the cell. We set the resolution in Faraday depth for the calculation of FDFs to $0.1~{\rm rad/m^2}$, which resolution is sufficiently high to follow the variation of the FDF, so that the results do not change significantly with higher resolution in Faraday depth.

In [@arkg13], the model parameters which determine global configurations of the galaxy such as the scale 