4 |
#include "GAD_OPTIONS.h" |
#include "GAD_OPTIONS.h" |
5 |
|
|
6 |
SUBROUTINE GAD_OS7MP_ADV_X( |
SUBROUTINE GAD_OS7MP_ADV_X( |
7 |
I bi,bj,k,deltaTloc, |
I bi,bj,k, calcCFL, deltaTloc, |
8 |
I uTrans, uFld, |
I uTrans, uFld, |
9 |
I maskLocW, Q, |
I maskLocW, Q, |
10 |
O uT, |
O uT, |
23 |
|
|
24 |
C == Routine arguments == |
C == Routine arguments == |
25 |
INTEGER bi,bj,k |
INTEGER bi,bj,k |
26 |
|
LOGICAL calcCFL |
27 |
_RL deltaTloc |
_RL deltaTloc |
28 |
_RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
29 |
_RL uFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL uFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
54 |
DO i=1-Olx+4,sNx+Olx-3 |
DO i=1-Olx+4,sNx+Olx-3 |
55 |
|
|
56 |
uLoc = uFld(i,j) |
uLoc = uFld(i,j) |
57 |
cfl = abs(uLoc*deltaTloc*recip_dxC(i,j,bi,bj)) |
cfl = uLoc |
58 |
|
IF ( calcCFL ) cfl = abs(uLoc*deltaTloc*recip_dxC(i,j,bi,bj)) |
59 |
|
|
60 |
IF (uLoc.gt.0.) THEN |
IF (uLoc.gt.0.) THEN |
61 |
Qippp = Q(i+2,j) |
Qippp = Q(i+2,j) |