anchor.diff:                  psum_int(1:nk) = psum_int(1:nk) + int_i(1:nk, j, i)
                   psum_int_sq(1:nk) = psum_int_sq(1:nk) + int_sq(1:nk, j, i)
               enddo
  +            if(i .eq. 37) then
  +                write(*,*) "Finished intensity calls on model", i
  +                write(*,*) int_sq(1:nk, 1, i)
  +            endif
   #ifdef DEBUG
               write(*,*) "Finished intensity calls on model", i
   #endif
  
positive.code:                     int_sq(1:nk, m, i) = int_i(1:nk, m, i)**2
                  endif
              enddo
          enddo
  #ifdef DEBUG
          write(*,*) "I am core", myid, "and I am past the int calls."
  
negative.code:         bin_max = int(r_max/fem_bin_width)+1
  
          allocate(gr_i(m_int%nelements,m_int%nelements, 0:bin_max), stat=istat)
          allocate(x1(size_pix_atoms), y1(size_pix_atoms),rr_a(size_pix_atoms), stat=istat)
          allocate(sum1(m_int%nelements,size_pix_atoms), stat=istat)
          allocate(znum_r(size_pix_atoms), stat=istat)
  
          gr_i = 0.0; int_i = 0.0; x1 = 0.0; y1 = 0.0; rr_a = 0.0; sum1 = 0.0; znum_r = 0.0
  
          do i=1, size_pix_atoms
              znum_r(i) = m_int%znum_r%ind(pix_atoms(i))
          enddo
  
          x2 = 0.0; y2 = 0.0
          const1 = twopi*(0.61/res)/fem_bin_width  ! (0.61/res = Q)
  