pkg/diagnostics/diagnostics_is_on.F

C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_is_on.F,v 1.3 2005/06/26 16:51:49 jmc Exp $
C $Name:  $

#include "DIAG_OPTIONS.h"

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

CBOP 0
C     !ROUTINE: DIAGNOSTICS_IS_ON

C     !INTERFACE:
      LOGICAL FUNCTION DIAGNOSTICS_IS_ON( diagName, myThid )

C     !DESCRIPTION:
C     *==========================================================*
C     | FUNCTION DIAGNOSTIC_IS_ON
C     | o Return TRUE if diagnostics "diagName" is Active
C     *==========================================================*

C     !USES:
      IMPLICIT NONE
#include "EEPARAMS.h"
#include "SIZE.h"
#include "DIAGNOSTICS_SIZE.h"
#include "DIAGNOSTICS.h"

C     !INPUT PARAMETERS:
C     diagName   ::  diagnostic identificator name (8 characters long)
C     myThid     ::  my thread Id number
      CHARACTER*8  diagName
      INTEGER      myThid
CEOP

C     !LOCAL VARIABLES:
      INTEGER n,m
      INTEGER ip,iSp
      INTEGER bi,bj

      DIAGNOSTICS_IS_ON = .FALSE.
      bi = myBxLo(myThid)
      bj = myByLo(myThid)

C-    search for this diag. in the active 2D/3D diagnostics list
      DO n=1,nlists
       DO m=1,nActive(n)
        IF ( diagName.EQ.flds(m,n) .AND. idiag(m,n).NE.0 ) THEN
          ip = ABS(idiag(m,n))
          IF ( ndiag(ip,bi,bj).GE.0 ) DIAGNOSTICS_IS_ON = .TRUE.
        ENDIF
       ENDDO
      ENDDO

      IF ( .NOT.DIAGNOSTICS_IS_ON ) THEN
C-    search for this diag. in the active statistics-diagnostics list
       DO n=1,diagSt_nbLists
        DO m=1,diagSt_nbActv(n)
         IF ( diagName.EQ.diagSt_Flds(m,n)
     &               .AND.iSdiag(m,n).NE.0 ) THEN
           iSp = ABS(iSdiag(m,n))
           IF (qSdiag(0,0,iSp,bi,bj).GE.0.) DIAGNOSTICS_IS_ON = .TRUE.
         ENDIF
        ENDDO
       ENDDO
      ENDIF

      RETURN
      END
1	jmc	1.4	C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_is_on.F,v 1.3 2005/06/26 16:51:49 jmc Exp $
2	jmc	1.1	C $Name: $
3
4			#include "DIAG_OPTIONS.h"
5
6			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7
8			CBOP 0
9			C !ROUTINE: DIAGNOSTICS_IS_ON
10
11			C !INTERFACE:
12			LOGICAL FUNCTION DIAGNOSTICS_IS_ON( diagName, myThid )
13
14			C !DESCRIPTION:
15			C ==========================================================
16			C \| FUNCTION DIAGNOSTIC_IS_ON
17			C \| o Return TRUE if diagnostics "diagName" is Active
18			C ==========================================================
19
20			C !USES:
21			IMPLICIT NONE
22			#include "EEPARAMS.h"
23			#include "SIZE.h"
24			#include "DIAGNOSTICS_SIZE.h"
25			#include "DIAGNOSTICS.h"
26
27			C !INPUT PARAMETERS:
28			C diagName :: diagnostic identificator name (8 characters long)
29			C myThid :: my thread Id number
30			CHARACTER*8 diagName
31			INTEGER myThid
32			CEOP
33
34			C !LOCAL VARIABLES:
35	jmc	1.3	INTEGER n,m
36			INTEGER ip,iSp
37	jmc	1.4	INTEGER bi,bj
38	jmc	1.1
39			DIAGNOSTICS_IS_ON = .FALSE.
40	jmc	1.4	bi = myBxLo(myThid)
41			bj = myByLo(myThid)
42	jmc	1.2
43			C- search for this diag. in the active 2D/3D diagnostics list
44	jmc	1.1	DO n=1,nlists
45			DO m=1,nActive(n)
46	jmc	1.3	IF ( diagName.EQ.flds(m,n) .AND. idiag(m,n).NE.0 ) THEN
47			ip = ABS(idiag(m,n))
48	jmc	1.4	IF ( ndiag(ip,bi,bj).GE.0 ) DIAGNOSTICS_IS_ON = .TRUE.
49	jmc	1.1	ENDIF
50			ENDDO
51			ENDDO
52
53	jmc	1.2	IF ( .NOT.DIAGNOSTICS_IS_ON ) THEN
54			C- search for this diag. in the active statistics-diagnostics list
55			DO n=1,diagSt_nbLists
56			DO m=1,diagSt_nbActv(n)
57	jmc	1.4	IF ( diagName.EQ.diagSt_Flds(m,n)
58	jmc	1.3	& .AND.iSdiag(m,n).NE.0 ) THEN
59			iSp = ABS(iSdiag(m,n))
60	jmc	1.4	IF (qSdiag(0,0,iSp,bi,bj).GE.0.) DIAGNOSTICS_IS_ON = .TRUE.
61	jmc	1.3	ENDIF
62	jmc	1.2	ENDDO
63			ENDDO
64			ENDIF
65
66	jmc	1.1	RETURN
67			END