external/io_netcdfpar/transpose.code

   1   ix=0
   2   jx=0
   3   kx=0
   4   call reorder(MemoryOrder,MemO)
   5   if(IACHAR(MemO(1:1)) > MaxUpperCase) ix=i2+i1
   6   if(IACHAR(MemO(2:2)) > MaxUpperCase) jx=j2+j1
   7   if(IACHAR(MemO(3:3)) > MaxUpperCase) kx=k2+k1
   8
   9 ! pjj/cray
  10   if(IO == 'write') then
  11 !dir$ concurrent
  12 !$OMP PARALLEL DO SCHEDULE(RUNTIME) PRIVATE(i,j,k)
  13   do k=k1,k2
  14     do j=j1,j2
  15 !dir$ prefervector
  16 !dir$ concurrent
  17 !cdir select(vector)
  18       do i=i1,i2
  19           DFIELD = Field(1:di,abs(ix-i),abs(jx-j),abs(kx-k))
  20       enddo
  21     enddo
  22   enddo
  23 !$OMP END PARALLEL DO
  24 else
  25 !dir$ concurrent
  26 !$OMP PARALLEL DO SCHEDULE(RUNTIME) PRIVATE(i,j,k)
  27   do k=k1,k2
  28     do j=j1,j2
  29 !dir$ prefervector
  30 !dir$ concurrent
  31 !cdir select(vector)
  32       do i=i1,i2
  33           Field(1:di,abs(ix-i),abs(jx-j),abs(kx-k)) = DFIELD
  34       enddo
  35     enddo
  36   enddo
  37 !$OMP END PARALLEL DO
  38 endif
  39
  40   return