In your pyxtra.f, add the following just before "call pyinit": C Leptonic decays of taus only do i=91,138 mdme(i,1)=0 enddo