esp restrain


Click here for full thread
Gets Around
Hi

It seems to me that I have implemented (method 1). File with the modified source, nwchem.exe and example linked.

By analogy with the "radius" I introduced in "ESP" module a command "q0":

      q0 <integer iatnum> <real q0>
where iatnum is the atomic number for which a q0 will be used.

However, there are some doubts.

1. I just followed the analogy without real understanding what is written in the code.

First, I do not understand what an undocumented variable "q0" do.
It does something strange... plus activates RESP3, whatever that meant.
So I did not used it, but introduced a new variable "q00" (array).

Secondly, I do not understand what is the b(j) variable. I simply ignored it.

(My changes are marked in the text by commentary "NIKI")

2. I do not know how global variables work in Fortran. So I made a file niki.fh. Although it looks like working, however, it is an obvious hack.


Could you be so kind as to look at my source code?

(Sorry I am working under Windows, so executable is a Windows exe file.)

Quote:Bert Oct 12th 12:09 pm
Looking at the code, it's a restraint on an atom (lines 235-265) but it's set the same for each atom, and it only is invoked for carbon and hydrogen. I would think this can easily be modified to suite your needs. If you develop a new functionality, you're more then welcome to contribute it back to the open-source development effort.

Yes, the grid file is x y z V.

Thanks,

Bert