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?