model {
for (j in 1:nsprings) {
 springVelMean[j] ~ dnorm(9,0.5)
 springVelVar[j] ~ dexp(20)
}

for (i in 1:40) {
 spring[i] ~ dcat(pi)
 vel[i] ~ dnorm(springVelMean[spring[i]], springVelVar[spring[i]])
}
springOut ~ dcat(pi)
velOut ~ dnorm(springVelMean[springOut], springVelVar[springOut])
}