From: Laura Joana (laurajoanalopes_at_gmail.com)
Date: Thu Jan 11 2018 - 10:36:31 CST
Dear all again,
Sorry for the two e-mails, but I decided to explore the issue a bit more.
Looking at the NAMD 2.12 release features I guess the problem is in this
sentence:
"The patch grid and other simulation paramters are *not* changed, so the
dimensions of the new and old structures must be compatible." (under Ability
to reload molecular structure without restarting).
Making a diff with my two logfiles I get this:
< Info: PATCH GRID IS 3 (PERIODIC) BY 3 (PERIODIC) BY 1 (PERIODIC)
< Info: PATCH GRID IS 2-AWAY BY 2-AWAY BY 1-AWAY
< Info: Reading from binary file
/home/silva-ll/cyclodextrine/ic/cycloA_A.vel
< Info: REMOVING COM VELOCITY -0.0432567 0.00241296 0.0339374
< Info: LARGEST PATCH (0) HAS 417 ATOMS
--- > Info: PATCH GRID IS 4 (PERIODIC) BY 2 (PERIODIC) BY 2 (PERIODIC) > Info: PATCH GRID IS 2-AWAY BY 1-AWAY BY 1-AWAY > Info: REMOVING COM VELOCITY 0.087986 -0.125532 0.0988174 > Info: LARGEST PATCH (3) HAS 237 ATOMS Does anyone know if I'm in the good road to find out what is the problem ? Best regards, Laura On Thu, Jan 11, 2018 at 3:26 PM, Laura Joana <laurajoanalopes_at_gmail.com> wrote: > Dear all, > > I'm having a problem involving the use of the command *reinitatoms file*. > I run two simulations, almost the same configuration file, the only > difference is that for the first one I call *bincoordinates, > binvelocities and extendedsystem* with some checkpoint and the other I > use the *temperature* command. For both I call *reinitatoms file* (some > other checkpoint) before any run command. I expected no difference between > those simulations. The problem is that the first one runs normally and the > second give me an error and do not even run ! The error is "Periodic cell > has become too small for original patch grid". > So the question is: what is the difference in using *reinitatoms file* > after *bincoordinates, binvelocities and extendedsystem* and after the > *temperature* command ? Shouldn't NAMD just read *file.coor, file.vel and > file.xsc* and ignore the other previous commands ? > Best regards, > > Laura >
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:47 CST