a small bug in inp_par.F

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
xiaozhu557
Posts: 62
Joined: Fri Sep 11, 2009 1:48 pm
Location: nmefc

a small bug in inp_par.F

#1 Post by xiaozhu557 »

Hello all,
I found a small bug in line 855 of the file inp_par.F.
The original version is

Code: Select all

854    #ifdef SOLVE3D
855            IF (LnudgeM3CLM(ng).or.                                         &
856         &      LnudgeM3CLM(ng).or.                                         &
857         &      ANY(LnudgeTCLM(:,ng))) THEN
858              Lnudging(ng)=.TRUE.
859            END IF
860    #else
861            IF (LnudgeM2CLM(ng)) THEN
862              Lnudging(ng)=.TRUE.
863            END IF
864    #endif
I think that the right version should be

Code: Select all

854    #ifdef SOLVE3D
855            IF (LnudgeM2CLM(ng).or.                                         &
856         &      LnudgeM3CLM(ng).or.                                         &
857         &      ANY(LnudgeTCLM(:,ng))) THEN
858              Lnudging(ng)=.TRUE.
859            END IF
860    #else
861            IF (LnudgeM2CLM(ng)) THEN
862              Lnudging(ng)=.TRUE.
863            END IF
864    #endif
Am I right?

User avatar
kate
Posts: 3893
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: a small bug in inp_par.F

#2 Post by kate »

Looks to me like you are right. Thanks!

User avatar
arango
Site Admin
Posts: 1138
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: a small bug in inp_par.F

#3 Post by arango »

Yes, good catch. Thank you.

Post Reply