[CP2K:4] Re: printing forces to __STD_OUT__

Nichols A. Romero naro... at gmail.com
Wed Apr 25 20:57:56 UTC 2007


Here you go. I've noticed some really odd behavior at times with the PRINT
section. For example, the output format of forces (dcd, xmol, atomic) being
affected by that Trajectory output format, or vice-versa (do not precisely
recollect).

&MOTION
  &MD
    ENSEMBLE NVE
    STEPS 2
    TIMESTEP 1.0
    PV_AVA T
  &END MD
  &PRINT
    &FORCES ON
      ADD_LAST NUMERIC
      FILENAME __STD_OUT__
    &END FORCES
    &TRAJECTORY ON
      ADD_LAST NUMERIC
      FILENAME ./Trajectory.dcd
      FORMAT DCD
    &END TRAJECTORY
    &VELOCITIES ON
      ADD_LAST NUMERIC
      FILENAME ./Velocities.dcd
      FORMAT DCD
    &END VELOCITIES
    &STRESS ON
      ADD_LAST NUMERIC
      FILENAME __STD_OUT__
    &END STRESS
  &END PRINT
&END MOTION

On 4/25/07, Teodoro Laino <teodor... at gmail.com> wrote:
>
> This works:&MOTION
>   &MD
>     ENSEMBLE NVE
>     STEPS 100
>     TIMESTEP 2.5
>     TEMPERATURE 298
>     &NOSE
>       LENGTH 3
>       YOSHIDA 3
>       TIMECON 1000
>       MTS 2
>     &END NOSE
>   &END MD
>   &PRINT
>     &FORCES ON
>       ADD_LAST NUMERIC
>       FILENAME __STD_OUT__
>     &END FORCES
>   &END
> &END MOTION
>
> Can you please post the full MOTION section?
>
> On 25 Apr 2007, at 22:46, Nichols A. Romero wrote:
>
> This is what I used. Have not tried it without the ADD_LAST NUMERIC.
>
>    &FORCES ON
>       ADD_LAST NUMERIC
>       FILENAME __STD_OUT__
>     &END FORCES
>
>
> On 4/25/07, Teodoro Laino <teodor... at gmail.com> wrote:
> >
> >
> > Ciao Nichols,
> >
> > If I use the following section it works:
> >
> > &MOTION
> >    &MD
> >      ENSEMBLE NVE
> >      STEPS 100
> >      TIMESTEP 2.5
> >      TEMPERATURE 298
> >      &NOSE
> >        LENGTH 3
> >        YOSHIDA 3
> >        TIMECON 1000
> >        MTS 2
> >      &END NOSE
> >    &END MD
> >    &PRINT
> >      &FORCES
> >        FILENAME __STD_OUT__
> >      &END
> >    &END
> > &END MOTION
> >
> > Can we know PRECISELY how your MOTION section looks like?
> > Beware if you use FORMAT DCD (that is implicitly binary) in the
> > FORCES section of MOTION you will get for sure the error you were
> > mentioning..
> >
> > Cheers
> > Teo
> > On 25 Apr 2007, at 22:29, Nichols A. Romero wrote:
> >
> > > This is a follow-up to an email that I posted yesterday and that
> > > Fawzi replied to.
> > >
> > > There are two places where one can print the forces to standard out.
> > > &FORCE_EVAL&PRINT&FORCES
> > > and
> > > &MOTION&PRINT&FORCES
> > >
> > > I've only manage to get the former one to print forces to std_out.
> > > The latter
> > > one always gives some I/O error. Just FYI.
> > >
> > > --
> > > Nichols A. Romero, Ph.D.
> > > 1613 Denise Dr. Apt. D
> > > Forest Hill, MD 21050
> > > 443-567-8328 (C)
> > > 410-306-0709 (O)
> > >
> > > >
> >
> >
> >
> >
> >
> >
> >
>
> >
>


-- 
Nichols A. Romero, Ph.D.
1613 Denise Dr. Apt. D
Forest Hill, MD 21050
443-567-8328 (C)
410-306-0709 (O)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20070425/d5eee27b/attachment.htm>


More information about the CP2K-user mailing list