Combining two distributions

You can save the rounded values in your model by wrapping them in a Deterministic:

indexes_continuous = pm.Normal(...)
indexes = pm.Deterministic('indexes', tt.round(indexes_continuous))

Then you can use indexes in your model and also check them in your trace without having to round a posteriori.