/[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.7 by adcroft, Mon Jun 25 20:35:23 2001 UTC revision 1.7.2.1 by adcroft, Tue Feb 26 16:04:49 2002 UTC
# Line 17  C     === Global data === Line 17  C     === Global data ===
17  #include "EEPARAMS.h"  #include "EEPARAMS.h"
18  #include "DYNVARS.h"  #include "DYNVARS.h"
19  #include "MONITOR.h"  #include "MONITOR.h"
20    #include "GRID.h"
21    
22  C     === Routine arguments ===  C     === Routine arguments ===
23        INTEGER myThid        INTEGER myThid
24    
25  C     === Local variables ====  C     === Local variables ====
26        INTEGER bi,bj,I,J,K        INTEGER bi,bj,I,J,K
27        _RL tmpVal,theMax,theMean        _RL tmpVal,theMax,theMean,theVolMean,theVol
28        INTEGER numPnts        INTEGER numPnts
29    
30        theMax=0.        theMax=0.
       theMean=0.  
31        numPnts=0        numPnts=0
32          theMean=0.
33          theVolMean=0.
34          theVol=0.
35    
36        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
37         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
38          DO K=1,Nr          DO K=1,Nr
39           DO J=1,sNy           DO J=1,sNy
40            DO I=1,sNx            DO I=1,sNx
41               theVol=theVol+ra(i,j,bi,bj)*drf(k)*hFacC(i,j,k,bi,bj)
42             tmpVal=0.25*( uVel( I , J ,K,bi,bj)*uVel( I , J ,K,bi,bj)             tmpVal=0.25*( uVel( I , J ,K,bi,bj)*uVel( I , J ,K,bi,bj)
43       &                  +uVel(I+1, J ,K,bi,bj)*uVel(I+1, J ,K,bi,bj)       &                  +uVel(I+1, J ,K,bi,bj)*uVel(I+1, J ,K,bi,bj)
44       &                  +vVel( I , J ,K,bi,bj)*vVel( I , J ,K,bi,bj)       &                  +vVel( I , J ,K,bi,bj)*vVel( I , J ,K,bi,bj)
# Line 44  C     === Local variables ==== Line 48  C     === Local variables ====
48              theMean=theMean+tmpVal              theMean=theMean+tmpVal
49              numPnts=numPnts+1              numPnts=numPnts+1
50             ENDIF             ENDIF
51               theVolMean=theVolMean+tmpVal
52         &           *ra(i,j,bi,bj)*drf(k)*hFacC(i,j,k,bi,bj)
53            ENDDO            ENDDO
54           ENDDO           ENDDO
55          ENDDO          ENDDO
# Line 54  C     === Local variables ==== Line 60  C     === Local variables ====
60        tmpVal=float(numPnts)        tmpVal=float(numPnts)
61        _GLOBAL_SUM_R8(tmpVal,myThid)        _GLOBAL_SUM_R8(tmpVal,myThid)
62        IF (tmpVal.NE.0.) theMean=theMean/tmpVal        IF (tmpVal.NE.0.) theMean=theMean/tmpVal
63          _GLOBAL_SUM_R8(theVol,myThid)
64          _GLOBAL_SUM_R8(theVolMean,myThid)
65    
66        CALL MON_OUT_RL( mon_string_none, theMax, mon_foot_max   ,myThid)        CALL MON_OUT_RL(mon_string_none,theMax,mon_foot_max,myThid)
67        CALL MON_OUT_RL( mon_string_none, theMean, mon_foot_mean ,myThid)        CALL MON_OUT_RL(mon_string_none,theMean,mon_foot_mean,myThid)
68          CALL MON_OUT_RL(mon_string_none,theVolMean,
69         &         mon_foot_volint,myThid)
70          IF (theVol.NE.0.) theVolMean=theVolMean/theVol
71          CALL MON_OUT_RL(mon_string_none,theVolMean,
72         &         mon_foot_volmean,myThid)
73    
74        RETURN        RETURN
75        END        END
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.7.2.1
  ViewVC Help
Powered by ViewVC 1.1.22