<div dir="ltr"><div dir="ltr">On Tue, Jan 19, 2021 at 10:15 AM Shivesh Sivakumar <<a href="mailto:shiveshsivakumar@gmail.com">shiveshsivakumar@gmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br>I tried something like 'gfortran v_of_rho.f90 <input> > <out>' but it didn't work out </div></blockquote><div><br></div><div>it's not that simple: that subroutine requires a bunch of stuff to be initialized and must be called with the needed input and output variables. You may have better luck by modifying one of the many post processing codes in PP/src that read final SCF data: they call 'read_file_new' that calls 'v_of_rho' to rebuild the potential from the charge. IIRC one of those codes (ppacf.f90) actually performs some splitting of exchange and correlation terms, in addition to many other operations.<br></div><div><br></div><div>Paolo<br></div><div><br></div></div>-- <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 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>