1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
#!/bin/bash
. ../MasterTest.sh
CleanFiles atomic.in fluct.*.dat dpdp.fluct.dat dpdp.adp.dat \
fluct.2.pdb occ.2.pdb scale.2.pdb fluct.1.pdb \
dpdp.adp.pdb myfluct.adp.dat heavy.adp.pdb
TESTNAME='Atomic fluctuations tests'
Requires netcdf
INPUT="atomic.in"
WriteInput() {
TOP="../tz2.parm7"
cat > $INPUT <<EOF
trajin ../tz2.nc
atomicfluct out fluct.$2.dat $1
EOF
RunCpptraj "Atomic fluctuations test [$1]."
DoTest fluct.$2.dat.save fluct.$2.dat
}
WriteInput "@C,CA,N byres bfactor" 1
WriteInput ":2-12 byatom" 2
WriteInput "bymask :3,4,5" 3
WriteInput "start 10 stop 30 offset 2 byres bfactor" 4
TOP=../DPDP.parm7
cat > $INPUT <<EOF
trajin ../DPDP.nc
rms first mass
atomicfluct MyFluct out dpdp.fluct.dat adpout dpdp.adp.dat
atomicfluct Heavy calcadp :2-21&!@/H
average crdset MyAvg
run
writedata myfluct.adp.dat MyFluct[ADP]
crdout MyAvg dpdp.adp.pdb adpdata MyFluct[ADP] bfacdata MyFluct
crdout MyAvg heavy.adp.pdb adpdata Heavy[ADP] bfacdata Heavy
EOF
RunCpptraj "Atomicfluct test with ADP output"
DoTest dpdp.fluct.dat.save dpdp.fluct.dat
DoTest dpdp.adp.dat.save dpdp.adp.dat
DoTest dpdp.adp.pdb.save dpdp.adp.pdb
DoTest myfluct.adp.dat.save myfluct.adp.dat
DoTest heavy.adp.pdb.save heavy.adp.pdb
TOP=../tz2.parm7
cat > $INPUT <<EOF
trajin ../tz2.nc
atomicfluct A0 :2-12
atomicfluct A1 @C,CA,N byres bfactor
average crdset MyAvg
run
crdout MyAvg fluct.2.pdb bfacdata A0 chainid " "
crdout MyAvg occ.2.pdb occdata A0 chainid " "
crdout MyAvg scale.2.pdb bfacdata A0 bfacscale chainid " "
crdout MyAvg fluct.1.pdb bfacdata A1 bfacbyres chainid " "
EOF
RunCpptraj "Atomicfluct test with PDB B-factor/occupancy output."
DoTest fluct.2.pdb.save fluct.2.pdb
DoTest occ.2.pdb.save occ.2.pdb
DoTest scale.2.pdb.save scale.2.pdb
DoTest fluct.1.pdb.save fluct.1.pdb
EndTest
exit 0
|
