anchor.diff:        end do
       end subroutine read_evops
   
  -    subroutine write_evops()
  -      integer (i4b)       :: i, file_ios
  -      integer(i4b)        :: Uelement, Uelement2,Utnemele,Utnemele2
  +    subroutine write_evops(code)
  +      integer (i4b)          :: i, file_ios
  +      integer(i4b)           :: Uelement, Uelement2,Utnemele,Utnemele2
  +      character, intent(in)  :: code
   
         do i=1,size(Evops,5)
         do Uelement=1,Nl1
  
positive.code:         !write(*,*) 'OUTPUTTING EVOPS'
          call flush()
          call open_files('w')
          call write_evops()
          call close_files('w')
  
          !write(*,*) 'OUTPUTTING EVOPS'
  
negative.code:         write(*,*) 'OUTPUTTING EVOPS'
          call flush()
          call open_files('w')
          call write_evops('w')
          call close_files('w')
  
      end subroutine only_convert_to_exciton
  