/[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.1 by utke, Tue Nov 20 15:20:57 2007 UTC revision 1.6 by utke, Sat May 12 13:26:48 2012 UTC
# Line 1  Line 1 
 cccC ###########################################################  
 ccc      subroutine active_read_xy_loc(  
 ccc     I                           active_var_file,  
 ccc     O                           active_var,  
 ccc     I                           irec,  
 ccc     I                           doglobalread,  
 ccc     I                           lAdInit,  
 ccc     I                           myOptimIter,  
 ccc     I                           mythid  
 ccc     I                         , dummy  
 ccc     &                         )  
 ccc      implicit none  
 ccc#include "EEPARAMS.h"  
 ccc#include "SIZE.h"  
 ccc      character*(*) active_var_file  
 ccc      _RL     active_var(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
 ccc      integer irec  
 ccc      integer myOptimIter  
 ccc      integer mythid  
 ccc      logical doglobalread  
 ccc      logical lAdInit  
 ccc      _RL     dummy (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)  
 ccc      dummy(1,1,1,1)=active_var(1,1,1,1)  
 ccc      end subroutine  
 cccC ###########################################################  
 ccc      subroutine active_read_xyz_loc(  
 ccc     I                            active_var_file,  
 ccc     O                            active_var,  
 ccc     I                            irec,  
 ccc     I                            doglobalread,  
 ccc     I                            lAdInit,  
 ccc     I                            myOptimIter,  
 ccc     I                            mythid  
 ccc     I                         , dummy  
 ccc     &                           )  
 ccc#include "EEPARAMS.h"  
 ccc#include "SIZE.h"  
 ccc      character*(*) active_var_file  
 ccc      _RL active_var(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 ccc      integer irec  
 ccc      integer myOptimIter  
 ccc      integer mythid  
 ccc      logical doglobalread  
 ccc      logical lAdInit  
 ccc      _RL     dummy(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 ccc      dummy(1,1,1,1,1)=active_var(1,1,1,1,1)  
 ccc      end subroutine  
 cccC ###########################################################  
 ccc      subroutine active_read_xyz(  
 ccc     I                            active_var_file,  
 ccc     O                            active_var,  
 ccc     I                            irec,  
 ccc     I                            doglobalread,  
 ccc     I                            lAdInit,  
 ccc     I                            myOptimIter,  
 ccc     I                            mythid  
 ccc     I                         , dummy  
 ccc     &                           )  
 ccc      implicit none  
 ccc#include "EEPARAMS.h"  
 ccc#include "SIZE.h"  
 ccc      character*(*) active_var_file  
 ccc      _RL active_var(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 ccc      integer irec  
 ccc      integer myOptimIter  
 ccc      integer mythid  
 ccc      logical doglobalread  
 ccc      logical lAdInit  
 ccc      _RL     dummy(1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
 ccc      dummy(1,1,1,1,1)=active_var(1,1,1,1,1)  
 ccc      end subroutine  
 cccC ###########################################################  
 ccc      SUBROUTINE EXCH_UV_AGRID_XY_RS( component1,  
 ccc     +                                component2,  
 ccc     +                                withSigns,  
 ccc     +                                myThid )  
 ccc      implicit none  
 ccc#include "SIZE.h"  
 ccc#include "EEPARAMS.h"  
 ccc#include "EESUPPORT.h"  
 ccc#include "EXCH.h"  
 ccc      _RS component1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
 ccc      _RS component2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
 ccc      LOGICAL withSigns  
 ccc      INTEGER myThid  
 ccc      component1(1,1,1,1)=component2(1,1,1,1)  
 ccc      component2(1,1,1,1)=component1(1,1,1,1)  
 ccc      end subroutine  
1  C ###########################################################  C ###########################################################
2        SUBROUTINE EXCH_UV_XY_RS(        SUBROUTINE EXCH1_RL(
3       U                       Uphi, Vphi,       U                 array,
4       I                       withSigns, myThid )       I                 myOLw, myOLe, myOLs, myOLn, myNz,
5        IMPLICIT NONE       I                 exchWidthX, exchWidthY,
6  #include "SIZE.h"       I                 cornerMode, myThid )
 #include "EEPARAMS.h"  
 #include "EESUPPORT.h"  
 #include "EXCH.h"  
7    
       _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
       _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  
       _RS temp  
       LOGICAL withSigns  
       INTEGER myThid  
       temp=Uphi(1,1,1,1)  
       Uphi(1,1,1,1)=Vphi(1,1,1,1)  
       Vphi(1,1,1,1)=temp  
       end subroutine  
 C ###########################################################  
       SUBROUTINE EXCH_UV_XYZ_RL(  
      U                       Uphi, Vphi,  
      I                       withSigns, myThid )  
8        IMPLICIT NONE        IMPLICIT NONE
9  #include "SIZE.h"  #include "SIZE.h"
10  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 #include "EESUPPORT.h"  
11  #include "EXCH.h"  #include "EXCH.h"
12        _RL Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        INTEGER myOLw, myOLe, myOLs, myOLn, myNz
13        _RL Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL     array( 1-myOLw:sNx+myOLe,
14        _RL temp       &               1-myOLs:sNy+myOLn,
15        LOGICAL withSigns       &               myNz, nSx, nSy )
16        INTEGER myThid        INTEGER exchWidthX
17        temp=Uphi(1,1,1,1,1)        INTEGER exchWidthY
18        Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1)        INTEGER cornerMode
19        Vphi(1,1,1,1,1)=temp        INTEGER myThid
20        end subroutine  C dummy self dependence (nontrivial so mfef90 doesn't kill it)
21          array(1,1,1,1,1)=2*array(1,1,1,1,1)
22          end subroutine
23    C ###########################################################
24    C      SUBROUTINE EXCH1_RS(
25    C     U                 array,
26    C     I                 myOLw, myOLe, myOLs, myOLn, myNz,
27    C     I                 exchWidthX, exchWidthY,
28    C     I                 cornerMode, myThid )
29    C
30    C      IMPLICIT NONE
31    C#include "SIZE.h"
32    C#include "EEPARAMS.h"
33    C#include "EXCH.h"
34    C      INTEGER myOLw, myOLe, myOLs, myOLn, myNz
35    C      _RS     array( 1-myOLw:sNx+myOLe,
36    C     &               1-myOLs:sNy+myOLn,
37    C     &               myNz, nSx, nSy )
38    C      INTEGER exchWidthX
39    C      INTEGER exchWidthY
40    C      INTEGER cornerMode
41    C      INTEGER myThid
42    C      end subroutine
43    C ###########################################################
44    C      SUBROUTINE GLOBAL_MAX_R8(
45    C     U                       maxphi,
46    C     I                       myThid )
47    C      IMPLICIT NONE
48    C#include "SIZE.h"
49    C#include "EEPARAMS.h"
50    C#include "EESUPPORT.h"
51    C#include "EXCH.h"
52    C      Real*8 maxPhi
53    C      INTEGER myThid
54    C      maxPhi=2*maxPhi
55    C      end subroutine
56  C ###########################################################  C ###########################################################
57        SUBROUTINE EXCH_UV_XYZ_RS(        SUBROUTINE GLOBAL_SUM_R8(
58       U                       Uphi, Vphi,       U                       sumphi,
      I                       withSigns, myThid )  
       IMPLICIT NONE  
 #include "SIZE.h"  
 #include "EEPARAMS.h"  
 #include "EESUPPORT.h"  
 #include "EXCH.h"  
       _RS Uphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RS Vphi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)  
       _RL temp  
       LOGICAL withSigns  
       INTEGER myThid  
       temp=Uphi(1,1,1,1,1)  
       Uphi(1,1,1,1,1)=Vphi(1,1,1,1,1)  
       Vphi(1,1,1,1,1)=temp  
       end subroutine  
 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  
 C ###########################################################  
       SUBROUTINE GLOBAL_MAX_R8(  
      U                       maxphi,  
59       I                       myThid )       I                       myThid )
60        IMPLICIT NONE        IMPLICIT NONE
61  #include "SIZE.h"  #include "SIZE.h"
62  #include "EEPARAMS.h"  #include "EEPARAMS.h"
63  #include "EESUPPORT.h"  #include "EESUPPORT.h"
64  #include "EXCH.h"  #include "EXCH.h"
65        Real*8 maxPhi        Real*8 sumPhi
66        INTEGER myThid        INTEGER myThid
67        maxPhi=2*maxPhi  C dummy self dependence (nontrivial so mfef90 doesn't kill it)
68          sumPhi=2*sumPhi
69        end subroutine        end subroutine
70  C ###########################################################  C ###########################################################
71        SUBROUTINE GLOBAL_SUM_R8(        SUBROUTINE CG2D(
72       U                       sumphi,       I                cg2d_b,
73       I                       myThid )       U                cg2d_x,
74         O                firstResidual,
75         O                minResidualSq,
76         O                lastResidual,
77         U                numIters,
78         O                nIterMin,
79         I                myThid )
80        IMPLICIT NONE        IMPLICIT NONE
81  #include "SIZE.h"  #include "SIZE.h"
82  #include "EEPARAMS.h"  #include "EEPARAMS.h"
83  #include "EESUPPORT.h"  #include "PARAMS.h"
84  #include "EXCH.h"  #include "CG2D.h"
85        Real*8 sumPhi        Real*8  cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86          Real*8  cg2d_x(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87          Real*8  firstResidual
88          Real*8  minResidualSq
89          Real*8  lastResidual
90          INTEGER numIters
91          INTEGER nIterMin
92        INTEGER myThid        INTEGER myThid
93        sumPhi=2*sumPhi  C dummy self dependence (nontrivial so mfef90 doesn't kill it)
94          cg2d_x(1,1,1,1)=2*cg2d_b(1,1,1,1)
95        end subroutine        end subroutine
 C ###########################################################  
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22