Hello,<br><br>I am seeing an internal error while trying to conduct joint x-ray/neutron refinement, and any help with the issue below would be greatly appreciated. Thank you very much in advance:<br>
<br>
================================ ADP refinement ===============================<br><div id=":5g">
<br>
Traceback (most recent call last):<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/command_line/refine.py&quot;,<br>
line 11, in &lt;module&gt;<br>
    command_line.run(command_name=&quot;phenix.refine&quot;, args=sys.argv[1:])<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/refinement/command_line.py&quot;,<br>
line 89, in run<br>
    call_back_handler=call_back_handler)<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/refinement/driver.py&quot;,<br>
line 1182, in run<br>
    call_back_handler       = call_back_handler)<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/phenix/phenix/refinement/strategies.py&quot;,<br>
line 735, in refinement_machine<br>
    h_params               = h_params)<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/cctbx_project/mmtbx/refinement/adp_refinement.py&quot;,<br>
line 192, in __init__<br>
    all_params = all_params)<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/cctbx_project/mmtbx/refinement/adp_refinement.py&quot;,<br>
line 216, in refine_adp<br>
    model.set_refine_individual_adp()<br>
  File &quot;/sware/phenix-1.6.4-486/linux64/phenix-1.6.4-486/cctbx_project/mmtbx/model.py&quot;,<br>
line 1035, in set_refine_individual_adp<br>
    iselection = selection_iso.iselection())<br>
RuntimeError: cctbx InternalError:<br>
/net/anaconda/scratch1/phenix/phenix-1.6.4-486/cctbx_project/cctbx/xray/scatterer_flags.h(417):<br>
CCTBX_ASSERT(f.use_u_iso()) failure.<br>
<br>
<br>
<br>
This error occurs during the 2nd macro cycle of refinement using this<br>
general command:<br>
<br>
phenix.refine model.pdb data_xray.hkl neutron_data.file_name=data_neutron.hkl<br>
input.xray_data.labels=FOBSx input.neutron_data.labels=FOBSn<br>
<br>
<br>
<br>
Thanks again.<br>
<br>
Regards,<br>
<br>
Louis</div>