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 |
IMPLICIT NONE |
C IMPLICIT NONE |
62 |
#include "SIZE.h" |
C#include "SIZE.h" |
63 |
#include "EEPARAMS.h" |
C#include "EEPARAMS.h" |
64 |
#include "EESUPPORT.h" |
C#include "EESUPPORT.h" |
65 |
#include "EXCH.h" |
C#include "EXCH.h" |
66 |
_RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
C _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
67 |
_RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
C _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
68 |
_RL temp |
C _RL temp |
69 |
LOGICAL withSigns |
C LOGICAL withSigns |
70 |
INTEGER myThid |
C INTEGER myThid |
71 |
temp=Uphi(1,1,1,1,1) |
C temp=Uphi(1,1,1,1,1) |
72 |
Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1) |
C Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1) |
73 |
Vphi(1,1,1,1,1)=temp |
C Vphi(1,1,1,1,1)=temp |
74 |
end subroutine |
C end subroutine |
75 |
C ########################################################### |
C ########################################################### |
76 |
SUBROUTINE EXCH_UV_AGRID_3D_RL( |
SUBROUTINE EXCH_UV_AGRID_3D_RL( |
77 |
U Uphi, Vphi, |
U Uphi, Vphi, |
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 ########################################################### |
C ########################################################### |