pkg/autodiff/adzero_adj.F

C $Header: /u/gcmpack/MITgcm/pkg/autodiff/adzero_adj.F,v 1.5 2010/02/15 20:03:56 gforget Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

C     !ROUTINE: ADZERO_ADJ_1D
C     !INTERFACE:
      SUBROUTINE ADZERO_ADJ_1D(
     I             NrIn, adFieldFwd, myThid)
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | o zero out a 1D adjoint field
C     | author: Gael Forget
C     *==========================================================*

      IMPLICIT NONE

C     == Global arguments ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
      integer NrIn,myThid
      _RL adFieldFwd(NrIn)

C     == Local arguments ==
      integer bi,bj,i,j,k

      DO k=1,NrIn
         adFieldFwd(k) = 0. _d 0
      ENDDO

      end


C     !ROUTINE: ADZERO_ADJ_LOC
C     !INTERFACE:
      SUBROUTINE ADZERO_ADJ_LOC(
     I             NrIn, adFieldFwd, myThid)
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | o zero out a local (no nsx,nsy dims.) adjoint field
C     | author: Gael Forget
C     *==========================================================*

      IMPLICIT NONE

C     == Global arguments ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
      integer NrIn,myThid
      _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)

C     == Local arguments ==
      integer bi,bj,i,j,k

      DO k=1,NrIn
       DO j=1-Oly,sNy+Oly
        DO i=1-Olx,sNx+Olx
         adFieldFwd(i,j,k) = 0. _d 0
        ENDDO
       ENDDO
      ENDDO

      end


C     !ROUTINE: ADZERO_ADJ
C     !INTERFACE:
      SUBROUTINE ADZERO_ADJ(
     I             NrIn, adFieldFwd, myThid)
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | o zero out a full (incl. nsx,nsy dims.) adjoint field
C     | author: Gael Forget
C     *==========================================================*

      IMPLICIT NONE

C     == Global arguments ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     == Routine arguments ==
C     myThid - Thread number for this instance of the routine.
      integer NrIn,myThid
      _RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)

C     == Local arguments ==
      integer bi,bj,i,j,k

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO k=1,NrIn
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
           adFieldFwd(i,j,k,bi,bj) = 0. _d 0
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO

      end

1	gforget	1.6	C $Header: /u/gcmpack/MITgcm/pkg/autodiff/adzero_adj.F,v 1.5 2010/02/15 20:03:56 gforget Exp $
2	jmc	1.2	C $Name: $
3	gforget	1.1
4			#include "CPP_OPTIONS.h"
5
6	gforget	1.6	C !ROUTINE: ADZERO_ADJ_1D
7			C !INTERFACE:
8			SUBROUTINE ADZERO_ADJ_1D(
9			I NrIn, adFieldFwd, myThid)
10			C !DESCRIPTION: \bv
11			C ==========================================================
12			C \| o zero out a 1D adjoint field
13			C \| author: Gael Forget
14			C ==========================================================
15
16			IMPLICIT NONE
17
18			C == Global arguments ==
19			#include "SIZE.h"
20			#include "EEPARAMS.h"
21			#include "PARAMS.h"
22
23			C == Routine arguments ==
24			C myThid - Thread number for this instance of the routine.
25			integer NrIn,myThid
26			_RL adFieldFwd(NrIn)
27
28			C == Local arguments ==
29			integer bi,bj,i,j,k
30
31			DO k=1,NrIn
32			adFieldFwd(k) = 0. _d 0
33			ENDDO
34
35			end
36
37
38	gforget	1.3	C !ROUTINE: ADZERO_ADJ_LOC
39	gforget	1.1	C !INTERFACE:
40	gforget	1.3	SUBROUTINE ADZERO_ADJ_LOC(
41	gforget	1.5	I NrIn, adFieldFwd, myThid)
42	gforget	1.1	C !DESCRIPTION: \bv
43			C ==========================================================
44	gforget	1.3	C \| o zero out a local (no nsx,nsy dims.) adjoint field
45	gforget	1.1	C \| author: Gael Forget
46			C ==========================================================
47
48	heimbach	1.4	IMPLICIT NONE
49
50			C == Global arguments ==
51	gforget	1.1	#include "SIZE.h"
52			#include "EEPARAMS.h"
53			#include "PARAMS.h"
54
55			C == Routine arguments ==
56			C myThid - Thread number for this instance of the routine.
57	heimbach	1.4	integer NrIn,myThid
58	gforget	1.3	_RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn)
59	gforget	1.1
60	heimbach	1.4	C == Local arguments ==
61			integer bi,bj,i,j,k
62
63	gforget	1.3	DO k=1,NrIn
64	gforget	1.1	DO j=1-Oly,sNy+Oly
65			DO i=1-Olx,sNx+Olx
66	heimbach	1.4	adFieldFwd(i,j,k) = 0. _d 0
67	gforget	1.1	ENDDO
68			ENDDO
69			ENDDO
70
71			end
72
73	gforget	1.3
74			C !ROUTINE: ADZERO_ADJ
75	gforget	1.1	C !INTERFACE:
76	gforget	1.3	SUBROUTINE ADZERO_ADJ(
77	gforget	1.5	I NrIn, adFieldFwd, myThid)
78	gforget	1.1	C !DESCRIPTION: \bv
79			C ==========================================================
80	gforget	1.3	C \| o zero out a full (incl. nsx,nsy dims.) adjoint field
81	gforget	1.1	C \| author: Gael Forget
82			C ==========================================================
83
84	heimbach	1.4	IMPLICIT NONE
85
86			C == Global arguments ==
87	gforget	1.1	#include "SIZE.h"
88			#include "EEPARAMS.h"
89			#include "PARAMS.h"
90
91			C == Routine arguments ==
92			C myThid - Thread number for this instance of the routine.
93	heimbach	1.4	integer NrIn,myThid
94	gforget	1.3	_RL adFieldFwd(1-Olx:sNx+Olx,1-Oly:sNy+Oly,NrIn,nSx,nSy)
95	gforget	1.1
96	heimbach	1.4	C == Local arguments ==
97			integer bi,bj,i,j,k
98
99	gforget	1.3	DO bj=myByLo(myThid),myByHi(myThid)
100			DO bi=myBxLo(myThid),myBxHi(myThid)
101			DO k=1,NrIn
102			DO j=1-Oly,sNy+Oly
103			DO i=1-Olx,sNx+Olx
104	heimbach	1.4	adFieldFwd(i,j,k,bi,bj) = 0. _d 0
105	gforget	1.3	ENDDO
106			ENDDO
107	gforget	1.1	ENDDO
108			ENDDO
109	gforget	1.3	ENDDO
110	gforget	1.1
111			end
112