anchor.diff:  
   function gauss_laguerre_5pt_capture(fpar) result(igral)
     !$acc routine seq
  +  !$acc data present(gauss_laguerre_5pt_xk, gauss_laguerre_5pt_wk)
     ! Do 5-pt Gauss Laguerre integration
     ! Weight function is w(x) = x^0 * exp(-x)
     use bl_types
  
positive.code: 
  function gauss_legendre_5pt_emission(fpar) result(igral)
    !$acc routine seq
    ! Do 5-pt Gauss Legendre integration
    ! Weight function is w(x) = 1
    use bl_types
  
negative.code: 
    igral = 0.0d0
    do j=1,N
      wkj = gauss_legendre_5pt_wk(j)
      xkj = gauss_legendre_5pt_xk(j)
      igral = igral + wkj*phase_emission_ne23(xkj, fpar)
    end do
    return
    !$acc end data
  end function gauss_legendre_5pt_emission
  
  function rate_capture_na23(rpar) result (lambda)
  