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 |
|
|
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 |
|
C since dpphys is RL but argument is RS. leave it like this for now |
|
|
C and will change it once the Regions are fully implemented. |
|
|
|
|
101 |
CALL DIAGSTATS_CALC( |
CALL DIAGSTATS_CALC( |
102 |
O statArr, |
O statArr, |
103 |
I inpArr, frcArr, scaleFact, power, useFract, |
I inpArr, frcArr, scaleFact, power, useFract, |
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 ) |
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 ) |
123 |
|
#endif /* REAL4_IS_SLOW */ |
124 |
|
|
125 |
c ENDIF |
c ENDIF |
126 |
#endif /* ALLOW_FIZHI */ |
#endif /* ALLOW_FIZHI */ |