100 |
& *maskS(i,j,k,bi,bj) |
& *maskS(i,j,k,bi,bj) |
101 |
& *maskC(i,j,k,bi,bj)*maskC(i,j-1,k,bi,bj) |
& *maskC(i,j,k,bi,bj)*maskC(i,j-1,k,bi,bj) |
102 |
theta_bar(k) = theta_bar(k) + |
theta_bar(k) = theta_bar(k) + |
|
cph & theta(i,j,k,bi,bj)*dxG(i,j,bi,bj) |
|
|
cph & *maskC(i,j,k,bi,bj) |
|
103 |
& 0.5*( cMeanTheta(i,j,k,bi,bj) |
& 0.5*( cMeanTheta(i,j,k,bi,bj) |
104 |
& +cMeanTheta(i,j-1,k,bi,bj) ) |
& +cMeanTheta(i,j-1,k,bi,bj) ) |
105 |
& *maskS(i,j,k,bi,bj)*dxG(i,j,bi,bj) |
& *maskS(i,j,k,bi,bj)*dxG(i,j,bi,bj) |
144 |
& *maskW(i,j,k,bi,bj) |
& *maskW(i,j,k,bi,bj) |
145 |
& *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj) |
& *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj) |
146 |
theta_bar(k) = theta_bar(k) + |
theta_bar(k) = theta_bar(k) + |
|
cph & + theta(i,j,k,bi,bj)*dyG(i,j,bi,bj) |
|
|
cph & *maskC(i,j,k,bi,bj) |
|
147 |
& 0.5*( cMeanTheta(i,j,k,bi,bj) |
& 0.5*( cMeanTheta(i,j,k,bi,bj) |
148 |
& +cMeanTheta(i-,j,k,bi,bj) ) |
& +cMeanTheta(i-,j,k,bi,bj) ) |
149 |
& *maskW(i,j,k,bi,bj)*dyG(i,j,bi,bj) |
& *maskW(i,j,k,bi,bj)*dyG(i,j,bi,bj) |