anchor.diff:      end subroutine read_evops
   
       subroutine write_evops(code)
  -      integer (i4b)          :: i, file_ios
  +      integer (i4b)          :: i, file_ios, LIMIT
         integer(i4b)           :: Uelement, Uelement2,Utnemele,Utnemele2
         character, intent(in)  :: code
   
  -      do i=1,MAXOUTIND
  +      LIMIT = MAXOUTIND
  +
  +      if(code == 'w') then
  +              LIMIT = size(Evops,5)
  +      end if
  +
  +      do i=1,LIMIT
         do Uelement=1,Nl1
         do Uelement2=1,Nl2
         do Utnemele=1,Nl1
  
positive.code:     subroutine renormalize_evops()
          integer(i4b) tind, i,j,k,l
          complex(dpc), allocatable, dimension(:,:,:,:) :: EvopsTmp
  
          if(type == 'E') then
  
  
negative.code:                 cycle
              end if
  
              do l=1, N1_from_type(type)
              do k=1, N1_from_type(type)
  
              if((i < j .and. k >= l) .or. (j < i .and. l >= k) ) then
                  cycle
              end if
  
  
              EvopsTmp(l,k,j,i) = EvopsTmp(l,k,j,i) + conjg(Evops(k,l,j,i,tind))
              EvopsTmp(k,l,j,i) = 0.0_dp
  
              if(tmp(k,l,j,i) == 0) then
                  write(*,*) 'removing Evops', k,l,j,i
                  tmp(k,l,j,i) = 1
              end if
  
              end do
              end do
  
              do k=1, N1_from_type(type)
                  EvopsTmp(k,k,j,i) = real(EvopsTmp(k,k,j,i))
  
                  if(tmp(k,k,j,i) == 0) then
                          write(*,*) 'realing Evops', k,k,j,i
                          tmp(k,k,j,i) = 1
                  end if
              end do
  
  
  