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

Diff of /MITgcm/pkg/monitor/mon_stats_rl.F

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

revision 1.12 by jmc, Mon Oct 15 00:18:40 2007 UTC revision 1.13 by jmc, Tue Apr 28 18:16:53 2009 UTC
# Line 96  c       theVar  = theVar  + tileVar (bi, Line 96  c       theVar  = theVar  + tileVar (bi,
96         ENDDO         ENDDO
97        ENDDO        ENDDO
98    
99  c     _GLOBAL_SUM_R8(theDel2,myThid)  c     _GLOBAL_SUM_RL(theDel2,myThid)
100  c     _GLOBAL_SUM_R8(theVol,myThid)  c     _GLOBAL_SUM_RL(theVol,myThid)
101  c     _GLOBAL_SUM_R8(theMean,myThid)  c     _GLOBAL_SUM_RL(theMean,myThid)
102  c     _GLOBAL_SUM_R8(theVar,myThid)  c     _GLOBAL_SUM_RL(theVar,myThid)
103        CALL GLOBAL_SUM_TILE_RL( tileDel2, theDel2, myThid )        CALL GLOBAL_SUM_TILE_RL( tileDel2, theDel2, myThid )
104        CALL GLOBAL_SUM_TILE_RL( tileVol , theVol , myThid )        CALL GLOBAL_SUM_TILE_RL( tileVol , theVol , myThid )
105        CALL GLOBAL_SUM_TILE_RL( tileMean, theMean, myThid )        CALL GLOBAL_SUM_TILE_RL( tileMean, theMean, myThid )
106        CALL GLOBAL_SUM_TILE_RL( tileVar , theVar , myThid )        CALL GLOBAL_SUM_TILE_RL( tileVar , theVar , myThid )
107        tmpVal=FLOAT(numPnts)        tmpVal=FLOAT(numPnts)
108        _GLOBAL_SUM_R8(tmpVal,myThid)        _GLOBAL_SUM_RL(tmpVal,myThid)
109        numPnts=NINT(tmpVal)        numPnts=NINT(tmpVal)
110    
111        IF (tmpVal.GT.0.) THEN        IF (tmpVal.GT.0.) THEN
# Line 118  c     _GLOBAL_SUM_R8(theVar,myThid) Line 118  c     _GLOBAL_SUM_R8(theVar,myThid)
118         theVar=theVar/theVol         theVar=theVar/theVol
119         IF ( noPnts ) theMin = theMean         IF ( noPnts ) theMin = theMean
120         theMin=-theMin         theMin=-theMin
121         _GLOBAL_MAX_R8(theMin,myThid)         _GLOBAL_MAX_RL(theMin,myThid)
122         theMin=-theMin         theMin=-theMin
123         IF ( noPnts ) theMax = theMean         IF ( noPnts ) theMax = theMean
124         _GLOBAL_MAX_R8(theMax,myThid)         _GLOBAL_MAX_RL(theMax,myThid)
125    
126         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
127          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 143  c        theSD = theSD + tileSD(bi,bj) Line 143  c        theSD = theSD + tileSD(bi,bj)
143          ENDDO          ENDDO
144         ENDDO         ENDDO
145    
146  c      _GLOBAL_SUM_R8(theSD,myThid)  c      _GLOBAL_SUM_RL(theSD,myThid)
147         CALL GLOBAL_SUM_TILE_RL( tileSD, theSD, myThid )         CALL GLOBAL_SUM_TILE_RL( tileSD, theSD, myThid )
148    
149         theSD = SQRT(theSD/theVol)         theSD = SQRT(theSD/theVol)
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22