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