anchor.diff:       endif
     enddo
   
  +  ierr=798
  +  call mpi_comm_free(CT_COMM_LOCAL,ierr)
  +  if (ierr.ne.0) then
  +     write(mpifileptr,*) "Error comm destroy simple_circ",ierr; call mpistop()
  +  endif
  +
  +  ierr=798
  +  call mpi_group_free(CT_GROUP_LOCAL,ierr)
  +  if (ierr.ne.0) then
  +     write(mpifileptr,*) "Error group destroy simple_circ",ierr; call mpistop()
  +  endif
  +
   end subroutine simple_circ
   
   
  
positive.code:   endif
  
    call myzfft1d_slowindex_mpi(in,tempout,pf(1),ctrank,ctset,dim1*howmany)
    call twiddlemult_mpi(tempout,outtemp,dim1,depth,newproclist,pf(1),ctrank,howmany)
    if (depth.eq.1) then
       call myzfft1d(outtemp,outtrans,dim1,howmany)
    else
  
negative.code:      out(:)=out(:)+work(:)*mat(ctrank,ibox)
    enddo
  
    ierr=798
    call mpi_comm_free(CT_COMM_LOCAL,ierr)
    if (ierr.ne.0) then
       write(mpifileptr,*) "Error comm destroy simple_circ",ierr; call mpistop()
    endif
  
    ierr=798
    call mpi_group_free(CT_GROUP_LOCAL,ierr)
    if (ierr.ne.0) then
       write(mpifileptr,*) "Error group destroy simple_circ",ierr; call mpistop()
    endif
  
  end subroutine simple_summa
  