/[MITgcm]/MITgcm/pkg/ebm/ebm_zonalmean.F
ViewVC logotype

Diff of /MITgcm/pkg/ebm/ebm_zonalmean.F

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

revision 1.1 by heimbach, Fri May 14 21:10:34 2004 UTC revision 1.2 by heimbach, Wed May 19 20:38:53 2004 UTC
# Line 43  C--   Top layer only Line 43  C--   Top layer only
43        k = 1        k = 1
44    
45  c--   Initialise  c--   Initialise
46        DO j=1-OLy,sNy+OLy        DO bj=myByLo(myThid),myByHi(myThid)
47           ZonalMeanSST(j) = 0.0         DO j=1-OLy,sNy+OLy
48           CountX(j) = 0.0           ZonalMeanSST(j,bj) = 0.0
49             CountX(j,bj) = 0.0
50           ENDDO
51        ENDDO        ENDDO
52    
53  C--   Calculate the zonal mean  C--   Calculate the zonal mean
# Line 60  C--   Calculate the zonal mean Line 62  C--   Calculate the zonal mean
62            CountX_tile(j,bi,bj) = CountX_tile(j,bi,bj) +              CountX_tile(j,bi,bj) = CountX_tile(j,bi,bj) +  
63       &         maskC(i,j,k,bi,bj)       &         maskC(i,j,k,bi,bj)
64           ENDDO           ENDDO
          ZonalMeanSST(j) = ZonalMeanSST(j) + ZonalMean_tile(j,bi,bj)  
          CountX(j) = CountX(j) + CountX_tile(j,bi,bj)  
65          ENDDO          ENDDO
66            ZonalMeanSST(j,bj) = ZonalMeanSST(j,bj) +
67         &   ZonalMean_tile(j,bi,bj)
68            CountX(j,bj) = CountX(j,bj) + CountX_tile(j,bi,bj)
69         ENDDO         ENDDO
70        ENDDO        ENDDO
71    
72        DO j=1-OLy,sNy+OLy        DO bj=myByLo(myThid),myByHi(myThid)
73         IF ( CountX(j) .GT. 0.0) THEN         DO j=1-OLy,sNy+OLy
74            ZonalMeanSST(j) = ZonalMeanSST(j)/CountX(j)          _GLOBAL_SUM_R8( CountX(j,bj), myThid )
75         ENDIF          _GLOBAL_SUM_R8( ZonalMeanSST(j,bj), myThid )
76            IF ( CountX(j,bj) .GT. 0.0) THEN
77              ZonalMeanSST(j,bj) = ZonalMeanSST(j,bj)/CountX(j,bj)
78            ENDIF
79           ENDDO
80        ENDDO        ENDDO
81    
82  #endif /* ALLOW_EBM */  #endif /* ALLOW_EBM */
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22