For recursive models, the tool in PyMC is a scan. It should be able to handle your model just fine, but there are some small details to pay attention to. Examples are here, here, and here, with some small discussion of the finer points here.
One question about the model: what is the law of motion for T_h?