<div dir="ltr"><div>It is a (very old) programming choice: "plotrho" was devised to be run interactively, so if it does not prompt for a variable, there is no need to supply it</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 14, 2020 at 4:41 PM <<a href="mailto:jqhuang16b@imr.ac.cn">jqhuang16b@imr.ac.cn</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><p>
        <span style="font-size:16px">Dear QE developers,</span> 
</p>
<p>
        <br>
</p>
<p>
        <span style="font-size:16px">There seems to be a problem in the </span><span style="font-size:16px">plotrho.f90 code(qe6.5, also in prev</span><span style="font-size:16px">ious version</span><span style="font-size:16px">), between line 93 and 99, as following:</span> 
</p>
<p>
        <br>
</p>
<p>
        <span style="font-size:16px">  IF (rhomin > 0.d0) THEN</span><br>
<span style="font-size:16px">     WRITE( stdout,'("Logarithmic scale (y/n)? > ")', advance="NO")</span><br>
<span style="font-size:16px">     READ (5, '(a)') ans</span><br>
<span style="font-size:16px">     logarithmic_scale = ans/='n'.and.ans/='N'</span><br>
<span style="font-size:16px">  ELSE</span><br>
<span style="font-size:16px">     logarithmic_scale = .false.</span><br>
<span style="font-size:16px">  ENDIF</span> 
</p>
<p>
        <br>
</p>
<p>
        <span style="font-size:16px">for the case rhomin < 0.d0 (may occur due to the error in</span><span style="font-size:16px"> the numerical calculation</span><span style="font-size:16px">), it won't read from the unit=5 file, </span> 
</p>
<p>
        <span style="font-size:16px">which leads to the program</span><span style="font-size:16px"> termination in the subsequent read</span><span style="font-size:16px"> because of the wrong data type.</span> 
</p>
<p>
        <span style="font-size:16px">So a tiny modification as follows can solve the problem</span> 
</p>
<p>
        <br>
</p>
<p>
        <span style="font-size:16px">  IF (rhomin > 0.d0) THEN</span><br>
<span style="font-size:16px">     WRITE( stdout,'("Logarithmic scale (y/n)? > ")', advance="NO")</span><br>
<span style="font-size:16px">     READ (5, '(a)') ans</span><br>
<span style="font-size:16px">     logarithmic_scale = ans/='n'.and.ans/='N'</span><br>
<span style="font-size:16px">  ELSE</span> 
</p>
<p>
        <span style="font-size:16px">      </span><span style="font-size:16px">READ (5, '(a)') ans</span><span style="font-size:16px"></span><br>
<span style="font-size:16px">     logarithmic_scale = .false.</span><br>
<span style="font-size:16px">  ENDIF</span>
</p>
<p>
        <span style="font-size:16px"></span>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <span style="font-size:16px">Best regards,</span>
</p>
<p>
        <span style="font-size:16px">Jian-qi Huang</span><span style="font-size:16px"> </span>
</p>
<p>
        <span style="font-size:16px"></span> 
</p>
<br>
<br>
<br>
<span>
<hr style="margin:0.5em 0px;width:10em;height:1px;background-color:rgb(153,153,153);border:medium none" align="left">
<p>
        Jian-qi Huang
</p>
<p>
        <span>Magnetism and Magnetic Materials Division</span><br>
<span>Institute of Metal Research </span><br>
<span>Chinese Academy of Sciences</span><br>
<span>72 Wenhua Road, Shenyang 110016, China</span> 
</p>
<p>
        <span>email:<a href="mailto:jqhuang16b@imr.ac.cn" target="_blank">jqhuang16b@imr.ac.cn</a></span> 
</p>
<p>
        <span><br>
</span> 
</p>
</span>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>