Yan Gao wrote:
Thank you Vandana and Justin, the cat command solve the problem.
Regarding the removement of the unnecessary lines in the middle, is
there a command that can do it automatically? or one has to make code by
self? Thanks.
That's a trivial exercise for a simple text editor.
-Justin
Yan
On Tue, May 25, 2010 at 3:25 PM, Justin A. Lemkul <[email protected]
<mailto:[email protected]>> wrote:
VANDANA KUMARI wrote:
Hi Yan,
You can do it by converting your both gro file into pdb files
and concatenating both pdb files into one file . and you can
convert pdb into gro using editconf
steps are as follows:
1. editconf -f 1.gro -o 1.pdb
2. editconf -f 2.gro -o 2.pdb
3. cat 1.pdb 2.pdb > both.pdb
you need to delete unnecessary line i.e. CRYST1 or END from you
both.pdb file
4. editconf -f both.pdb -o both.gro
Conversion back and forth between .pdb and .gro is not required (but
it does work). The main point of creating a concatenated file is
that both molecules are placed correctly within the same box, so for
each:
editconf -f 1.gro -o 1_newbox.gro -center x y z -box x y z
(and the same for molecule 2)
Getting the molecule positions and box vectors right is the real trick.
Then concatenate, remove unnecessary box and header lines in the
middle of the file, correct the atom count and you're done.
-Justin
--
Vandana Kumari
------------------------------------------------------------------------
*From:* [email protected]
<mailto:[email protected]>
[[email protected]
<mailto:[email protected]>] on behalf of Yan Gao
[[email protected] <mailto:[email protected]>]
*Sent:* Tuesday, May 25, 2010 5:33 PM
*To:* Discussion list for GROMACS users
*Subject:* [gmx-users] How to combine two gro files in one big file
Hi there,
I want to simulate two molecules which are identical, and one is
rotated/translated from another. I have the gro file for one
molecule.
I tried editconf, which helped manipulate the molecule. Then I
want to combine it with the original gro file, so that the two
identical molecules are placed in one file. I searched the
manual but could not find a proper command. Could anyone help?
Thanks.
--
Yan
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu <http://vt.edu> | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
gmx-users mailing list [email protected]
<mailto:[email protected]>
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before
posting!
Please don't post (un)subscribe requests to the list. Use the www
interface or send it to [email protected]
<mailto:[email protected]>.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
Yan Gao
Jacobs School of Engineering
University of California, San Diego
Tel: 858-952-2308
Email: [email protected] <mailto:[email protected]>
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [email protected].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php