/[MITgcm]/MITgcm/pkg/diagnostics/diagstats_local.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagstats_local.F

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

revision 1.10 by jmc, Mon Sep 3 20:29:47 2012 UTC revision 1.11 by jmc, Fri Aug 8 19:29:48 2014 UTC
# Line 12  C     !INTERFACE: Line 12  C     !INTERFACE:
12       I                  scaleFact, power, useFract, sizF,       I                  scaleFact, power, useFract, sizF,
13       I                  sizI1,sizI2,sizJ1,sizJ2,sizK,sizTx,sizTy,       I                  sizI1,sizI2,sizJ1,sizJ2,sizK,sizTx,sizTy,
14       I                  iRun,jRun,kIn,biIn,bjIn,       I                  iRun,jRun,kIn,biIn,bjIn,
15       I                  k,bi,bj, region2fill, ndId, parsFld,       I                  k,bi,bj, bibjFlg, region2fill,
16       I                  myThid )       I                  ndId, parsFld, myThid )
17    
18  C     !DESCRIPTION:  C     !DESCRIPTION:
19  C     Update array statFld  C     Update array statFld
# Line 54  C     iRun,jRun   :: range of 1rst & 2nd Line 54  C     iRun,jRun   :: range of 1rst & 2nd
54  C     kIn         :: level index of inpFld array to process  C     kIn         :: level index of inpFld array to process
55  C     biIn,bjIn   :: tile indices of inpFld array to process  C     biIn,bjIn   :: tile indices of inpFld array to process
56  C     k,bi,bj     :: level and tile indices used for weighting (mask,area ...)  C     k,bi,bj     :: level and tile indices used for weighting (mask,area ...)
57    C     bibjFlg     :: passed from calling S/R (see diagstats_fill.F)
58  C     region2fill :: indicates whether to compute statistics over this region  C     region2fill :: indicates whether to compute statistics over this region
59  C     ndId        :: Diagnostics Id Number (in available diag. list)  C     ndId        :: Diagnostics Id Number (in available diag. list)
60  C     parsFld     :: parser field with characteristics of the diagnostics  C     parsFld     :: parser field with characteristics of the diagnostics
# Line 67  C     myThid      :: my Thread Id number Line 68  C     myThid      :: my Thread Id number
68        INTEGER power        INTEGER power
69        LOGICAL useFract        LOGICAL useFract
70        INTEGER iRun, jRun, kIn, biIn, bjIn        INTEGER iRun, jRun, kIn, biIn, bjIn
71        INTEGER k, bi, bj, ndId        INTEGER k, bi, bj, bibjFlg
72        INTEGER region2fill(0:nRegions)        INTEGER region2fill(0:nRegions)
73          INTEGER ndId
74        CHARACTER*16 parsFld        CHARACTER*16 parsFld
75        INTEGER myThid        INTEGER myThid
76  CEOP  CEOP
# Line 222  c        ELSEIF ( parsFld(10:10) .EQ. '1 Line 224  c        ELSEIF ( parsFld(10:10) .EQ. '1
224    
225  C     Update cumulative statistics array  C     Update cumulative statistics array
226           IF ( statLoc(0).GT.0. ) THEN           IF ( statLoc(0).GT.0. ) THEN
227              IF ( statFld(0,n).LE.0. ) THEN             IF ( statFld(0,n).LE.0. ) THEN
228                statFld(im,n) = statLoc(im)               statFld(im,n) = statLoc(im)
229                statFld(ix,n) = statLoc(ix)               statFld(ix,n) = statLoc(ix)
230              ELSE             ELSE
231                statFld(im,n) = MIN( statFld(im,n), statLoc(im) )               statFld(im,n) = MIN( statFld(im,n), statLoc(im) )
232                statFld(ix,n) = MAX( statFld(ix,n), statLoc(ix) )               statFld(ix,n) = MAX( statFld(ix,n), statLoc(ix) )
233              ENDIF             ENDIF
234               IF ( bibjFlg.GE.0 ) THEN
235              DO i=0,iv              DO i=0,iv
236                statFld(i,n) = statFld(i,n) + statLoc(i)               statFld(i,n) = statFld(i,n) + statLoc(i)
237              ENDDO              ENDDO
238               ELSE
239                DO i=1,iv
240                 statFld(i,n) = statFld(i,n) + statLoc(i)
241                ENDDO
242               ENDIF
243           ENDIF           ENDIF
244    
245  C---   processing region "n" ends here.  C---   processing region "n" ends here.
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22