Hi,
I'm trying to build a system with an Fe2+ center, and using OPC water. The
manual states that, in AmberTools 22, MCPB.py should be compatible with OPC
water.
However, when trying to build the system with OPC water, I get
the following message:
<--
$ MCPB.py -i MCPB.in -s 1
[...]
Totally there are 105 atoms in the large model.
Totally there are 450 electrons in the large model.
Traceback (most recent call last):
File "/opt/miniconda3/envs/AmberTools22/bin/MCPB.py", line 651, in
<module>
gene_model_files(orpdbf, ionids, addres, addbpairs, gname, ff_choice,
File
"/opt/amber/current/lib/python3.7/site-packages/pymsmt/mcpb/gene_model_files.py",
line 1909, in gene_model_files
build_large_model(mol, reslist, lmsresids, lmsresace, lmsresnme,
lmsresgly,
File
"/opt/amber/current/lib/python3.7/site-packages/pymsmt/mcpb/gene_model_files.py",
line 1600, in build_large_model
IonLJParaDict = get_ionljparadict(watermodel)
File
"/opt/amber/current/lib/python3.7/site-packages/pymsmt/mol/element.py",
line 325, in get_ionljparadict
IonLJParaDict.update(diljpara2)
UnboundLocalError: local variable 'diljpara2' referenced before assignment
-->
If I just replace the OPC model by TIP3P in the water_model line, it works
just fine.
--
Gustavo Seabra.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Oct 13 2022 - 15:00:02 PDT
