Dear Andrew, your contribution has been included in the espresso repository. Could you please check that the results are consistent with your tests? I had to make small changes in order to fit the programming requirements for flib (no USE statements that involve calculation parameters. This is to avoid cross dependencies).<br>
<br>for downloading the repository follow the instruction at <a href="http://www.qe-forge.org/gf/project/q-e/scmsvn/?action=AccessInfo">http://www.qe-forge.org/gf/project/q-e/scmsvn/?action=AccessInfo</a><br><br>cheers<br>
<br>Layla<br><br><div class="gmail_quote">2012/8/31 Richard Andrew <span dir="ltr"><<a href="mailto:Richard.Andrew@up.ac.za" target="_blank">Richard.Andrew@up.ac.za</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="font-family:Tahoma,sans-serif;font-size:13px">
Dear QE team,<br><br>You requested we send any code additions we may have.   I work in 2D systems and I noticed most DFT codes have no way for you to relax the 2D cell shape while keeping the area constant.  The cell_dofree='shape' option in QE relaxes all degrees of freedom for the cell shape at constant volume which for 2D systems means also changing the vacuum layer.  I have modified the cell_base.f90 module to include a cell_dofree='2Dshape' option which will relax the in-plane (xy plane) cell shape while keeping the area (therefore the vacuum layer) constant.  This also required a change in the deviatoric.f90 code in the flib.  Please find attached these modified files for your consideration.  All modifications are marked between '! RICHARD' comment lines.<br>
<br>I also included in cell_base.f90 a cell_dofree='2Dxy' option where the xy  plane is truly free as you had suggested in your code comment for cell_dofree='xy'.<br><br>Regards,<br><br>Richard Charles Andrew<br>
<br><div style="clear:both">Postgraduate Student<br>Physics Department<br>University of Pretoria<br>South Africa</div><br></div>
<br>_______________________________________________<br>
Q-e-developers mailing list<br>
<a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
<a href="http://qe-forge.org/mailman/listinfo/q-e-developers" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
<br></blockquote></div><br>