1 |
mmeineke |
1049 |
#include <stdio.h> |
2 |
|
|
#include <stdlib.h> |
3 |
|
|
|
4 |
|
|
|
5 |
|
|
#include "params.h" |
6 |
|
|
#include "tcProps.h" |
7 |
|
|
|
8 |
|
|
|
9 |
|
|
int main( int argC, char *argV[] ){ |
10 |
|
|
|
11 |
|
|
// list of 'a priori' constants |
12 |
|
|
|
13 |
|
|
const int nLipAtoms = 19; |
14 |
|
|
const int nBonds = 18; |
15 |
|
|
const int nLipids = NLIPIDS; |
16 |
|
|
const int nSSD = NSSD; |
17 |
|
|
const int nAtoms = nLipAtoms * nLipids + nSSD; |
18 |
|
|
|
19 |
|
|
struct atomCoord atoms[nAtoms]; |
20 |
|
|
|
21 |
|
|
|
22 |
|
|
|
23 |
|
|
int i,j,k; |
24 |
|
|
|
25 |
|
|
// initialize the arrays |
26 |
|
|
|
27 |
|
|
for(i=0;i<nLipids;i++){ |
28 |
|
|
|
29 |
|
|
atoms[nLipAtoms*i+0].type = HEAD; |
30 |
|
|
|
31 |
|
|
atoms[nLipAtoms*i+1].type = CH2; |
32 |
|
|
|
33 |
|
|
atoms[nLipAtoms*i+2].type = CH; |
34 |
|
|
|
35 |
|
|
atoms[nLipAtoms*i+3].type = CH2; |
36 |
|
|
|
37 |
|
|
atoms[nLipAtoms*i+4].type = CH2; |
38 |
|
|
|
39 |
|
|
atoms[nLipAtoms*i+5].type = CH2; |
40 |
|
|
|
41 |
|
|
atoms[nLipAtoms*i+6].type = CH2; |
42 |
|
|
|
43 |
|
|
atoms[nLipAtoms*i+7].type = CH2; |
44 |
|
|
|
45 |
|
|
atoms[nLipAtoms*i+8].type = CH2; |
46 |
|
|
|
47 |
|
|
atoms[nLipAtoms*i+9].type = CH2; |
48 |
|
|
|
49 |
|
|
atoms[nLipAtoms*i+10].type = CH3; |
50 |
|
|
|
51 |
|
|
atoms[nLipAtoms*i+11].type = CH2; |
52 |
|
|
|
53 |
|
|
atoms[nLipAtoms*i+12].type = CH2; |
54 |
|
|
|
55 |
|
|
atoms[nLipAtoms*i+13].type = CH2; |
56 |
|
|
|
57 |
|
|
atoms[nLipAtoms*i+14].type = CH2; |
58 |
|
|
|
59 |
|
|
atoms[nLipAtoms*i+15].type = CH2; |
60 |
|
|
|
61 |
|
|
atoms[nLipAtoms*i+16].type = CH2; |
62 |
|
|
|
63 |
|
|
atoms[nLipAtoms*i+17].type = CH2; |
64 |
|
|
|
65 |
|
|
atoms[nLipAtoms*i+18].type = CH3; |
66 |
|
|
} |
67 |
|
|
|
68 |
|
|
for(i=;i<nAtoms;i++){ |
69 |
|
|
atoms[i].type = SSD; |
70 |
|
|
} |