AMBER Archive (2003)

Subject: RE: AMBER: Carnal Input file

From: Ming Lei (MILei_at_clarku.edu)
Date: Mon Sep 22 2003 - 22:34:05 CDT


In my opinion, you should output a series of pdb files along your trajectory
using carnal. Then you could easily edit one shell script to grep the
certain line including the atom information in the pdb file( in Unix/linux
OS).

I don't think CARNAL could help more in this case. Just FYI.

Best,
Ming

-----Original Message-----
From: Bill Ross
To: amber_at_scripps.edu
Sent: 9/22/2003 7:23 PM
Subject: Re: AMBER: Carnal Input file

> FILES_IN
> PARM p1 pidimer.parm;
> STREAM s1 pidimer.mdcrd;
> FILES_OUT
> COORD c1 pidimer.W120_carbon_ALFA.carnal;
> DECLARE
> GROUP carbon (ATOM 1939);
> OUTPUT
> COORD c1 s1;
> END

If you want the coordinates of that atom in the coordinate
file, you need to put the group id in the final COORD statement.

> ... what I was expecting is a file with a time
> parameter, showing the time step of the simulation, in one column and
3
> columns with x,y,z positions.

You could post-process the mdcrd file with e.g. an awk script
(or a perl script, or visual basic, etc) to put it in the form
you want. Possibly adding a group with one atom to an OUTPUT
TABLE statement would cause that atom's coordinates to appear
in the table. However there is no way for carnal to know what
time steps you are using and how often you print to mdcrd, so
unless carnal's input language were extended to give a way to
specify this, there is no way it will appear in carnal output.

Hence you must postprocess the carnal output (no matter how you
get it) to add the times in the 1st column.

Bill Ross

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu