model{
for (i in 1:T){
	mu[i] <- m + alpha[halfcent[i] ]
	temp[i] ~ dnorm(mu[i],tau) 
}
alpha[1] <- -sum(alpha[2:H])

#priors
m~dnorm(0,0.001)
for (i in 2:H){alpha[i]~dnorm(0,0.001)}
tau~dgamma(0.01,0.01)
s<-sqrt(1/tau)
}