model;
{
for (i in 1:k) 
  {
    t[i] ~ dnorm(0, tau)  
  }

  tau<-0.0001*phi

  phi ~ dgamma(0.0001, 0.0001)  

  for (i in 1:n)
  {
    mean[i] <- t[location[i]]  
    weight[i] ~ dnorm(mean[i], phi)  
  }
}