Parent Directory | Revision Log | Revision Graph | Patch
revision 1.27 by jmc, Sun Feb 9 02:58:39 2003 UTC | revision 1.28 by jmc, Tue Feb 11 02:31:32 2003 UTC | |
---|---|---|
# | Line 173 C conserve KE+PE exactly even | Line 173 C conserve KE+PE exactly even |
173 | C | C |
174 | IF (k.LT.Nr) phiHyd(i,j,k+1)=phiHyd(i,j,k) | IF (k.LT.Nr) phiHyd(i,j,k+1)=phiHyd(i,j,k) |
175 | & + drF(K)*gravity*alphaRho(i,j)*recip_rhoConst | & + drF(K)*gravity*alphaRho(i,j)*recip_rhoConst |
176 | phiHyd(i,j,k)=phiHyd(i,j,k)+ | phiHyd(i,j,k)=phiHyd(i,j,k) |
177 | & + half*drF(K)*gravity*alphaRho(i,j)*recip_rhoConst | & + half*drF(K)*gravity*alphaRho(i,j)*recip_rhoConst |
178 | ||
179 | ENDDO | ENDDO |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |