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 no_save dihedral  P O5' C5' C4' angle2 angle22
secondary select active end no_save dihedral O5' C5' C4' C3' angle3 angle32
secondary select active end no_save dihedral C5' C4' C3' O3' angle4 angle42
secondary select active end no_save dihedral C4' C3' O3' +P angle5 angle52
secondary select active end no_save dihedral C3' O3' +P +O5' angle6 angle62
calc bond sele active .a atom name C4' O4' end dist = 2.0
image over false
return