Fix OEP full - spin polarized case

Description

The OEP full was not taking the right occupations for the spin-polarized case.

News snippet

The OEP full was not taking the right occupations for the spin-polarized case.

Checklist

• I have checked that my code follows the Octopus coding standards
• I have added tests for all the new features added in this request.

src/electrons/v_ks.F90

@@ -317,14 +317,6 @@ contains
     !%End
     call messages_obsolete_variable(namespace, 'OEPPtX', 'EnablePhotons')
     call parse_variable(namespace, 'EnablePhotons', .false., ks%has_photons)
-    if (ks%has_photons) then
-      call messages_experimental('EnablePhotons = yes')
-      call photon_mode_init(ks%pt, namespace, gr, space%dim, st%qtot)
-      write(message(1), '(a,i5,a)') 'Happy to have ', ks%pt%nmodes, ' photon modes with us.'
-      call messages_info(1, namespace=namespace)
-      call mf_init(ks%pt_mx, ks%gr, st, ions, ks%pt)
-    end if
-
 
     ! initialize XC modules
 
@@ -378,7 +370,7 @@ contains
     call messages_obsolete_variable(namespace, 'NonInteractingElectrons', 'TheoryLevel')
     call messages_obsolete_variable(namespace, 'HartreeFock', 'TheoryLevel')
 
-    ! Due to how the code is make, we need to set this to have theory level other than DFT
+    ! Due to how the code is made, we need to set this to have theory level other than DFT
     ! correct...
     ks%sic%amaldi_factor = M_ONE
 
@@ -587,6 +579,15 @@ contains
       ks%vdw_self_consistent = .false.
     end if
 
+    if (ks%has_photons) then
+      call messages_experimental('EnablePhotons = yes')
+      call photon_mode_init(ks%pt, namespace, gr, space%dim, st%qtot)
+      write(message(1), '(a,i5,a)') 'Happy to have ', ks%pt%nmodes, ' photon modes with us.'
+      call messages_info(1, namespace=namespace)
+      call mf_init(ks%pt_mx, ks%gr, st, ions, ks%pt)
+    end if
+
+
     POP_SUB(v_ks_init)
 
   contains