| 48 |
|
const double waterVol = 4.0 / waterRho; // volume occupied by 4 waters |
| 49 |
|
const double waterCell = 4.929; // fcc unit cell length |
| 50 |
|
|
| 51 |
< |
const double water_padding = 2.5; |
| 52 |
< |
const double lipid_spaceing = 5.0; |
| 51 |
> |
const double water_padding = 6.0; |
| 52 |
> |
const double lipid_spaceing = 8.0; |
| 53 |
|
|
| 54 |
|
|
| 55 |
|
int i,j,k, l; |
| 224 |
|
|
| 225 |
|
int targetWaters = nWaters + n_deleted * nLipids; |
| 226 |
|
|
| 227 |
+ |
targetWaters = (int) ( targetWaters * 1.2 ); |
| 228 |
+ |
|
| 229 |
|
// find the best box size for the sim |
| 230 |
|
|
| 231 |
|
int testTot; |
| 470 |
|
bsInfo.boxY = box_y; |
| 471 |
|
bsInfo.boxZ = box_z; |
| 472 |
|
|
| 473 |
+ |
double boxVector[3]; |
| 474 |
|
simnfo = new SimInfo(); |
| 475 |
|
simnfo->n_atoms = nAtoms; |
| 476 |
< |
simnfo->box_x = bsInfo.boxX; |
| 477 |
< |
simnfo->box_y = bsInfo.boxY; |
| 478 |
< |
simnfo->box_z = bsInfo.boxZ; |
| 479 |
< |
|
| 476 |
> |
boxVector[0] = bsInfo.boxX; |
| 477 |
> |
boxVector[1] = bsInfo.boxY; |
| 478 |
> |
boxVector[2] = bsInfo.boxZ; |
| 479 |
> |
simnfo->setBox( boxVector ); |
| 480 |
> |
|
| 481 |
|
sprintf( simnfo->sampleName, "%s.dump", bsInfo.outPrefix ); |
| 482 |
|
sprintf( simnfo->finalName, "%s.init", bsInfo.outPrefix ); |
| 483 |
|
|
| 486 |
|
// set up the writer and write out |
| 487 |
|
|
| 488 |
|
writer = new DumpWriter( simnfo ); |
| 489 |
< |
writer->writeFinal(); |
| 489 |
> |
writer->writeFinal( 0.0 ); |
| 490 |
|
|
| 491 |
|
// clean up the memory |
| 492 |
|
|
| 674 |
|
bsInfo.boxY = nCells * cell; |
| 675 |
|
bsInfo.boxZ = nCells * cell; |
| 676 |
|
|
| 677 |
+ |
double boxVector[3]; |
| 678 |
|
simnfo = new SimInfo(); |
| 679 |
|
simnfo->n_atoms = nAtoms; |
| 680 |
< |
simnfo->box_x = bsInfo.boxX; |
| 681 |
< |
simnfo->box_y = bsInfo.boxY; |
| 682 |
< |
simnfo->box_z = bsInfo.boxZ; |
| 683 |
< |
|
| 680 |
> |
boxVector[0] = bsInfo.boxX; |
| 681 |
> |
boxVector[1] = bsInfo.boxY; |
| 682 |
> |
boxVector[2] = bsInfo.boxZ; |
| 683 |
> |
simnfo->setBox( boxVector ); |
| 684 |
> |
|
| 685 |
|
sprintf( simnfo->sampleName, "%s.dump", bsInfo.outPrefix ); |
| 686 |
|
sprintf( simnfo->finalName, "%s.init", bsInfo.outPrefix ); |
| 687 |
|
|
| 690 |
|
// set up the writer and write out |
| 691 |
|
|
| 692 |
|
writer = new DumpWriter( simnfo ); |
| 693 |
< |
writer->writeFinal(); |
| 693 |
> |
writer->writeFinal(0.0); |
| 694 |
|
|
| 695 |
|
// clean up the memory |
| 696 |
|
|