[gmx-users] how to use g_angle to get the average dihedral angle ?
maggin
maggin.chen at gmail.com
Wed Apr 16 09:03:40 CEST 2014
Hi,
I use g_angle and want to get the average dihedral angle, so operation as
follow:
(1) mk_angndx -s md_0_1.tpr -n angle.ndx -type dihedral
and get angle.ndx :
[ Phi=180.0_2_43.93 ]
4 17 15 16 9 13 12 14 19 33 31 32
35 49 47 48 50 63 61 62 65 80 78 79
82 102 100 101 104 124 122 123 126 141 139 140
143 158 156 157 160 180 178 179 182 192 190 191
184 188 187 189 194 212 210 211 214 231 229 230
233 246 244 245 238 242 241 243 248 268 266 267
270 287 285 286 289 294 292 293 296 301 299 300
303 319 318 320
[ Phi=180.0_2_4.18 ]
15 19 17 18 31 35 33 34 47 50 49 58
61 65 63 64 74 72 76 77 78 82 80 81
100 104 102 103 122 126 124 125 139 143 141 142
152 150 154 155 156 160 158 159 178 182 180 181
190 194 192 193 210 214 212 213 229 233 231 232
244 248 246 247 266 270 268 269 285 289 287 288
292 296 294 295 299 303 301 302
[ Phi=180.0_2_4.60 ]
67 70 72 71 70 76 72 73 71 76 74 75
145 148 150 149 148 154 150 151 149 154 152 153
196 199 202 200 199 204 200 201 199 206 202 203
200 208 204 205 202 208 206 207 204 206 208 209
(2) g_angle -n angle.ndx -type dihedral -ov -all
and get:
Group 0 (Phi=180.0_2_43.93) has 88 elements
Group 1 (Phi=180.0_2_4.18) has 80 elements
Group 2 (Phi=180.0_2_4.60) has 48 elements
Select a group:
At this step, I want to get the average dihedral angle of a peptide, but it
appear Group 0, Group 1, Group 2. So, how to get average dihedral angle as a
whole , not seperate in groups?
Thank you very much!
maggin
--
View this message in context: http://gromacs.5086.x6.nabble.com/how-to-use-g-angle-to-get-the-average-dihedral-angle-tp5015873.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
More information about the gromacs.org_gmx-users
mailing list