52 |
#endif |
#endif |
53 |
_RL hFacCtmp |
_RL hFacCtmp |
54 |
_RL hFacMnSz |
_RL hFacMnSz |
55 |
|
_RL tileArea |
56 |
CEOP |
CEOP |
57 |
|
|
58 |
C- Calculate lopping factor hFacC : over-estimate the part inside of the domain |
C- Calculate lopping factor hFacC : over-estimate the part inside of the domain |
161 |
& 'Model Ro_surf (ini_masks_etc)', 1, myThid) |
& 'Model Ro_surf (ini_masks_etc)', 1, myThid) |
162 |
|
|
163 |
C Calculate quantities derived from XY depth map |
C Calculate quantities derived from XY depth map |
164 |
|
globalArea = 0. _d 0 |
165 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
166 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
167 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
177 |
ENDIF |
ENDIF |
178 |
ENDDO |
ENDDO |
179 |
ENDDO |
ENDDO |
180 |
|
C- Compute the domain Area: |
181 |
|
tileArea = 0. _d 0 |
182 |
|
DO j=1,sNy |
183 |
|
DO i=1,sNx |
184 |
|
tileArea = tileArea + rA(i,j,bi,bj)*maskH(i,j,bi,bj) |
185 |
|
ENDDO |
186 |
|
ENDDO |
187 |
|
globalArea = globalArea + tileArea |
188 |
ENDDO |
ENDDO |
189 |
ENDDO |
ENDDO |
190 |
C _EXCH_XY_R4( recip_Rcol, myThid ) |
C _EXCH_XY_R4( recip_Rcol, myThid ) |
191 |
|
_GLOBAL_SUM_R8( globalArea, myThid ) |
192 |
|
|
193 |
C hFacW and hFacS (at U and V points) |
C hFacW and hFacS (at U and V points) |
194 |
DO bj=myByLo(myThid), myByHi(myThid) |
DO bj=myByLo(myThid), myByHi(myThid) |