[Pw_forum] pw.x not able to specity initial ion velocities
Tim Teatro
timtro at gmail.com
Wed Aug 23 17:59:50 CEST 2006
Hello all.
This is probably an easy one. I'm sorry to bother you with it.
When I try to specify starting velocities for ions in a pw.x input file
(obviously for dynamics), I set ion_velocities = 'from_input' and specify an
ATOMIC_VELOCITIES card. When it runs, I get a warning that a card is
ignored. The ignored card is (for whatever reason) the last entry in my
ATOMIC_VELOCITIES card ( Species(n) Vx(n) Vy(n) Vz(n) ), and the
temperature calculated at the end of the first step certainly does not
reflect the temperature that the system should have given the velocities I
specified.
I looked in the source code, and according to read_cards.f90, the program
should hit the velocity card and read of the next nat lines, so I have no
idea why it would try to read off my last velocity entry as a card. There
also doesn't seem to be any comments that this bit of code is restricted to
cp.x or something. But is it?
Thanks!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: /pipermail/attachments/20060823/7c0650cd/attachment.htm
More information about the Pw_forum
mailing list