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 ########################################################### |
233 |
Real*8 lastResidual |
Real*8 lastResidual |
234 |
INTEGER numIters |
INTEGER numIters |
235 |
INTEGER myThid |
INTEGER myThid |
236 |
cg2d_x(1,1,1,1)=cg2d_b(1,1,1,1) |
C dummy self dependence (nontrivial so mfef90 doesn't kill it) |
237 |
|
cg2d_x(1,1,1,1)=2*cg2d_b(1,1,1,1) |
238 |
end subroutine |
end subroutine |
239 |
C ########################################################### |
C ########################################################### |