anchor.diff:         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_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(nwfsx,info%nspin))
  +       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)=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(:)
  
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:       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,ocs,oc
        use ld1_parameters, only: nwfx,nwfsx
  
        type ps_io_type
  