62 |
ab05 = -0.5 - abEps |
ab05 = -0.5 - abEps |
63 |
|
|
64 |
C-- stagger time stepping: grad Phi_Hyp is not in gU,gV and needs to be added: |
C-- stagger time stepping: grad Phi_Hyp is not in gU,gV and needs to be added: |
65 |
IF (.FALSE.) THEN |
c IF (.FALSE.) THEN |
66 |
c IF (staggerTimeStep) THEN |
IF (staggerTimeStep) THEN |
67 |
phxFac = pfFacMom |
phxFac = pfFacMom |
68 |
phyFac = pfFacMom |
phyFac = pfFacMom |
69 |
ELSE |
ELSE |
79 |
& +ab05*(etaNm1(i,j,bi,bj)*Bo_surf(i,j,bi,bj) ) |
& +ab05*(etaNm1(i,j,bi,bj)*Bo_surf(i,j,bi,bj) ) |
80 |
ENDDO |
ENDDO |
81 |
ENDDO |
ENDDO |
82 |
IF (staggerTimeStep) THEN |
IF (.FALSE.) THEN |
83 |
|
c IF (staggerTimeStep) THEN |
84 |
DO j=jMin,jMax |
DO j=jMin,jMax |
85 |
DO i=iMin,iMax |
DO i=iMin,iMax |
86 |
pf(i,j) = pf(i,j)+phi_hyd(i,j,k) |
pf(i,j) = pf(i,j)+phi_hyd(i,j,k) |