pkg/monitor/mon_stats_rs.F

C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rs.F,v 1.3 2003/05/13 18:18:05 adcroft Exp $
C $Name:  $

#include "MONITOR_OPTIONS.h"

      SUBROUTINE MON_STATS_RS(
     I                myNr, arr,
     O                theMin,theMax,theMean,theSD,
     I                myThid )
C     /==========================================================\
C     | SUBROUTINE MON_STATS_RS                                  |
C     | o Calculate bare statistics of global array "_RS arr"    |
C     |==========================================================|
C     \==========================================================/
      IMPLICIT NONE

C     === Global data ===
#include "SIZE.h"
#include "EEPARAMS.h"

C     === Routine arguments ===
      INTEGER myNr
      _RS arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
      _RL theMin
      _RL theMax
      _RL theMean,theMeanTile
      _RL theSD,theSDTile
      INTEGER myThid

C     === Local variables ====
      INTEGER bi,bj,I,J,K
      INTEGER numPnts
      LOGICAL noPnts
      _RL tmpVal,rNumPnts
      _RL theVar,theVarTile

      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)
        theVarTile=0.
        theMeanTile=0.
        DO K=1,myNr
         DO J=1,sNy
          DO I=1,sNx
           tmpVal=arr(I,J,K,bi,bj)
           IF (tmpVal.NE.0. .AND. noPnts) THEN
            theMin=tmpVal
            theMax=tmpVal
            noPnts=.FALSE.
           ENDIF
           IF (tmpVal.NE.0.) THEN
            theMin=min(theMin,tmpVal)
            theMax=max(theMax,tmpVal)
            theMeanTile=theMeanTile+tmpVal
            theVarTile=theVarTile+tmpVal**2
            numPnts=numPnts+1
           ENDIF
          ENDDO
         ENDDO
        ENDDO
        theMean=theMean+theMeanTile
        theVar=theVar+theVarTile
       ENDDO
      ENDDO

      theMin=-theMin
      _GLOBAL_MAX_R8(theMin,myThid)
      theMin=-theMin
      _GLOBAL_MAX_R8(theMax,myThid)
      _GLOBAL_SUM_R8(theMean,myThid)
      _GLOBAL_SUM_R8(theVar,myThid)
      tmpVal=FLOAT(numPnts)
      _GLOBAL_SUM_R8(tmpVal,myThid)
      numPnts=INT(tmpVal+0.5)

      IF (tmpVal.GT.0.) THEN
       rNumPnts=1./tmpVal
       theMean=theMean*rNumPnts
       theVar=theVar*rNumPnts

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

       _GLOBAL_SUM_R8(theSD,myThid)

       theSD=sqrt(theSD*rNumPnts)
c      theSD=sqrt(theVar-theMean**2)
      ENDIF

      RETURN
      END
1	adcroft	1.4	C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rs.F,v 1.3 2003/05/13 18:18:05 adcroft Exp $
2	adcroft	1.2	C $Name: $
3	adcroft	1.1
4	adcroft	1.3	#include "MONITOR_OPTIONS.h"
5	adcroft	1.1
6			SUBROUTINE MON_STATS_RS(
7			I myNr, arr,
8			O theMin,theMax,theMean,theSD,
9			I myThid )
10			C /==========================================================\
11			C \| SUBROUTINE MON_STATS_RS \|
12			C \| o Calculate bare statistics of global array "_RS arr" \|
13			C \|==========================================================\|
14			C \==========================================================/
15			IMPLICIT NONE
16
17			C === Global data ===
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20
21			C === Routine arguments ===
22			INTEGER myNr
23			_RS arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
24			_RL theMin
25			_RL theMax
26	adcroft	1.4	_RL theMean,theMeanTile
27			_RL theSD,theSDTile
28	adcroft	1.1	INTEGER myThid
29
30			C === Local variables ====
31			INTEGER bi,bj,I,J,K
32			INTEGER numPnts
33			LOGICAL noPnts
34			_RL tmpVal,rNumPnts
35	adcroft	1.4	_RL theVar,theVarTile
36	adcroft	1.1
37			theMin=0.
38			theMax=0.
39			theMean=0.
40			theSD=0.
41			theVar=0.
42			numPnts=0
43			noPnts=.TRUE.
44
45			DO bj=myByLo(myThid),myByHi(myThid)
46			DO bi=myBxLo(myThid),myBxHi(myThid)
47	adcroft	1.4	theVarTile=0.
48			theMeanTile=0.
49	adcroft	1.1	DO K=1,myNr
50			DO J=1,sNy
51			DO I=1,sNx
52			tmpVal=arr(I,J,K,bi,bj)
53			IF (tmpVal.NE.0. .AND. noPnts) THEN
54			theMin=tmpVal
55			theMax=tmpVal
56			noPnts=.FALSE.
57			ENDIF
58			IF (tmpVal.NE.0.) THEN
59			theMin=min(theMin,tmpVal)
60			theMax=max(theMax,tmpVal)
61	adcroft	1.4	theMeanTile=theMeanTile+tmpVal
62			theVarTile=theVarTile+tmpVal**2
63	adcroft	1.1	numPnts=numPnts+1
64			ENDIF
65			ENDDO
66			ENDDO
67			ENDDO
68	adcroft	1.4	theMean=theMean+theMeanTile
69			theVar=theVar+theVarTile
70	adcroft	1.1	ENDDO
71			ENDDO
72
73			theMin=-theMin
74			_GLOBAL_MAX_R8(theMin,myThid)
75			theMin=-theMin
76			_GLOBAL_MAX_R8(theMax,myThid)
77			_GLOBAL_SUM_R8(theMean,myThid)
78			_GLOBAL_SUM_R8(theVar,myThid)
79			tmpVal=FLOAT(numPnts)
80			_GLOBAL_SUM_R8(tmpVal,myThid)
81			numPnts=INT(tmpVal+0.5)
82
83			IF (tmpVal.GT.0.) THEN
84			rNumPnts=1./tmpVal
85			theMean=theMean*rNumPnts
86			theVar=theVar*rNumPnts
87
88			DO bj=myByLo(myThid),myByHi(myThid)
89			DO bi=myBxLo(myThid),myBxHi(myThid)
90	adcroft	1.4	theSDtile=0.
91	adcroft	1.1	DO K=1,myNr
92			DO J=1,sNy
93			DO I=1,sNx
94			tmpVal=arr(I,J,K,bi,bj)
95			IF (tmpVal.NE.0.) THEN
96	adcroft	1.4	theSDtile=theSDtile+(tmpVal-theMean)**2
97	adcroft	1.1	ENDIF
98			ENDDO
99			ENDDO
100			ENDDO
101	adcroft	1.4	theSD=theSD+theSDtile
102	adcroft	1.1	ENDDO
103			ENDDO
104
105			_GLOBAL_SUM_R8(theSD,myThid)
106
107	adcroft	1.2	theSD=sqrt(theSD*rNumPnts)
108			c theSD=sqrt(theVar-theMean**2)
109	adcroft	1.1	ENDIF
110
111			RETURN
112			END