source leaprc.protein.ff19SB
source leaprc.phosaa19SB
loadamberparams frcmod.ff19SB
source leaprc.gaff2

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
