anchor.diff:        use ld1inc, only: zed,enne,zval,nwf,nn,ll,jj,el,rcut,octsc,&
         enltsc,phis,lloc,nbeta,rhos,pawsetup,enls,vpsloc,betas,rho,&
         lls,isws,nlcc,&
  -      rhoc,rcore,nns,nwfs,vnl,ikk
  +      rhoc,rcore,nns,nwfs,vnl,ikk,ocs,oc
         use ld1_parameters, only: nwfx,nwfsx
   
         type ps_io_type
  
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:     
          !
          title = ""
  
      
  !        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
          print *, title
          
      
          !Mesh
          call mesh_null(new_m)
  