pkg/gmredi/gmredi_diagnostics_impl.F

C $Header:  $
C $Name:  $

#include "GMREDI_OPTIONS.h"

CBOP
C     !ROUTINE: GMREDI_DIAGNOSTICS_IMPL
C     !INTERFACE:
      SUBROUTINE GMREDI_DIAGNOSTICS_IMPL(
     I                  myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE GMREDI_DIAGNOSTICS_IMPL
C     | o fill GM-Redi diagnostics for backward (implicit) terms
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global data ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"
#include "GMREDI.h"

C     !INPUT/OUTPUT PARAMETERS:
      _RL     myTime
      INTEGER myIter
      INTEGER myThid

#ifdef ALLOW_DIAGNOSTICS
C     !FUNCTIONS:
      LOGICAL     DIAGNOSTICS_IS_ON
      EXTERNAL    DIAGNOSTICS_IS_ON
#endif /* ALLOW_DIAGNOSTICS */
CEOP

#ifdef ALLOW_GMREDI
#ifdef ALLOW_DIAGNOSTICS

C     !LOCAL VARIABLES:
      INTEGER bi,bj
      INTEGER i,j,k
      _RL tmpVar(1-OLx:sNx+OLx,1-OLy:sNy+OLy)

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      IF ( useDiagnostics .AND. implicitDiffusion .AND.
     &     Nr.GE.2 .AND. tempStepping ) THEN

       IF ( DIAGNOSTICS_IS_ON('GM_KwzTz', myThid ) ) THEN
         DO bj=myByLo(myThid),myByHi(myThid)
          DO bi=myBxLo(myThid),myBxHi(myThid)
           CALL DIAGNOSTICS_COUNT( 'GM_KwzTz', bi, bj, myThid )
           DO k=2,Nr
             DO j = 1, sNy
              DO i = 1, sNx
                tmpVar(i,j) = -recip_drC(k)*rkSign
     &            *(theta(i,j,k,bi,bj)-theta(i,j,k-1,bi,bj))
     &            * maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
                tmpVar(i,j) = tmpVar(i,j)*Kwz(i,j,k,bi,bj)
     &                                   *rA(i,j,bi,bj)
              ENDDO
             ENDDO
             CALL DIAGNOSTICS_FILL( tmpVar, 'GM_KwzTz',
     &                              k, 1, 2, bi, bj, myThid )
           ENDDO
          ENDDO
         ENDDO
       ENDIF

      ENDIF

#endif /* ALLOW_DIAGNOSTICS */
#endif /* ALLOW_GMREDI */

      RETURN
      END
1	C $Header: $
2	C $Name: $
3
4	#include "GMREDI_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: GMREDI_DIAGNOSTICS_IMPL
8	C !INTERFACE:
9	SUBROUTINE GMREDI_DIAGNOSTICS_IMPL(
10	I myTime, myIter, myThid )
11	C !DESCRIPTION: \bv
12	C ==========================================================
13	C \| SUBROUTINE GMREDI_DIAGNOSTICS_IMPL
14	C \| o fill GM-Redi diagnostics for backward (implicit) terms
15	C ==========================================================
16	C \ev
17
18	C !USES:
19	IMPLICIT NONE
20	C === Global data ===
21	#include "SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PARAMS.h"
24	#include "GRID.h"
25	#include "DYNVARS.h"
26	#include "GMREDI.h"
27
28	C !INPUT/OUTPUT PARAMETERS:
29	_RL myTime
30	INTEGER myIter
31	INTEGER myThid
32
33	#ifdef ALLOW_DIAGNOSTICS
34	C !FUNCTIONS:
35	LOGICAL DIAGNOSTICS_IS_ON
36	EXTERNAL DIAGNOSTICS_IS_ON
37	#endif /* ALLOW_DIAGNOSTICS */
38	CEOP
39
40	#ifdef ALLOW_GMREDI
41	#ifdef ALLOW_DIAGNOSTICS
42
43	C !LOCAL VARIABLES:
44	INTEGER bi,bj
45	INTEGER i,j,k
46	_RL tmpVar(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
47
48	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
49
50	IF ( useDiagnostics .AND. implicitDiffusion .AND.
51	& Nr.GE.2 .AND. tempStepping ) THEN
52
53	IF ( DIAGNOSTICS_IS_ON('GM_KwzTz', myThid ) ) THEN
54	DO bj=myByLo(myThid),myByHi(myThid)
55	DO bi=myBxLo(myThid),myBxHi(myThid)
56	CALL DIAGNOSTICS_COUNT( 'GM_KwzTz', bi, bj, myThid )
57	DO k=2,Nr
58	DO j = 1, sNy
59	DO i = 1, sNx
60	tmpVar(i,j) = -recip_drC(k)*rkSign
61	& *(theta(i,j,k,bi,bj)-theta(i,j,k-1,bi,bj))
62	& * maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
63	tmpVar(i,j) = tmpVar(i,j)*Kwz(i,j,k,bi,bj)
64	& *rA(i,j,bi,bj)
65	ENDDO
66	ENDDO
67	CALL DIAGNOSTICS_FILL( tmpVar, 'GM_KwzTz',
68	& k, 1, 2, bi, bj, myThid )
69	ENDDO
70	ENDDO
71	ENDDO
72	ENDIF
73
74	ENDIF
75
76	#endif /* ALLOW_DIAGNOSTICS */
77	#endif /* ALLOW_GMREDI */
78
79	RETURN
80	END