18 |
Vphi(1,1,1,1)=temp |
Vphi(1,1,1,1)=temp |
19 |
end subroutine |
end subroutine |
20 |
C ########################################################### |
C ########################################################### |
21 |
|
SUBROUTINE EXCH_UV_XY_RL( |
22 |
|
U Uphi, Vphi, |
23 |
|
I withSigns, myThid ) |
24 |
|
IMPLICIT NONE |
25 |
|
#include "SIZE.h" |
26 |
|
#include "EEPARAMS.h" |
27 |
|
#include "EESUPPORT.h" |
28 |
|
#include "EXCH.h" |
29 |
|
|
30 |
|
_RL Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
31 |
|
_RL Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
32 |
|
_RL temp |
33 |
|
LOGICAL withSigns |
34 |
|
INTEGER myThid |
35 |
|
temp=Uphi(1,1,1,1) |
36 |
|
Uphi(1,1,1,1)=Vphi(1,1,1,1) |
37 |
|
Vphi(1,1,1,1)=temp |
38 |
|
end subroutine |
39 |
|
C ########################################################### |
40 |
SUBROUTINE EXCH_UV_XYZ_RL( |
SUBROUTINE EXCH_UV_XYZ_RL( |
41 |
U Uphi, Vphi, |
U Uphi, Vphi, |
42 |
I withSigns, myThid ) |
I withSigns, myThid ) |
55 |
Vphi(1,1,1,1,1)=temp |
Vphi(1,1,1,1,1)=temp |
56 |
end subroutine |
end subroutine |
57 |
C ########################################################### |
C ########################################################### |
58 |
SUBROUTINE EXCH_UV_XYZ_RS( |
C SUBROUTINE EXCH_UV_XYZ_RS( |
59 |
U Uphi, Vphi, |
C U Uphi, Vphi, |
60 |
I withSigns, myThid ) |
C I withSigns, myThid ) |
61 |
|
C IMPLICIT NONE |
62 |
|
C#include "SIZE.h" |
63 |
|
C#include "EEPARAMS.h" |
64 |
|
C#include "EESUPPORT.h" |
65 |
|
C#include "EXCH.h" |
66 |
|
C _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
67 |
|
C _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
68 |
|
C _RL temp |
69 |
|
C LOGICAL withSigns |
70 |
|
C INTEGER myThid |
71 |
|
C temp=Uphi(1,1,1,1,1) |
72 |
|
C Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1) |
73 |
|
C Vphi(1,1,1,1,1)=temp |
74 |
|
C end subroutine |
75 |
|
C ########################################################### |
76 |
|
SUBROUTINE EXCH_UV_AGRID_3D_RL( |
77 |
|
U Uphi, Vphi, |
78 |
|
I withSigns, myNz, myThid ) |
79 |
IMPLICIT NONE |
IMPLICIT NONE |
80 |
#include "SIZE.h" |
#include "SIZE.h" |
81 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
82 |
#include "EESUPPORT.h" |
INTEGER myNz |
83 |
#include "EXCH.h" |
_RL Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
84 |
_RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
_RL Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
|
_RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
|
|
_RL temp |
|
85 |
LOGICAL withSigns |
LOGICAL withSigns |
86 |
INTEGER myThid |
INTEGER myThid |
87 |
temp=Uphi(1,1,1,1,1) |
temp=Uphi(1,1,1,1,1) |
89 |
Vphi(1,1,1,1,1)=temp |
Vphi(1,1,1,1,1)=temp |
90 |
end subroutine |
end subroutine |
91 |
C ########################################################### |
C ########################################################### |
92 |
|
C SUBROUTINE EXCH_UV_AGRID_3D_RS( |
93 |
|
C U Uphi, Vphi, |
94 |
|
C I withSigns, myNz, myThid ) |
95 |
|
C IMPLICIT NONE |
96 |
|
C#include "SIZE.h" |
97 |
|
C#include "EEPARAMS.h" |
98 |
|
C INTEGER myNz |
99 |
|
C _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
100 |
|
C _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
101 |
|
C LOGICAL withSigns |
102 |
|
C INTEGER myThid |
103 |
|
C temp=Uphi(1,1,1,1,1) |
104 |
|
C Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1) |
105 |
|
C Vphi(1,1,1,1,1)=temp |
106 |
|
C end subroutine |
107 |
|
C ########################################################### |
108 |
|
SUBROUTINE EXCH_UV_3D_RL( |
109 |
|
U Uphi, Vphi, |
110 |
|
I withSigns, myNz, myThid ) |
111 |
|
IMPLICIT NONE |
112 |
|
#include "SIZE.h" |
113 |
|
#include "EEPARAMS.h" |
114 |
|
INTEGER myNz |
115 |
|
_RL Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
116 |
|
_RL Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
117 |
|
LOGICAL withSigns |
118 |
|
INTEGER myThid |
119 |
|
temp=Uphi(1,1,1,1,1) |
120 |
|
Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1) |
121 |
|
Vphi(1,1,1,1,1)=temp |
122 |
|
end subroutine |
123 |
|
C ########################################################### |
124 |
|
SUBROUTINE EXCH_XY_RS( |
125 |
|
U phi, |
126 |
|
I myThid ) |
127 |
|
IMPLICIT NONE |
128 |
|
#include "SIZE.h" |
129 |
|
#include "EEPARAMS.h" |
130 |
|
#include "EESUPPORT.h" |
131 |
|
#include "EXCH.h" |
132 |
|
_RS phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
133 |
|
INTEGER myThid |
134 |
|
phi(1,1,1,1)=2*phi(1,1,1,1) |
135 |
|
end subroutine |
136 |
|
C ########################################################### |
137 |
SUBROUTINE EXCH_XY_RL( |
SUBROUTINE EXCH_XY_RL( |
138 |
U phi, |
U phi, |
139 |
I myThid ) |
I myThid ) |
159 |
INTEGER myThid |
INTEGER myThid |
160 |
phi(1,1,1,1,1)=2*phi(1,1,1,1,1) |
phi(1,1,1,1,1)=2*phi(1,1,1,1,1) |
161 |
end subroutine |
end subroutine |
162 |
CC ########################################################### |
C ########################################################### |
163 |
C SUBROUTINE GLOBAL_MAX_R8( |
C SUBROUTINE GLOBAL_MAX_R8( |
164 |
C U maxphi, |
C U maxphi, |
165 |
C I myThid ) |
C I myThid ) |
173 |
C maxPhi=2*maxPhi |
C maxPhi=2*maxPhi |
174 |
C end subroutine |
C end subroutine |
175 |
C ########################################################### |
C ########################################################### |
176 |
|
SUBROUTINE EXCH_3D_RL( |
177 |
|
U phi, |
178 |
|
I myNz, myThid ) |
179 |
|
IMPLICIT NONE |
180 |
|
#include "SIZE.h" |
181 |
|
#include "EEPARAMS.h" |
182 |
|
#include "EESUPPORT.h" |
183 |
|
#include "EXCH.h" |
184 |
|
INTEGER myNz |
185 |
|
_RL phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,1:myNz,nSx,nSy) |
186 |
|
INTEGER myThid |
187 |
|
phi(1,1,1,1,1)=2*phi(1,1,1,1,1) |
188 |
|
end subroutine |
189 |
|
C ########################################################### |
190 |
|
C SUBROUTINE EXCH_Z_3D_RS( |
191 |
|
C U phi, |
192 |
|
C I myNz, myThid ) |
193 |
|
C IMPLICIT NONE |
194 |
|
C#include "SIZE.h" |
195 |
|
C#include "EEPARAMS.h" |
196 |
|
C#include "EESUPPORT.h" |
197 |
|
C#include "EXCH.h" |
198 |
|
C INTEGER myNz |
199 |
|
C _RS phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,myNz,nSx,nSy) |
200 |
|
C INTEGER myThid |
201 |
|
C phi(1,1,1,1,1)=2*phi(1,1,1,1,1) |
202 |
|
C end subroutine |
203 |
|
C ########################################################### |
204 |
SUBROUTINE GLOBAL_SUM_R8( |
SUBROUTINE GLOBAL_SUM_R8( |
205 |
U sumphi, |
U sumphi, |
206 |
I myThid ) |
I myThid ) |
211 |
#include "EXCH.h" |
#include "EXCH.h" |
212 |
Real*8 sumPhi |
Real*8 sumPhi |
213 |
INTEGER myThid |
INTEGER myThid |
214 |
|
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
215 |
sumPhi=2*sumPhi |
sumPhi=2*sumPhi |
216 |
end subroutine |
end subroutine |
217 |
C ########################################################### |
C ########################################################### |
219 |
I cg2d_b, |
I cg2d_b, |
220 |
U cg2d_x, |
U cg2d_x, |
221 |
O firstResidual, |
O firstResidual, |
222 |
|
O minResidualSq, |
223 |
O lastResidual, |
O lastResidual, |
224 |
U numIters, |
U numIters, |
225 |
|
O nIterMin, |
226 |
I myThid ) |
I myThid ) |
227 |
IMPLICIT NONE |
IMPLICIT NONE |
228 |
#include "SIZE.h" |
#include "SIZE.h" |
232 |
Real*8 cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
Real*8 cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
233 |
Real*8 cg2d_x(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
Real*8 cg2d_x(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
234 |
Real*8 firstResidual |
Real*8 firstResidual |
235 |
|
Real*8 minResidualSq |
236 |
Real*8 lastResidual |
Real*8 lastResidual |
237 |
INTEGER numIters |
INTEGER numIters |
238 |
|
INTEGER nIterMin |
239 |
INTEGER myThid |
INTEGER myThid |
240 |
cg2d_x(1,1,1,1)=cg2d_b(1,1,1,1) |
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
241 |
|
cg2d_x(1,1,1,1)=2*cg2d_b(1,1,1,1) |
242 |
end subroutine |
end subroutine |
243 |
|
C ########################################################### |