117 |
|
|
118 |
C- Re-calculate Reference surface position, taking into account hFacC |
C- Re-calculate Reference surface position, taking into account hFacC |
119 |
C initialize Total column fluid thickness and surface k index |
C initialize Total column fluid thickness and surface k index |
120 |
|
C Note: if no fluid (continent) ==> ksurf = Nr+1 |
121 |
DO J=1-Oly,sNy+Oly |
DO J=1-Oly,sNy+Oly |
122 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
123 |
tmpfld(I,J,bi,bj) = 0. |
tmpfld(I,J,bi,bj) = 0. |
124 |
k_surf(I,J,bi,bj) = Nr |
ksurfC(I,J,bi,bj) = Nr+1 |
125 |
Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj) |
Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj) |
126 |
DO K=Nr,1,-1 |
DO K=Nr,1,-1 |
127 |
Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj) |
Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj) |
128 |
& + drF(k)*hFacC(I,J,K,bi,bj) |
& + drF(k)*hFacC(I,J,K,bi,bj) |
129 |
IF (hFacC(I,J,K,bi,bj).NE.0.) THEN |
IF (hFacC(I,J,K,bi,bj).NE.0.) THEN |
130 |
k_surf(I,J,bi,bj) = k |
ksurfC(I,J,bi,bj) = k |
131 |
tmpfld(i,j,bi,bj) = tmpfld(i,j,bi,bj) + 1. |
tmpfld(i,j,bi,bj) = tmpfld(i,j,bi,bj) + 1. |
132 |
ENDIF |
ENDIF |
133 |
ENDDO |
ENDDO |
241 |
ENDDO |
ENDDO |
242 |
ENDDO |
ENDDO |
243 |
ENDDO |
ENDDO |
244 |
|
C- Calculate surface k index for interface W & S (U & V points) |
245 |
|
DO J=1-Oly,sNy+Oly |
246 |
|
DO I=1-Olx,sNx+Olx |
247 |
|
ksurfW(I,J,bi,bj) = Nr+1 |
248 |
|
ksurfS(I,J,bi,bj) = Nr+1 |
249 |
|
DO k=Nr,1,-1 |
250 |
|
IF (hFacW(I,J,K,bi,bj).NE.0.) ksurfW(I,J,bi,bj) = k |
251 |
|
IF (hFacS(I,J,K,bi,bj).NE.0.) ksurfS(I,J,bi,bj) = k |
252 |
|
ENDDO |
253 |
|
ENDDO |
254 |
|
ENDDO |
255 |
|
C - end bi,bj loops. |
256 |
ENDDO |
ENDDO |
257 |
ENDDO |
ENDDO |
258 |
C _EXCH_XYZ_R4(recip_HFacC , myThid ) |
C _EXCH_XYZ_R4(recip_HFacC , myThid ) |