160 |
tv = tv+p5*(theta(i,j,k,bi,bj) |
tv = tv+p5*(theta(i,j,k,bi,bj) |
161 |
$ + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj) |
$ + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj) |
162 |
$ * _dxG(i,j,bi,bj) |
$ * _dxG(i,j,bi,bj) |
163 |
& * drF(k)*_hFacS(i,j,k,bi,bj) |
& * drF(k)*maskS(i,j,k,bi,bj) |
164 |
$ *HeatCapacity_Cp*rhoNil |
$ *HeatCapacity_Cp*rhoNil |
165 |
endif |
endif |
166 |
enddo |
enddo |
204 |
$ +p5*(theta(i,j,k,bi,bj) |
$ +p5*(theta(i,j,k,bi,bj) |
205 |
$ + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj) |
$ + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj) |
206 |
$ * _dxG(i,j,bi,bj) |
$ * _dxG(i,j,bi,bj) |
207 |
& * drF(k)*_hFacS(i,j,k,bi,bj) |
& * drF(k)*maskS(i,j,k,bi,bj) |
208 |
$ *HeatCapacity_Cp*rhoNil |
$ *HeatCapacity_Cp*rhoNil |
209 |
endif |
endif |
210 |
enddo |
enddo |