model;
{
   psi0 <- 1/(1.1*1.1)
   theta0 ~ dnorm(82.4,psi0)
   theta ~ dnorm( 82.4,psi0)
   x ~ dnorm(theta,psi)
   psi <- 1/(1.7*1.7)
   pi <- step(83.0-theta0)
   p <- step(83.0-theta)
}

data;
list(x=82.1)