anchor.diff:  
           !my siesta title thing
           !do i=info%n_wfs-nwfts+1, info%n_wfs
  +        print *, info%wfs_occ(:,1)
           l=1
           do i=info%n_wfs-info%n_pwfs+1, info%n_wfs
             if (info%wfs_occ(i, 1) .le. 0) then
  
positive.code:        call ps_io_save(info)
        endsubroutine siesta_output
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
        subroutine ps_io_type_fill(grid_in,info)
         type(radial_grid_type), intent(in) :: grid_in
         type(ps_io_type) ,intent(out) :: info
         !mesh parameters
         type(mesh_type) :: m
         integer :: no_mesh_points,i,j,l,ir,ns,lam,n
         real(R8) :: rn,r1
         !generalities
         character, external :: atom_name*2
  
negative.code:          end if
         enddo
  
         do l=1,nwfsx!define the proper cutoff values of the wfs
           do n=1,grid_in%mesh
             if ( info%m%r(n) .lt. rcut(l) ) then
               info%wfs_rc(l) = info%m%r(n)
             endif
           enddo
         enddo
  
        endsubroutine ps_io_type_fill
  
        subroutine ps_io_save(info)!this seems only necessary to setup the filename and assign io units
  