Adding 06simulate_mcl_pbs.R which reads PBS variable from the environment to set cores, and a PBS script