128 |
IF (k.EQ.1) THEN |
IF (k.EQ.1) THEN |
129 |
DO j=jMin,jMax |
DO j=jMin,jMax |
130 |
DO i=iMin,iMax |
DO i=iMin,iMax |
131 |
phiHyd(i,j,k) = phi0surf(i,j,bi,bj) |
c phiHyd(i,j,k) = phi0surf(i,j,bi,bj) |
132 |
c phiHyd(i,j,k) = 0. |
phiHyd(i,j,k) = 0. |
133 |
ENDDO |
ENDDO |
134 |
ENDDO |
ENDDO |
135 |
ENDIF |
ENDIF |
227 |
IF (k.EQ.1) THEN |
IF (k.EQ.1) THEN |
228 |
DO j=jMin,jMax |
DO j=jMin,jMax |
229 |
DO i=iMin,iMax |
DO i=iMin,iMax |
230 |
phiHyd(i,j,k) = phi0surf(i,j,bi,bj) |
c phiHyd(i,j,k) = phi0surf(i,j,bi,bj) |
231 |
c phiHyd(i,j,k) = 0. |
phiHyd(i,j,k) = 0. |
232 |
ENDDO |
ENDDO |
233 |
ENDDO |
ENDDO |
234 |
ENDIF |
ENDIF |
331 |
& -((rC(K)/atm_Po)**atm_kappa) ) |
& -((rC(K)/atm_Po)**atm_kappa) ) |
332 |
DO j=jMin,jMax |
DO j=jMin,jMax |
333 |
DO i=iMin,iMax |
DO i=iMin,iMax |
334 |
c phiHyd(i,j,K)= |
c phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
335 |
phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
phiHyd(i,j,K)= |
336 |
& ddPIp*maskC(i,j,K,bi,bj) |
& ddPIp*maskC(i,j,K,bi,bj) |
337 |
& *(tFld(I,J,K,bi,bj)-tRef(K)) |
& *(tFld(I,J,K,bi,bj)-tRef(K)) |
338 |
ENDDO |
ENDDO |
372 |
& -((rC(K)/atm_Po)**atm_kappa) ) |
& -((rC(K)/atm_Po)**atm_kappa) ) |
373 |
DO j=jMin,jMax |
DO j=jMin,jMax |
374 |
DO i=iMin,iMax |
DO i=iMin,iMax |
375 |
c phiHyd(i,j,K)= |
c phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
376 |
phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
phiHyd(i,j,K)= |
377 |
& ddPIp*_hFacC(I,J, K ,bi,bj) |
& ddPIp*_hFacC(I,J, K ,bi,bj) |
378 |
& *(tFld(I,J, K ,bi,bj)-tRef( K )) |
& *(tFld(I,J, K ,bi,bj)-tRef( K )) |
379 |
ENDDO |
ENDDO |
411 |
& -((rC(Kp1)/atm_Po)**atm_kappa) ) |
& -((rC(Kp1)/atm_Po)**atm_kappa) ) |
412 |
DO j=jMin,jMax |
DO j=jMin,jMax |
413 |
DO i=iMin,iMax |
DO i=iMin,iMax |
414 |
c phiHyd(i,j,K)= |
c phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
415 |
phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
phiHyd(i,j,K)= |
416 |
& ( ddPIm*max(zero, _hFacC(i,j,K,bi,bj)-half) |
& ( ddPIm*max(zero, _hFacC(i,j,K,bi,bj)-half) |
417 |
& +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)-half) ) |
& +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)-half) ) |
418 |
& *(tFld(i,j, K ,bi,bj)-tRef( K )) |
& *(tFld(i,j, K ,bi,bj)-tRef( K )) |
457 |
& -((rC(Kp1)/atm_Po)**atm_kappa) ) |
& -((rC(Kp1)/atm_Po)**atm_kappa) ) |
458 |
DO j=jMin,jMax |
DO j=jMin,jMax |
459 |
DO i=iMin,iMax |
DO i=iMin,iMax |
460 |
c phiHyd(i,j,K)= |
c phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
461 |
phiHyd(i,j,K)= phi0surf(i,j,bi,bj)+ |
phiHyd(i,j,K)= |
462 |
& ( ddPIm*max(zero,(_hFacC(i,j,K,bi,bj)-one)*ratioRm+half) |
& ( ddPIm*max(zero,(_hFacC(i,j,K,bi,bj)-one)*ratioRm+half) |
463 |
& +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)*ratioRp -half) ) |
& +ddPIp*min(zero, _hFacC(i,j,K,bi,bj)*ratioRp -half) ) |
464 |
& *(tFld(i,j, K ,bi,bj)-tRef( K )) |
& *(tFld(i,j, K ,bi,bj)-tRef( K )) |