source leaprc.water.$SOLVENT_MODEL
mol = loadpdb $INPUT_PDB_FILENAME
addionsrand mol Na+ $NUM_NA_IONS Cl- $NUM_CL_IONS $SEPARATION
savepdb mol $OUTPUT_PDB_FILENAME
quit
