本帖最后由 yaonaite 于 2026-6-17 12:38 编辑
相较于上一版shelx2019,shelx2025版本能够更好地处理孪晶,使得一些原本无法使用的单晶数据有了抢救的可能。
然而,当测得的单晶数据中有孪晶时,或存在重原子以至于旁边有鬼峰时,总有一些Q峰需要使用OMIT命令处理,一大堆坐标实在是难以手敲。
不过这个过程非常机械,非常适合用脚本自动化。以下提供convert_omit.py和convert_multiple_omit.py,分别处理一个和多个数据。
把需要处理的Q相关的行写进xxxxomit.txt或xxxxOMIT.txt中即可处理。设置了工作目录为D:\forrings,可以很方便地在脚本最前一段修改。
示例如下:
输入:
Q1 35 (-4,8,10), (4,-8,-10)
输出:
Q1 35 (-4,8,10), (4,-8,-10)
OMIT -4 8 10
OMIT 4 -8 -10
输入:
Q1 35 (-4,8,10), (4,-8,-10)
Q2 40 (-2,-1,1)
Q3 47 (13,10,-4), (-4,10,-12)
输出:
Q1 35 (-4,8,10), (4,-8,-10)
OMIT -4 8 10
OMIT 4 -8 -10
Q2 40 (-2,-1,1)
OMIT -2 -1 1
Q3 47 (13,10,-4), (-4,10,-12)
OMIT 13 10 -4
OMIT -4 10 -12
|