Bugfix for slice plots and HDF5

io/cdn_io.F90

@@ -230,7 +230,14 @@ MODULE m_cdn_io
                            1,-1.0,fermiEnergy,l_qfix,den)
             END IF
          ELSE
-            WRITE(*,*) TRIM(ADJUSTL(filename))//' file or relevant density entry not found.'
+            INQUIRE(FILE=TRIM(ADJUSTL(filename)),EXIST=l_exist)
+            IF(l_exist) THEN
+               WRITE(*,*) 'densityType is ', densityType
+               WRITE(*,*) 'Relevant density entry is '//TRIM(ADJUSTL(archiveName))//'.'
+               WRITE(*,*) 'Entry not found in '//TRIM(ADJUSTL(filename))//'.'
+            ELSE
+               WRITE(*,*) TRIM(ADJUSTL(filename))//' file not found.'
+            END IF
             WRITE(*,*) 'Falling back to stream access.'
             mode = CDN_STREAM_MODE
          END IF

main/cdngen.F90

@@ -120,7 +120,7 @@ SUBROUTINE cdngen(eig_id,mpi,input,banddos,sliceplot,vacuum,&
 
    IF (sliceplot%slice) THEN
       IF (mpi%irank.EQ.0) THEN
-         CALL writeDensity(stars,vacuum,atoms,cell,sphhar,input,sym,oneD,archiveType,CDN_INPUT_DEN_const,&
+         CALL writeDensity(stars,vacuum,atoms,cell,sphhar,input,sym,oneD,CDN_ARCHIVE_TYPE_CDN_const,CDN_INPUT_DEN_const,&
                            0,-1.0,0.0,.FALSE.,outDen,'cdn_slice')
       END IF
       CALL juDFT_end("slice OK",mpi%irank)