56 |
C-- Initialise arrays (NLFS using r-coordinate): |
C-- Initialise arrays (NLFS using r-coordinate): |
57 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
58 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
59 |
hFac_surfC(i,j,bi,bj) = _hFacC(i,j,1,bi,bj) |
hFac_surfC(i,j,bi,bj) = 0. |
60 |
hFac_surfW(i,j,bi,bj) = _hFacW(i,j,1,bi,bj) |
hFac_surfW(i,j,bi,bj) = 0. |
61 |
hFac_surfS(i,j,bi,bj) = _hFacS(i,j,1,bi,bj) |
hFac_surfS(i,j,bi,bj) = 0. |
62 |
hFac_surfNm1C(i,j,bi,bj) = _hFacC(i,j,1,bi,bj) |
hFac_surfNm1C(i,j,bi,bj) = 0. |
63 |
hFac_surfNm1W(i,j,bi,bj) = _hFacW(i,j,1,bi,bj) |
hFac_surfNm1W(i,j,bi,bj) = 0. |
64 |
hFac_surfNm1S(i,j,bi,bj) = _hFacS(i,j,1,bi,bj) |
hFac_surfNm1S(i,j,bi,bj) = 0. |
65 |
PmEpR(i,j,bi,bj) = 0. |
PmEpR(i,j,bi,bj) = 0. |
66 |
Rmin_surf(i,j,bi,bj) = Ro_surf(i,j,bi,bj) |
Rmin_surf(i,j,bi,bj) = Ro_surf(i,j,bi,bj) |
67 |
ENDDO |
ENDDO |