pkg/monitor/mon_stats_rl.F

C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rl.F,v 1.11 2005/01/27 16:36:24 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, arrMask,arrhFac, arrArea, arrDr,
     O     theMin,theMax,theMean,theSD,theDel2,theVol,
     I     myThid )

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)
      _RS arrMask(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
      _RS arrhFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
      _RS arrArea(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RS arrDr(myNr)
      _RL theMin, theMax, theMean, theSD, theDel2, theVol
      INTEGER myThid
CEOP

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

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

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        tileDel2(bi,bj) = 0.
        tileVol (bi,bj) = 0.
        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)
           IF (arrMask(I,J,K,bi,bj).NE.0. .AND. noPnts) THEN
            theMin=tmpVal
            theMax=tmpVal
            noPnts=.FALSE.
           ENDIF
           IF (arrMask(I,J,K,bi,bj).NE.0.) THEN
            theMin=min(theMin,tmpVal)
            theMax=max(theMax,tmpVal)
            tileDel2(bi,bj) = tileDel2(bi,bj)
     &      + 0.25*ABS(
     &         (arr(I+1,J,K,bi,bj)-tmpVal)*arrMask(I+1,J,K,bi,bj)
     &        +(arr(I-1,J,K,bi,bj)-tmpVal)*arrMask(I-1,J,K,bi,bj)
     &        +(arr(I,J+1,K,bi,bj)-tmpVal)*arrMask(I,J+1,K,bi,bj)
     &        +(arr(I,J-1,K,bi,bj)-tmpVal)*arrMask(I,J-1,K,bi,bj)
     &                )
            numPnts=numPnts+1
            tmpVol = arrArea(I,J,bi,bj)*arrhFac(I,J,K,bi,bj)*arrDr(K)
     &                                 *arrMask(I,J,K,bi,bj)
            tileVol (bi,bj) = tileVol (bi,bj) + tmpVol
            tileMean(bi,bj) = tileMean(bi,bj) + tmpVol*tmpVal
            tileVar (bi,bj) = tileVar (bi,bj) + tmpVol*tmpVal*tmpVal
           ENDIF
          ENDDO
         ENDDO
        ENDDO
c       theDel2 = theDel2 + tileDel2(bi,bj)
c       theVol  = theVol  + tileVol(bi,bj)
c       theMean = theMean + tileMean(bi,bj)
c       theVar  = theVar  + tileVar (bi,bj)
       ENDDO
      ENDDO

c     _GLOBAL_SUM_R8(theDel2,myThid)
c     _GLOBAL_SUM_R8(theVol,myThid)
c     _GLOBAL_SUM_R8(theMean,myThid)
c     _GLOBAL_SUM_R8(theVar,myThid)
      CALL GLOBAL_SUM_TILE_RL( tileDel2, theDel2, myThid )
      CALL GLOBAL_SUM_TILE_RL( tileVol , theVol , myThid )
      CALL GLOBAL_SUM_TILE_RL( tileMean, theMean, myThid )
      CALL GLOBAL_SUM_TILE_RL( tileVar , theVar , myThid )
      tmpVal=FLOAT(numPnts)
      _GLOBAL_SUM_R8(tmpVal,myThid)
      numPnts=NINT(tmpVal)

      IF (tmpVal.GT.0.) THEN
       rNumPnts=1. _d 0/tmpVal
       theDel2=theDel2*rNumPnts
      ENDIF

      IF (theVol.GT.0.) THEN
       theMean=theMean/theVol
       theVar=theVar/theVol
       IF ( noPnts ) theMin = theMean
       theMin=-theMin
       _GLOBAL_MAX_R8(theMin,myThid)
       theMin=-theMin
       IF ( noPnts ) theMax = theMean
       _GLOBAL_MAX_R8(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)
            IF (arrMask(I,J,K,bi,bj).NE.0.) THEN
             tmpVol=arrArea(I,J,bi,bj)*arrhFac(I,J,K,bi,bj)*arrDr(K)
     &                                *arrMask(I,J,K,bi,bj)
             tileSD(bi,bj) = tileSD(bi,bj)
     &                     + tmpVol*(tmpVal-theMean)*(tmpVal-theMean)
            ENDIF
           ENDDO
          ENDDO
         ENDDO
c        theSD = theSD + tileSD(bi,bj)
        ENDDO
       ENDDO

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

       theSD = SQRT(theSD/theVol)
c      theSD = SQRT(theVar-theMean*theMean)
      ENDIF

      RETURN
      END
1	jmc	1.12	C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rl.F,v 1.11 2005/01/27 16:36:24 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	edhill	1.10	I myNr, arr, arrMask,arrhFac, arrArea, arrDr,
13			O theMin,theMax,theMean,theSD,theDel2,theVol,
14			I myThid )
15
16			C Calculate bare statistics of global array ``\_RL arr''.
17
18			C !USES:
19	adcroft	1.1	IMPLICIT NONE
20			#include "SIZE.h"
21			#include "EEPARAMS.h"
22
23	edhill	1.10	C !INPUT PARAMETERS:
24	adcroft	1.1	INTEGER myNr
25			_RL arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
26	jmc	1.3	_RS arrMask(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
27			_RS arrhFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
28			_RS arrArea(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
29			_RS arrDr(myNr)
30	edhill	1.10	_RL theMin, theMax, theMean, theSD, theDel2, theVol
31	adcroft	1.1	INTEGER myThid
32	edhill	1.10	CEOP
33	adcroft	1.1
34	edhill	1.10	C !LOCAL VARIABLES:
35	adcroft	1.1	INTEGER bi,bj,I,J,K
36			INTEGER numPnts
37			LOGICAL noPnts
38			_RL tmpVal,rNumPnts
39	jmc	1.12	_RL theVar
40	jmc	1.4	_RL tmpVol
41	jmc	1.12	_RL tileMean(nSx,nSy)
42			_RL tileVar (nSx,nSy)
43			_RL tileSD (nSx,nSy)
44			_RL tileDel2(nSx,nSy)
45			_RL tileVol (nSx,nSy)
46	adcroft	1.1
47			theMin=0.
48			theMax=0.
49			theMean=0.
50			theSD=0.
51			theVar=0.
52	jmc	1.3	theDel2=0.
53			theVol=0.
54	adcroft	1.1	numPnts=0
55			noPnts=.TRUE.
56
57			DO bj=myByLo(myThid),myByHi(myThid)
58			DO bi=myBxLo(myThid),myBxHi(myThid)
59	jmc	1.12	tileDel2(bi,bj) = 0.
60			tileVol (bi,bj) = 0.
61			tileMean(bi,bj) = 0.
62			tileVar (bi,bj) = 0.
63	adcroft	1.1	DO K=1,myNr
64			DO J=1,sNy
65			DO I=1,sNx
66			tmpVal=arr(I,J,K,bi,bj)
67	jmc	1.3	IF (arrMask(I,J,K,bi,bj).NE.0. .AND. noPnts) THEN
68	adcroft	1.1	theMin=tmpVal
69			theMax=tmpVal
70			noPnts=.FALSE.
71			ENDIF
72	jmc	1.3	IF (arrMask(I,J,K,bi,bj).NE.0.) THEN
73	adcroft	1.1	theMin=min(theMin,tmpVal)
74			theMax=max(theMax,tmpVal)
75	jmc	1.12	tileDel2(bi,bj) = tileDel2(bi,bj)
76			& + 0.25*ABS(
77	jmc	1.3	& (arr(I+1,J,K,bi,bj)-tmpVal)*arrMask(I+1,J,K,bi,bj)
78			& +(arr(I-1,J,K,bi,bj)-tmpVal)*arrMask(I-1,J,K,bi,bj)
79			& +(arr(I,J+1,K,bi,bj)-tmpVal)*arrMask(I,J+1,K,bi,bj)
80			& +(arr(I,J-1,K,bi,bj)-tmpVal)*arrMask(I,J-1,K,bi,bj)
81	jmc	1.12	& )
82	adcroft	1.1	numPnts=numPnts+1
83	jmc	1.3	tmpVol = arrArea(I,J,bi,bj)arrhFac(I,J,K,bi,bj)arrDr(K)
84	jmc	1.6	& *arrMask(I,J,K,bi,bj)
85	jmc	1.12	tileVol (bi,bj) = tileVol (bi,bj) + tmpVol
86			tileMean(bi,bj) = tileMean(bi,bj) + tmpVol*tmpVal
87			tileVar (bi,bj) = tileVar (bi,bj) + tmpVoltmpValtmpVal
88	adcroft	1.1	ENDIF
89			ENDDO
90			ENDDO
91			ENDDO
92	jmc	1.12	c theDel2 = theDel2 + tileDel2(bi,bj)
93			c theVol = theVol + tileVol(bi,bj)
94			c theMean = theMean + tileMean(bi,bj)
95			c theVar = theVar + tileVar (bi,bj)
96	adcroft	1.1	ENDDO
97			ENDDO
98
99	jmc	1.12	c _GLOBAL_SUM_R8(theDel2,myThid)
100			c _GLOBAL_SUM_R8(theVol,myThid)
101			c _GLOBAL_SUM_R8(theMean,myThid)
102			c _GLOBAL_SUM_R8(theVar,myThid)
103			CALL GLOBAL_SUM_TILE_RL( tileDel2, theDel2, myThid )
104			CALL GLOBAL_SUM_TILE_RL( tileVol , theVol , myThid )
105			CALL GLOBAL_SUM_TILE_RL( tileMean, theMean, myThid )
106			CALL GLOBAL_SUM_TILE_RL( tileVar , theVar , myThid )
107	adcroft	1.1	tmpVal=FLOAT(numPnts)
108			_GLOBAL_SUM_R8(tmpVal,myThid)
109	jmc	1.9	numPnts=NINT(tmpVal)
110	adcroft	1.1
111			IF (tmpVal.GT.0.) THEN
112	jmc	1.11	rNumPnts=1. _d 0/tmpVal
113	jmc	1.3	theDel2=theDel2*rNumPnts
114			ENDIF
115
116			IF (theVol.GT.0.) THEN
117			theMean=theMean/theVol
118			theVar=theVar/theVol
119	jmc	1.9	IF ( noPnts ) theMin = theMean
120			theMin=-theMin
121			_GLOBAL_MAX_R8(theMin,myThid)
122			theMin=-theMin
123			IF ( noPnts ) theMax = theMean
124			_GLOBAL_MAX_R8(theMax,myThid)
125	adcroft	1.1
126			DO bj=myByLo(myThid),myByHi(myThid)
127			DO bi=myBxLo(myThid),myBxHi(myThid)
128	jmc	1.12	tileSD(bi,bj)=0.
129	adcroft	1.1	DO K=1,myNr
130			DO J=1,sNy
131			DO I=1,sNx
132			tmpVal=arr(I,J,K,bi,bj)
133	jmc	1.3	IF (arrMask(I,J,K,bi,bj).NE.0.) THEN
134			tmpVol=arrArea(I,J,bi,bj)arrhFac(I,J,K,bi,bj)arrDr(K)
135	jmc	1.6	& *arrMask(I,J,K,bi,bj)
136	jmc	1.12	tileSD(bi,bj) = tileSD(bi,bj)
137			& + tmpVol(tmpVal-theMean)(tmpVal-theMean)
138	adcroft	1.1	ENDIF
139			ENDDO
140			ENDDO
141			ENDDO
142	jmc	1.12	c theSD = theSD + tileSD(bi,bj)
143	adcroft	1.1	ENDDO
144			ENDDO
145
146	jmc	1.12	c _GLOBAL_SUM_R8(theSD,myThid)
147			CALL GLOBAL_SUM_TILE_RL( tileSD, theSD, myThid )
148	adcroft	1.1
149	jmc	1.12	theSD = SQRT(theSD/theVol)
150			c theSD = SQRT(theVar-theMean*theMean)
151	adcroft	1.1	ENDIF
152
153			RETURN
154			END