model 

{for (i in 1:n)
   {c[i]~dcat(q[])
    t[i]~dgamma(alpha[c[i]],beta[c[i]])
    }

 for (j in 1:2)
   {alpha[j]~dgamma(3,0.1)
    beta[j]<-alpha[j]/lambda[j]
    lambda[j]<-exp(mu[j])
    }

 mumean[1]<-mu0-0.2
 mumean[2]<-mu0+0.2
 for (j in 1:2)
     {mudash[j]~dnorm(mumean[j],3.3) 
      }
 mu[1:2]<-sort(mudash) # This imposes the order constraint.

 mu0~dnorm(4.0,p.mu)
 p.mu<-1/0.3

 pi~dbeta(3,3)
 q[1]<-pi
 q[2]<-1-pi
 }