Minor fixes at I/O. HDF is now working with FullyFullyNoco but not functioning...

Minor fixes at I/O. HDF is now working with FullyFullyNoco but not functioning yet since the offdiagonal elements of the density are not stroed in hdf yet.

io/cdnpot_io_hdf.F90

@@ -1690,7 +1690,7 @@ MODULE m_cdnpot_io_hdf
 
             dimsInt(:4)=(/jmtd,nlhd+1,ntype,input%jspins/)
             CALL h5dopen_f(groupID, 'fr', frSetID, hdfError)
-            CALL io_write_real4(frSetID,(/1,1,1,1/),dimsInt(:4),den%mt)
+            CALL io_write_real4(frSetID,(/1,1,1,1/),dimsInt(:4),den%mt(:,:,:,:input%jspins))
             CALL h5dclose_f(frSetID, hdfError)
 
             dimsInt(:3)=(/2,ng3,input%jspins/)
@@ -1757,7 +1757,7 @@ MODULE m_cdnpot_io_hdf
             CALL h5screate_simple_f(4,dims(:4),frSpaceID,hdfError)
             CALL h5dcreate_f(groupID, "fr", H5T_NATIVE_DOUBLE, frSpaceID, frSetID, hdfError)
             CALL h5sclose_f(frSpaceID,hdfError)
-            CALL io_write_real4(frSetID,(/1,1,1,1/),dimsInt(:4),den%mt)
+            CALL io_write_real4(frSetID,(/1,1,1,1/),dimsInt(:4),den%mt(:,:,:,:input%jspins))
             CALL h5dclose_f(frSetID, hdfError)
 
             dims(:3)=(/2,ng3,input%jspins/)
@@ -1861,7 +1861,7 @@ MODULE m_cdnpot_io_hdf
          CALL h5screate_simple_f(4,dims(:4),frSpaceID,hdfError)
          CALL h5dcreate_f(groupID, "fr", H5T_NATIVE_DOUBLE, frSpaceID, frSetID, hdfError)
          CALL h5sclose_f(frSpaceID,hdfError)
-         CALL io_write_real4(frSetID,(/1,1,1,1/),dimsInt(:4),den%mt)
+         CALL io_write_real4(frSetID,(/1,1,1,1/),dimsInt(:4),den%mt(:,:,:,:input%jspins))
          CALL h5dclose_f(frSetID, hdfError)
 
          dims(:3)=(/2,ng3,input%jspins/)

optional/flipcdn.f90

@@ -116,8 +116,7 @@ IF (input%l_removeMagnetisationFromInterstitial) THEN
       den%pw(:,2)=(den%pw(:,1)+den%pw(:,2))*0.5 !mean value 
       den%pw(:,1)=den%pw(:,2)
       IF (noco%l_noco) THEN
-         den%pw(:,3)=0.0
-         den%pw(:,4)=0.0
+         den%pw(:,3)=CMPLX(0.0,0.0)
       END IF
    END IF