anchor.diff:      
           !
           title = ""
  -    
  -        do i=1, info%n_wfs
  -        print *, info%wfs_n(i),info%wfs_l(i)
  +
  +    
  +!        do l = 0, 3
  +!          n = minval(info%wfs_n, mask=info%wfs_l == l)
  +!          !j = minval(info%wfs_j, mask=(info%wfs_l == l .and. info%wfs_n == n))
  +!    
  +!          do i = 1, info%n_wfs
  +!            !if (info%wfs_n(i) == n .and. info%wfs_l(i) == l .and. info%wfs_j(i) == j) then
  +!            if (info%wfs_n(i) == n .and. info%wfs_l(i) == l) then
  +!              !select case (irel)
  +!              !case('nrl')
  +!                write(title,'(A)') (title)
  +!                write(title,'(A)') trim(title)
  +!                !here
  +!                write(title,'(A,A2,F5.2,"  r=",F5.2,"/")') trim(title),&
  +!                                info%wfs_label(i), info%wfs_occ(i, 1),&
  +!                                 info%wfs_rc(i)
  +!     !!          case('isp')
  +!     !!            write(title,'(A,A2,F4.2,1X,F4.2,1X,F4.2,"/")') &
  +!      trim(title), &
  +!     !!                     info%wfs_label(i), info%wfs_occ(i, 1),&
  +!       info%wfs_occ(i, 2), &
  +!     !!                     info%wfs_rc(i)
  +!     !!          case ('rel')
  +!     !!            occ = sum(info%wfs_occ(:, 1), mask=(info%wfs_l == l&
  +!      .and. &
  +!     !!                 (info%wfs_j == l + M_HALF .or. info%wfs_j == l&
  +!       - M_HALF)))
  +!     !!            write(title,'(A,A2,F5.2,"  r=",F5.2,"/")') &
  +!      trim(title), &
  +!     !!                                           info%wfs_label(i), &
  +!      occ, info%wfs_rc(i)
  +!     !!          end select
  +!            end if
  +!          end do
  +!        end do
  +
  +        !my siesta title thing
  +        !do i=info%n_wfs-nwfts+1, info%n_wfs
  +        l=1
  +        do i=info%n_wfs-4+1, info%n_wfs
  +          if (info%wfs_occ(i, 1) .le. 0) then
  +            info%wfs_occ(i, 1)=0.0
  +          endif
  +          write(title,'(A,A2,F5.2,"  r=",F5.2,"/")') trim(title),&
  +                                info%wfs_label(i),&
  +                                info%wfs_occ(i, 1),&
  +                                info%wfs_rc(l)
  +          l=l+1
           enddo
  -    
  -        do l = 0, 3
  -          n = minval(info%wfs_n, mask=info%wfs_l == l)
  -          !j = minval(info%wfs_j, mask=(info%wfs_l == l .and. info%wfs_n == n))
  -    
  -          do i = 1, info%n_wfs
  -            !if (info%wfs_n(i) == n .and. info%wfs_l(i) == l .and. info%wfs_j(i) == j) then
  -            if (info%wfs_n(i) == n .and. info%wfs_l(i) == l) then
  -              !select case (irel)
  -              !case('nrl')
  -                print *, "Here is what we want"
  -                print *, info%wfs_rc
  -                write(title,'(A)') (title)
  -                write(title,'(A)') trim(title)
  -                write(title,'(A,A2,F5.2,"  r=",F5.2,"/")') trim(title), &
  -                                info%wfs_label(i), info%wfs_occ(i, 1), info%wfs_rc(i)
  -     !!          case('isp')
  -     !!            write(title,'(A,A2,F4.2,1X,F4.2,1X,F4.2,"/")') trim(title), &
  -     !!                     info%wfs_label(i), info%wfs_occ(i, 1), info%wfs_occ(i, 2), &
  -     !!                     info%wfs_rc(i)
  -     !!          case ('rel')
  -     !!            occ = sum(info%wfs_occ(:, 1), mask=(info%wfs_l == l .and. &
  -     !!                 (info%wfs_j == l + M_HALF .or. info%wfs_j == l - M_HALF)))
  -     !!            write(title,'(A,A2,F5.2,"  r=",F5.2,"/")') trim(title), &
  -     !!                                           info%wfs_label(i), occ, info%wfs_rc(i)
  -     !!          end select
  -            end if
  -          end do
  -    
  -        end do
  +        print *, title
  +        
       
           !Mesh
           call mesh_null(new_m)
  
positive.code:         end select
            
          call date_and_time(values=values)
          write(header,'(A,A,3X,I4.4,"/",I2.2,"/",I2.2)') "APE-QE Version-", PACKAGE_VERSION, values(1), values(2), values(3)
          select case (info%scheme)
          case (HAM)
            header = trim(header)//"   Hamann"
  
negative.code:        info%wfs_n(:)=nn(:)
         allocate(info%wfs_l(nwfx))
         info%wfs_l(:)=ll(:)
         !allocate(info%wfs_j(nwfx))
         !info%wfs_j(:)=jj(:)
         allocate(info%wfs_label(nwfx))
         info%wfs_label(:)=el(:)
         allocate(info%wfs_rc(nwfsx))!cut off for pseudos (NO US allowed)
         info%wfs_rc(:)=rcut(:)
         allocate(info%wfs_occ(nwfx,info%nspin))
         !here we are using the variale configuration occupation listing, just the first one though
         !info%wfs_occ(:,1)=octsc(:,1)!pretend there's only a single spin channel
         info%wfs_occ(:,1)=oc(:)!pretend there's only a single spin channel
         allocate(info%wfs_ev(nwfsx))
         !info%wfs_ev(:)=enltsc(:,1)
         info%wfs_ev(:)=enls(:)
  