Hi everyone,

                    I would like to merge (not concatenating using catdcd) 2 DCD files (not concatenating). Each DCD file has a protein chain. The total frames are 20001. I used the following script (using topotools http://www.ks.uiuc.edu/Research/vmd/plugins/topotools/ ), but it is giving segmentation fault (works fine for a single pdb file insted of dcd):

mol new psf1.psf

mol addfile dcd1.dcd first 0 last -1 step 1 waitfor all

mol new psf2.psf

mol addfile dcd2.dcd first 0 last -1 step 1 waitfor all

set sel1 [atomselect 0 all]

set sel2 [atomselect 1 all]

package require topotools

set out [::TopoTools::selections2mol "$sel1 $sel2"]

animate write psf merged.psf $out

animate write dcd merged.dcd beg 0 end 20000 $out

------

Can someone help me?

Similar questions and discussions