* merge: merge trajectories
*

! { parameters passed }
set 0 JOB: merge trajectorie

set struct wt_10ns
set inp wt_lz.dyna_desolv
set dcd wt_lz.dyna_prod1

set stop     10000000
set skip     500
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
set begin    2500000

bomb -2

open unit 1 card read name ../lib/top_all27_prot_na.rtf
read RTF card unit 1
close unit 1

open unit 1 card read name ../lib/par_all27_prot_na.prm
read PARA card unit 1
close unit 1

stream ../lib/aco.rtf

open unit 1 read formatted name ../scr/@inp.psf
read psf card unit 1
close unit 1
open read formatted unit 17 name ../scr/@inp.crd
read coor card unit 17
coor copy comp

open read unit 14 unfor NAME ../scr/@dcd.cor

open write unformatted unit 99 name ../scr/@inp.for_quasi.dcd
merge skip @skip firstu 14 begin @begin stop @stop nunit 1 outp 99 -
      sele all end -
      orient sele type ca end

open unit 2 write form name ../scr/@{struct}.for_quasi.pdb
write coor pdb unit 2
close unit 2
open unit 2 write form name ../scr/@{struct}.for_quasi.psf
write psf card unit 2
close unit 2
open unit 2 write form name ../scr/@{struct}.for_quasi.crd
write coor card unit 2
close unit 2