on a qrsh ~18G it takes about 20 minutes. When I can't get a node with that kind of memory its a few hours. I've been dividing up the iterations as separate qsubs which has been helping. Simplifying it to a right fischer model also greatly helps with speed..( constant population size) but this assumption has drawbacks in modeling true variation seen in pandemic outbreaks where you obviously don't have a constant population size.