anchor.diff:                          ((coupling_Array(13,tau).ne.0).and.(coupling_Array(13,taup).ne.0)).or. &
                           ((coupling_Array(14,tau).ne.0).and.(coupling_Array(11,taup).ne.0)).or. &
                           ((coupling_Array(11,tau).ne.0).and.(coupling_Array(14,taup).ne.0))) then
  -                    p_tot = p_tot + RS1(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array)  * w**2 * sumW(w,vesc,i,tau,taup,3,1)&
  -                                - RS1(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array)  * 1./(4.*mu_N**2) * sumW(w,vesc,i,tau,taup,3,2)
  +                    p_tot = p_tot + RS1(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array)  * w**2 * sumW(w,vesc,i,tau,taup,3,1) * 2/3
                   end if
                   ! c15,c15
                   if ((coupling_Array(14,tau).ne.0).and.(coupling_Array(14,taup).ne.0)) then
  -                    p_tot = p_tot + RS1(mnuc,c0,tau,taup,v2q4,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,3,2) &
  -                                - RS1(mnuc,c0,tau,taup,v2q4,j_chi,coupling_Array) * 1./(4.*mu_N**2) * sumW(w,vesc,i,tau,taup,3,3) 
  +                    p_tot = p_tot + RS1(mnuc,c0,tau,taup,v2q4,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,3,2) * 2/3
                   end if
                   
                   ! RP2 (c, q2, q4)
  
positive.code:                 end if
                  ! c8,c8
                  if ((coupling_Array(7,tau).ne.0).and.(coupling_Array(7,taup).ne.0)) then
                      p_tot = p_tot + RM(mnuc,c0,tau,taup,v2,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,1,0) &
                                  - RM(mnuc,c0,tau,taup,v2,j_chi,coupling_Array) * 1./(4.*mu_N**2) * sumW(w,vesc,i,tau,taup,1,1) 
                  end if
                  ! c11,c11
                  if ((coupling_Array(10,tau).ne.0).and.(coupling_Array(10,taup).ne.0)) then
  
negative.code:                 ! c7,c7 or c12,c12
                  if (((coupling_Array(6,tau).ne.0).and.(coupling_Array(6,taup).ne.0)).or. &
                          ((coupling_Array(11,tau).ne.0).and.(coupling_Array(11,taup).ne.0))) then
                      p_tot = p_tot + RS1(mnuc,c0,tau,taup,v2,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,3,0) * 2/3
                  end if
                  ! c9,c9
                  if ((coupling_Array(8,tau).ne.0).and.(coupling_Array(8,taup).ne.0)) then
  