36 |
|
|
37 |
#include "SIZE.h" |
#include "SIZE.h" |
38 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
|
c#include "EESUPPORT.h" |
|
39 |
#include "W2_EXCH2_SIZE.h" |
#include "W2_EXCH2_SIZE.h" |
40 |
#include "W2_EXCH2_TOPOLOGY.h" |
#include "W2_EXCH2_TOPOLOGY.h" |
41 |
|
|
68 |
negOne = 1. |
negOne = 1. |
69 |
IF (withSigns) negOne = -1. |
IF (withSigns) negOne = -1. |
70 |
|
|
|
IF ( useCubedSphereExchange ) THEN |
|
|
C--- using CubedSphereExchange: |
|
|
|
|
71 |
C-- First call the exchanges for the two components, ignoring the Sign |
C-- First call the exchanges for the two components, ignoring the Sign |
72 |
C note the order: vPhi,uPhi on D-grid are co-located with (u,v)_Cgrid |
C note the order: vPhi,uPhi on D-grid are co-located with (u,v)_Cgrid |
73 |
|
|
95 |
I .FALSE., myNz, myThid ) |
I .FALSE., myNz, myThid ) |
96 |
#endif |
#endif |
97 |
|
|
98 |
|
IF ( useCubedSphereExchange ) THEN |
99 |
|
C--- using CubedSphereExchange: |
100 |
|
|
101 |
C-- Then we may need to switch the signs depending on which cube face |
C-- Then we may need to switch the signs depending on which cube face |
102 |
C we are located. |
C we are located. |
103 |
|
|
107 |
|
|
108 |
C- Now choose what to do at each edge of the halo based on which face |
C- Now choose what to do at each edge of the halo based on which face |
109 |
C (we assume that bj is always=1) |
C (we assume that bj is always=1) |
110 |
myTile = W2_myTileList(bi) |
myTile = W2_myTileList(bi,bj) |
111 |
myFace = exch2_myFace(myTile) |
myFace = exch2_myFace(myTile) |
112 |
|
|
113 |
C-- Loops on level index: |
C-- Loops on level index: |
200 |
ENDDO |
ENDDO |
201 |
ENDDO |
ENDDO |
202 |
|
|
|
ELSE |
|
|
C--- not using CubedSphereExchange: |
|
|
|
|
|
#ifndef AUTODIFF_EXCH2 |
|
|
CALL EXCH_RX( uPhi, |
|
|
I OLw, OLe, OLs, OLn, myNz, |
|
|
I exchWidthX, exchWidthY, |
|
|
I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid ) |
|
|
CALL EXCH_RX( vPhi, |
|
|
I OLw, OLe, OLs, OLn, myNz, |
|
|
I exchWidthX, exchWidthY, |
|
|
I FORWARD_SIMULATION, EXCH_UPDATE_CORNERS, myThid ) |
|
|
#endif |
|
|
|
|
203 |
C--- using or not using CubedSphereExchange: end |
C--- using or not using CubedSphereExchange: end |
204 |
ENDIF |
ENDIF |
205 |
|
|