pkg/monitor/mon_stats_rl.F

C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rs.F,v 1.9 2009/04/28 18:16:53 jmc Exp $
C $Name:  $

#include "MONITOR_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP
C     !ROUTINE: MON_STATS_RL

C     !INTERFACE:
      SUBROUTINE MON_STATS_RL(
     I     myNr, arr,
     O     theMin,theMax,theMean,theSD,
     I     myThid )

C     !DESCRIPTION:
C     Calculate bare statistics of global array ``\_RL arr''.

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

C     !INPUT PARAMETERS:
      INTEGER myNr
      _RL arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
      _RL theMin, theMax, theMean, theSD
      INTEGER myThid
CEOP

C     !LOCAL VARIABLES:
      INTEGER bi,bj,I,J,K
      INTEGER numPnts
      LOGICAL noPnts
      _RL tmpVal,rNumPnts
      _RL theVar
      _RL tileMean(nSx,nSy)
      _RL tileVar (nSx,nSy)
      _RL tileSD  (nSx,nSy)

      theMin = 0.
      theMax = 0.
      theMean= 0.
      theSD  = 0.
      theVar = 0.
      numPnts= 0
      noPnts = .TRUE.

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        tileMean(bi,bj) = 0.
        tileVar (bi,bj) = 0.
        DO K=1,myNr
         DO J=1,sNy
          DO I=1,sNx
           tmpVal=arr(I,J,K,bi,bj)
c          IF (tmpVal.NE.0. .AND. noPnts) THEN
           IF ( noPnts ) THEN
            theMin = tmpVal
            theMax = tmpVal
            noPnts = .FALSE.
           ENDIF
c          IF (tmpVal.NE.0.) THEN
            theMin = MIN(theMin,tmpVal)
            theMax = MAX(theMax,tmpVal)
            tileMean(bi,bj)=tileMean(bi,bj)+tmpVal
            tileVar (bi,bj)=tileVar (bi,bj)+tmpVal*tmpVal
            numPnts=numPnts+1
c          ENDIF
          ENDDO
         ENDDO
        ENDDO
c       theMean=theMean+tileMean(bi,bj)
c       theVar =theVar +tileVar (bi,bj)
       ENDDO
      ENDDO

c     _GLOBAL_SUM_RL(theMean,myThid)
c     _GLOBAL_SUM_RL(theVar,myThid)
      CALL GLOBAL_SUM_TILE_RL( tileMean, theMean, myThid )
      CALL GLOBAL_SUM_TILE_RL( tileVar , theVar , myThid )
      tmpVal=FLOAT(numPnts)
      _GLOBAL_SUM_RL(tmpVal,myThid)
      numPnts=NINT(tmpVal)

      IF (tmpVal.GT.0.) THEN
       rNumPnts=1. _d 0/tmpVal
       theMean=theMean*rNumPnts
       theVar=theVar*rNumPnts
       IF ( noPnts ) theMin = theMean
       theMin=-theMin
       _GLOBAL_MAX_RL(theMin,myThid)
       theMin=-theMin
       IF ( noPnts ) theMax = theMean
       _GLOBAL_MAX_RL(theMax,myThid)

       DO bj=myByLo(myThid),myByHi(myThid)
        DO bi=myBxLo(myThid),myBxHi(myThid)
         tileSD(bi,bj)=0.
         DO K=1,myNr
          DO J=1,sNy
           DO I=1,sNx
            tmpVal=arr(I,J,K,bi,bj)
c           IF (tmpVal.NE.0.) THEN
             tileSD(bi,bj) = tileSD(bi,bj)
     &                     + (tmpVal-theMean)*(tmpVal-theMean)
c           ENDIF
           ENDDO
          ENDDO
         ENDDO
c        theSD = theSD + tileSD(bi,bj)
        ENDDO
       ENDDO

c      _GLOBAL_SUM_RL(theSD,myThid)
       CALL GLOBAL_SUM_TILE_RL( tileSD, theSD, myThid )

       theSD = SQRT(theSD*rNumPnts)
c      theSD = SQRT(theVar-theMean*theMean)
      ENDIF

      RETURN
      END
1	jmc	1.14	C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rs.F,v 1.9 2009/04/28 18:16:53 jmc Exp $
2	adcroft	1.2	C $Name: $
3	adcroft	1.1
4	adcroft	1.5	#include "MONITOR_OPTIONS.h"
5	adcroft	1.1
6	edhill	1.10	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7			CBOP
8			C !ROUTINE: MON_STATS_RL
9
10			C !INTERFACE:
11	adcroft	1.1	SUBROUTINE MON_STATS_RL(
12	jmc	1.14	I myNr, arr,
13			O theMin,theMax,theMean,theSD,
14	edhill	1.10	I myThid )
15
16	jmc	1.14	C !DESCRIPTION:
17	edhill	1.10	C Calculate bare statistics of global array ``\_RL arr''.
18
19			C !USES:
20	adcroft	1.1	IMPLICIT NONE
21			#include "SIZE.h"
22			#include "EEPARAMS.h"
23
24	edhill	1.10	C !INPUT PARAMETERS:
25	adcroft	1.1	INTEGER myNr
26			_RL arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
27	jmc	1.14	_RL theMin, theMax, theMean, theSD
28	adcroft	1.1	INTEGER myThid
29	edhill	1.10	CEOP
30	adcroft	1.1
31	edhill	1.10	C !LOCAL VARIABLES:
32	adcroft	1.1	INTEGER bi,bj,I,J,K
33			INTEGER numPnts
34			LOGICAL noPnts
35			_RL tmpVal,rNumPnts
36	jmc	1.12	_RL theVar
37			_RL tileMean(nSx,nSy)
38			_RL tileVar (nSx,nSy)
39			_RL tileSD (nSx,nSy)
40	adcroft	1.1
41	jmc	1.14	theMin = 0.
42			theMax = 0.
43			theMean= 0.
44			theSD = 0.
45			theVar = 0.
46			numPnts= 0
47			noPnts = .TRUE.
48	adcroft	1.1
49			DO bj=myByLo(myThid),myByHi(myThid)
50			DO bi=myBxLo(myThid),myBxHi(myThid)
51	jmc	1.12	tileMean(bi,bj) = 0.
52			tileVar (bi,bj) = 0.
53	adcroft	1.1	DO K=1,myNr
54			DO J=1,sNy
55			DO I=1,sNx
56			tmpVal=arr(I,J,K,bi,bj)
57	jmc	1.14	c IF (tmpVal.NE.0. .AND. noPnts) THEN
58			IF ( noPnts ) THEN
59			theMin = tmpVal
60			theMax = tmpVal
61			noPnts = .FALSE.
62	adcroft	1.1	ENDIF
63	jmc	1.14	c IF (tmpVal.NE.0.) THEN
64			theMin = MIN(theMin,tmpVal)
65			theMax = MAX(theMax,tmpVal)
66			tileMean(bi,bj)=tileMean(bi,bj)+tmpVal
67			tileVar (bi,bj)=tileVar (bi,bj)+tmpVal*tmpVal
68	adcroft	1.1	numPnts=numPnts+1
69	jmc	1.14	c ENDIF
70	adcroft	1.1	ENDDO
71			ENDDO
72			ENDDO
73	jmc	1.14	c theMean=theMean+tileMean(bi,bj)
74			c theVar =theVar +tileVar (bi,bj)
75	adcroft	1.1	ENDDO
76			ENDDO
77
78	jmc	1.13	c _GLOBAL_SUM_RL(theMean,myThid)
79			c _GLOBAL_SUM_RL(theVar,myThid)
80	jmc	1.12	CALL GLOBAL_SUM_TILE_RL( tileMean, theMean, myThid )
81			CALL GLOBAL_SUM_TILE_RL( tileVar , theVar , myThid )
82	adcroft	1.1	tmpVal=FLOAT(numPnts)
83	jmc	1.13	_GLOBAL_SUM_RL(tmpVal,myThid)
84	jmc	1.9	numPnts=NINT(tmpVal)
85	adcroft	1.1
86			IF (tmpVal.GT.0.) THEN
87	jmc	1.11	rNumPnts=1. _d 0/tmpVal
88	jmc	1.14	theMean=theMean*rNumPnts
89			theVar=theVar*rNumPnts
90	jmc	1.9	IF ( noPnts ) theMin = theMean
91			theMin=-theMin
92	jmc	1.13	_GLOBAL_MAX_RL(theMin,myThid)
93	jmc	1.9	theMin=-theMin
94			IF ( noPnts ) theMax = theMean
95	jmc	1.13	_GLOBAL_MAX_RL(theMax,myThid)
96	adcroft	1.1
97			DO bj=myByLo(myThid),myByHi(myThid)
98			DO bi=myBxLo(myThid),myBxHi(myThid)
99	jmc	1.12	tileSD(bi,bj)=0.
100	adcroft	1.1	DO K=1,myNr
101			DO J=1,sNy
102			DO I=1,sNx
103			tmpVal=arr(I,J,K,bi,bj)
104	jmc	1.14	c IF (tmpVal.NE.0.) THEN
105	jmc	1.12	tileSD(bi,bj) = tileSD(bi,bj)
106	jmc	1.14	& + (tmpVal-theMean)*(tmpVal-theMean)
107			c ENDIF
108	adcroft	1.1	ENDDO
109			ENDDO
110			ENDDO
111	jmc	1.12	c theSD = theSD + tileSD(bi,bj)
112	adcroft	1.1	ENDDO
113			ENDDO
114
115	jmc	1.13	c _GLOBAL_SUM_RL(theSD,myThid)
116	jmc	1.12	CALL GLOBAL_SUM_TILE_RL( tileSD, theSD, myThid )
117	adcroft	1.1
118	jmc	1.14	theSD = SQRT(theSD*rNumPnts)
119	jmc	1.12	c theSD = SQRT(theVar-theMean*theMean)
120	adcroft	1.1	ENDIF
121
122			RETURN
123			END