/[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.18 by jmc, Sat Nov 28 20:52:54 2009 UTC revision 1.19 by jmc, Mon Nov 30 03:58:22 2009 UTC
# Line 30  CEOP Line 30  CEOP
30    
31  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
32        INTEGER bi,bj,i,j,k,kp1        INTEGER bi,bj,i,j,k,kp1
33        _RL numPnts,theVol,tmpVal,mskp1        _RL numPnts,theVol,tmpVal, mskp1, msk_1
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)
# Line 56  C     !LOCAL VARIABLES: Line 56  C     !LOCAL VARIABLES:
56          DO k=1,Nr          DO k=1,Nr
57           kp1 = MIN(k+1,Nr)           kp1 = MIN(k+1,Nr)
58           mskp1 = 1.           mskp1 = 1.
59           IF ( k.GE.Nr) mskp1 = 0.           IF ( k.GE.Nr ) mskp1 = 0.
60    C- Note: Present NH implementation does not account for D.w/dt at k=1.
61    C        Consequently, wVel(k=1) does not contribute to NH KE (msk_1=0).
62             msk_1 = 1.
63             IF ( k.EQ. 1 ) msk_1 = 0.
64           DO j=1,sNy           DO j=1,sNy
65            DO i=1,sNx            DO i=1,sNx
66             tileVol(bi,bj) = tileVol(bi,bj)             tileVol(bi,bj) = tileVol(bi,bj)
# Line 91  C    with no assumption on how grid spac Line 95  C    with no assumption on how grid spac
95  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
96             IF ( nonHydrostatic ) THEN             IF ( nonHydrostatic ) THEN
97              tmpWke = 0.25*              tmpWke = 0.25*
98       &        ( wVel(i,j, k, bi,bj)*wVel(i,j, k, bi,bj)       &        ( wVel(i,j, k, bi,bj)*wVel(i,j, k, bi,bj)*msk_1
99       &                             *deepFac2F( k )*rhoFacF( k )       &                             *deepFac2F( k )*rhoFacF( k )
100       &         +wVel(i,j,kp1,bi,bj)*wVel(i,j,kp1,bi,bj)*mskp1       &         +wVel(i,j,kp1,bi,bj)*wVel(i,j,kp1,bi,bj)*mskp1
101       &                             *deepFac2F(kp1)*rhoFacF(kp1)       &                             *deepFac2F(kp1)*rhoFacF(kp1)
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22