79 |
DO j=jMin,jMax |
DO j=jMin,jMax |
80 |
DO i=iMin,iMax |
DO i=iMin,iMax |
81 |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
82 |
& +foFacMom*surfaceTendencyU(i,j,bi,bj) |
& +foFacMom*surfaceForcingU(i,j,bi,bj) |
83 |
& *_maskW(i,j,kLev,bi,bj) |
& *recip_drF(kLev)*recip_hFacW(i,j,kLev,bi,bj) |
84 |
ENDDO |
ENDDO |
85 |
ENDDO |
ENDDO |
86 |
ENDIF |
ENDIF |
167 |
DO j=jMin,jMax |
DO j=jMin,jMax |
168 |
DO i=iMin,iMax |
DO i=iMin,iMax |
169 |
gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) |
gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) |
170 |
& +foFacMom*surfaceTendencyV(i,j,bi,bj) |
& +foFacMom*surfaceForcingV(i,j,bi,bj) |
171 |
& *_maskS(i,j,kLev,bi,bj) |
& *recip_drF(kLev)*recip_hFacS(i,j,kLev,bi,bj) |
172 |
ENDDO |
ENDDO |
173 |
ENDDO |
ENDDO |
174 |
ENDIF |
ENDIF |
264 |
DO j=jMin,jMax |
DO j=jMin,jMax |
265 |
DO i=iMin,iMax |
DO i=iMin,iMax |
266 |
gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj) |
gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj) |
267 |
& +maskC(i,j,kLev,bi,bj)*surfaceTendencyT(i,j,bi,bj) |
& +surfaceForcingT(i,j,bi,bj) |
268 |
|
& *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj) |
269 |
ENDDO |
ENDDO |
270 |
ENDDO |
ENDDO |
271 |
ENDIF |
ENDIF |
289 |
& -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj) |
& -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj) |
290 |
& -swfracb(2)*maskC(i,j,kp1, bi,bj)) |
& -swfracb(2)*maskC(i,j,kp1, bi,bj)) |
291 |
& *recip_Cp*recip_rhoConst*recip_drF(klev) |
& *recip_Cp*recip_rhoConst*recip_drF(klev) |
292 |
|
c & *recip_hFacC(i,j,kLev,bi,bj) |
293 |
ENDDO |
ENDDO |
294 |
ENDDO |
ENDDO |
295 |
#endif |
#endif |
379 |
DO j=jMin,jMax |
DO j=jMin,jMax |
380 |
DO i=iMin,iMax |
DO i=iMin,iMax |
381 |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
382 |
& +maskC(i,j,kLev,bi,bj)*surfaceTendencyS(i,j,bi,bj) |
& +surfaceForcingS(i,j,bi,bj) |
383 |
|
& *recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj) |
384 |
ENDDO |
ENDDO |
385 |
ENDDO |
ENDDO |
386 |
ENDIF |
ENDIF |