64 |
C-- Calculate tracer cell face open areas |
C-- Calculate tracer cell face open areas |
65 |
DO j=jMin,jMax |
DO j=jMin,jMax |
66 |
DO i=iMin,iMax |
DO i=iMin,iMax |
67 |
xA(i,j) = dyG(i,j,bi,bj)*dzF(k)*hFacW(i,j,k,bi,bj) |
xA(i,j) = _dyG(i,j,bi,bj)*dzF(k)*_hFacW(i,j,k,bi,bj) |
68 |
yA(i,j) = _dxG(i,j,bi,bj)*dzF(k)*hFacS(i,j,k,bi,bj) |
yA(i,j) = _dxG(i,j,bi,bj)*dzF(k)*_hFacS(i,j,k,bi,bj) |
69 |
ENDDO |
ENDDO |
70 |
ENDDO |
ENDDO |
71 |
|
|
92 |
DO j=jMin,jMax |
DO j=jMin,jMax |
93 |
DO i=iMin,iMax |
DO i=iMin,iMax |
94 |
maskC(i,j) = 1. |
maskC(i,j) = 1. |
95 |
if (hFacC(i,j,k,bi,bj).eq.0.) maskC(i,j)=0. |
if (_hFacC(i,j,k,bi,bj).eq.0.) maskC(i,j)=0. |
96 |
maskUp(i,j) = 1. |
maskUp(i,j) = 1. |
97 |
if (hFacC(i,j,kM1,bi,bj).eq.0.) maskUp(i,j)=0. |
if (_hFacC(i,j,kM1,bi,bj).eq.0.) maskUp(i,j)=0. |
98 |
ENDDO |
ENDDO |
99 |
ENDDO |
ENDDO |
100 |
|
|