anchor.diff:      ! declarations
       real(dp), dimension(:,:), allocatable  :: Ham
       complex(dpc), allocatable:: rho1(:,:), prhox1(:,:), prhodx1(:,:)
  -    complex(dpc), dimension(:,:,:,:,:), allocatable      :: Ueg, Uee
  +    complex(dpc), dimension(:,:,:,:,:), allocatable      :: Evops
       complex(dpc), dimension(:,:), private, allocatable   :: A, U, VT
       complex(dpc), dimension(:), private, allocatable     :: B
       real(dp), dimension(:), private, allocatable         :: EIGVAL
  
positive.code:         write(*,*) 'READING EXTERNAL EVOPS'
          call flush()
          call open_files('r')
          call read_evops(Uee)
          call close_files('r')
  
          write(*,*) 'CALCULATION OF THE DESIGN MATRIX'
  
negative.code:       end do
  
        deallocate(calc)
  
        write(*,*) 'ACHIEVED CHI-SQUARE:', chisq
  
      end subroutine write_fitted_evops
  
  