/[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.19 by jmc, Mon Nov 30 03:58:22 2009 UTC revision 1.20 by jmc, Mon Dec 21 00:06:03 2009 UTC
# Line 126  C       not sure for atmos/ocean in P ; Line 126  C       not sure for atmos/ocean in P ;
126             tmpVal = tmpVal             tmpVal = tmpVal
127       &            + phi0surf(i,j,bi,bj)*etaN(i,j,bi,bj)       &            + phi0surf(i,j,bi,bj)*etaN(i,j,bi,bj)
128             tilePEav(bi,bj) = tilePEav(bi,bj)             tilePEav(bi,bj) = tilePEav(bi,bj)
129       &            + tmpVal*rA(i,j,bi,bj)*deepFac2F(1)*maskH(i,j,bi,bj)       &            + tmpVal*rA(i,j,bi,bj)*deepFac2F(1)
130         &                    *maskInC(i,j,bi,bj)
131  c          tmpVal = etaN(i,j,bi,bj)  c          tmpVal = etaN(i,j,bi,bj)
132  c    &            + phi0surf(i,j,bi,bj)*recip_Bo(i,j,bi,bj)  c    &            + phi0surf(i,j,bi,bj)*recip_Bo(i,j,bi,bj)
133  c          tilePEav(bi,bj) = tilePEav(bi,bj)  c          tilePEav(bi,bj) = tilePEav(bi,bj)
134  c    &        + 0.5 _d 0*Bo_surf(i,j,bi,bj)*tmpVal*tmpVal  c    &        + 0.5 _d 0*Bo_surf(i,j,bi,bj)*tmpVal*tmpVal
135  c    &                  *rA(i,j,bi,bj)*maskH(i,j,bi,bj)  c    &                  *rA(i,j,bi,bj)*maskInC(i,j,bi,bj)
136            ENDDO            ENDDO
137           ENDDO           ENDDO
 c        theMean    = theMean    + tileMean(bi,bj)  
 c        theVol     = theVol     + tileVol(bi,bj)  
 c        theVolMean = theVolMean + tileVlAv(bi,bj)  
 c        potEnMean  = potEnMean  + tilePEav(bi,bj)  
138  C- end bi,bj loops  C- end bi,bj loops
139         ENDDO         ENDDO
140        ENDDO        ENDDO
141        _GLOBAL_SUM_RL(numPnts,myThid)        _GLOBAL_SUM_RL(numPnts,myThid)
142        _GLOBAL_MAX_RL(theMax,myThid)        _GLOBAL_MAX_RL(theMax,myThid)
 c     _GLOBAL_SUM_RL(theMean,myThid)  
 c     _GLOBAL_SUM_RL(theVol,myThid)  
 c     _GLOBAL_SUM_RL(theVolMean,myThid)  
 c     _GLOBAL_SUM_RL(potEnMean, myThid)  
143        CALL GLOBAL_SUM_TILE_RL( tileMean, theMean   , myThid )        CALL GLOBAL_SUM_TILE_RL( tileMean, theMean   , myThid )
144        CALL GLOBAL_SUM_TILE_RL( tileVol , theVol    , myThid )        CALL GLOBAL_SUM_TILE_RL( tileVol , theVol    , myThid )
145        CALL GLOBAL_SUM_TILE_RL( tileVlAv, theVolMean, myThid )        CALL GLOBAL_SUM_TILE_RL( tileVlAv, theVolMean, myThid )

Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20

  ViewVC Help
Powered by ViewVC 1.1.22