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

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

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

revision 1.17 by jmc, Tue Apr 28 18:16:53 2009 UTC revision 1.18 by jmc, Sat Nov 28 20:52:54 2009 UTC
# Line 29  C     !INPUT PARAMETERS: Line 29  C     !INPUT PARAMETERS:
29  CEOP  CEOP
30    
31  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
32        INTEGER bi,bj,i,j,k        INTEGER bi,bj,i,j,k,kp1
33        _RL numPnts,theVol,tmpVal        _RL numPnts,theVol,tmpVal,mskp1
34        _RL theMax,theMean,theVolMean,potEnMean        _RL theMax,theMean,theVolMean,potEnMean
35        _RL tileMean(nSx,nSy)        _RL tileMean(nSx,nSy)
36        _RL tileVlAv(nSx,nSy)        _RL tileVlAv(nSx,nSy)
37        _RL tilePEav(nSx,nSy)        _RL tilePEav(nSx,nSy)
38        _RL tileVol (nSx,nSy)        _RL tileVol (nSx,nSy)
39    #ifdef ALLOW_NONHYDROSTATIC
40          _RL tmpWke
41    #endif
42    
43        numPnts=0.        numPnts=0.
44        theVol=0.        theVol=0.
# Line 51  C     !LOCAL VARIABLES: Line 54  C     !LOCAL VARIABLES:
54          tileVlAv(bi,bj) = 0. _d 0          tileVlAv(bi,bj) = 0. _d 0
55          tilePEav(bi,bj) = 0. _d 0          tilePEav(bi,bj) = 0. _d 0
56          DO k=1,Nr          DO k=1,Nr
57             kp1 = MIN(k+1,Nr)
58             mskp1 = 1.
59             IF ( k.GE.Nr) mskp1 = 0.
60           DO j=1,sNy           DO j=1,sNy
61            DO i=1,sNx            DO i=1,sNx
62             tileVol(bi,bj) = tileVol(bi,bj)             tileVol(bi,bj) = tileVol(bi,bj)
# Line 82  C    with no assumption on how grid spac Line 88  C    with no assumption on how grid spac
88       &                     + tmpVal*deepFac2C(k)*rhoFacC(k)*drF(k)       &                     + tmpVal*deepFac2C(k)*rhoFacC(k)*drF(k)
89             tmpVal= tmpVal*_recip_hFacC(i,j,k,bi,bj)*recip_rA(i,j,bi,bj)             tmpVal= tmpVal*_recip_hFacC(i,j,k,bi,bj)*recip_rA(i,j,bi,bj)
90    
91    #ifdef ALLOW_NONHYDROSTATIC
92               IF ( nonHydrostatic ) THEN
93                tmpWke = 0.25*
94         &        ( wVel(i,j, k, bi,bj)*wVel(i,j, k, bi,bj)
95         &                             *deepFac2F( k )*rhoFacF( k )
96         &         +wVel(i,j,kp1,bi,bj)*wVel(i,j,kp1,bi,bj)*mskp1
97         &                             *deepFac2F(kp1)*rhoFacF(kp1)
98         &        )*maskC(i,j,k,bi,bj)
99                tileVlAv(bi,bj) = tileVlAv(bi,bj)
100         &             + tmpWke*rA(i,j,bi,bj)*drF(k)*_hFacC(i,j,k,bi,bj)
101                tmpVal = tmpVal
102         &             + tmpWke*recip_deepFac2C(k)*recip_rhoFacC(k)
103               ENDIF
104    #endif
105    
106             theMax=MAX(theMax,tmpVal)             theMax=MAX(theMax,tmpVal)
107             IF (tmpVal.NE.0.) THEN             IF (tmpVal.NE.0.) THEN
108              tileMean(bi,bj)=tileMean(bi,bj)+tmpVal              tileMean(bi,bj)=tileMean(bi,bj)+tmpVal
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22