subroutine \
	real angle1 real angle12 real angle2 real angle22 \
	real angle3 real angle32 real angle4 real angle42 \
	real angle5 real angle52 real angle6 real angle62 

dele bond sele active .a atom name C4' O4' end
secondary select active end dihedral -O3' P O5' C5' angle1 angle12
secondary select active end false dihedral  P O5' C5' C4' angle2 angle22
secondary select active end false dihedral O5' C5' C4' C3' angle3 angle32
secondary select active end false dihedral C5' C4' C3' O3' angle4 angle42
secondary select active end false dihedral C4' C3' O3' +P angle5 angle52
secondary select active end false dihedral C3' O3' +P +O5' angle6 angle62
calc bond sele active .a atom name C4' O4' end dist = 2.0
image over false
return