anchor.diff:          ! I'll need mu_T to include in the prefactor when there is a v^2 term
           a = AtomicNumber_oper(eli)
           mu_T = (mnuc*a*mdm)/(mnuc*a+mdm)
  -        yConverse = (1.38*10.**-27) * ( ((mnuc*a)**(1./3) + 0.33) / (hbar*c0) )**2
  +        ! yConverse = (1.38*10.**-27) * ( ((mnuc*a)**(1./3) + 0.33) / (hbar*c0) )**2
   
           ! the current response function type in order: M, S2, S1, P2, MP2, P1, D, S1D
           do funcType = 1,8
  
positive.code:     ! specific to captn_oper
      integer :: funcType, tau, taup, term_R, term_W, q_pow, w_pow ! loop indicies
      integer :: q_functype, q_index
      double precision :: J, j_chi, RFuncConst, WFuncConst, yConverse, mu_T, prefactor_functype
      double precision :: RD, RM, RMP2, RP1, RP2, RS1, RS1D, RS2 !R functions stored in their own source files
      double precision :: prefactor_array(niso,11,2)
  
  
negative.code: 
                                      ! calculates the total number of q^2
                                      q_index = 1 + q_functype + term_W - 1 + floor((term_R-1.)/2.)
                                      prefactor_current = prefactor_functype*RFuncConst*WFuncConst*yConverse_array(eli)**(term_W-1)
  
                                      ! check if term_R is even (in my index convention this corresponds to it having a v^2 in the term)
                                      ! v^2 = w^2 - q^2/(2mu_T)^2
                                      if ( mod(term_R,2).eq.0 ) then
                                          ! this is the -q^2/(2mu_T)^2 contribution
                                          prefactor_array(eli,q_index+1,1) = prefactor_array(eli,q_index+1,1) - prefactor_current * &
                                              (c0**2/(4.*mu_T**2)) ! The Rfunctions are programmed with the 1/c0^2 in their v_perp^2 term (so I need to un-correct it for the- q^2/(2*mu_T)^2, and leave it be for the w^2/c^2)
                                          ! this is the +w^2 contribution
                                          prefactor_array(eli,q_index,2) = prefactor_array(eli,q_index,2) + prefactor_current
                                          
                                      else
                                          prefactor_array(eli,q_index,1) = prefactor_array(eli,q_index,1) + prefactor_current
  
                                      end if
                                  end if
  