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	C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_is_on.F,v 1.3 2005/06/26 16:51:49 jmc Exp $
2	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	INTEGER n,m
36	INTEGER ip,iSp
37	INTEGER bi,bj
38
39	DIAGNOSTICS_IS_ON = .FALSE.
40	bi = myBxLo(myThid)
41	bj = myByLo(myThid)
42
43	C- search for this diag. in the active 2D/3D diagnostics list
44	DO n=1,nlists
45	DO m=1,nActive(n)
46	IF ( diagName.EQ.flds(m,n) .AND. idiag(m,n).NE.0 ) THEN
47	ip = ABS(idiag(m,n))
48	IF ( ndiag(ip,bi,bj).GE.0 ) DIAGNOSTICS_IS_ON = .TRUE.
49	ENDIF
50	ENDDO
51	ENDDO
52
53	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	IF ( diagName.EQ.diagSt_Flds(m,n)
58	& .AND.iSdiag(m,n).NE.0 ) THEN
59	iSp = ABS(iSdiag(m,n))
60	IF (qSdiag(0,0,iSp,bi,bj).GE.0.) DIAGNOSTICS_IS_ON = .TRUE.
61	ENDIF
62	ENDDO
63	ENDDO
64	ENDIF
65
66	RETURN
67	END