model/src/do_statevars_diags.F

C $Header: /u/gcmpack/MITgcm/model/src/do_statevars_diags.F,v 1.1 2004/07/06 21:12:51 jmc Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"
#ifdef ALLOW_GCHEM
#include "GCHEM_OPTIONS.h"
#endif


CBOP
C     !ROUTINE: DO_STATEVARS_DIAGS
C     !INTERFACE:
      SUBROUTINE DO_STATEVARS_DIAGS( myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE DO_STATEVARS_DIAGS                                
C     | o Controlling routine for state variables diagnostics
C     *==========================================================*
C     | Computing statistics of the model state (state-variables) 
C     | is done at this level (after updating the halo region),
C     | as opposed to other diagnostic calls (fluxes, tendencies)
C     | that remain within the computation sequence.  
C     | Note: IO are not supposed to be done at this level 
C     |       but later (in DO_THE_MODEL_IO)
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myTime - Current time of simulation ( s )
C     myIter - Iteration number
C     myThid - Thread number for this instance of the routine.
      _RL     myTime
      INTEGER myIter
      INTEGER myThid

CEOP

C     !LOCAL VARIABLES:
C     == Local variables
      INTEGER bi,bj

#ifdef ALLOW_TIMEAVE
C--   Fill-in TimeAve pkg diagnostics (for state-variables)

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)

        IF ( taveFreq.GT.0. ) THEN
          CALL TIMEAVE_STATVARS(myTime, myIter, bi, bj, myThid)
# ifdef ALLOW_PTRACERS
          IF ( usePTRACERS ) THEN       
           CALL PTRACERS_STATVARS(myTime, myIter, bi, bj, myThid)
          ENDIF 
# endif /* ALLOW_PTRACERS */
        ENDIF

C--    End of bi,bj loop
       ENDDO
      ENDDO
#endif /* ALLOW_TIMEAVE */

#ifdef ALLOW_DIAGNOSTICS
C--   Fill-in Diagnostics pkg storage array (for state-variables)

      IF ( usediagnostics ) THEN
        CALL DIAGNOSTICS_FILL_STATE(myThid)
# ifdef ALLOW_PTRACERS
       IF (useptracers) THEN
        CALL PTRACERS_FILL_DIAGNOSTICS(myThid)
       ENDIF
# endif
      ENDIf

#endif /* ALLOW_DIAGNOSTICS */

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/do_statevars_diags.F,v 1.1 2004/07/06 21:12:51 jmc Exp $
2	C $Name: $
3
4	#include "PACKAGES_CONFIG.h"
5	#include "CPP_OPTIONS.h"
6	#ifdef ALLOW_GCHEM
7	#include "GCHEM_OPTIONS.h"
8	#endif
9
10
11	CBOP
12	C !ROUTINE: DO_STATEVARS_DIAGS
13	C !INTERFACE:
14	SUBROUTINE DO_STATEVARS_DIAGS( myTime, myIter, myThid )
15	C !DESCRIPTION: \bv
16	C ==========================================================
17	C \| SUBROUTINE DO_STATEVARS_DIAGS
18	C \| o Controlling routine for state variables diagnostics
19	C ==========================================================
20	C \| Computing statistics of the model state (state-variables)
21	C \| is done at this level (after updating the halo region),
22	C \| as opposed to other diagnostic calls (fluxes, tendencies)
23	C \| that remain within the computation sequence.
24	C \| Note: IO are not supposed to be done at this level
25	C \| but later (in DO_THE_MODEL_IO)
26	C ==========================================================
27	C \ev
28
29	C !USES:
30	IMPLICIT NONE
31	C == Global variables ===
32	#include "SIZE.h"
33	#include "EEPARAMS.h"
34	#include "PARAMS.h"
35
36	C !INPUT/OUTPUT PARAMETERS:
37	C == Routine arguments ==
38	C myTime - Current time of simulation ( s )
39	C myIter - Iteration number
40	C myThid - Thread number for this instance of the routine.
41	_RL myTime
42	INTEGER myIter
43	INTEGER myThid
44
45	CEOP
46
47	C !LOCAL VARIABLES:
48	C == Local variables
49	INTEGER bi,bj
50
51	#ifdef ALLOW_TIMEAVE
52	C-- Fill-in TimeAve pkg diagnostics (for state-variables)
53
54	DO bj=myByLo(myThid),myByHi(myThid)
55	DO bi=myBxLo(myThid),myBxHi(myThid)
56
57	IF ( taveFreq.GT.0. ) THEN
58	CALL TIMEAVE_STATVARS(myTime, myIter, bi, bj, myThid)
59	# ifdef ALLOW_PTRACERS
60	IF ( usePTRACERS ) THEN
61	CALL PTRACERS_STATVARS(myTime, myIter, bi, bj, myThid)
62	ENDIF
63	# endif /* ALLOW_PTRACERS */
64	ENDIF
65
66	C-- End of bi,bj loop
67	ENDDO
68	ENDDO
69	#endif /* ALLOW_TIMEAVE */
70
71	#ifdef ALLOW_DIAGNOSTICS
72	C-- Fill-in Diagnostics pkg storage array (for state-variables)
73
74	IF ( usediagnostics ) THEN
75	CALL DIAGNOSTICS_FILL_STATE(myThid)
76	# ifdef ALLOW_PTRACERS
77	IF (useptracers) THEN
78	CALL PTRACERS_FILL_DIAGNOSTICS(myThid)
79	ENDIF
80	# endif
81	ENDIf
82
83	#endif /* ALLOW_DIAGNOSTICS */
84
85	RETURN
86	END