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	mlosch	1.2	C $Header: /u/gcmpack/MITgcm/model/src/do_statevars_diags.F,v 1.1 2004/07/06 21:12:51 jmc Exp $
2	jmc	1.1	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	mlosch	1.2	# ifdef ALLOW_PTRACERS
60			IF ( usePTRACERS ) THEN
61			CALL PTRACERS_STATVARS(myTime, myIter, bi, bj, myThid)
62			ENDIF
63			# endif /* ALLOW_PTRACERS */
64	jmc	1.1	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