112 |
|
|
113 |
C-- First call the exchanges for the two components |
C-- First call the exchanges for the two components |
114 |
|
|
115 |
CALL EXCH2_RX1_CUBE( uPhi, 'T ', |
CALL EXCH2_RX1_CUBE( uPhi, .FALSE., 'T ', |
116 |
I OLw, OLe, OLs, OLn, myNz, |
I OLw, OLe, OLs, OLn, myNz, |
117 |
I exchWidthX, exchWidthY, |
I exchWidthX, exchWidthY, |
118 |
I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid ) |
I EXCH_UPDATE_CORNERS, myThid ) |
119 |
CALL EXCH2_RX1_CUBE( uPhi, 'T ', |
CALL EXCH2_RX1_CUBE( uPhi, .FALSE., 'T ', |
120 |
I OLw, OLe, OLs, OLn, myNz, |
I OLw, OLe, OLs, OLn, myNz, |
121 |
I exchWidthX, exchWidthY, |
I exchWidthX, exchWidthY, |
122 |
I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid ) |
I EXCH_UPDATE_CORNERS, myThid ) |
123 |
|
|
124 |
CALL EXCH2_RX1_CUBE( vPhi, 'T ', |
CALL EXCH2_RX1_CUBE( vPhi, .FALSE., 'T ', |
125 |
I OLw, OLe, OLs, OLn, myNz, |
I OLw, OLe, OLs, OLn, myNz, |
126 |
I exchWidthX, exchWidthY, |
I exchWidthX, exchWidthY, |
127 |
I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid ) |
I EXCH_UPDATE_CORNERS, myThid ) |
128 |
CALL EXCH2_RX1_CUBE( vPhi, 'T ', |
CALL EXCH2_RX1_CUBE( vPhi, .FALSE., 'T ', |
129 |
I OLw, OLe, OLs, OLn, myNz, |
I OLw, OLe, OLs, OLn, myNz, |
130 |
I exchWidthX, exchWidthY, |
I exchWidthX, exchWidthY, |
131 |
I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid ) |
I EXCH_UPDATE_CORNERS, myThid ) |
132 |
|
|
133 |
C- note: can substitute the low-level S/R calls above with: |
C- note: can substitute the low-level S/R calls above with: |
134 |
c CALL EXCH2_3D_RX( uPhi, myNz, myThid ) |
c CALL EXCH2_3D_RX( uPhi, myNz, myThid ) |
149 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
150 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
151 |
|
|
152 |
C- Now choose what to do at each edge of the halo based on which face |
C- Choose what to do at each edge of the halo based on which face we are |
|
C (we assume that bj is always=1) |
|
153 |
myTile = W2_myTileList(bi,bj) |
myTile = W2_myTileList(bi,bj) |
154 |
myFace = exch2_myFace(myTile) |
myFace = exch2_myFace(myTile) |
155 |
|
|