79 |
DO j = 1,sNy |
DO j = 1,sNy |
80 |
DO i = 1,sNx+1 |
DO i = 1,sNx+1 |
81 |
#ifdef ALLOW_LONGSTEP |
#ifdef ALLOW_LONGSTEP |
82 |
C at first timestep we don't have averaged velocities yet - |
C at first timestep we do not have averaged velocities yet - |
83 |
C use initial velocities instead |
C use initial velocities instead |
84 |
IF ( trIter.GT.nIter0 ) THEN |
IF ( trIter.GT.nIter0 ) THEN |
85 |
dummy(i,j,k,bi,bj) = |
dummy(i,j,k,bi,bj) = |
112 |
DO j = 1,sNy+1 |
DO j = 1,sNy+1 |
113 |
DO i = 1,sNx |
DO i = 1,sNx |
114 |
#ifdef ALLOW_LONGSTEP |
#ifdef ALLOW_LONGSTEP |
115 |
C at first timestep we don't have averaged velocities yet - |
C at first timestep we do not have averaged velocities yet - |
116 |
C use initial velocities instead |
C use initial velocities instead |
117 |
IF ( trIter.GT.nIter0 ) THEN |
IF ( trIter.GT.nIter0 ) THEN |
118 |
dummy(i,j,k,bi,bj) = |
dummy(i,j,k,bi,bj) = |
146 |
DO j = 1,sNy |
DO j = 1,sNy |
147 |
DO i = 1,sNx |
DO i = 1,sNx |
148 |
#ifdef ALLOW_LONGSTEP |
#ifdef ALLOW_LONGSTEP |
149 |
C at first timestep we don't have averaged velocities yet - |
C at first timestep we do not have averaged velocities yet - |
150 |
C use initial velocities instead |
C use initial velocities instead |
151 |
IF ( trIter.GT.nIter0 ) THEN |
IF ( trIter.GT.nIter0 ) THEN |
152 |
dummy(i,j,k,bi,bj) = LS_wVel(i,j,k,bi,bj) * |
dummy(i,j,k,bi,bj) = LS_wVel(i,j,k,bi,bj) * |