anchor.diff:      ! populate_array will place the non-zero value into a chosen slot at runtime
       coupling_Array = reshape((/0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, &
                                   0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0, 0d0/), (/14, 2/))
  +
  +    ! calculate all the yConversion factors out now (only depends oon the isotope)
  +    do i = 1, 16
  +        yConverse_array(i) = (1.38*10.**-27) * ( ((mnuc*AtomicNumber_oper(i))**(1./3) + 0.33) / (hbar*c0) )**2
  +    end do
   end subroutine captn_init_oper
   
   !   this is the integral over R in eqn 2.3 in 1501.03729
  
positive.code:         ! 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
  
          ! the current response function type in order: M, S2, S1, P2, MP2, P1, D, S1D
          do funcType = 1,8
  
negative.code:                                                         0., 0., 0., 0., 0., 0./) !spins pulled from https://physics.nist.gov/PhysRefData/Handbook/element_name.htm
      double precision :: coupling_Array(14,2)
      double precision :: W_array(8,16,2,2,7)
      double precision :: yConverse_array(16)
  
      integer :: q_shared
      logical :: w_shared
  