Dear CCP4BB,
I'm trying to genereate some self rotations with polarrfn to decide on the NCS
and number of molecules in my AU.
using polarrfn script (copied from the examples on the web):
#
# Self-rotation function
#
set resolution = 15 3.5
set radius = 20
# Calculate whole map, & search for peaks
polarrfn hklin p2221_33A mapout p2221_self << eof-1
title p2221 Self-Rotation function, resolution ${resolution} radius ${radius}
self ${radius}
resolution ${resolution}
crystal file 1 bfac -20
labin file 1 F=F SIGF=SIGF
map
find 5 100
eof-1
# Now plot just the kappa = 180 deg section
plot:
polarrfn hklin p2221_33A mapin p2221_self plot p2221_self << eof-2
title p2221 Self-Rotation function, resolution ${resolution} radius ${radius}
self ${radius}
resolution ${resolution}
crystal file 1 bfac -20
labin file 1 F=F SIGF=SIGF
read 180 180
plot 10 10
eof-2
As you can guess, I think I'm in p2221... The script seems to run fine and
spits out map and .plo files.
However, when I try to run pltdev to generate a postscript file (this used to
work in Unix in a previous lab...)
> pltdev -i p2221_self.plo
I get a file (plot84.ps) with lots of nan records (see below).
I'm running Suse 10.3, CCP4 6.0.2
Any pointers in the right direction much appreciated.
Thanks
Ed
%!PS-Adobe-3.0
%%Creator: pltdev
%%Pages: 1
%%BoundingBox: -2147483648 -2147483648 -2147483648 -2147483648
%%EndComments
%%BeginProlog
/L {lineto} bind def
/M {moveto} bind def
/S {stroke} bind def
/N {newpath} bind def
%%EndProlog
%%Page: 1 1
%%PageBoundingBox: -2147483648 -2147483648 -2147483648 -2147483648
0.8 setlinewidth 1 setlinejoin N 0 0 M
S N nan nan M
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
nan nan L
