Re: [AMBER] MMGBSA bad atom type error not fixed by source code edits

From: Matthew Guberman-Pfeffer <matthew.guberman-pfeffer.uconn.edu>
Date: Sat, 3 Jul 2021 08:12:09 -0400

Hi Dr. Case,

Thanks, I was doing what no one should—copying and pasting the code without learning what the syntax meant. I don’t know Fortran (yet).
I’ll try removing the (1:1) so it checks against the two characters.

Best,
Matthew

> On Jul 3, 2021, at 8:04 AM, David A Case <david.case.rutgers.edu> wrote:
>
> *Message sent from a system outside of UConn.*
>
>
> On Sat, Jul 03, 2021, Matthew Guberman-Pfeffer wrote:
>>
>> I edited amber20_src/AmberTools/src/sander/mdread2.F90 to contain the
>> following blocks for Fe, as recommended in a previous post. I still get the
>> “Bad Atom Type” message after recompiling. What can I do to fix this issue?
>>
>> else if (atype(1:1) == 'FE' .or. atype(1:1) == 'Fe' .or. atype(1:1) == 'Fe') then
>
> This if branch will never be taken. "atype(1:1)" gives the first character
> of the atype string. It can never match a two-character constant like "FE"
> or "Fe".
>
> ....dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Flists.ambermd.org%2Fmailman%2Flistinfo%2Famber&amp;data=04%7C01%7C%7C1451f77741f74807cf0908d93e1abe76%7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C637609106846452121%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&amp;sdata=Rsqn1%2FWUppk8WupibM3selOgTMR7lUCPo%2B828aWwo%2BI%3D&amp;reserved=0


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Jul 03 2021 - 05:30:03 PDT
Custom Search