small bug in reporting

[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.

[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.