/[MITgcm]/MITgcm/pkg/gmredi/gmredi_calc_psi_b.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_calc_psi_b.F

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

revision 1.1 by jmc, Sun Dec 16 18:54:49 2001 UTC revision 1.2 by heimbach, Thu Nov 14 22:43:49 2002 UTC
# Line 24  C     == Global variables == Line 24  C     == Global variables ==
24  #include "GMREDI.h"  #include "GMREDI.h"
25  #include "GMREDI_DIAGS.h"  #include "GMREDI_DIAGS.h"
26    
27    #ifdef ALLOW_AUTODIFF_TAMC
28    #include "tamc.h"
29    #include "tamc_keys.h"
30    #endif /* ALLOW_AUTODIFF_TAMC */
31    
32  C     == Routine arguments ==  C     == Routine arguments ==
33  C  C
34        _RL sigmaX(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL sigmaX(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
# Line 47  C     == Local variables == Line 52  C     == Local variables ==
52    
53  C-    Initialization : <= done in S/R gmredi_init  C-    Initialization : <= done in S/R gmredi_init
54    
55    #ifdef ALLOW_AUTODIFF_TAMC
56              act1 = bi - myBxLo(myThid)
57              max1 = myBxHi(myThid) - myBxLo(myThid) + 1
58              act2 = bj - myByLo(myThid)
59              max2 = myByHi(myThid) - myByLo(myThid) + 1
60              act3 = myThid - 1
61              max3 = nTx*nTy
62              act4 = ikey_dynamics - 1
63              ikey = (act1 + 1) + act2*max1
64         &                      + act3*max1*max2
65         &                      + act4*max1*max2*max3
66    #endif /* ALLOW_AUTODIFF_TAMC */
67    
68        IF (GM_AdvForm) THEN        IF (GM_AdvForm) THEN
69         DO k=2,Nr         DO k=2,Nr
70  C-----         km1 = k-1
71    
72    #ifdef ALLOW_AUTODIFF_TAMC
73           kkey = (ikey-1)*Nr + k
74           DO j=1-Oly,sNy+Oly
75            DO i=1-Olx,sNx+Olx
76             SlopeX(i,j)       = 0. _d 0
77             SlopeY(i,j)       = 0. _d 0
78             dSigmaDrW(i,j)    = 0. _d 0
79             dSigmaDrS(i,j)    = 0. _d 0
80            ENDDO
81           ENDDO
82    #endif
83    
       km1 = k-1  
84        DO j=1-Oly+1,sNy+Oly-1        DO j=1-Oly+1,sNy+Oly-1
85         DO i=1-Olx+1,sNx+Olx-1         DO i=1-Olx+1,sNx+Olx-1
86    
# Line 71  C      Gradient of Sigma below U and V p Line 100  C      Gradient of Sigma below U and V p
100  C     Calculate slopes , taper and/or clip  C     Calculate slopes , taper and/or clip
101        CALL GMREDI_SLOPE_PSI_B(        CALL GMREDI_SLOPE_PSI_B(
102       I             dSigmaDrW, dSigmaDrS,       I             dSigmaDrW, dSigmaDrS,
103       I             rF(K),       I             rF(K),K,
104       U             SlopeX, SlopeY,       U             SlopeX, SlopeY,
105       O             taperX, taperY,       O             taperX, taperY,
106       I             bi, bj, myThid )       I             bi, bj, myThid )
107    
108    #ifdef ALLOW_AUTODIFF_TAMC
109    CADJ STORE SlopeX(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
110    CADJ STORE SlopeY(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
111    CADJ STORE taperX(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
112    CADJ STORE taperY(:,:)       = comlev1_bibj_k, key=kkey, byte=isbyte
113    #endif /* ALLOW_AUTODIFF_TAMC */
114    
115        DO j=1-Oly+1,sNy+Oly-1        DO j=1-Oly+1,sNy+Oly-1
116         DO i=1-Olx+1,sNx+Olx-1         DO i=1-Olx+1,sNx+Olx-1
117    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22