71 |
_RL oneSixth |
_RL oneSixth |
72 |
PARAMETER(oneSixth=1.D0/6.D0) |
PARAMETER(oneSixth=1.D0/6.D0) |
73 |
|
|
74 |
|
C loop range for computing vertical advection tendency |
75 |
|
C iMinAdvR,iMaxAdvR :: 1rst index (X-dir) loop range for vertical advection |
76 |
|
C jMinAdvR,jMaxAdvR :: 2nd index (Y-dir) loop range for vertical advection |
77 |
|
INTEGER iMinAdvR, iMaxAdvR, jMinAdvR, jMaxAdvR |
78 |
|
c PARAMETER ( iMinAdvR = 1-OLx , iMaxAdvR = sNx+OLx ) |
79 |
|
c PARAMETER ( jMinAdvR = 1-OLy , jMaxAdvR = sNy+OLy ) |
80 |
|
C- note: we use to compute vertical advection tracer tendency everywhere |
81 |
|
C (overlap included) as above, but really needs valid tracer tendency |
82 |
|
C in interior only (as below): |
83 |
|
PARAMETER ( iMinAdvR = 1 , iMaxAdvR = sNx ) |
84 |
|
PARAMETER ( jMinAdvR = 1 , jMaxAdvR = sNy ) |
85 |
|
|
86 |
C Differentiate between tracers (needed for KPP - arrgh!!!) |
C Differentiate between tracers (needed for KPP - arrgh!!!) |
87 |
cph and GMRedi arrgh*arrgh!!!) |
cph and GMRedi arrgh*arrgh!!!) |
88 |
cph indices are used for TAF key computations, so need to |
cph indices are used for TAF key computations, so need to |