[cctbxbb] Extracting data from MTZ

[David Waterman]

Hi all,

I have a script that extracts Fobs and Fcalc from an MTZ file produced by
Refmac for some analysis. I was just doing this in a trivial way:

      fobs = m.extract_reals("F").data
      fc = m.extract_reals("FC_ALL").data

and this worked fine, until I came across a file with missing observations,
so that the lengths of fobs and fc do not match.

Can anybody advise a better way to do this please? I'm inexperienced with
iotbx.mtz. I don't need anything other than the values of the structure
factors, but I need each element of fobs and fc to contain the equivalent
values.

Thanks

-- David
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://phenix-online.org/pipermail/cctbxbb/attachments/20180814/7d858ef8/attachment.htm>