170 |
phiHyd(i,j,K)=phiHyd(i,j,K-1) |
phiHyd(i,j,K)=phiHyd(i,j,K-1) |
171 |
& -( ddRm*(theta(I,J,K-1,bi,bj)-tRef(K-1)) |
& -( ddRm*(theta(I,J,K-1,bi,bj)-tRef(K-1)) |
172 |
& +ddRp*(theta(I,J, K ,bi,bj)-tRef( K )) ) |
& +ddRp*(theta(I,J, K ,bi,bj)-tRef( K )) ) |
173 |
C Old code bug looked like this |
C Old code (atmos-exact) looked like this |
174 |
Cold phiHyd(i,j,K)=phiHyd(i,j,K-1)-(ddRm1* |
Cold phiHyd(i,j,K)=phiHyd(i,j,K-1) - ddRm1* |
175 |
Cold & (theta(I,J,K-1,bi,bj)+theta(I,J,K,bi,bj))-tRef(K)) |
Cold & (theta(I,J,K-1,bi,bj)+theta(I,J,K,bi,bj)-2.*tRef(K)) |
176 |
ENDDO |
ENDDO |
177 |
ENDDO |
ENDDO |
178 |
ENDIF |
ENDIF |