small bug in reporting

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
johnluick

small bug in reporting

#1 Unread post by johnluick »

I think that in read_phypar.F:

Code: Select all

# ifdef T_PASSIVE
          DO i=1,NPT
            itrc=inert(i)
            WRITE (out,185) LtracerSrc(i,ng), 'LtracerSrc', itrc,       &
     &            'Processing point sources/Sink on tracer ', itrc,     &
     &            TRIM(Vname(1,idTvar(itrc)))
          END DO
# endif 
LtracerSrc(i,ng) should be changed to LtracerSrc(NAT+i,ng)

It only affects reporting.

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

Re: small bug in reporting

#2 Unread post by arango »

Yes, but it is better to use itrc=inert(i) for those specific inert tracers. So we just need:

Code: Select all

 ifdef T_PASSIVE
          DO i=1,NPT
            itrc=inert(i)
            WRITE (out,185) LtracerSrc(itrc,ng), 'LtracerSrc', itrc,    &
     &            'Processing point sources/Sink on tracer ', itrc,     &
     &            TRIM(Vname(1,idTvar(itrc)))
          END DO
# endif
Thank you for reporting this typo.

Post Reply