33 |
#include "PTRACERS_SIZE.h" |
#include "PTRACERS_SIZE.h" |
34 |
#include "PTRACERS_PARAMS.h" |
#include "PTRACERS_PARAMS.h" |
35 |
c #include "PTRACERS_FIELDS.h" |
c #include "PTRACERS_FIELDS.h" |
36 |
c #include "GAD.h" |
#include "GAD.h" |
37 |
|
|
38 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
39 |
C == Routine arguments == |
C == Routine arguments == |
52 |
C !LOCAL VARIABLES: |
C !LOCAL VARIABLES: |
53 |
C == Local variables == |
C == Local variables == |
54 |
C I, J :: Loop counters |
C I, J :: Loop counters |
55 |
INTEGER i,j,iTr |
INTEGER i,j,iTr,trIdentity |
56 |
_RL KbryanLewis79 |
_RL KbryanLewis79 |
57 |
#ifdef ALLOW_BL79_LAT_VARY |
#ifdef ALLOW_BL79_LAT_VARY |
58 |
_RL KbryanLewisEQ |
_RL KbryanLewisEQ |
69 |
C Loop over tracers |
C Loop over tracers |
70 |
DO iTr=1,PTRACERS_numInUse |
DO iTr=1,PTRACERS_numInUse |
71 |
|
|
72 |
|
trIdentity = iTr + GAD_TR1 - 1 |
73 |
|
|
74 |
IF ( .NOT. PTRACERS_useKPP(iTr) ) THEN |
IF ( .NOT. PTRACERS_useKPP(iTr) ) THEN |
75 |
DO j = 1-Oly, sNy+Oly |
DO j = 1-Oly, sNy+Oly |
76 |
DO i = 1-Olx, sNx+Olx |
DO i = 1-Olx, sNx+Olx |
105 |
CALL GMREDI_CALC_DIFF( |
CALL GMREDI_CALC_DIFF( |
106 |
I bi,bj,iMin,iMax,jMin,jMax,k,1, |
I bi,bj,iMin,iMax,jMin,jMax,k,1, |
107 |
U KappaRtr(1-Olx,1-Oly,iTr), |
U KappaRtr(1-Olx,1-Oly,iTr), |
108 |
I myThid) |
I trIdentity,myThid) |
109 |
ENDIF |
ENDIF |
110 |
#endif |
#endif |
111 |
|
|