| 611 |
|
int targetNlipids, targetNwaters; |
| 612 |
|
double targetWaterLipidRatio; |
| 613 |
|
double maxZ, minZ, zHeight; |
| 614 |
– |
double maxY, minY; |
| 615 |
– |
double maxX, minX; |
| 614 |
|
|
| 615 |
|
molStart = NULL; |
| 616 |
|
|
| 785 |
|
} |
| 786 |
|
|
| 787 |
|
atoms[0]->getPos( myPos ); |
| 790 |
– |
|
| 791 |
– |
maxX = myPos[0]; |
| 792 |
– |
minX = myPos[0]; |
| 793 |
– |
|
| 794 |
– |
maxY = myPos[1]; |
| 795 |
– |
minY = myPos[1]; |
| 788 |
|
|
| 789 |
|
maxZ = myPos[2]; |
| 790 |
|
minZ = myPos[2]; |
| 791 |
|
|
| 792 |
|
for(i=0;i<nAtoms;i++){ |
| 793 |
|
atoms[i]->getPos( myPos ); |
| 802 |
– |
minX = (minX > myPos[0]) ? myPos[0] : minX; |
| 803 |
– |
maxX = (maxX < myPos[0]) ? myPos[0] : maxX; |
| 794 |
|
|
| 805 |
– |
minY = (minY > myPos[1]) ? myPos[1] : minY; |
| 806 |
– |
maxY = (maxY < myPos[1]) ? myPos[1] : maxY; |
| 807 |
– |
|
| 795 |
|
minZ = (minZ > myPos[2]) ? myPos[2] : minZ; |
| 796 |
|
maxZ = (maxZ < myPos[2]) ? myPos[2] : maxZ; |
| 797 |
|
} |
| 798 |
|
|
| 812 |
– |
boxX = (maxX - minX)+2.0; |
| 813 |
– |
boxY = (maxY - minY)+2.0; |
| 799 |
|
boxZ = (maxZ - minZ)+2.0; |
| 800 |
|
|
| 801 |
|
double centerX, centerY, centerZ; |
| 802 |
|
|
| 803 |
< |
centerX = ((maxX - minX) / 2.0) + minX; |
| 804 |
< |
centerY = ((maxY - minY) / 2.0) + minY; |
| 803 |
> |
centerX = (boxX / 2.0); |
| 804 |
> |
centerY = (boxY / 2.0); |
| 805 |
|
centerZ = ((maxZ - minZ) / 2.0) + minZ; |
| 806 |
|
|
| 807 |
|
// set up water coordinates |