Previous Thread
Next Thread
Print Thread
Page 2 of 2 1 2
Re: rmsd-rgyr.inp
r_hz #370 11/22/06 01:06 AM
Joined: Oct 2005
Posts: 206
P
Forum Member
Offline
Forum Member
P
Joined: Oct 2005
Posts: 206
from the example
open unit 51 read unform name traj1.cor
open unit 52 read unform name traj2.cor
correl maxtime 10000
enter v1 rms
enter v2 gyra

traj firstu 51 nunit 2 begin 10500 skip 500 end 55000

why there is two *.cor files , both these files have same values repeated twice or values of nsavc starting from 10500 to 55000 partly filled in each other.

Re: rmsd-rgyr.inp
prakash #371 11/22/06 01:49 AM
Joined: Sep 2003
Posts: 8,483
rmv Offline
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 8,483
The latter is the usual case; perhaps steps 500 thru 30000 are in the first file, and steps 30500 thru 60000 are in the second, so both are needed to analyze 10500 thru 55000.


Rick Venable
computational chemist

Re: rmsd-rgyr.inp
rmv #372 11/22/06 02:14 AM
Joined: Oct 2005
Posts: 206
P
Forum Member
Offline
Forum Member
P
Joined: Oct 2005
Posts: 206
thank you

Re: rmsd-rgyr.inp
lennart #32003 05/29/13 06:27 PM
Joined: Oct 2012
Posts: 24
S
Forum Member
Offline
Forum Member
S
Joined: Oct 2012
Posts: 24
I ran the script using psf, pdb, and traj made from your merge script and deleted water.
The error I get running the rmsd-rgyr script is:
At line 959 of file dynio.f
Fortran runtime error: End of file

The output file is attached.
This is the input file:
Code:
!*FILENAME: rmsd-rgyr.inp
!*PURPOSE: compute rmsd vs initial structure and radius of gyration from
!trajectory
!*AUTHOR: Lennart Nilsson, Karolinska Institutet, October 2003
!*
!unix environment variable CHM_HOME points to CHARMM installation
!directory
stream toppar.str

read psf card name "new_psf.psf"
open unit 2 card read name -
"test.pdb"
read coor pdb resid unit 2

! put reference coordinates into comparison set as well
coor copy comp
open unit 11 write form name rmsd.dat
write title unit 11
* time rmsd
*
open unit 12 write form name rgyr.dat
write title unit 12
        * time rgyr
*
! assume this is a trajectory where overall translation/rotation has been
!removed and
! without water or other unnecessary things

open unit 51 read unform name -
"traj_all.cor"

correl maxtime 10000000
enter v1 rms
enter v2 gyra

traj firstu 51

write v1 dumb time unit 11
*hi
*

write v2 dumb time unit 12
*hi
*

end

Attached Files
test.out.txt (24.59 KB, 529 downloads)
Re: rmsd-rgyr.inp
lennart #35170 08/04/15 11:38 PM
Joined: Jul 2015
Posts: 9
Forum Member
Offline
Forum Member
Joined: Jul 2015
Posts: 9
im doing this for a .pdb but show me this

The following time series will be filled:
V2


VOPEN> Attempting to open::1_wt.pdb::

***** LEVEL -2 WARNING FROM *****
***** Apparent problem with endianness
******************************************
BOMLEV ( 0) IS REACHED - TERMINATING. WRNLEV IS 2

Re: rmsd-rgyr.inp
lennart #35171 08/05/15 06:41 AM
Joined: Sep 2003
Posts: 4,783
Likes: 2
lennart Offline OP
Forum Member
OP Offline
Forum Member
Joined: Sep 2003
Posts: 4,783
Likes: 2
This script analyzes the time-dependent evolution of RMSD and RGYR from a trajectory. To analyze a single coordinate set contained in a PDB file:

read coor pdb name 1_wt.pdb
read coor pdb name ref.pdb comp
coor rgyr
coor orient


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Page 2 of 2 1 2

Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.4
(Release build 20200307)
Responsive Width:

PHP: 5.6.33-0+deb8u1 Page Time: 0.011s Queries: 27 (0.005s) Memory: 0.9473 MB (Peak: 1.0522 MB) Data Comp: Off Server Time: 2020-08-12 01:15:38 UTC
Valid HTML 5 and Valid CSS