RE: is possible to remove the waterbox before running VMD

From: Wolbach, Jeffrey P \(WOLBACH\) (WOLBACH_at_juniata.edu)
Date: Thu Mar 09 2006 - 13:45:20 CST

Lily,
 
It is possible, and actually fairly easy, to create a .dcd file
containing only the protein and not the water molecules. Below is a tcl
script to do this.
 
# Start tcl script
# Initially, load the solvated system psf and dcd files.
mol load psf Swelled_300_Neutral.psf dcd 300_Simulation_2.dcd
 
# Find the number of frames.
set nf [molinfo top get numframes]
 
# Loop over each frame, extracting the solute from the solvated system.
# The exact syntax for the atomselect command will vary on a
system-by-system basis.
for {set i 0 } {$i < $nf } {incr i } {
  set sel [atomselect top "protein" frame $i]
  $sel writepdb $i.pdb
}
 
# One enters the path to the protein psf file. This psf file should be
for the protein alone,
# not the solvated system!
mol load psf PHK_Trimer.psf
 
for {set i 0 } {$i < $nf } {incr i } {
  animate read pdb $i.pdb }
 
# One has to create a unique name for the extracted protein dcd file.
animate write dcd {Tetramer_ReCent_Trio.dcd} waitfor all top
 
for {set i 0 } {$i < $nf } {incr i } {
    puts "Removing file $i.pdb"
    rm $i.pdb}
 
exit
exit
# End tcl script
 
The script is then executed in VMD, in text mode:
    vmd -dispdev text -e scriptname.tcl
 
After the new .dcd is created you can load and view in VMD like you did
for the solvated system.
 
JPW

________________________________

From: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] On
Behalf Of lily jin
Sent: Thursday, March 09, 2006 11:34 AM
To: NAMD-L; VMD-L
Subject: namd-l: is possible to remove the waterbox before running VMD

I need the waterbox to do simulation. But when I show the result in VMD,
the waterbox is not necessary, which takes a lot of memory. Is it
possible to do simulation with waterbox and remove it from the
trajectory before displaying in VMD?

Thanks!

Lily

________________________________

Yahoo! Mail
Bring photos to life! New PhotoMail
<http://pa.yahoo.com/*http://us.rd.yahoo.com/evt=39174/*http://photomail
.mail.yahoo.com> makes sharing a breeze.

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:41:43 CST