hi jesseb,
i've tried your script in the new 1.50.2 release and it seems to work fine.
you could add an option to change the analysis method just before the line of the compute() command. these are the current options (you have first to deactivate the bishop method as actually it is the default method):
set("Method","BishopSim","off","Method","GLE/M-P","on")
or
set("Method","BishopSim","off","Method","JanbuSim","on")
or
set("Method","BishopSim","off","Method","Spencer","on")