parm ../../build/full.topo.hmr [traj]

strip :WAT,Na+,K+,Cl-

autoimage origin

rmsd :1,2&!@H= first

nativecontacts :1@N1,C6,C5,C4,N4,N3,C2,O2 :2@N1,N2,C2,N3,C4,C5,C6,O6,N7,C8,N9 mindist
distance COM :1@N1,C6,H6,C5,H5,C4,N4,H41,H42,N3,C2,O2 :2@N1,C2,H2,N3,C4,C5,C6,N2,H21,H22,N7,C8,H8,N9,O6

vector v_base1 :1@N1,C6,C5,C4,N4,N3,C2,O2 corrplane
vector v_base2 :2@N1,C2,N3,C4,C5,C6,N2,N7,C8,N9,O6 corrplane
vectormath name normalangle vec1 v_base1 vec2 v_base2 dotangle

trajin ../../traj.1.01 remdtraj remdtrajtemp 396.40 trajnames ../../traj.1.02,../../traj.1.03,../../traj.1.04,../../traj.1.05,../../traj.1.06,../../traj.1.07,../../traj.1.08,../../traj.1.09,../../traj.1.10,../../traj.1.11,../../traj.1.12,../../traj.1.13,../../traj.1.14,../../traj.1.15,../../traj.1.16,../../traj.1.17,../../traj.1.18,


run

writedata criteria-raw_396.40.dat Contacts_00002[mindist] COM normalangle
