* Timeseries template input
*

set dcd wt_lz.dyna_prod1
set inp wt_lz.dyna_desolv

!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Topology, parameters and PSF
!!!!!!!!!!!!!!!!!!!!!!!!!!!
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

!  Read in psf and coordinates

open read formatted unit 20 name ../scr/@inp.psf
read psf card unit 20
open read formatted unit 21 name ../scr/@inp.crd
read coor card unit 21
coor copy comp

OPEN UNIT 14 READ unformatted NAME ../scr/@dcd.cor

open unit 98 write form name z1n_aco.ts
correl maxseries 12 maxtimesteps 10000000  maxatom 15
enter ts1 bond z 1 n l 156 s1p
enter ts2 bond z 1 n l 156 o
enter ts3 bond z 1 n l 156 c
  trajectory firstunit 14 nunits 1 skip 500 begin 0 stop 10000000
  edit ts1 veccode 3! veccode as maxseries indicates the number of timeseries here two...
  write ts1 unit 98 card
*# JOBNAME: h-bonds
*
stop