anchor.diff:      double precision :: ier,alist,blist,rlist,elist,iord,last!for integrator
       double precision, allocatable :: u_int_res(:)
   
  -    ! ! stuff for openmp testing
  -    ! integer :: threadnums, threadid
  -    ! integer, external :: omp_get_num_threads, omp_get_thread_num
  -
       dimension alist(1000),blist(1000),elist(1000),iord(1000),rlist(1000)!for integrator
       !external gausstest !this is just for testing
       external integrand_oper
  
positive.code:             end if
              tally = tally + W_array(Wtype,iso,tau,tauprime,k) * yConverse**(k-1) * G
          end do
          sumW = tally
      end function sumW
      
  
negative.code:                     p_tot(w,vesc,pickIsotope)
              end if
          else ! if all the isotopes are being run
              !!$OMP parallel default(none) shared(niso, mdm, u, w, tab_starrho, rindex, tab_mfr_oper, vesc) &
              !!$OMP private(i, J, mu, muplus)
              !!$OMP do
              do i = 1,niso
                  J = AtomicSpin_oper(i)
                  mu = mdm/mnuc/AtomicNumber_oper(i)
  