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 |