1      subroutine tce_sortacc_2(unsorted,sorted,a,b,i,j,factor)
2      implicit none
3      integer a,b
4      integer i,j
5      integer id(2),jd(2),ia,ib,j1,j2
6      double precision sorted(a*b)
7      double precision unsorted(a*b)
8      double precision factor
9      ia = 0
10      jd(1) = a
11      jd(2) = b
12      do j1 = 1,a
13       id(1) = j1
14       do j2 = 1,b
15        id(2) = j2
16        ia = ia + 1
17        ib = id(j)+jd(j)*(id(i)-1)
18        sorted(ib) = sorted(ib) + unsorted(ia) * factor
19       enddo
20      enddo
21      return
22      end
23