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 |
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) |
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) |