readdvscf : error # 80 ../phonons/save/pb.dvscf_q2 too short, check ecut

Post here questions linked with issue while running the EPW code

Moderator: stiwari

Post Reply
zega14
Posts: 15
Joined: Wed Oct 21, 2020 8:05 pm
Affiliation: University of Arkansas

readdvscf : error # 80 ../phonons/save/pb.dvscf_q2 too short, check ecut

Post by zega14 »

Dear All

When I run epw, I come up with the readdvscf error as indicated on the subject. I tried to look for readdvsc.f90 file in EPW/src/ folder but there is no file in that name. On this forum before 4 years ago, the same error had been discussed by Jiachen and Prof. Samuel, and the solution was to remove the checks. Does anyone know in which .f90 file the following code exist or how can I solve this error? I am using the latest QE (6.7). Thank you!

INTEGER(kind=8) :: mult_unit

CALL set_ndnmbr ( 0, iq, 1, nqc, filelab)
tempfile = trim(dvscf_dir) // trim(prefix) // '.dvscf_q' // filelab
unf_recl = DIRECT_IO_FACTOR * lrdrho

mult_unit = unf_recl
mult_unit = recn * mult_unit

open (iudvscf, file = tempfile, form = 'unformatted', &
access = 'direct', iostat=ios,recl = unf_recl,status='old')
IF (ios /= 0) call errore ('readdvscf','error opening ' // tempfile, iudvscf)

#if defined(__NAG)
CALL fstat( iudvscf, statb, errno=ios)
IF (mult_unit .gt. statb%st_size) call errore('readdvscf', &
trim(tempfile)//' too short, check ecut', iudvscf)
#endif
#if ! defined(__NAG)
ios = fstat ( iudvscf, statb)
!IF (recn * unf_recl .gt. statb(8)) call errore('readdvscf', &
IF (mult_unit .gt. statb(8)) call errore('readdvscf', &
trim(tempfile)//' too short, check ecut', iudvscf)
#endif
!
read (iudvscf, rec = recn) dvscf
close (iudvscf, status = 'keep')
!
!
end subroutine readdvscf



Sincerely, Zegnet

hpaudya1
Posts: 190
Joined: Tue Mar 21, 2017 7:11 pm
Affiliation:

Re: readdvscf : error # 80 ../phonons/save/pb.dvscf_q2 too short, check ecut

Post by hpaudya1 »

Hi Zegnet,

I think here;
https://gitlab.com/QEF/q-e/-/blob/devel ... .f90#L1485

Please make sure you read the discussion and comments by Samuel.

Best,
Hari Paudyal

zega14
Posts: 15
Joined: Wed Oct 21, 2020 8:05 pm
Affiliation: University of Arkansas

Re: readdvscf : error # 80 ../phonons/save/pb.dvscf_q2 too short, check ecut

Post by zega14 »

Dear Paudyal

Thank you so much! I miss to open that link.

Sincerely, Zegnet

Post Reply