subroutine char key_1 char key_2 real DENSITY
! set rad sele segm name SEGMENTS end vdw
 poin surf all sele key_1 .a .not atom name H* end \
       sele  key_2 .a .not atom name H* end dens DENSITY init
return