Fix library-path projection overflow when num_proj > num_wann

The library interface called get_projections(lcount=.false.) directly,
allocating input_proj at size num_wann.  When the projections block
defines more functions than num_wann (valid for disentanglement), the
fill loop writes beyond the array bounds.

Add the counting pass (lcount=.true.) before the fill pass so that
num_proj is updated to the actual projection count before input_proj
is allocated.

Co-Authored-By: Claude Sonnet 4.6 <noreply@anthropic.com>

Author: 2 people

src/wannier90_readwrite.F90

@@ -1606,6 +1606,12 @@ subroutine w90_wannier90_readwrite_read_projections(settings, proj, proj_input,
 
       if (atom_data%num_species > 0) then
         ! only read projections if atom positions are known
+        ! Count pass first: num_proj may exceed num_wann when there are more initial
+        ! projections than target Wannier functions (valid for disentanglement).
+        ! Without this, input_proj is allocated at size num_wann and overflows.
+        call w90_readwrite_get_projections(settings, num_proj, atom_data, num_wann, proj_input, &
+                                           recip_lattice, .true., spinors, bohr, stdout, error, comm)
+        if (allocated(error)) return
         call w90_readwrite_get_projections(settings, num_proj, atom_data, num_wann, proj_input, &
                                            recip_lattice, .false., spinors, bohr, stdout, error, comm)
         if (allocated(error)) return