anchor.diff:  #endif
     implicit none
     integer,intent(in)              :: ik
  -  complex(8),intent(in)           :: tpsi(0:PNLz-1,0:PNLy-1,0:PNLx-1)
  +  complex(8),intent(in)           ::  tpsi(0:PNLz-1,0:PNLy-1,0:PNLx-1)
     complex(8),intent(out)          :: htpsi(0:PNLz-1,0:PNLy-1,0:PNLx-1)
     complex(8),intent(out),optional :: ttpsi(0:PNLz-1,0:PNLy-1,0:PNLx-1)
     real(8) :: k2,k2lap0_2
  
positive.code:   end subroutine
  end subroutine
  
  
  subroutine hpsi_omp_KB_RT(ik,tpsi,htpsi,ttpsi)
    use timer
    use Global_Variables, only: NLx,NLy,NLz,kAc,lapx,lapy,lapz,nabx,naby,nabz,Vloc,Mps,uV,iuV,Hxyz,ekr_omp,Nlma,a_tbl
    use opt_variables, only: lapt,PNLx,PNLy,PNLz,PNL
  
negative.code: 
    tid = omp_get_thread_num()
    call init(tpsi,zhtpsi(:,1,tid))
    call hpsi_omp_KB_base(ik,zhtpsi(:,1,tid),zhtpsi(:,2,tid),zttpsi(:,tid))
    call copyout(zhtpsi(:,2,tid),zttpsi(:,tid),htpsi,ttpsi)
  
  contains
  