4 |
#include "GAD_OPTIONS.h" |
#include "GAD_OPTIONS.h" |
5 |
|
|
6 |
SUBROUTINE GAD_OS7MP_ADV_Y( |
SUBROUTINE GAD_OS7MP_ADV_Y( |
7 |
I bi,bj,k,deltaTloc, |
I bi,bj,k, calcCFL, deltaTloc, |
8 |
I vTrans, vFld, |
I vTrans, vFld, |
9 |
I maskLocS, Q, |
I maskLocS, Q, |
10 |
O vT, |
O vT, |
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 vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
29 |
_RL vFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL vFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
56 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
57 |
|
|
58 |
vLoc = vFld(i,j) |
vLoc = vFld(i,j) |
59 |
cfl = abs(vLoc*deltaTloc*recip_dyC(i,j,bi,bj)) |
cfl = vLoc |
60 |
|
IF ( calcCFL ) cfl = abs(vLoc*deltaTloc*recip_dyC(i,j,bi,bj)) |
61 |
|
|
62 |
IF (vLoc.gt.0.) THEN |
IF (vLoc.gt.0.) THEN |
63 |
Qippp = Q(i,j+2) |
Qippp = Q(i,j+2) |