anchor.diff:                  end if
                   ! c5,c5
                   if ((coupling_Array(4,tau).ne.0).and.(coupling_Array(4,taup).ne.0)) then
  -                    p_tot = p_tot + RM(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,1,1) * 2/3
  +                    p_tot = p_tot + RM(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,1,1) &
  +                                - RM(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array) * 1./(4.*mu_N**2) * sumW(w,vesc,i,tau,taup,1,2) 
                   end if
                   
                   ! RS2 (c, v2, q2, v2q2, q4)
  
positive.code:                         ((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) * 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) * 2/3
                  end if
                  
                  ! RP2 (c, q2, q4)
  
negative.code:                 end if
                  ! c13,c13
                  if ((coupling_Array(12,tau).ne.0).and.(coupling_Array(12,taup).ne.0)) then
                      p_tot = p_tot + RS2(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array) * w**2 * sumW(w,vesc,i,tau,taup,2,1) &
                                  - RS2(mnuc,c0,tau,taup,v2q2,j_chi,coupling_Array) * 1./(4.*mu_N**2) * sumW(w,vesc,i,tau,taup,2,2)
                  end if
                  ! c6,c6
                  if ((coupling_Array(5,tau).ne.0).and.(coupling_Array(5,taup).ne.0)) then
  