<div class="iw_mail" dir="auto" style="font-size: 13px;"><div>Dear Prof. Dal Corso,</div><div><br></div><div>In order to do a phonon-disp calculation within DFT+U+J0, I am convinced that "thermo_pw-1.7.1+QE-7.1" is not ready to do that and as you mentioned in previous post, it is not supported yet.</div><div>So, I used "thermo_pw-1.7.0+QE-7.0" which uses old-fashined inputs: lda_plus_u=.true., . . . </div><div>It successfully computes the scf before using "ph_control". Unfortunately, the code does not go to "$TEMPDIR/g1" folder to initialize the phonon calculation and issues the "davcio" error. However, when I use  "$TEMPDIR/g1" as outdir in ph_control file, the code "ph.x" successfully starts calculation. Can I somehow (by slight modification) force thermo_pw.x to search in "$TEMPDIR/g1" ?</div><div>I would like to mention that the elegance of thermo_pw makes me attach to it as much as possible.<br></div><div>Thank you in advance.</div><div><br></div><div>Best regards,</div><div>Mahmoud Payami</div><div><br></div></div>