/[MITgcm]/MITgcm_contrib/heimbach/OpenAD/code_regress/externalDummies.F
ViewVC logotype

Diff of /MITgcm_contrib/heimbach/OpenAD/code_regress/externalDummies.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by utke, Thu Sep 18 15:23:43 2008 UTC revision 1.8 by utke, Fri Nov 8 17:55:05 2013 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3  C ###########################################################  C ###########################################################
4  C      SUBROUTINE EXCH_UV_XY_RS(        SUBROUTINE EXCH1_RL(
5  C     U                       Uphi, Vphi,       U                 array,
6  C     I                       withSigns, myThid )       I                 myOLw, myOLe, myOLs, myOLn, myNz,
7  C      IMPLICIT NONE       I                 exchWidthX, exchWidthY,
8  C#include "SIZE.h"       I                 cornerMode, myThid )
9  C#include "EEPARAMS.h"  
 C#include "EESUPPORT.h"  
 C#include "EXCH.h"  
 C  
 C      _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
 C      _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
 C      _RS temp  
 C      LOGICAL withSigns  
 C      INTEGER myThid  
 C      temp=Uphi(1,1,1,1)  
 C      Uphi(1,1,1,1)=Vphi(1,1,1,1)  
 C      Vphi(1,1,1,1)=temp  
 C      end subroutine  
 C ###########################################################  
       SUBROUTINE EXCH_UV_XYZ_RL(  
      U                       Uphi, Vphi,  
      I                       withSigns, myThid )  
10        IMPLICIT NONE        IMPLICIT NONE
11  #include "SIZE.h"  #include "SIZE.h"
12  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 #include "EESUPPORT.h"  
13  #include "EXCH.h"  #include "EXCH.h"
14        _RL Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        INTEGER myOLw, myOLe, myOLs, myOLn, myNz
15        _RL Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL     array( 1-myOLw:sNx+myOLe,
16        _RL temp       &               1-myOLs:sNy+myOLn,
17        LOGICAL withSigns       &               myNz, nSx, nSy )
18          INTEGER exchWidthX
19          INTEGER exchWidthY
20          INTEGER cornerMode
21        INTEGER myThid        INTEGER myThid
22        temp=Uphi(1,1,1,1,1)  C dummy self dependence (nontrivial so mfef90 doesn't kill it)
23        Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1)        array(1,1,1,1,1)=2*array(1,1,1,1,1)
       Vphi(1,1,1,1,1)=temp  
24        end subroutine        end subroutine
25  C ###########################################################  C ###########################################################
26  C      SUBROUTINE EXCH_UV_XYZ_RS(  C      SUBROUTINE EXCH1_RS(
27  C     U                       Uphi, Vphi,  C     U                 array,
28  C     I                       withSigns, myThid )  C     I                 myOLw, myOLe, myOLs, myOLn, myNz,
29    C     I                 exchWidthX, exchWidthY,
30    C     I                 cornerMode, myThid )
31    C
32  C      IMPLICIT NONE  C      IMPLICIT NONE
33  C#include "SIZE.h"  C#include "SIZE.h"
34  C#include "EEPARAMS.h"  C#include "EEPARAMS.h"
 C#include "EESUPPORT.h"  
35  C#include "EXCH.h"  C#include "EXCH.h"
36  C      _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  C      INTEGER myOLw, myOLe, myOLs, myOLn, myNz
37  C      _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  C      _RS     array( 1-myOLw:sNx+myOLe,
38  C      _RL temp  C     &               1-myOLs:sNy+myOLn,
39  C      LOGICAL withSigns  C     &               myNz, nSx, nSy )
40    C      INTEGER exchWidthX
41    C      INTEGER exchWidthY
42    C      INTEGER cornerMode
43  C      INTEGER myThid  C      INTEGER myThid
 C      temp=Uphi(1,1,1,1,1)  
 C      Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1)  
 C      Vphi(1,1,1,1,1)=temp  
44  C      end subroutine  C      end subroutine
45  C ###########################################################  C ###########################################################
       SUBROUTINE EXCH_XY_RL(  
      U                       phi,  
      I                       myThid )  
       IMPLICIT NONE  
 #include "SIZE.h"  
 #include "EEPARAMS.h"  
 #include "EESUPPORT.h"  
 #include "EXCH.h"  
       _RL phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
       INTEGER myThid  
       phi(1,1,1,1)=2*phi(1,1,1,1)  
       end subroutine  
 C ###########################################################  
       SUBROUTINE EXCH_XYZ_RL(  
      U                       phi,  
      I                       myThid )  
       IMPLICIT NONE  
 #include "SIZE.h"  
 #include "EEPARAMS.h"  
 #include "EESUPPORT.h"  
 #include "EXCH.h"  
       _RL phi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       INTEGER myThid  
       phi(1,1,1,1,1)=2*phi(1,1,1,1,1)  
       end subroutine  
 CC ###########################################################  
46  C      SUBROUTINE GLOBAL_MAX_R8(  C      SUBROUTINE GLOBAL_MAX_R8(
47  C     U                       maxphi,  C     U                       maxphi,
48  C     I                       myThid )  C     I                       myThid )
# Line 93  C      INTEGER myThid Line 56  C      INTEGER myThid
56  C      maxPhi=2*maxPhi  C      maxPhi=2*maxPhi
57  C      end subroutine  C      end subroutine
58  C ###########################################################  C ###########################################################
59        SUBROUTINE GLOBAL_SUM_R8(  C      SUBROUTINE GLOBAL_SUM_R8(
60    C     U                       sumphi,
61    C     I                       myThid )
62    C      IMPLICIT NONE
63    C#include "SIZE.h"
64    C#include "EEPARAMS.h"
65    C#include "EESUPPORT.h"
66    C#include "EXCH.h"
67    C      Real*8 sumPhi
68    C      INTEGER myThid
69    CC dummy self dependence (nontrivial so mfef90 doesn't kill it)
70    C      sumPhi=2*sumPhi
71    C      end subroutine
72    C ###########################################################
73          SUBROUTINE GLOBAL_SUM_TILE_RL(
74         U                       phiTile,
75       U                       sumphi,       U                       sumphi,
76       I                       myThid )       I                       myThid )
77        IMPLICIT NONE        IMPLICIT NONE
# Line 101  C ###################################### Line 79  C ######################################
79  #include "EEPARAMS.h"  #include "EEPARAMS.h"
80  #include "EESUPPORT.h"  #include "EESUPPORT.h"
81  #include "EXCH.h"  #include "EXCH.h"
82        Real*8 sumPhi        _RL phiTile(nSx,nSy)
83          _RL sumPhi
84        INTEGER myThid        INTEGER myThid
85        sumPhi=2*sumPhi  C dummy self dependence (nontrivial so mfef90 doesn't kill it)
86          sumPhi=2*phiTile(1,1)
87        end subroutine        end subroutine
88  C ###########################################################  C ###########################################################
89        SUBROUTINE CG2D(        SUBROUTINE CG2D(
90       I                cg2d_b,       I                cg2d_b,
91       U                cg2d_x,       U                cg2d_x,
92       O                firstResidual,       O                firstResidual,
93         O                minResidualSq,
94       O                lastResidual,       O                lastResidual,
95       U                numIters,       U                numIters,
96         O                nIterMin,
97       I                myThid )       I                myThid )
98        IMPLICIT NONE        IMPLICIT NONE
99  #include "SIZE.h"  #include "SIZE.h"
# Line 121  C ###################################### Line 103  C ######################################
103        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)
104        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)
105        Real*8  firstResidual        Real*8  firstResidual
106          Real*8  minResidualSq
107        Real*8  lastResidual        Real*8  lastResidual
108        INTEGER numIters        INTEGER numIters
109          INTEGER nIterMin
110        INTEGER myThid        INTEGER myThid
111        cg2d_x(1,1,1,1)=cg2d_b(1,1,1,1)  C dummy self dependence (nontrivial so mfef90 doesn't kill it)
112          cg2d_x(1,1,1,1)=2*cg2d_b(1,1,1,1)
113        end subroutine        end subroutine
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22