<div dir="ltr"><div><div><div>Hello Andrea, <br><br></div></div>I agree with you. However, I guess, for a simple <br>tetragonal (ibrav=6) cell, it will not be very complicated.  <br>I tried it with ibrav=0 as well. For ibrav=0 no such error as for ibrav=6,  <br>but I got arbitrary numbers and inappropriate results, in accord with expectation. <br><br></div><div>Thanks & regards, <br></div><div>Suza.    <br></div><div><div><div><div><br><br><div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 7, 2015 at 2:24 PM, Andrea Dal Corso <span dir="ltr"><<a href="mailto:dalcorso@sissa.it" target="_blank">dalcorso@sissa.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">If something is not implemented it cannot work. This is a present<br>
limitation of the thermo_pw code. All the software needed to deal with<br>
anharmonic properties of anisotropic systems has not been written yet.<br>
It will require a substantial amount of work, so do not expect it soon.<br>
<br>
Andrea<br>
<span class="im HOEnZb"><br>
<br>
<br>
On Wed, 2015-01-07 at 13:33 +0100, Suza W wrote:<br>
</span><div class="HOEnZb"><div class="h5">> Dear Prof. Giannozzi,<br>
><br>
> grep noncubic -B1 thermo_readin.f90<br>
><br>
><br>
> gives :<br>
><br>
> IF (ibrav > 3 .AND. what=='mur_lc_t') CALL errore('thermo_pw','option not &<br>
>               & available for noncubic systems',1)<br>
><br>
> Thanks and regards,<br>
> Suza.<br>
><br>
><br>
> On Wed, Jan 7, 2015 at 12:12 PM, Paolo Giannozzi <<a href="mailto:paolo.giannozzi@uniud.it">paolo.giannozzi@uniud.it</a>><br>
> wrote:<br>
><br>
> > On Wed, 2014-12-24 at 15:07 +0100, Suza W wrote:<br>
> ><br>
> > >      Error in routine thermo_pw (1):<br>
> > >      option not  available for noncubic systems<br>
> ><br>
> > $ grep  thermo_pw */*.f90 */*/*.f90<br>
> ><br>
> > nothing found<br>
> ><br>
> > --<br>
> >  Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
> >  Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
> >  Phone +39-0432-558216, fax +39-0432-558222<br>
> ><br>
> > _______________________________________________<br>
> > Pw_forum mailing list<br>
> > <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> > <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
> ><br>
> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<br>
</div></div><span class="HOEnZb"><font color="#888888">--<br>
Andrea Dal Corso                    Tel. 0039-040-3787428<br>
SISSA, Via Bonomea 265              Fax. 0039-040-3787249<br>
I-34136 Trieste (Italy)             e-mail: <a href="mailto:dalcorso@sissa.it">dalcorso@sissa.it</a><br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div></div></div></div></div></div></div>