This code was put together to automate the PAML-codeml package (Yang, 2007). It requires BioPython.
This code works well if you have many (i.e. tens-thousands) of genes/alignments that you want to test for positive selection or fit any codeml model to.