[MITgcm-devel] monitor for ptracers

[Martin Losch]

Hi,
currently the monitor statistics for the ptracers package (for which Ed 
heroically wrote a test in testreport) are computed and written in a 
routine called GCHEM_MONITOR that is called from GCHEM_DIAGS. Because 
the ptracers-statistics are not only interesting for the GCHEM package, 
I would like to move their computation to a place that is more general 
(and rename GCHEM_MONITOR to PTRACERS_MONITOR). Where would that place 
be. Basically two options:
1. The simple, but not very clean option is to put the calls to 
mon_printstats_rl inside MONITOR inside #ifdef, such as:
> #ifdef ALLOW_PTRACERS
> C       Print stats for passive tracers
>         IF ( usePTRACERS ) THEN
>          CALL MON_SET_PREF('trcstat',myThid)
>          DO iTracer = 1, PTRACERS_numInUse
>           WRITE(suff,'(A8,I2.2)') '_ptracer',iTracer
>           CALL MON_PRINTSTATS_RL(
>      &         Nr,pTracer(1-Olx,1-Oly,1,1,1,iTracer),suff,
>      &         maskC,hFacC,rA ,drF,myThid)
>          ENDDO
>         ENDIF
> #endif /* ALLOW_PTRACERS */
This works, and I wouldn't have to worry about "mnc-ing" the ptracers 
statistics, but you then have trcstat_ptracer??_XXX within the "dynamic 
field statistics".
2. More cleanly, call a separate PTRACERS_MONITOR somewhere, but where? 
I would put this call into the MONITOR routine (at the end), but is 
this the right place?

Martin