Added season subset parameter and put failsafe in so probability cannot be 0 or 1.