/[MITgcm]/MITgcm/pkg/autodiff/adzero_adj.F
ViewVC logotype

Diff of /MITgcm/pkg/autodiff/adzero_adj.F

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

revision 1.3 by gforget, Wed Jun 11 18:29:11 2008 UTC revision 1.4 by heimbach, Wed Jul 23 20:26:12 2008 UTC
# Line 13  C     | o zero out a local (no nsx,nsy d Line 13  C     | o zero out a local (no nsx,nsy d
13  C     | author: Gael Forget  C     | author: Gael Forget
14  C     *==========================================================*  C     *==========================================================*
15    
16          IMPLICIT NONE
17    
18    C     == Global arguments ==
19  #include "SIZE.h"  #include "SIZE.h"
20  #include "EEPARAMS.h"  #include "EEPARAMS.h"
21  #include "PARAMS.h"  #include "PARAMS.h"
22    
23  C     == Routine arguments ==  C     == Routine arguments ==
24  C     myThid - Thread number for this instance of the routine.  C     myThid - Thread number for this instance of the routine.
25        integer myThid,bi,bj,i,j,k        integer NrIn,myThid
26        _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)        _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)
27    
28    C     == Local arguments ==
29          integer bi,bj,i,j,k
30    
31        DO k=1,NrIn        DO k=1,NrIn
32         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
33          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
34           adFieldFwd(i,j,k) = 0.           adFieldFwd(i,j,k) = 0. _d 0
35          ENDDO          ENDDO
36         ENDDO         ENDDO
37        ENDDO        ENDDO
# Line 43  C     | o zero out a full (incl. nsx,nsy Line 49  C     | o zero out a full (incl. nsx,nsy
49  C     | author: Gael Forget  C     | author: Gael Forget
50  C     *==========================================================*  C     *==========================================================*
51    
52          IMPLICIT NONE
53    
54    C     == Global arguments ==
55  #include "SIZE.h"  #include "SIZE.h"
56  #include "EEPARAMS.h"  #include "EEPARAMS.h"
57  #include "PARAMS.h"  #include "PARAMS.h"
58    
59  C     == Routine arguments ==  C     == Routine arguments ==
60  C     myThid - Thread number for this instance of the routine.  C     myThid - Thread number for this instance of the routine.
61        integer myThid,bi,bj,i,j,k        integer NrIn,myThid
62        _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)        _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)
63    
64    C     == Local arguments ==
65          integer bi,bj,i,j,k
66    
67        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
68         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
69          DO k=1,NrIn          DO k=1,NrIn
70           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
71            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
72             adFieldFwd(i,j,k,bi,bj) = 0.             adFieldFwd(i,j,k,bi,bj) = 0. _d 0
73            ENDDO            ENDDO
74           ENDDO           ENDDO
75          ENDDO          ENDDO
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22