19 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
20 |
#include "PARAMS.h" |
#include "PARAMS.h" |
21 |
#include "GRID.h" |
#include "GRID.h" |
22 |
|
#include "GAD.h" |
23 |
|
#ifdef ALLOW_LONGSTEP |
24 |
|
#include "LONGSTEP.h" |
25 |
|
#endif |
26 |
#ifdef ALLOW_GMREDI |
#ifdef ALLOW_GMREDI |
27 |
#include "GMREDI.h" |
#include "GMREDI.h" |
28 |
#endif |
#endif |
54 |
DO k=1,MIN(Nr,kSize) |
DO k=1,MIN(Nr,kSize) |
55 |
DO j=jMin,jMax |
DO j=jMin,jMax |
56 |
DO i=iMin,iMax |
DO i=iMin,iMax |
57 |
|
IF (tracerIdentity .LT. GAD_TR1) THEN |
58 |
|
KappaRx(i,j,k) = KappaRx(i,j,k)+Kwz(i,j,k,bi,bj) |
59 |
|
ELSE |
60 |
|
#ifdef ALLOW_LONGSTEP |
61 |
|
KappaRx(i,j,k) = KappaRx(i,j,k)+LS_Kwz(i,j,k,bi,bj) |
62 |
|
#else |
63 |
KappaRx(i,j,k) = KappaRx(i,j,k)+Kwz(i,j,k,bi,bj) |
KappaRx(i,j,k) = KappaRx(i,j,k)+Kwz(i,j,k,bi,bj) |
64 |
|
#endif |
65 |
|
ENDIF |
66 |
ENDDO |
ENDDO |
67 |
ENDDO |
ENDDO |
68 |
ENDDO |
ENDDO |
71 |
k = MIN(kArg,kSize) |
k = MIN(kArg,kSize) |
72 |
DO j=jMin,jMax |
DO j=jMin,jMax |
73 |
DO i=iMin,iMax |
DO i=iMin,iMax |
74 |
|
IF (tracerIdentity .LT. GAD_TR1) THEN |
75 |
KappaRx(i,j,k) = KappaRx(i,j,k)+Kwz(i,j,kArg,bi,bj) |
KappaRx(i,j,k) = KappaRx(i,j,k)+Kwz(i,j,kArg,bi,bj) |
76 |
|
ELSE |
77 |
|
#ifdef ALLOW_LONGSTEP |
78 |
|
KappaRx(i,j,k) = KappaRx(i,j,k)+LS_Kwz(i,j,kArg,bi,bj) |
79 |
|
#else |
80 |
|
KappaRx(i,j,k) = KappaRx(i,j,k)+Kwz(i,j,kArg,bi,bj) |
81 |
|
#endif |
82 |
|
ENDIF |
83 |
ENDDO |
ENDDO |
84 |
ENDDO |
ENDDO |
85 |
ENDIF |
ENDIF |