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

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

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

revision 1.1 by jmc, Mon Jan 23 22:31:11 2006 UTC revision 1.3 by jmc, Mon Sep 3 20:29:47 2012 UTC
# Line 81  C     !LOCAL VARIABLES: Line 81  C     !LOCAL VARIABLES:
81        LOGICAL useWeight        LOGICAL useWeight
82        INTEGER kl        INTEGER kl
83        _RL drLoc        _RL drLoc
84    #ifndef REAL4_IS_SLOW
85          INTEGER i,j
86          _RS tmp_hFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
87    #endif
88    
89  c     IF ( useFIZHI ) THEN  c     IF ( useFIZHI ) THEN
90    
# Line 93  c     IF ( useFIZHI ) THEN Line 97  c     IF ( useFIZHI ) THEN
97          ENDIF          ENDIF
98          drLoc = 1. _d 0          drLoc = 1. _d 0
99    
100  C- jmc: here we have a Problem if RL & RS are not the same:  #ifdef REAL4_IS_SLOW
101  C    since dpphys is RL but argument is RS. leave it like this for now          CALL DIAGSTATS_CALC(
 C    and will change it once the Regions are fully implemented.  
   
         CALL  DIAGSTATS_CALC(  
102       O            statArr,       O            statArr,
103       I            inpArr, frcArr, scaleFact, power, useFract,       I            inpArr, frcArr, scaleFact, power, useFract,
104       I            regId, regMskVal,       I            regId, regMskVal,
105       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
106       I            regMask, arrMask,       I            regMask, arrMask,
107       I            dpphys(1-Olx,1-Oly,kl,bi,bj), arrArea,       I            dpphys(1-OLx,1-OLy,kl,bi,bj), arrArea,
108         I            drLoc, specialVal, exclSpVal, useWeight, myThid )
109    #else /* REAL4_IS_SLOW */
110    C    make local copy of dpphys (RL type) into RS array tmp_hFac
111            DO j=1-OLy,sNy+OLy
112             DO i=1-OLx,sNx+OLx
113              tmp_hFac(i,j) = dpphys(i,j,kl,bi,bj)
114             ENDDO
115            ENDDO
116            CALL DIAGSTATS_CALC(
117         O            statArr,
118         I            inpArr, frcArr, scaleFact, power, useFract,
119         I            regId, regMskVal,
120         I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
121         I            regMask, arrMask, tmp_hFac, arrArea,
122       I            drLoc, specialVal, exclSpVal, useWeight, myThid )       I            drLoc, specialVal, exclSpVal, useWeight, myThid )
123    #endif /* REAL4_IS_SLOW */
124    
125  c     ENDIF  c     ENDIF
126  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
# Line 200  c     IF ( useLand ) THEN Line 216  c     IF ( useLand ) THEN
216          ENDIF          ENDIF
217          useWeight = .TRUE.          useWeight = .TRUE.
218    
219          CALL  DIAGSTATS_CALC(          CALL DIAGSTATS_CALC(
220       O            statArr,       O            statArr,
221       I            inpArr, frcArr, scaleFact, power, useFract,       I            inpArr, frcArr, scaleFact, power, useFract,
222       I            regId, regMskVal,       I            regId, regMskVal,
223       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,       I            nStats,sizI1,sizI2,sizJ1,sizJ2,iRun,jRun,
224       I            regMask, aim_landFr(1-Olx,1-Oly,bi,bj),       I            regMask, aim_landFr(1-OLx,1-OLy,bi,bj),
225       I            aim_landFr(1-Olx,1-Oly,bi,bj), arrArea,       I            aim_landFr(1-OLx,1-OLy,bi,bj), arrArea,
226       I            drLoc, specialVal, exclSpVal, useWeight, myThid )       I            drLoc, specialVal, exclSpVal, useWeight, myThid )
227    
228  c     ENDIF  c     ENDIF
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22