/[MITgcm]/MITgcm/pkg/monitor/mon_stats_rs.F
ViewVC logotype

Annotation of /MITgcm/pkg/monitor/mon_stats_rs.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.9 - (hide annotations) (download)
Tue Apr 28 18:16:53 2009 UTC (15 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint62u, checkpoint62t, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint63g, checkpoint64, checkpoint65, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint65o, checkpoint62b, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint61n, checkpoint61q, checkpoint61o, checkpoint61m, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61z, checkpoint61x, checkpoint61y, HEAD
Changes since 1.8: +7 -7 lines 
change macros (EXCH & GLOBAL_SUM/MAX) sufix _R4/_R8 to _RS/_RL
 when applied to _RS/_RL variable
1 jmc 1.9 C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_stats_rs.F,v 1.8 2007/10/15 00:19:38 jmc Exp $
2 adcroft 1.2 C $Name: $
3 adcroft 1.1
4 adcroft 1.3 #include "MONITOR_OPTIONS.h"
5 adcroft 1.1
6 edhill 1.6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP
8     C !ROUTINE: MON_STATS_RS
9    
10     C !INTERFACE:
11 adcroft 1.1 SUBROUTINE MON_STATS_RS(
12 edhill 1.6 I myNr, arr,
13     O theMin,theMax,theMean,theSD,
14     I myThid )
15    
16     C !DESCRIPTION:
17     C Calculate bare statistics of global array ``\_RS arr''.
18    
19     C !USES:
20 adcroft 1.1 IMPLICIT NONE
21     #include "SIZE.h"
22     #include "EEPARAMS.h"
23    
24 edhill 1.6 C !INPUT PARAMETERS:
25 adcroft 1.1 INTEGER myNr
26     _RS arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNr,nSx,nSy)
27 edhill 1.6 _RL theMin, theMax, theMean, theSD
28 adcroft 1.1 INTEGER myThid
29 edhill 1.6 CEOP
30 adcroft 1.1
31 edhill 1.6 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.8 _RL theVar
37     _RL tileMean(nSx,nSy)
38     _RL tileVar (nSx,nSy)
39     _RL tileSD (nSx,nSy)
40    
41     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.8 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.8 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.8 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.8 c ENDIF
70 adcroft 1.1 ENDDO
71     ENDDO
72     ENDDO
73 jmc 1.8 c theMean=theMean+tileMean(bi,bj)
74     c theVar =theVar +tileVar (bi,bj)
75 adcroft 1.1 ENDDO
76     ENDDO
77    
78 jmc 1.9 c _GLOBAL_SUM_RL(theMean,myThid)
79     c _GLOBAL_SUM_RL(theVar,myThid)
80 jmc 1.8 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.9 _GLOBAL_SUM_RL(tmpVal,myThid)
84 jmc 1.7 numPnts=NINT(tmpVal)
85 adcroft 1.1
86     IF (tmpVal.GT.0.) THEN
87 jmc 1.7 rNumPnts=1. _d 0/tmpVal
88 adcroft 1.1 theMean=theMean*rNumPnts
89     theVar=theVar*rNumPnts
90 jmc 1.5 IF ( noPnts ) theMin = theMean
91     theMin=-theMin
92 jmc 1.9 _GLOBAL_MAX_RL(theMin,myThid)
93 jmc 1.5 theMin=-theMin
94     IF ( noPnts ) theMax = theMean
95 jmc 1.9 _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.8 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.8 c IF (tmpVal.NE.0.) THEN
105     tileSD(bi,bj) = tileSD(bi,bj)
106     & + (tmpVal-theMean)*(tmpVal-theMean)
107     c ENDIF
108 adcroft 1.1 ENDDO
109     ENDDO
110     ENDDO
111 jmc 1.8 c theSD = theSD + tileSD(bi,bj)
112 adcroft 1.1 ENDDO
113     ENDDO
114    
115 jmc 1.9 c _GLOBAL_SUM_RL(theSD,myThid)
116 jmc 1.8 CALL GLOBAL_SUM_TILE_RL( tileSD, theSD, myThid )
117 adcroft 1.1
118 jmc 1.8 theSD = SQRT(theSD*rNumPnts)
119     c theSD = SQRT(theVar-theMean*theMean)
120 adcroft 1.1 ENDIF
121    
122     RETURN
123     END
  ViewVC Help
Powered by ViewVC 1.1.22