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/grid/vtk_inc.F90

@@ -48,11 +48,9 @@ subroutine X(vtk_out_cf)(filename, namespace, fieldname, ierr, cf_in, cube, spac
     trim(conf%git_commit), " configuration: ",  trim(conf%config_time)
   write(iunit, '(1a)') 'BINARY'
   write(iunit, '(1a)') 'DATASET STRUCTURED_POINTS'
-  write(iunit, '(1a,3i5)') 'DIMENSIONS ', cube%rs_n_global(1), cube%rs_n_global(2), cube%rs_n_global(3)
-  write(iunit, '(1a,3f12.6)') 'ORIGIN ', - int(cube%rs_n_global(1)/2)* spacing(1),&
-    - int(cube%rs_n_global(2)/2)* spacing(2),&
-    - int(cube%rs_n_global(3)/2)* spacing(3)
-  write(iunit, '(1a,3f12.6)') 'SPACING ', spacing(1), spacing(2), spacing(3)
+  write(iunit, '(1a,3i5)') 'DIMENSIONS ', cube%rs_n_global(1:3)
+  write(iunit, '(1a,3(e16.6, x))') 'ORIGIN ', - int(cube%rs_n_global/2)* spacing
+  write(iunit, '(1a,3e16.6)') 'SPACING ', spacing(1:3)
   write(iunit, '(1a,1i9)') 'POINT_DATA ', np
 #ifdef R_TCOMPLEX
   write(iunit, '(3a)') 'SCALARS Re_'//trim(fieldname), ' double 1'
@@ -132,11 +130,9 @@ subroutine X(vtk_out_cf_vector)(filename, namespace, fieldname, ierr, cf_in, vec
     trim(conf%git_commit), " configuration: ",  trim(conf%config_time)
   write(iunit, '(1a)') 'BINARY'
   write(iunit, '(1a)') 'DATASET STRUCTURED_POINTS'
-  write(iunit, '(1a,3i5)') 'DIMENSIONS ', cube%rs_n_global(1), cube%rs_n_global(2), cube%rs_n_global(3)
-  write(iunit, '(1a,3f12.6)') 'ORIGIN ', - int(cube%rs_n_global(1)/2)* spacing(1),&
-    - int(cube%rs_n_global(2)/2)* spacing(2),&
-    - int(cube%rs_n_global(3)/2)* spacing(3)
-  write(iunit, '(1a,3f12.6)') 'SPACING ', spacing(1), spacing(2), spacing(3)
+  write(iunit, '(1a,3i5)') 'DIMENSIONS ', cube%rs_n_global(1:3)
+  write(iunit, '(1a,3(e16.6, x))') 'ORIGIN ', - int(cube%rs_n_global/2) * spacing
+  write(iunit, '(1a,3e16.6)') 'SPACING ', spacing(1:3)
   write(iunit, '(1a,1i9)') 'POINT_DATA ', np
 
 #ifdef R_TCOMPLEX