Yes, but the correct solution is to have:
( defined SALINITY && !defined ANA_SSFLUX && \
(defined BULK_FLUXES && !defined EMINUSP)) || \
since the freshwater flux is only computed in bulk_flux.F when EMINUSP is activated.
Also corrected a typo in ana_passive.h. We need to have instead:
# ifdef AGE_MEAN
DO ip=1,NPT,2
itrc=inert(ip)
iage=inert(ip+1)
DO k=1,N(ng)
DO j=JstrT,JendT
DO i=IstrT,IendT
t(i,j,k,1,itrc)=???
t(i,j,k,2,itrc)=t(i,j,k,1,itrc)
t(i,j,k,1,iage)=0.0_r8
t(i,j,k,2,iage)=t(i,j,k,1,iage)
END DO
END DO
END DO
END DO
# else
...
# endif
The first DO-loop need to have a step-size of 2. Otherwise, the age concentration is overwritten. Many thanks to John Wilkin for bringing this to my attention.