read 
 mol pqr comp13.pqr #System 
 mol pqr Protein13.pqr #Protein 
 mol pqr BEC13.pqr #BEC 
end
APOLAR name Hcomp13 
 grid 0.1 0.1 0.1 #System 
 mol 1 #System 
 srfm sacc #System 
 swin 0.3 #System 
 srad 1.4 #System 
 press 0.000 #System
 gamma 0.0227 #System 
 bconc 0.00 #System 
 sdens 10.0 #System 
 dpos 0.20 #System 
 temp 300 #System 
 calcenergy total #System 
 calcforce total #System 
end 

APOLAR name HProtein13 
 grid 0.1 0.1 0.1 #Protein 
 mol 2 #Protein 
 srfm sacc #Protein 
 swin 0.3 #Protein 
 srad 1.4 #Protein 
 press 0.000 #Protein
 gamma 0.0227 #Protein 
 bconc 0.00 #Protein 
 sdens 10.0 #Protein 
 dpos 0.20 #Protein 
 temp 300 #Protein 
 calcenergy total #Protein 
 calcforce total #Protein 
end 

APOLAR name HBEC13 
 grid 0.1 0.1 0.1 #BEC 
 mol 3 #BEC 
 srfm sacc #BEC 
 swin 0.3 #BEC 
 srad 1.4 #BEC 
 press 0.000 #BEC 
 gamma 0.0227 #BEC 
 bconc 0.00 #BEC 
 sdens 10.0 #BEC 
 dpos 0.20 #BEC 
 temp 300 #BEC 
 calcenergy total #BEC 
 calcforce total #BEC 
end 

print apolEnergy Hcomp13 - HProtein13 - HBEC13 end 
quit
