************************************************************************ * * SetEWCouplings.f: * * This subroutine sets the values of the vector and axial couplings * of the up- and down-type quarks. * ************************************************************************ subroutine SetEWCouplings(vd,vu,ad,au) * implicit none * include "../commons/EWCouplings.h" * * Variables * double precision vd,vu,ad,au * VectorD = vd VectorU = vu AxialD = ad AxialU = au * * If they are all zero, tell the code to use the standard couplings * ExtCoup = .true. if(vd.eq.0d0.and.vu.eq.0d0.and. 1 ad.eq.0d0.and.au.eq.0d0) ExtCoup = .false. * InEWCouplings = "done" * return end