*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
! 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
! psf and coordinates in MAIN coordinate set
! selection SOLUTE defined
! 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
!scalar wmain statistics select "NONPOLAR" end


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden