#!/bin/csh -f
#      
# self rotation function
#
polarrfn hklin ../gstp1.mtz \
mapout self.map \
PLOT self.plt \
>self.out <<EOF
SELF 18.0 		!Patterson radius 
LABIN FILE 1 F=FP SIGF=SIGFP 
RESOLUTION 20.0 2.5 
LIMITS 0 180 2.5   0 180 2.5   0 180 10.0 
NOPRINT 
MAP 
FIND 20 40 OUTPUT selfpeak.list
PLOT 30.0 10.0 !contour start and increment 10 10 is OK
CRYSTAL FILE 1 ORTH 1 BFAC 26.4 FLIMITS 20 1700 
EOF
#
#
pltdev -lni 0.2 -i self.plt -o self_all.ps
#
#pltdev -i self.plt -o self_all.ps
#
#rm self.plt
#
