220 |
DO bi=myBxLo(myThid), myBxHi(myThid) |
DO bi=myBxLo(myThid), myBxHi(myThid) |
221 |
DO K=1, Nr |
DO K=1, Nr |
222 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
223 |
|
hFacW(1-OLx,J,K,bi,bj)= 0. |
224 |
DO I=2-Olx,sNx+Olx |
DO I=2-Olx,sNx+Olx |
225 |
hFacW(I,J,K,bi,bj)= |
hFacW(I,J,K,bi,bj)= |
226 |
& MIN(hFacC(I,J,K,bi,bj),hFacC(I-1,J,K,bi,bj)) |
& MIN(hFacC(I,J,K,bi,bj),hFacC(I-1,J,K,bi,bj)) |
227 |
ENDDO |
ENDDO |
228 |
ENDDO |
ENDDO |
229 |
|
DO I=1-Olx,sNx+Olx |
230 |
|
hFacS(I,1-OLy,K,bi,bj)= 0. |
231 |
|
ENDDO |
232 |
DO J=2-Oly,sNy+oly |
DO J=2-Oly,sNy+oly |
233 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
234 |
hFacS(I,J,K,bi,bj)= |
hFacS(I,J,K,bi,bj)= |
352 |
|
|
353 |
c #ifdef ALLOW_NONHYDROSTATIC |
c #ifdef ALLOW_NONHYDROSTATIC |
354 |
C-- Calculate "recip_hFacU" = reciprocal hfac distance/volume for W cells |
C-- Calculate "recip_hFacU" = reciprocal hfac distance/volume for W cells |
355 |
C not used ; computed locally in CALC_GW |
C NOTE: not used ; computed locally in CALC_GW |
356 |
c #endif |
c #endif |
357 |
|
|
358 |
RETURN |
RETURN |