128 |
DO I=1-Olx,sNx+Olx |
DO I=1-Olx,sNx+Olx |
129 |
tmpfld(I,J,bi,bj) = 0. |
tmpfld(I,J,bi,bj) = 0. |
130 |
ksurfC(I,J,bi,bj) = Nr+1 |
ksurfC(I,J,bi,bj) = Nr+1 |
131 |
|
maskH(i,j,bi,bj) = 0. |
132 |
Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj) |
Ro_surf(I,J,bi,bj) = R_low(I,J,bi,bj) |
133 |
DO K=Nr,1,-1 |
DO K=Nr,1,-1 |
134 |
Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj) |
Ro_surf(I,J,bi,bj) = Ro_surf(I,J,bi,bj) |
135 |
& + drF(k)*hFacC(I,J,K,bi,bj) |
& + drF(k)*hFacC(I,J,K,bi,bj) |
136 |
IF (hFacC(I,J,K,bi,bj).NE.0.) THEN |
IF (hFacC(I,J,K,bi,bj).NE.0.) THEN |
137 |
ksurfC(I,J,bi,bj) = k |
ksurfC(I,J,bi,bj) = k |
138 |
|
maskH(i,j,bi,bj) = 1. |
139 |
tmpfld(i,j,bi,bj) = tmpfld(i,j,bi,bj) + 1. |
tmpfld(i,j,bi,bj) = tmpfld(i,j,bi,bj) + 1. |
140 |
ENDIF |
ENDIF |
141 |
ENDDO |
ENDDO |