  In this paper we are concerned with $N$-urn susceptible-infected-removed
epidemics, where each urn is in one of three states, namely `susceptible',
`infected' and `removed'. We assume that recovery rates of infected urns and
infection rates between infected and susceptible urns are all
coordinate-dependent. We show that the hydrodynamic limit of our model is
driven by a deterministic $(C[0, 1])^\prime$-valued process with density which
is the solution to a nonlinear $C[0, 1]$-valued ordinary differential equation
consistent with a mean-field analysis. We further show that the fluctuation of
our process is driven by a generalized Ornstein-Uhlenbeck process. A key step
in proofs of above main results is to show that sates of different urns are
approximately independent as $N\rightarrow+\infty$.
