You are absolutely right! The difficulty in getting from MTZ to any
other format or back is unacceptable. Expecting working
crystallographers to write Fortran format statements is ridiculous. I've
been trying to address this by adding support for other formats to
clipper, but my pace has been glacial, owing to the limited academic
rewards for such work. Even then, it needs someone to put together an
application and GUI to use it with.
I will try and remember to raise this issue at the CCP4 developers
meeting this year. If you prod me in April and I'll report back to the BB.
Kevin
George M. Sheldrick wrote:
Dear Ian,
It is all part of a diabolical CCP4 plot to make it as inconvenient as
possible to move from a REFMAC refinement to SHELXL! I hope that i do
not get excommunicated like DVD for this comment.
To summarize your and Martin's suggestions, I know of only two ways to
move from mtz to hkl:
=====================================================================
1] run mtz2various with the keywords:
LABIN FP=FP SIGFP=SIGFP FREE=FreeR_flag
OUTPUT SHELX
FSQUAR
then edit the resulting .hkl file with any text editor to remove the
header (which SHELXL cannot read) and all occurences of the word 'FREE'.
I have from time to time suggested that mtz2various be changed so that
it no longer outputs the extra junk, but it is still there (at least in
CCP4-6.0.2).
This produces an 'inferior' SHELXL HKLF4 format file because the
intensities and their standard deviations have been converted to F and
back to I, which requires the assumption of an intensity distribution
function that may not be completely valid (e.g. if NCS is present) and
so produces inferior sigma(I) values. In practice one can live with
this, but it is not very scientific.
=====================================================================
2] Use Tim Gruene's mtz2sca to convert from .mtz to .sca.
mtz2sca name
reads name.mtz and writes name.sca. This program is avaliable from Tim
or from the SHELX download area (Linux only). This has the advantage
that it uses intensities if it can find them in the .mtz (e.g. if they
are still there from SCALA) and otherwise uses F's. mtz2sca is fine if
you want to do exprimental phasing (e.g. read the data into SHELXC
either directly or via Thomas Schneider's hkl2map) but unfortunately
.sca format does not know about Free-R flags. So, as Martin suggested,
one can read the .sca file into xprep (start xprep without a file name
and give the full name of the .sca file when prompted) and then use the
xprep option to transfer the Free-R flags from the 'inferior' .hkl file
from mtz2various and write it out in SHELXL HKLF4 format. This produces
a slightly better file for refinement when the .mtz file contains
intensities, but requires xprep.
====================================================================
George
Martin Hallberg wrote:
I would probably prioritize keeping the same R-free set (thus using
the F^2 output by mtz2various) over
going through the scalepack format and loosing track of it.
You can however use XPREP to transfer the R-free set from the
"inferior" F^2 HKLF4 file output by mtz2various to the proper
intensities read from the scalepack file generated by mtz2sca. Then
you can write out an HKLF4 file for refinement in SHELXL-97. That
HKLF4 file will have proper intensities and the same Rfree set that
you have used in Refmac5 (or perhaps Restrain?). But then one would
wish that mtz2various did this from the beginning...
Best regards,
Martin
On Mar 1, 2007, at 7:58 PM, Ethan Merritt wrote:
On Thursday 01 March 2007 10:17, Ian Tickle wrote:
All, I thought this would be a simple task, but for the life of me I
can't see how to do it! All I want to do is convert an MTZ file to
Shel-X format for refinement. I thought it would take me 2 secs, but
it's taken me at least 5 attempts, and it's still not right!
Do the conversion on the shelx side, rather than the CCP4 side.
http://shelx.uni-ac.gwdg.de/~tg/mtz2sca/mtz2sca.html
.
B. Martin Hallberg, PhD
Molecular Cell Biology Program
Department of Cell and Molecular Biology
Karolinska Institute
Nobels väg 3
SE-171 77 Stockholm