95 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
96 |
|
|
97 |
C-- Calculate gradient in X direction: |
C-- Calculate gradient in X direction: |
98 |
|
#ifndef ALLOW_AUTODIFF_TAMC |
99 |
IF ( .NOT.Shap_alwaysExchTr |
IF ( .NOT.Shap_alwaysExchTr |
100 |
& .AND. useCubedSphereExchange ) THEN |
& .AND. useCubedSphereExchange ) THEN |
101 |
C to compute d/dx(tmpFld), fill corners with appropriate values: |
C to compute d/dx(tmpFld), fill corners with appropriate values: |
103 |
& tmpFld(1-Olx,1-Oly,k,bi,bj), |
& tmpFld(1-Olx,1-Oly,k,bi,bj), |
104 |
& bi,bj, myThid ) |
& bi,bj, myThid ) |
105 |
ENDIF |
ENDIF |
106 |
|
#endif |
107 |
IF ( n.LE.nShapComput ) THEN |
IF ( n.LE.nShapComput ) THEN |
108 |
C- Computational space: del_i |
C- Computational space: del_i |
109 |
DO j=0,sNy+1 |
DO j=0,sNy+1 |
126 |
ENDIF |
ENDIF |
127 |
|
|
128 |
C-- Calculate gradient in Y direction: |
C-- Calculate gradient in Y direction: |
129 |
|
#ifndef ALLOW_AUTODIFF_TAMC |
130 |
IF ( .NOT.Shap_alwaysExchTr |
IF ( .NOT.Shap_alwaysExchTr |
131 |
& .AND. useCubedSphereExchange ) THEN |
& .AND. useCubedSphereExchange ) THEN |
132 |
C to compute d/dy(tmpFld), fill corners with appropriate values: |
C to compute d/dy(tmpFld), fill corners with appropriate values: |
134 |
& tmpFld(1-Olx,1-Oly,k,bi,bj), |
& tmpFld(1-Olx,1-Oly,k,bi,bj), |
135 |
& bi,bj, myThid ) |
& bi,bj, myThid ) |
136 |
ENDIF |
ENDIF |
137 |
|
#endif |
138 |
IF ( n.LE.nShapComput ) THEN |
IF ( n.LE.nShapComput ) THEN |
139 |
C- Computational space: del_j |
C- Computational space: del_j |
140 |
DO j=0,sNy+2 |
DO j=0,sNy+2 |