*FILENAME: asa-protein.str
*PURPOSE: compute Lee-Richards type accessible surface area of a protein probe radius 1.4A
*AUTHOR: Lennart Nilsson, Karolinska Institutet, October 2003
*
!RETURNS:
! ASA computed for the selection defined as SOLUTE (which also gives the context to be used)
! is returned in variable ASA
! One could extend this to subdivide the calculation based on atom or aa type
!ASSUMPTIONS:
! psf and coordinates in MAIN coordinate set
! selection SOLUTE defined
!USAGE EXAMPLE
! read rtf,para,psf, coor
! define SOLUTE sele segid prot end
! stream asa-protein.str

coor surf rprobe 1.4 select SOLUTE end
! ASA per atom is now in WMAIN array, sum up over selected atoms
scalar wmain statistics select SOLUTE end
set ASA ?STOT
!scalar wmain statistics select "NONPOLAR" end
!set ASANP ?STOT

return


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden