#Load extrema1.pqr, extrema2.pqr ... produced by ESPpt_pqr.bat to show ESP extrema on vdW surface

set nsystem 1
for {set i 1} {$i<=$nsystem} {incr i} {
mol new extrema$i.pqr
mol modstyle 0 top VDW 0.07 20
mol modselect 0 top name C
mol modcolor 0 top ColorID 32
mol addrep top
mol modstyle 1 top VDW 0.07 20
mol modselect 1 top name O
mol modcolor 1 top ColorID 21
}