Re: [AMBER] MMPSA: PB Bomb in pb_aaradi(): ...

From: Jason Swails <>
Date: Tue, 7 Aug 2012 10:21:39 -0400

On Tue, Aug 7, 2012 at 9:53 AM, Vlad Cojocaru <> wrote:

> Jason,
> However, one mistery remains. Why using inp=2 and radiopt=0 still gives
> the same error but inp=1 and radiopt=0 does not ? As far as I got it,
> radiopt=0 should force to read radii from the topology file and
> therefore in should work with any inp value. Or maybe when setting inp=2
> the program overwrites automatically the radiopt=0 from the input file ?

It seems that your latter guess is right. From the PBSA source code:

   if ( inp == 2 .and. radiopt /= 1 ) then
      write(6, *) 'PB Warning in pb_read(): radiopt has been reset to 1
because inp=2.'
      radiopt = 1
   end if

Do you see that warning in the output file? In the long run, I think
modifying the source code to recognize ILDN atom types is the best
solution, but the fix I mentioned previously is a good
workaround for the time being.


Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Candidate
AMBER mailing list
Received on Tue Aug 07 2012 - 07:30:02 PDT
Custom Search