/[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.11 by jmc, Fri Aug 8 19:29:48 2014 UTC revision 1.12 by jmc, Mon Aug 25 21:59:18 2014 UTC
# Line 83  C     !FUNCTIONS: Line 83  C     !FUNCTIONS:
83    
84  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
85  C     i,j    :: loop indices  C     i,j    :: loop indices
86        INTEGER i, n, kFr, kRegMsk        INTEGER i, n, kFr, kRegMsk, lReg
87        INTEGER im, ix, iv        INTEGER im, ix, iv
88        PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats )        PARAMETER ( iv = nStats - 2 , im = nStats - 1 , ix = nStats )
89        LOGICAL exclSpVal        LOGICAL exclSpVal
# Line 110  C---+----1----+----2----+----3----+----4 Line 110  C---+----1----+----2----+----3----+----4
110  C---   Compute statistics for this tile, level and region:  C---   Compute statistics for this tile, level and region:
111    
112           kRegMsk = diagSt_kRegMsk(n)           kRegMsk = diagSt_kRegMsk(n)
113             lReg = 0
114             IF ( n.GE.1 ) THEN
115               lReg = 1
116               IF ( parsFld(2:2).EQ.'U' ) lReg = 2
117               IF ( parsFld(2:2).EQ.'V' ) lReg = 3
118             ENDIF
119    
120           IF ( parsFld(10:10) .EQ. 'R' ) THEN           IF ( parsFld(10:10) .EQ. 'R' ) THEN
121    
# Line 123  C---   Compute statistics for this tile, Line 129  C---   Compute statistics for this tile,
129       O            statLoc,       O            statLoc,
130       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
131       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
132       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
133       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
134       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
135       I            maskW(1-OLx,1-OLy,k,bi,bj),       I            maskInW(1-OLx,1-OLy,bi,bj),
136       I            hFacW(1-OLx,1-OLy,k,bi,bj), rAw(1-OLx,1-OLy,bi,bj),       I            hFacW(1-OLx,1-OLy,k,bi,bj), rAw(1-OLx,1-OLy,bi,bj),
137       I            drLoc, specialVal, exclSpVal, useWeight, myThid )       I            drLoc, specialVal, exclSpVal, useWeight, myThid )
138  c    I            drLoc, k,bi,bj, parsFld, myThid )  c    I            drLoc, k,bi,bj, parsFld, myThid )
# Line 135  c    I            drLoc, k,bi,bj, parsFl Line 141  c    I            drLoc, k,bi,bj, parsFl
141       O            statLoc,       O            statLoc,
142       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
143       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
144       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
145       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
146       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
147       I            maskS(1-OLx,1-OLy,k,bi,bj),       I            maskInS(1-OLx,1-OLy,bi,bj),
148       I            hFacS(1-OLx,1-OLy,k,bi,bj), rAs(1-OLx,1-OLy,bi,bj),       I            hFacS(1-OLx,1-OLy,k,bi,bj), rAs(1-OLx,1-OLy,bi,bj),
149       I            drLoc, specialVal, exclSpVal, useWeight, myThid )       I            drLoc, specialVal, exclSpVal, useWeight, myThid )
150            ELSE            ELSE
# Line 146  c    I            drLoc, k,bi,bj, parsFl Line 152  c    I            drLoc, k,bi,bj, parsFl
152       O            statLoc,       O            statLoc,
153       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
154       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
155       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
156       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
157       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
158       I            maskC(1-OLx,1-OLy,k,bi,bj),       I            maskInC(1-OLx,1-OLy,bi,bj),
159       I            hFacC(1-OLx,1-OLy,k,bi,bj), rA(1-OLx,1-OLy,bi,bj),       I            hFacC(1-OLx,1-OLy,k,bi,bj), rA(1-OLx,1-OLy,bi,bj),
160       I            drLoc, specialVal, exclSpVal, useWeight, myThid )       I            drLoc, specialVal, exclSpVal, useWeight, myThid )
161            ENDIF            ENDIF
# Line 161  c    I            drLoc, k,bi,bj, parsFl Line 167  c    I            drLoc, k,bi,bj, parsFl
167       O            statLoc,       O            statLoc,
168       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
169       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
170       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
171       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
172       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
173       I            maskInC(1-OLx,1-OLy,bi,bj), rA(1-OLx,1-OLy,bi,bj),       I            maskInC(1-OLx,1-OLy,bi,bj), rA(1-OLx,1-OLy,bi,bj),
# Line 174  c    I            drLoc, k,bi,bj, parsFl Line 180  c    I            drLoc, k,bi,bj, parsFl
180       O            statLoc,       O            statLoc,
181       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
182       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
183       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
184       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
185       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
186       I            rA(1-OLx,1-OLy,bi,bj),       I            rA(1-OLx,1-OLy,bi,bj),
# Line 190  c        ELSEIF ( parsFld(10:10) .EQ. '1 Line 196  c        ELSEIF ( parsFld(10:10) .EQ. '1
196       O            statLoc,       O            statLoc,
197       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
198       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
199       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
200       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
201       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
202       I            maskInW(1-OLx,1-OLy,bi,bj),       I            maskInW(1-OLx,1-OLy,bi,bj),
# Line 201  c        ELSEIF ( parsFld(10:10) .EQ. '1 Line 207  c        ELSEIF ( parsFld(10:10) .EQ. '1
207       O            statLoc,       O            statLoc,
208       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
209       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
210       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
211       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
212       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
213       I            maskInS(1-OLx,1-OLy,bi,bj),       I            maskInS(1-OLx,1-OLy,bi,bj),
# Line 212  c        ELSEIF ( parsFld(10:10) .EQ. '1 Line 218  c        ELSEIF ( parsFld(10:10) .EQ. '1
218       O            statLoc,       O            statLoc,
219       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),       I            inpFld(sizI1,sizJ1,kIn,biIn,bjIn),
220       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),       I            frcFld(sizI1,sizJ1,kFr,biIn,bjIn),
221       I            scaleFact, power, useFract, n, diagSt_vRegMsk(n),       I            scaleFact, power, useFract, lReg, diagSt_vRegMsk(n),
222       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
223       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),       I            diagSt_regMask(1-OLx,1-OLy,kRegMsk,bi,bj),
224       I            maskInC(1-OLx,1-OLy,bi,bj),       I            maskInC(1-OLx,1-OLy,bi,bj),
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22