/[MITgcm]/MITgcm/pkg/ctrl/ctrl_map_gentim2d.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_map_gentim2d.F

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

revision 1.2 by jmc, Sat Feb 2 23:26:57 2013 UTC revision 1.3 by gforget, Wed Mar 27 00:10:26 2013 UTC
# Line 59  C     == Local variables == Line 59  C     == Local variables ==
59    
60        _RL fac        _RL fac
61        _RL xx_gentim2d_loc(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)        _RL xx_gentim2d_loc(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
62          CHARACTER*(MAX_LEN_MBUF) msgBuf
63          _RL LOCsumTile(nSx,nSy), LOCsumGlob
64    
65  c     == external ==  c     == external ==
66        integer  ilnblnk        integer  ilnblnk
# Line 99  C Line 101  C
101       I      zeroRL, zeroRL,       I      zeroRL, zeroRL,
102       I      mytime, myiter, mythid )       I      mytime, myiter, mythid )
103  C  C
104          DO bj=myByLo(myThid),myByHi(myThid)
105           DO bi=myBxLo(myThid),myBxHi(myThid)
106                do j = 1,sNy
107                  do i =  1,sNx
108                  if (xx_gentim2d_cumsum(iarr)) then
109                    xx_gentim2d(i,j,bi,bj,iarr)=xx_gentim2d(i,j,bi,bj,iarr)
110         &            +xx_gentim2d_loc(i,j,bi,bj)
111                  else
112                    xx_gentim2d(i,j,bi,bj,iarr)=xx_gentim2d_loc(i,j,bi,bj)
113                  endif
114                  enddo
115                enddo
116            enddo
117          enddo
118    C
119          if (xx_gentim2d_glosum(iarr)) then
120    
121          LOCsumGlob=0. _d 0
122          DO bj=myByLo(myThid),myByHi(myThid)
123           DO bi=myBxLo(myThid),myBxHi(myThid)
124              LOCsumTile(bi,bj)=0. _d 0
125                do j = 1,sNy
126                  do i =  1,sNx
127                    LOCsumTile(bi,bj)=LOCsumTile(bi,bj)+
128         &            maskC(i,j,1,bi,bj)*rA(i,j,bi,bj)
129         &            *xx_gentim2d(i,j,bi,bj,iarr)
130                  enddo
131                enddo
132            enddo
133          enddo
134    
135          CALL GLOBAL_SUM_TILE_RL( LOCsumTile, LOCsumGlob, myThid )
136    
137           DO bj = myByLo(myThid), myByHi(myThid)
138            DO bi = myBxLo(myThid), myBxHi(myThid)
139             DO J = 1-Oly,sNy+Oly
140              DO I = 1-Olx,sNx+Olx
141                 xx_gentim2d(I,J,bi,bj,iarr) =
142         &            LOCsumGlob/globalArea*maskC(i,j,1,bi,bj)
143              ENDDO
144             ENDDO
145            ENDDO
146           ENDDO
147    
148             WRITE(msgBuf,'(A,I6,A,I6,A,1PE21.14)') ' xx_gentim2d ',
149         &    iarr,' : iter=', myiter, ' ; global sum = ', LOCsumGlob
150             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
151         &                       SQUEEZE_RIGHT, myThid )
152    
153           endif
154    
155        ENDDO        ENDDO
156    
157  #endif /* ALLOW_GENTIM2D_CONTROL */  #endif /* ALLOW_GENTIM2D_CONTROL */
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22