Dear all,
I want to read the source code for GB/SA simulations and
find file egb.f. In the subroutine egb(x,f,...,...), I find the "if(gbsa ==
2) then " and "if (gbsa ==1) then" . I think they are use to control
different models. Now what I want to do is to add another "if" sentence to
use another model, e.g. "if (gbsa == 3) then call XXX"(XXX is the
subroutine for my model). But when I change the code it seems nothing
happens. For example, I add "esurf = 1.0d0" after "# include 'gbsa.h'",but
when I run "minimization", the result shown in "mdinfo" is the same as
before.(Actually, no matter what I add to change the value of esurf and
totsasa after "# include 'gbsa.h'" , the value of "esurf" shown in "mdinfo"
is the same.) My input file "min.in"for minimization is as follows:
cat
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Aug 14 2012 - 16:30:03 PDT
