Parent Directory | Revision Log | Revision Graph | Patch
revision 1.8 by jmc, Tue Oct 18 16:03:55 2005 UTC | revision 1.9 by heimbach, Wed Jun 7 01:55:14 2006 UTC | |
---|---|---|
# | Line 54 C with no need to compute thetaM ( | Line 54 C with no need to compute thetaM ( |
54 | ||
55 | c uFld = uVel(i,j,k,bi,bj) | c uFld = uVel(i,j,k,bi,bj) |
56 | uFld = uTrans(i,j)*recip_dyG(i,j,bi,bj) | uFld = uTrans(i,j)*recip_dyG(i,j,bi,bj) |
57 | & *recip_drF(k)*recip_hFacW(i,j,k,bi,bj) | & *recip_drF(k)*_recip_hFacW(i,j,k,bi,bj) |
58 | cfl=abs(uFld*deltaTloc*recip_dxC(i,j,bi,bj)) | cfl=abs(uFld*deltaTloc*recip_dxC(i,j,bi,bj)) |
59 | d0=(2. _d 0 -cfl)*(1. _d 0 -cfl)*oneSixth | d0=(2. _d 0 -cfl)*(1. _d 0 -cfl)*oneSixth |
60 | d1=(1. _d 0 -cfl*cfl)*oneSixth | d1=(1. _d 0 -cfl*cfl)*oneSixth |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |