Talk:REFERENCE DATA SET

If the mtz file to be used for reference data set creation contains three digit reflection indices with negative signs (e.g., -102), the output from the original script will not be recognizable to XDS due to lack of space between reflection indices. An alternative method using mtz2various (below) avoids this issue.

#!/bin/bash

#starting mtz needs to be named refine.mtz with an F column called FP and a SIGF column called SIGFP
rm head.hkl temp.hkl I-model.mtz
echo square model amplitudes:
sftools <<EOF
read refine.mtz
select col FP = present
calc col I-model = col FP col FP *
write I-model.mtz col I-model SIGFP
quit
EOF
echo convert to ASCII format:
mtz2various hklin I-model.mtz hklout temp.hkl <<eof
LABIN FP=I-model SIGFP=SIGFP
RESOLUTION 50 1.8
OUTPUT USER '(3I5,7X,F9.2,5X,F6.2)'
END
eof
echo prepare I-model.hkl:
echo \!FORMAT=XDS_ASCII > head.hkl
echo \!END_OF_HEADER   >> head.hkl
cat head.hkl temp.hkl > I-model.hkl
rm head.hkl temp.hkl I-model.mtz