[AMBER] too many titratable groups

From: <piotrek.vega.umcs.lublin.pl>
Date: Fri, 20 Aug 2010 20:26:42 +0200 (CEST)

Dear Amber users,


 I'm still learing how to use Amber, so sorry if the question is too silly.

 I'm running a const-pH calculations on a large protein, and it looks like
I reached the limit of titratable groups, crash with error:

 forrtl: severe (18): too many values for NAMELIST variable, unit 18, file
/home/d3x582/PIOTR/test_const_ph_even_more_titratable_groups/ph_5/cpin,
line 40, position 76


and lines 40,41 in cpin are:

'Residue: GL4 341','Residue: AS4 366','Residue: AS4 368','Residue: LYS 381',
 'Residue: AS4 382','Residue: AS4 387','Residue: GL4 391','Residue: GL4 394',


I am able to run calculations with smaller number of titrable groups (e.g.
only HIS residues, or HIS,GLU). I did not see any fixed-size arrays in
constph.f file, the value for the titratable gropus counter (trescnt) is
set by cpin, so ... how should I modify amber code to increase the number
of protonatable-residues.

 Thank you for help,
 Piotrek


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Aug 20 2010 - 11:30:03 PDT
Custom Search